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

Lx00591

matlab求多项式的系数:
$ m$ R4 L% Z# C' K- i9 `: \
单变量：

syms x y
4 V5 v+ u- ?) Z# V: _f1=x^4+2*x+1;
5 }' S2 g* g/ l# cf2=y^6+5*y^3+3;
f3=x^5+2*x^3*y^4+x*y^2+4;
) r1 y# k1 a* d5 @0 Uc = sym2poly(f1)
c =
    1     0     0     2     1
% c' l! a1 i6 V1 I  @>> c = sym2poly(f2)
c =
    1     0     0     5     0     0     3
. |& e4 F+ ^( ]7 o* z3 v# [
0 z- c2 }! T, `& l* y& Y. Q2 Z" c多变量：
0 |2 v) k' x+ c+ U! Y/ ~
6 K9 R' J2 A/ M$ ], a' A  d9 nunction coef=poly_coef(f,var)
%提取多项式f中指定变量var的系数，将结果赋给数组coef
5 K" B; M* v2 F! ]& ^1 ~%f可以是含多变量的多项式
%var是多项式中指定的变量，可选，默认是x
%要用到函数poly_degree()来获得f中指定变量var的最高次幂
IF nargin==1
8 g* G# k( D# X% [5 b3 Cvar=sym('x');
( d8 A  U7 `3 G( T" R+ jend
9 l7 P5 H9 }5 ~1 T, G9 \degree=poly_degree(f,var);
temp_f=f;
8 w3 v0 M& v  O2 A+ s. ~coef(degree+1)=subs(temp_f,var,0);
for n=1:degree
temp_f=simple((temp_f-coef(degree+2-n))/var);
1 U: _. ~% g. C# {3 @$ Wcoef(degree+1-n)=subs(temp_f,var,0);
end
) M4 c) V' K4 I5 u% B7 {9 q. r. Nend
) a3 r$ q% }; E举几个例子：
复制内容到剪贴板
/ X0 C- F% w9 ]: A. S2 v/ J代码:
& z( L: O9 R# B, t/ J/ ^
1 j" w/ @- x/ }8 ^- [9 j$ l) n>> syms x y
>> f1=x^4+2*x+1;
$ B% z# j- r3 }) w>> f2=y^6+5*y^3+3;
3 X- U6 f( ~) \>> f3=x^5+2*x^3*y^4+x*y^2+4;
2 }/ {9 B, x2 a6 f2 p8 d4 N>> poly_coef(f1)

$ I# N) \1 ^+ a' W( J' d& Pans =

                               [1, 0, 0, 2, 1]
+ V9 C' U* z1 q  ~! m>> poly_coef(f1,y)

* f7 ?0 X$ m2 Y: Eans =

+ L( B7 k# P/ x1 d% P3 m                               [ 4       ]
                               [x   + 2 x + 1]
# U2 B* s* w) Z0 z1 ?2 Y, i>> poly_coef(f2)

ans =
0 i2 h. Y3 d7 A& d1 O4 U2 X$ |; v
" U" ^" q4 K% m' `                               [ 6    3 ]
6 m8 `7 d" `5 N' B4 a. F& q8 g                               [y   + 5 y   + 3]
>> poly_coef(f2,y)
8 |3 {6 Y  d$ r' e" s; B
ans =

; }9 T7 a  K2 t! }  R                            [1, 0, 0, 5, 0, 0, 3]
) t3 U! W7 _4 S7 K- I2 {>> poly_coef(f3)
0 ~# j2 U# `2 c1 P% W9 N) b
ans =

                        [       4    2 ]
; e5 P) l0 |, j  S: w4 [                        [1, 0, 2 y , 0, y , 4]
  Y# D7 j' }, X  {/ _, y>> poly_coef(f3,y)

& b8 M: w7 s: C5 Kans =
) i1 m, p6 B5 K, \
8 ^+ `) u8 R3 d# i5 F2 K                        [ 3          5 ]
4 X: h$ l5 ~& T- Z( X' H8 e8 N                        [2 x , 0, x, 0, x   + 4]

Demyar

matlab求多项式的系数。