第5讲数值总结.pptVIP

MATLAB工程应用 MATLAB工程应用 MATLAB工程应用 1. MATLAB的多项式表示 在MATLAB中，对多项式 用其系数的行向量 来表示： 其中，该向量的元素按幂指数降序排列，变量x的幂次已隐含在系数元素的排序中。 注意：如果x的某次幂的系数为零，这个零必须列入系数向量中。例如一个一元三次多项式： 可表示成行向量： MATLAB工程应用 MATLAB工程应用 MATLAB工程应用 2.多项式的加减运算 MATLAB规定，只有长度相同的向量才能相加。因此如果两个多项式向量大小相同，则可利用标准的数组加减法来进行多项式的加减运算。 例：把多项式 与多项式 相加。 a=[1,2,3,4]; b=[1,4,9,16]; d=a+b 结果：d = 2 6 12 20 注意：若两多项式系数向量长度不等，必须在短的向量中补以若干个零元素，才能进行加减。 MATLAB工程应用 3.多项式相乘 多项式相乘就是两个多项式系数向量的卷积，利用函数conv来实现。格式： 　　　　　　 w=conv(a,b) 功能：返回a、b两向量的卷积，也就是a、b代表的两个多项式的乘积。 例： a=[1 2 3 4]; b=[1 4 9 16]; w=conv(a,b) MATLAB工程应用 4.多项式相除 　　函数[q,r]=deconv(u,v)用于对多项式u和v作除法运算。其中q返回多项式u除以v的商式，r返回u除以v的余式。 　　这里，q和r仍是多项式系数向量。 deconv是conv的逆函数，即有 　　　　　　u=conv(v,q)+r 见下页例题。 MATLAB工程应用 例：a=[1 2 3 4]; b=[1 4 9 16]; w=conv(a,b) [q,r]=deconv(w,b) 如令除数a1=a+1,则： a1=a+1; [q1,r1]=deconv(w,a1) 　　可以用商式与除式相乘，再加上余式的方法来检验： w1=conv(q1,a1)+r1 MATLAB工程应用 MATLAB工程应用 MATLAB工程应用 2. 多项式的根 　　求解多项式的根，即a(x)=0的解。MATLAB中求解多项式的根由roots函数命令来完成。 格式： r=roots(p) ； 返回多项式a(x) 的根。 　 注意： MATLAB按惯例规定，多项式是行向量， 根是列向量。 MATLAB工程应用 MATLAB工程应用 MATLAB工程应用 MATLAB工程应用 MATLAB工程应用 MATLAB工程应用 MATLAB工程应用 例：曲线拟合程序范例 t = [0 .3 .8 1.1 1.6 2.3]; y = [0.5 0.82 1.14 1.25 1.35 1.40]; %给出拟合曲线 p=polyfit(t,y,2) %计算拟合曲线函数在更多点上的值 t2 = 0:.1:2.8; y2=polyval(p,t2); %画出原来数据点(t,y)，并利用计算出的数据(t2,y2)给出拟合函数曲线 plot(t,y,o,t2,y2) MATLAB工程应用 MATLAB工程应用 MATLAB工程应用 MATLAB工程应用 例 不同阶次拟合曲线的拟合效果对比。 clear x=0:0.5:20; y=polyval([3,5,1,2],x)+randn(size(x)); p1=polyfit(x,y,1); y1=polyval(p1,x); p2=polyfit(x,y,2); y2=polyval(p2,x); p3=polyfit(x,y,3); y3=polyval(p3,x); plot(x,y,.,x,y1,-.,x,y2,--,x,y3,-) legend(data,linear,square,cubic) MATLAB工程应用 MATLAB工程应用 MATLAB工程应用 MATLAB工程应用 MATLAB工程应用 MATLAB工程应用 MATLAB工程应用 MATLAB工程应用 MATLAB工程应用 MATLAB工程应用 MATLAB工程应用 例：用四种方法对