4065400彭显春matlab2013实用matlab课件6章.pptx

MATLAB 2013 第6章MATLAB数值处理2012.10.12本章要点： 多项式运算 曲线插值 曲线拟合 数值微分 数值积分 极值与过零点 线性代数运算 Matlab与概率论 6.1数值运算6.1.1多项式1、多项式的构造最简单的方法为直接输入向量 P=[1 7 14 8] %以向量的形式显示多项式P = 1 7 14 8poly2sym(P) %以符号表达式的形式显示多项式ans = x^3+7*x^2+14*x+8用多项式的根逆推多项式 roots=[-0.5,-0.3+0.4*i,-0.3-0.4*i]; P=poly (roots) %以向量的形式显示多项式P =1.0000 1.1000 0.5500 0.1250创建方阵A的特征多项式 A= [0 1 0 0; 0 0 1 0; 0 0 0 1;1 0 0 0]; P=poly (A) %以向量的形式显示多项式P = 1.0000 -0.0000 -0.0000 0 -1.00002、多项式的运算函数调用方式说明convC=conv(a, b)计算多项式a和b的乘积deconv[q, r]=deconv(a, b)除法polyP=poly(r)用根构造多项式polyderpolyder(a)对多项式或有理多项式求导PolyintPolyint(a)对多项式a积分polyfitpolyfit(x, y, n)多项式数据拟合polyvalY=polyval(p, s)以数组规则计算s点中多项式值polyvalmY=polyvalm(p, x)以矩阵规则计算X点中多项式值residue[r, p, k]=residue(a, b)部分分式展开式residue[a, b]=residue(r, p, k)部分分式组合rootsroots(a)求多项式的根多项式的求值运算Matlab中，用函数polyval和polyvalm来计算多项式在给定点的值 Y=polyval(p, s)，说明： 若s为一数值，则求多项式在该点的值；若s为向量或矩阵，则对向量或矩阵中的每个元素求其多项式的值。Y=polyvalm(p, x)，说明：用来求矩阵多项式的值，其调用格式与polyval相同，但含义不同，p为多项式，x为方阵，它以方阵为自变量求多项式的值 设A为方阵，P代表多项式x3-5x2+8，那么polyvalm(P,A)的含义是：A*A*A-5*A*A+8*eye(size(A))而polyval(P,A)的含义是：A.*A.*A-5*A.*A+8*ones(size(A))前者以矩阵规则计算多项式的值，后者以数组规则计算多项式的值。 P=[1,2,3];A=[1,2;3,4]; polyval(P,A)%A可为矩阵或向量，采用数组运算（点运算）ans = 6 11 18 27 polyvalm(P,A) % A必须为方阵，并以X为自变量，采用矩阵运算ans = 12 14 21 33多项式的四则运算 Matlab中，采用+，-运算符来实现多项式的加减运算，采用conv函数和deconv函数来实现多项式的乘除运算。 p1=[1 -2 -1 2];p2=[0 1 0 -1]; %输入多项式，当相加减的多项式的阶次不同时，低阶多项式用首零填补 p3=p1+p2p3 = 1 -1 -1 1 p4=p1-p2p4 = 1 -3 -1 3 p5=conv(p1,p2) %用conv函数计算多项式的乘法p5 = 0 1 -2 -2 4 1 -2 p2=[ 1 0 -1]; p5=conv(p1,p2)%多项式的乘或除p1和p2的长度可以不等p5 = 1 -2 -2 4 1 -2 p2=[ 1 0 -1]; %除法运算中多项式的首系数不能为0.重新输入p2 [p6,r]=deconv(p1,p2) %用deconv函数计算多项式的除法p6 = 1 -2%该返回值为p1除以p2的商r = 0 0 0 0%该返回值为 p1除以p2的余子式求多项式的根Matlab中使用roots函数来计算多项式的根，根以列向量的形式表示。反之，已知根，也可通过poly函数获得多项式 p=[2 1 4 5] ; r=roots(p) r = 0.2500 + 1.5612i 0.2500 - 1.5612i -1.0000部分分式展开与组合 用residue函数实现多项式的部分分式展开和组合，其调用方式分别为 [R,P,K] =