仿真工具教程_04数值计算.pptVIP

第四章 数值计算 4.1 数值微积分 4.1.2数值求和与近似数值积分 4.1.3计算精度可控的数值积分 4.1.4 函数极值的数值求解 4.1.5常微分方程的数值解 4.2矩阵和代数方程 4.2.2矩阵的变换和特征值分解 4.2.3 线性方程的解 3.矩阵逆 4.2.4一般代数方程的解 4.4 多项式运算和卷积 4.4.3两个有限长序列的卷积 指令：inv(A) 对于方程AX=b的求解方法： 左除法：x=A\b 逆矩阵法：x=inv(A)*b 注意： 方法（1）速度快，且精度比方法（2）高得多 对于函数f(x)=0求零点： 步骤(1):作图观察获得近似解 先确定一个零点可能存在的自变量空间，通过 plot指令作图，观察f(x)与横轴的交点坐标，用 zoom指令局部放大，借助ginput指令获得交点 坐标值。 步骤（2）指令求精确解： [x,favl]=fzero(fun,x0) %求一元函数零点最简格式 [x,favl]x=fsolve(fun,x0) %求解多元函数零点 注意：x0为零点初始猜测值，可以通过作图先得到。 例：求f(x)=x3-2x+5的零点 %第一步，绘制一元曲线图 x=-5:0.1:5; f=x.^3-2*x-5; xlabel(x); ylabel(f(x)); plot(x,f,x,0) %x轴划0刻度 %通过观察零点大致在2附近 f=@(x)x.^3-2*x-5 % 计算2附近的零点和该点的函数值 [x,f] = fzero(f,2) x = 2.0946 f = -8.8818e-016 例：求方程f(x,y,z)=2x+3y2+xz在（2,2,2）附 近的零点 y=@(x)(2*x(1)+3*x(2)^2+x(1)*x(3)); %变量不要写成x,y,z x0 = [2 2 2]; % 在同一工作目录下，给出初始值 [x,fval] = fsolve(y,x0) % 计算零点和该点的函数值 x = -0.0040 -0.0650 1.1829 fval = 4.9537e-012 4.4.1多项式的运算函数 1.多项式表达方式的约定 降幂多项式P(x)=a1xn+a2xn-1+…+anx+an+1用 以下系数行向量表示： P=[a1 a2…an an+1],所缺少的系数项为0。 2.多项式运算函数 见表4.4-1(P,177) （1）求商及余多项式 format rat %避免浮点显示 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）验证 qp2=conv(q,p2); pp1=qp2+r; pp1==p1 ans = 1 1 1 1 1 【例 4.4-2】矩阵和特征多项式，特征值和多项式根。 （1）特征多项式 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 + 1.8303e-014 （2）特征值和多项式根。 s=eig(A) r=roots(PA) s = 45.3965 -0.3965 0.0000 r = 45.3965 -0.3965 0.0000 【例 4.4-4】多项式求值指令polyval与polyvalm 的本质差别。 （1）给定多项式和（2×2）数组 clear p=[1,2,3]; poly2str(p,x) X=[1,2;3,4] ans = x^2 + 2 x + 3 X = 1 2 3 4 * MATLAB语言 电子教案 ?通信与电子工程学院 * 4.1 数值微积分 4.2 矩阵和代数方程 4.4 多项式运算和卷积 1.近似数值极限及导数 Matlab数值计算中，没有提供专