内容简介作者简介 - 科学与工程计算系

P=poly(R)
PR=real(P)
PPR=poly2str(PR,'x')
P =
1.0000 1.1000 0.5500 0.1250
PR =
1.0000 1.1000 0.5500 0.1250
PPR =
x^3 + 1.1 x^2 + 0.55 x + 0.125
2.8.2 多项式运算函数
2
( s + 2)(
s + 4)(
s + 1)
【例 2.8.2-1】求
的“商”及“余”多项式。
3
s + s + 1
p1=conv([1,0,2],conv([1,4],[1,1]));
p2=[1 0 1 1];
[q,r]=deconv(p1,p2);
cq='商多项式为 '; cr='余多项式为 ';
disp([cq,poly2str(q,'s')]),disp([cr,poly2str(r,'s')])
商多项式为 s + 5
余多项式为 5 s^2 + 4 s + 3
【例 2.8.2-2】两种多项式求值指令的差别。
S=pascal(4)
P=poly(S);
PP=poly2str(P,'s')
PA=polyval(P,S)
PM=polyvalm(P,S)
S =
1 1 1 1
1 2 3 4
1 3 6 10
1 4 10 20
PP =
s^4 - 29 s^3 + 72 s^2 - 29 s + 1
PA =
1.0e+004 *
0.0016 0.0016 0.0016 0.0016
0.0016 0.0015 -0.0140 -0.0563
0.0016 -0.0140 -0.2549 -1.2089
0.0016 -0.0563 -1.2089 -4.3779
PM =
1.0e-010 *
0.0016 0.0033 0.0090 0.0205
0.0045 0.0101 0.0286 0.0697
0.0095 0.0210 0.0653 0.1596
0.0163 0.0387 0.1226 0.3019
【例 2.8.2-3】部分分式展开。
a=[1,3,4,2,7,2];
b=[3,2,5,4,6];
[r,s,k]=residue(b,a)
r =
1.1274 + 1.1513i
1.1274 - 1.1513i
-0.0232 - 0.0722i
-0.0232 + 0.0722i
0.7916
s =
6