Matlab数值计算.pptVIP

天津大学数学系 MATLAB与数值计算 刘东毅 数学与应用数学系 天津大学理学院 MATLAB在数值分析中的应用举例 线性代数方程组的数值解法 Gauss消取法解方程组； 多项式和矩阵的特征系统(eig) 函数的插值 Lagrange多项式插值 曲线拟合与函数的数值逼近 构造Legendre正交多项式 MATLAB在数值分析中的应用举例 数值积分与数值微分 解算子quad, quadl。 计算椭圆积分 常微分方程（组）数值解 解初值问题的解算子 ODE23,ODE45,ODE113 非线性方程和方程组的数值解法 解算子FZERO，FSOLVE 1.线性代数方程组的数值解法 【例A1.1】求解方程组 线性代数方程组的数值解法 (1) 在键盘上输入下列内容 A = [1,2,3; 4,5,6; 7,8,0]%节尾没有分号； b = [37;85;61]; %节尾有分号； x=A\b %节尾没有分号； (2) 每按一次【Enter】键，指令就被马上执行（逐行执行）。由于第二条指令节尾有分号，其结果不被显示出来，其它两条指令的结果被马上显示出来。最后在指令窗中将显示以下结果： 线性代数方程组的数值解法 A = x = 1 2 3 3.0000 4 5 6 5.0000 7 8 0 8.0000 ▲BACK 2.多项式和矩阵的特征系统 2.1多项式 MATLAB约定：用系数行向量 P=[a0,a1,…,an-1,an] 来表示多项式 多项式的四则运算 乘除法：MATLAB提供了卷积和解卷函数 乘法：p = conv(p1,p2)，它表示多项式p为多项式p1与多项式p2的积。 除法（带余除法）： [q, r] = deconv(p1, p2)，它表示多项式p1被p2除的商为多项式q而余项是多项式r，即满足p1 = q*p2+r。 【例A.2.1】求下列多项式的“商”及“余” 2.2 多项式求根及其逆问题 n次多项式具有n个根（实根或成对的共轭复根）。MATLAB提供的roots函数用于求多项式的全部根，其调用格式为： x=roots(P) 其中P为多项式的系数向量，根赋给向量x，即x(1),x(2),…,x(n)分别代表多项式的n个根。反之，P = poly( x ) 函数poly生成以向量x为根的多项式. 【例A.2.2】多项式求根及其逆问题 R=[1,-2, -0.3+0.5*i,-0.3-0.5*i]; P=poly(R) PR=real(P) PPR=poly2str(PR,x) Rr=roots(P) 结果依次为 P = 1.0000 1.6000 -1.0600 -0.8600 -0.6800 PR = 1.0000 1.6000 -1.0600 -0.8600 -0.6800 PPR=x^4+1.6 x^3 - 1.06 x^2 - 0.86 x - 0.68 Rr = -2.0000 1.0000 -0.3000 + 0.5000i -0.3000 - 0.5000i 2.3 矩阵的特征系统 【例A.2.3】计算的矩阵A全部特征值，特征向量与特征多项式。 矩阵的特征系统 解：在指令窗输入： A=[3,-4,3;-4,6,3;3,3,1] [V,D]=eig(A) P = poly(A) 便可得到结果： V = -0.5818 0.6312 0.5130 -0.4534 0.2719 -0.8488 0.6752 0.7264 -0.1280 每一列为A的一个特征向量。 矩阵的特征系统 D = -3.5995 0 0 0 4.7296 0 0 0 8.8699 D的对角线上的元素为对应的特征值，即对应V的每一个列向量。