matlab第四章数值运算基础.pptVIP

第4章 数值运算基础;第1节　多项式polynomial 　MATLAB用行向量表示多项式。将多项式的系数按降幂次序存放在行向量中。如：;一、创建多项式 1、系数矢量直接输入法 　在命令窗口直接输入多式的系数向量;2、方阵特征多项式输入法 p=poly(A)　 　若A为n×n的矩阵，则返回值P将是一个含有n+1个元素的行向量，也就是该矩阵特征多项式的系数 ;3、根矢量创建 　 p=poly(A) A为待求多项式的根矢量，则返回值将是对应多项式的系数行矢量，该多项式的根为矢量A。此时p=poly(A)与A=roots(p)互逆。系统定义P0=1。;【例4-3】根据根矢量[-0.5 -0.3+0.4i -0.3-0.4i]，创建多项式 r=[-0.5 -0.3+0.4i -0.3-0.4i]; p=poly(r) pr=real(p) ppr=poly2sym(pr);【例4-4】求多项式3x^2+2x+1在5、7和9处的值。 p = [3 2 1]; polyval(p,[5 7 9]);2、求多项式的根 　格式：C=roots(p) p为多项式的系数矢量，C为函数返回多项式的根矢量 如果C为复数，则必成对出现。;3、多项式的乘除运算 多项式的乘法conv格式： c=conv(a,b) 多项式的乘法运算,也是矢量的卷积运算 向量a长度为m，向量b长度为n，a和b的卷积定义为：;;4、多项式微积分 polyder(p) 返回多项式系数向量p 的导数;5、多项式的部分分式展开 MATLAB提供了residue命令来执行部分分式展开或多项式系数之间的转换。该命令通常用于信号与控制领域中。格式如下： [r,p,k]=residue(b,a) 该命令是求多项式之比b(s)/a(s)的部分分式展开，返回留数r、极点p和直项向量k。a和b分别是分母和分子多项式的系数向量;; b(x) 5x^3+3x^2-2x+7 【例3-11】两多项式的比为—— = ————————, a(x) -4x^3+8x+3 求部分分式展开。 a = [-4 0 8 3]; b = [ 5 3 -2 7]; [r, p, k] = residue(b,a) [b1,a1] = residue(r,p,k) %分母最高次项归1 [r2,p2,k2] = residue([1 1],[1 -2 1]) %出现重根;6、多项式拟合 对于给定的一组数据{(xi ,yi)，i=1,2,…,n},如果要采用多项式模型对数据组进行描述，形成如多项式y(x)=f(x,p)=p1 xn+ p2 xn-1+ p3 xn-2+…+ pn+1的形式，求取参数p使得量值χ2(p)的值最小的过程，称为对数据组进行多项式拟合，其中;p=polyfit(x,y,n)　其中，x,y分别表示横、纵坐标向量；n是给定的拟合??项式的最高阶数，返回一个多项式系数向量p。 如n=3,若p=[1 0.5 1 2],则 　 y=1*x^3+0.5*x^2+1*x^1+2*x^0;【例3-13】 求误差函数的6阶拟合多项式。 x = (0: 0.1: 2.5); % 生成0至2.5间隔为0.1的自变量 y = erf(x); % 计算误差函数 p = polyfit(x,y,6) % 求6阶拟合多项式 x = (0: 0.1: 5); % 生成0至5间隔为0.1的自变量 y1 = erf(x); % 计算误差函数 f = polyval(p,x); % 计算拟合函数的值 plot(x,y1,o,x,f,-) % 绘图函数 [p0,s0,mu0] = polyfit(x,y,6) %x=(x-mean(x))/std(x) [p1,s1,mu1] = polyfit((x-mu(1))/mu(2),y,6);第2节　线性代数 　　给定两个矩阵A和B,求X的解,使得: 　　　AX=B 　　　XA=B; 根据矩阵A的结构(m, n），可以将方程分为以下3类： m=n 方阵系统，可偿试求精确解 mn 超定系统，可偿试求最小二乘解 mn 欠定系统，可偿试求少m个非零解; ;;二、超定系统 　　方程组Ax=b，A为m×n矩阵,如果A列满秩,且mn,则方程没有精确解,此时方程组为超定方程组。一般采用最小二乘法。;一组实验数据，时间t和测