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

Lx00591

matlab求多项式的系数:

/ d7 y" U& R# a( X) b6 X单变量：
5 O& z3 w% U: `: J/ ?/ m$ H
4 W" ]# U: ^: [* S: }syms x y
7 l" j4 R" H" t) l8 d- ~3 D2 [f1=x^4+2*x+1;
' a) B) S+ v5 }! A9 A0 Uf2=y^6+5*y^3+3;
f3=x^5+2*x^3*y^4+x*y^2+4;
c = sym2poly(f1)
c =
( H+ m/ }% E8 l6 m0 k1 l    1     0     0     2     1
>> c = sym2poly(f2)
/ j$ B" E  e* t/ I( C( K4 v; b$ cc =
    1     0     0     5     0     0     3
) ]. `: u. P7 I
多变量：

$ \8 t/ M4 Q$ n- L) n1 J% Zunction coef=poly_coef(f,var)
%提取多项式f中指定变量var的系数，将结果赋给数组coef
%f可以是含多变量的多项式
) M* A6 e7 a4 z%var是多项式中指定的变量，可选，默认是x
& u& o9 ]; @2 i' N%要用到函数poly_degree()来获得f中指定变量var的最高次幂
IF nargin==1
var=sym('x');
1 K7 F  n6 ~2 g# W/ aend
degree=poly_degree(f,var);
: U5 k. \* P# |5 U% t/ Qtemp_f=f;
' P' D$ v+ }- T$ ~, ~coef(degree+1)=subs(temp_f,var,0);
for n=1:degree
+ {. b" S6 V1 O( ctemp_f=simple((temp_f-coef(degree+2-n))/var);
7 n" F/ n1 Z2 `coef(degree+1-n)=subs(temp_f,var,0);
  ]4 k5 \: P* E; N5 x$ I+ _end
: [) p" }. O! vend
举几个例子：
复制内容到剪贴板
代码:

: b( ]% l0 j( [8 T; c>> syms x y
! Z( Q/ [+ n4 u, G. F$ p>> f1=x^4+2*x+1;
>> f2=y^6+5*y^3+3;
>> f3=x^5+2*x^3*y^4+x*y^2+4;
, U1 [8 ^& O9 D5 X' [>> poly_coef(f1)
. D- w0 F, {) m: I' P2 l: f
ans =
. q7 K! ^1 d. J* ^' N
                               [1, 0, 0, 2, 1]
>> poly_coef(f1,y)

ans =
" ]$ V  ?9 v0 h2 g, C
                               [ 4       ]
* V1 X4 A' w% M3 k1 Q$ ?, u                               [x   + 2 x + 1]
3 ?5 g  i: X* W2 j9 t% D>> poly_coef(f2)
' C* L" `/ i3 g! d* C: f$ O+ S* e
ans =

                               [ 6    3 ]
( q7 ^5 k2 G3 X  w* n                               [y   + 5 y   + 3]
0 m, O9 r' P! p2 w8 w>> poly_coef(f2,y)
8 ^3 [4 l1 m1 b# i" G+ u& ?9 x, S
ans =
$ U5 ~+ j) {$ k6 Q# b$ J" D
3 Y% I9 V% L# J0 R" ]. N                            [1, 0, 0, 5, 0, 0, 3]
>> poly_coef(f3)

' Y4 v0 I& Q0 ^) M& Xans =
$ `6 L# a8 U( o
+ ~9 |/ C4 f/ b' }/ e: {                        [       4    2 ]
( u5 g# j1 x$ Q6 p. n, {" p+ F                        [1, 0, 2 y , 0, y , 4]
, x  C( n( J# W# K& {. a>> poly_coef(f3,y)
. ?6 F2 D" R9 W7 \' O" I
ans =

' P9 ~' {, ?( Q* ~/ u! [                        [ 3          5 ]
                        [2 x , 0, x, 0, x   + 4]

Demyar

matlab求多项式的系数。