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

Lx00591

matlab求多项式的系数:

单变量：

- t5 H. K. a; Lsyms x y
f1=x^4+2*x+1;
f2=y^6+5*y^3+3;
1 I- l& h9 O, v4 p* vf3=x^5+2*x^3*y^4+x*y^2+4;
c = sym2poly(f1)
" a9 A5 O& K" J. n4 h: Xc =
, ^* P* k; d$ m    1     0     0     2     1
>> c = sym2poly(f2)
1 ~" W1 c  p- j4 rc =
+ ^0 c* t# t. Z2 d    1     0     0     5     0     0     3

2 o. @& D- [) C$ i  H+ \2 C" Q+ V2 |多变量：
5 V4 R7 Z) v. x# M: m
unction coef=poly_coef(f,var)
9 j: P  w; b1 [9 @; W, x8 |5 Y$ w8 k%提取多项式f中指定变量var的系数，将结果赋给数组coef
9 ^( G/ C8 `7 i: K+ b%f可以是含多变量的多项式
%var是多项式中指定的变量，可选，默认是x
%要用到函数poly_degree()来获得f中指定变量var的最高次幂
: t! B4 P6 Q+ f! z# ]$ R8 HIF nargin==1
var=sym('x');
end
& A! B* n( Y. [; f* z& D9 Wdegree=poly_degree(f,var);
5 a  U* ], p+ f- |. etemp_f=f;
coef(degree+1)=subs(temp_f,var,0);
for n=1:degree
temp_f=simple((temp_f-coef(degree+2-n))/var);
coef(degree+1-n)=subs(temp_f,var,0);
, K6 |; v( `: S% F- L$ rend
end
+ w$ D  k0 c- _; B6 m$ G$ d- n1 O举几个例子：
复制内容到剪贴板
代码:
$ K7 O9 l9 }  F  [. ^, }
>> syms x y
) u% @5 ^% b" s6 I$ e; }>> f1=x^4+2*x+1;
6 H  V0 ^5 t- r. T- @! ~" E>> f2=y^6+5*y^3+3;
- x" p- g8 [- w- m- g) k1 o8 I>> f3=x^5+2*x^3*y^4+x*y^2+4;
8 |% c2 j% M0 j8 s# ~7 b" h>> poly_coef(f1)
" O  {  m  K! L% T; x( B0 L
! l9 a1 w/ e3 H. Fans =

, Y5 k9 R3 o7 D6 [                               [1, 0, 0, 2, 1]
>> poly_coef(f1,y)

: O4 D) {9 _4 g8 V- Hans =
) e; P  V3 P% G3 N
                               [ 4       ]
                               [x   + 2 x + 1]
>> poly_coef(f2)

) p4 G) |% k1 P( H" z6 d- bans =
$ o- o  K3 r  p/ R  r6 q9 x
7 r2 f1 X+ }6 S5 v. _                               [ 6    3 ]
$ Z/ l/ T7 I% N5 `  r  t                               [y   + 5 y   + 3]
>> poly_coef(f2,y)

ans =
9 Y# R, V! I3 f4 J% @7 K
, z- v& Q( A4 w) c& t6 a# g                            [1, 0, 0, 5, 0, 0, 3]
>> poly_coef(f3)

ans =
8 G8 d3 G" C; L! H! F% [
                        [       4    2 ]
                        [1, 0, 2 y , 0, y , 4]
$ J3 j6 n& X# u# O& B8 o/ E- H>> poly_coef(f3,y)
; C2 H; ?- M( ?7 V/ J# H
1 P4 D' J4 t$ s0 xans =
  Q- ^2 D1 m5 q: @2 m
                        [ 3          5 ]
                        [2 x , 0, x, 0, x   + 4]

Demyar

matlab求多项式的系数。