matlab 多项式表达,MATLAB 多项式的表达方式及其操作

3.8 多项式的表达方式及其操作

3.8.1.2 多项式行向量的创建方法

【 * 例 3.8.1 .2-1 】求 3 阶方阵 A 的特征多项式。

A=[11 12 13;14 15 16;17 18 19];

PA=poly(A) %A 的特征多项式

PPA=poly2str(PA,'s') % 以较习惯的方式显示多项式

PA =

1.0000 -45.0000 -18.0000 -0.0000

PPA =

s^3 - 45 s^2 - 18 s - 2.8387e-015

【 * 例 3.8.1 .2-2 】由给定根向量求多项式系数向量。

R=[-0.5,-0.3+0.4*i,-0.3-0.4*i]; % 根向量

P=poly(R) %R 的特征多项式

PR=real(P) % 求 PR 的实部

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

3.8.2 多项式运算函数

【 * 例 3.8.2 -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

【 * 例 3.8.2 -2 】两种多项式求值指令的差别。

S=pascal(4) % 生成一个 4 阶方阵

P=poly(S);PP=poly2str(P,'s')

PA=polyval(P,S) % 独立变量取数组 S 元素时的多项式值

PM=polyvalm(P,S) % 独立变量取矩阵 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-011 *

-0.0077 0.0053 -0.0096 0.0430

-0.0068 0.0481 -0.0110 0.1222

0.0075 0.1400 -0.0095 0.2608

0.0430 0.2920 -0.0007 0.4737

【 * 例 3.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 =

-1.7680 + 1.2673i

-1.7680 - 1.2673i

0.4176 + 1.1130i

0.4176 - 1.1130i

-0.2991

k =

[]

[1]