第3章MATLAB在高等数学中的应用.pptVIP

格式二：[Q,R,E]=qr(A) 功能：产生一个置换矩阵E，一个上三角矩阵R(其对角线元素降序排列)和一个归一化矩阵Q，满足A*E=Q*R； 格式二：X=diag(A) 功能：当A为n元向量时，等同于k=0时的X=diag(A,k)，即产生A的元素处于主对角线的对角方阵。当A为矩阵时，X=diag(A)相当于k=0。 3．多项式相乘运算 格式：w=conv(u,v) 功能：返回u和v两向量的卷积，也就是u和v代表的两多项式的乘积。 4．多项式相除 格式：[q ,r]=deconv(u , v) 功能：给出商多项式q和余数多项式r ，u为被除多项式 2．多项式的加减运算 格式：A=B+C 1．多项式求导函数polyder 格式一：k=polyder(p) 功能：返回多项式p的一阶导数。 格式二：k=polyder(u,v) 功能：返回多项式u与v乘积的导数。 格式三：[q,d]=polyer(u,v) 功能：返回多项式商u/v的导数, 返回的格式为：q为分子, d为分母。 3.2.2 多项式求导、求根和求值 w=[1 6 20 5] w = 1 6 20 5 polyder(w) ans = 3 12 20 a=[1 2 3]; b=[1 4 9]; k=polyder(a,b) k = 4 18 40 30 [q,d]=polyder(a,b) q = 2 12 6 d = 1 8 34 72 81 2．多项式的根 求解多项式的根，即p(x)=0的解。 格式：r=roots(p) 功能：返回多项式p(x)的根。注意，MATLAB按惯例规定，多项式是行向量，根是列向量。 3．多项式求值函数polyval( ) 利用函数polyval可以求得多项式在某一点的值。 格式：y=polyval(p,x) 功能：返回多项式p在x处的值。其中x可以是复数，也可以是数组。 当多项式的变量是矩阵时，构成的矩阵多项式可以利用polyvalm函数求值。 格式：Y=polyvalm(p,X) 功能：返回矩阵多项式p在X处的值。 4．部分分式展开函数residue( ) 格式一：[r,p,k]=residue(b,a) 功能：把b(s)/a(s)展开成： 格式二：[b, a]=residue(r, p, k) 功能：格式一的逆作用。 其中r代表余数数组，p代表极点数组，ks代表部分分式展开的常数项。当分母多项式的阶次数高于分子多项式的阶次数时ks=0。 1．多项式拟合函数polyfit( ) 格式：p=polyfit(x,y,n) 功能：利用已知的数据向量x和y所确定的数据点，采用最小二乘法构造出n阶多项式去逼近已知的离散数据，实现多项式曲线的拟合。其中p是求出的多项式系数，n阶多项式应该有n+1个系数，故p的长度为n+1。 3.3.3 多项式拟合与多项式插值 2．多项式插值 插值和拟合的不同点在于：①插值函数通常是分段的，人们关心的不是函数的表达式，而是插值出的数据点；②插值函数应通过给定的数据点。 (1)一维插值函数interpl( ) 格式：yi＝interpl (x, y, xi, method) 功能：为给定的数据对(x,y)以及x坐标上的插值范围向量xi，用指定所使用的插值方法method实现插值。yi是插值后的对应数据点集的y坐标。插值的方法method有以下6种可供选择： nearest(最邻近插值法)、linear(线性插值)、 spline(三次样条插值)、cubic(立方插值)、pchip(三次Hermite插值，该方法与cubic类似) 。 例 ： 一维插值函数示例。 001 %INTERP_EX1 一维插值计算示例 002 % 准备数据 003 x = 0:10; 004 y = cos(x); 005 % 插值点 006 xi = 0:0.2:10; 007 % 进行插值运算 013 title(一维插值计算示例) 008 yin = interp1(x,y,xi,nearest); 009 yic = interp1(x,y,xi,cubic); 010 % 绘制结果 011 plot(x,y,o,xi,yin,*,xi,yic) 012 legend(origin,nearest,cubic) 结