第三讲Matlab的高级计算3.10.pptVIP

（三）Matlab的高级数值计算;matlab语言把多项式表达成一个行向量， 该向量中的元素是按多项式降幂排列的。 f(x)=anxn+an-1xn-1+……+a0 可用行向量 p=[an an-1 …… a1 +a0]表示 (1) poly —— 产生特征多项式系数向量 特征多项式一定是n+1维的;例:a=[1 2 3;4 5 6;7 8 0]; p=poly(a) p =1.00 -6.00 -72.00 -27.00 这是多项式p(x)=x3-6x2-72x-27的matlab描述方法， 可用： p1=poly2str(p,‘x’)— 函数文件，显示数学多项式的形式 p1 =x^3 - 6 x^2 - 72 x - 27;利用roots求多项式的根 r=roots(p) r = 12.1229 -5.7345 -0.3884;求多项式零点: roots(p);;(3) Polyval 计算系数为p的多项式在标量或向量x处的值 X = pascal(4) X =? ? 1? ? 1? ? 1? ? 1? ? 1? ? 2? ? 3? ? 4? ? 1? ? 3? ? 6? ?10? ? 1? ? 4? ?10? ?20 p = poly(X)p = ? ? 1? ? -29? ? 72? ? -29? ? 1polyval(p,X) ans =? ? 16? ?? ?16? ?? ?16? ?? ?16? ? 16? ?? ?15? ? -140? ? -563? ? 16? ? -140? ?-2549??-12089? ? 16? ? -563??-12089??-43779;(4).conv多项式乘运算;(5).deconv多项式除运算;(6).多项式微分; k = conv(p,q) [k,r] = deconv(p,q);4.数值微积分;被积函数由一个表格定义 在MATLAB中，对由表格形式定义的函数关系的求定积分问题用trapz(X,Y)函数。其中向量X,Y定义函数关系Y=f(X)。 例 用trapz函数计算定积分。 命令如下： X=1:0.01:2.5; Y=exp(-X); %生成函数关系数据向量 trapz(X,Y) ans = 0.28579682416393;例：用数值积分法求 在x=0到x=10之间所围面积，并讨论步长和积分方法对精度的影响 ;MATLAB程序 for dx=[2,1,0.5,0.1] % 设不同步长 x=0:.1:10; y=-x.*x+115; plot(x,y,g),hold % 画出被积曲线 x1=0:dx:10; %根据不同步长取样点 y1=-x1.*x1+115; % 求取样点上的y1 n=length(x1); s=sum(y1(1:n-1))*dx; % 用欧拉法求积分 q=trapz(y1)*dx; % 用梯形法求积分 stairs(x1,y1), plot(x1,y1) % 画出欧拉法及梯形法的积分区域 [dx,s,q],pause(1),hold off %给出不同步长下的结果 end ;运行结果： 步长 欧拉法解 梯形法解 2 910 810 1 865 815 0.5 841.25 816.25 0.1 821.65 816.65;变步长辛普生法 基于变步长辛普生法，MATLAB给出了quad函数来求定积分。该函数的调用格式为： [I,n]=quad(‘fname’,a,b,tol,trace) 其中fname是被积函数名。a和b分别是定积分的下限和上限。 tol用来控制积分精度，缺省时取tol=0.001。trace控制是否展现积分过程，若取非0则展现积分过程，取0则不展现，缺省时取trace=0。 需要事先建立函数文件fname;牛顿－柯特斯法 基于牛顿－柯特斯法，MATLAB给出了quad8函数来求定积分。该函数的调用格式为： [I,n]=qua