分享一个MATLAB求多项式的系数

Lx00591

matlab求多项式的系数:

- y/ }, c8 q/ `7 j: h单变量：

; y3 O, k' T* G$ o# Bsyms x y
f1=x^4+2*x+1;
) K: M% i9 c1 |! B' Q% Mf2=y^6+5*y^3+3;
f3=x^5+2*x^3*y^4+x*y^2+4;
- R. h# y3 N4 v; A( Y: Zc = sym2poly(f1)
c =
    1     0     0     2     1
>> c = sym2poly(f2)
c =
4 ^8 ^( H- B1 Z    1     0     0     5     0     0     3

多变量：
3 P. r( F, ?, O" f& q
unction coef=poly_coef(f,var)
%提取多项式f中指定变量var的系数，将结果赋给数组coef
%f可以是含多变量的多项式
%var是多项式中指定的变量，可选，默认是x
/ J! M9 v& f4 g8 r. z* K7 U5 g%要用到函数poly_degree()来获得f中指定变量var的最高次幂
IF nargin==1
var=sym('x');
end
degree=poly_degree(f,var);
temp_f=f;
coef(degree+1)=subs(temp_f,var,0);
9 G8 U0 A  h7 @$ Y! n+ v6 Rfor n=1:degree
temp_f=simple((temp_f-coef(degree+2-n))/var);
6 |; l! a- N2 T3 A, Ccoef(degree+1-n)=subs(temp_f,var,0);
end
end
0 i" [( g8 K& E7 b. T) m举几个例子：
复制内容到剪贴板
代码:

$ a$ ~$ d- G  H: ^. e9 B>> syms x y
>> f1=x^4+2*x+1;
) t8 q: u2 W/ x>> f2=y^6+5*y^3+3;
>> f3=x^5+2*x^3*y^4+x*y^2+4;
>> poly_coef(f1)
6 T/ P5 h- r6 e1 |
1 y. B0 u' d# _& {3 Z) D' v) G' V5 D) |ans =
9 z, C; \. ~* |/ p  }  Z6 B
* c8 k; h6 V+ T; r' `0 Q( i                               [1, 0, 0, 2, 1]
" `2 z4 W# I6 _4 ?( D: B>> poly_coef(f1,y)
" I4 G6 k, v7 n) i! h
* Z* M. B1 G# h7 |ans =
, g0 |" d% o! X+ H& {; l$ g
$ c+ `0 ^# g9 X  {( ~* t) m- a                               [ 4       ]
& T  {9 f( V1 M' t6 v" k                               [x   + 2 x + 1]
/ z0 F. t- b1 x( ~9 Y>> poly_coef(f2)

ans =
0 _. f% J) |1 Y9 O
                               [ 6    3 ]
                               [y   + 5 y   + 3]
' D3 X* H6 G, j) `) Y2 G  G>> poly_coef(f2,y)

ans =

                            [1, 0, 0, 5, 0, 0, 3]
1 l; o/ M: y# j1 V) V>> poly_coef(f3)
3 h$ a' m+ ?& c
( p8 D, q8 R, X: A9 V8 |ans =

1 q" K1 ^  w- Z5 ]: N1 r                        [       4    2 ]
  t4 m; W5 q8 [) X0 N4 G) H                        [1, 0, 2 y , 0, y , 4]
. t$ Y1 c5 q) m>> poly_coef(f3,y)

ans =

* S/ G7 d5 p* n- p' J                        [ 3          5 ]
0 X/ L4 \7 R8 p                        [2 x , 0, x, 0, x   + 4]

Demyar

matlab求多项式的系数。