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

Lx00591

matlab求多项式的系数:

! k: P  h4 x) [单变量：

+ H2 m- z  s& ^1 ~$ usyms x y
f1=x^4+2*x+1;
; W; i% P& l' s; ^/ pf2=y^6+5*y^3+3;
# x) ^& y3 Z1 W4 Cf3=x^5+2*x^3*y^4+x*y^2+4;
/ _  M0 Z# I8 J5 V1 Pc = sym2poly(f1)
( V/ a9 @. R3 E; r% N  Xc =
    1     0     0     2     1
>> c = sym2poly(f2)
c =
    1     0     0     5     0     0     3
- w+ ^' z9 @; |* `
多变量：
4 r& o+ Q8 X- e' T- l
3 v$ D( |4 n1 h8 q! d: W0 Gunction coef=poly_coef(f,var)
%提取多项式f中指定变量var的系数，将结果赋给数组coef
%f可以是含多变量的多项式
%var是多项式中指定的变量，可选，默认是x
%要用到函数poly_degree()来获得f中指定变量var的最高次幂
  e4 N9 R: m$ U; p: Y# x- o2 zIF nargin==1
var=sym('x');
end
degree=poly_degree(f,var);
temp_f=f;
coef(degree+1)=subs(temp_f,var,0);
for n=1:degree
( r# h% \+ \. G+ k$ d' Ztemp_f=simple((temp_f-coef(degree+2-n))/var);
coef(degree+1-n)=subs(temp_f,var,0);
  `0 |& @# P9 `0 r8 S2 o  oend
2 c! W% Q  I! \$ c9 X0 [7 d0 send
& u& e+ b& j8 ~举几个例子：
复制内容到剪贴板
0 e" `6 j4 F  J% y8 P. T1 t代码:

# J# `/ G6 A6 s, J6 }5 R2 ^8 b1 ~>> syms x y
# F( c3 ^/ G& |1 K; K>> 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;
>> poly_coef(f1)

* T$ f! U* [9 ]0 |ans =
  x0 ]2 Y$ d# E9 M
                               [1, 0, 0, 2, 1]
>> poly_coef(f1,y)
( \9 w( I) T5 [7 k
ans =

                               [ 4       ]
7 f% X: d: j! O: z! O% r6 D8 l                               [x   + 2 x + 1]
. b4 d: G5 E, T>> poly_coef(f2)

ans =

                               [ 6    3 ]
                               [y   + 5 y   + 3]
4 C& \6 v. o2 x. R2 O2 B>> poly_coef(f2,y)

5 M/ R$ h& X# D( I0 q2 F4 U. r4 Eans =

                            [1, 0, 0, 5, 0, 0, 3]
>> poly_coef(f3)

ans =

7 x7 g4 u0 E1 `# B5 j2 G2 B3 @                        [       4    2 ]
                        [1, 0, 2 y , 0, y , 4]
>> poly_coef(f3,y)

ans =

                        [ 3          5 ]
, U" s5 e! w, N7 h  ^) e                        [2 x , 0, x, 0, x   + 4]

Demyar

matlab求多项式的系数。