第六章matlab数值计算.docVIP

第六章MATLAB数值计算 ———————————————————————————————— 作者： ———————————————————————————————— 日期： 第六章 MATLAB 数值计算 6-1 多项式的运算 6-1-1 多项式的生成和表达 多项式的表达 在MATLAB环境下多项式是用向量的形式表达的. 向量最右边的元素表示多项式的0阶，向左数依次表示多项式的第1阶、第2阶、第3阶…。 例如多项式表示为：[5 0 3 2 1]。 多项式的生成 语法： P=ploy(MA) 说明： 假设MA为方阵，那么生成的多项式P为方阵MA的特征多项式。 假设MA为向量，那么向量和多项式满足这样一种关系：，生成的多项式为： 直接输入的方式生成多项式。 例6-1 利用方阵M=[5 6 7;8 9 1;11 12 13]生成一个多项式〔为方阵M的特征多项式〕。 程序设计： clear M=[5 6 7;8 9 1;11 12 13]; P=poly(M); %产生多项式的向量表达式 Px=poly2str(P, x); %生成常见的多项式表示形式 P,Px 运行结果： P = Px = x^3 - 27 x^2 + 90 x + 54 例6-2 利用向量A=[2 3 4 5]生成一个多项式。 程序设计： clear A=[2 3 4 5]; P=poly(A); %产生多项式的向量表达式 Px=poly2str(P, x); %生成常见的多项式表示形式 P,Px 运行结果： P = 1 -14 71 -154 120 Px = x^4 - 14 x^3 + 71 x^2 - 154 x + 120 6-1-2 多项式的乘除 语法： c=conv(a,b) [q,r]=decony(c,a) 说明： a、b和c分别是多项式的向量表示形式。A表示两个多项式的乘积运算，B表示两个多项式的除法运算。 q表示除运算的商，r表示除运算的余数。 例6-3 求多项式和的乘积M(x)。 程序设计： clear a=[1 5 0]; %第一个多项式F(x) b=[2 1]; %第二个多项式G(x) c=conv(a,b); %求两个多项式的乘积 Mx=poly2str(c, x); %用常用的方式表示多项式的积 c,Mx %end 运行结果： c = 2 11 5 0 Mx = 2 x^3 + 11 x^2 + 5 x 例6-4 求多项式和的除运算D(x)。 程序设计： clear c=[1 5 0]; %第一个多项式F(x) a=[2 1]; %第二个多项式G(x) [q,r]=deconv(c,a); %求F(x)/ G(x) Dx=poly2str(q, x); %用常用的方式表示多项式的积 q,r,Dx %end 运行结果： q = r = Dx = 程序说明： 在运行结果中变量q是F(x)除以G(x)的商，而r那么是除不尽的余数。 运行结果变量Dx表示的商没有加上余数。 6-1-3多项式的求导 语法： Dp=polyder(p) 说明： p为向量表示的多项式。 例6-5 求多项式和的一阶导数。 我们容易知道以上两个方程的导数手工验算结果为：F (x) =2x+5和G (x)=2 我们看MATLAB的计算结果。 程序设计： clear f=[1 5 0]; %第一个多项式F(x) g=[2 1]; %第二个多项式G(x) Df=polyder(f); %求F(x)的导数 Dg= polyder(g); %求G (x)的导数 Dfx=poly2str(Df, x); Dgx=poly2str(Dg, x); D