用Matlab进行现代科学运算答题.pptVIP

第五章 用Matlab进行现代科学计算;5.1 解析解与数值解;2、解析解存在但不实用 如求解n元一次代数方程组的问题，从理论上讲，总可以把多元一次方程组简化成解析可解的形式。 然而当n较大时，需要的基本运算次数是非常惊人的，对于求解实际问题所需要的计算成本是根本无法接受的，只能采用数值解法。;5.2 数值线性代数问题及求解;5.2.2 矩阵的特征参数运算;矩阵的范数 norm(A,p);矩阵的特征多项式： c=poly(A);矩阵的特征值与特征向量; 矩阵求逆;多项式及多项式矩阵求值;面向矩阵各个元素的函数;对矩阵进行数值分析的函数;5.2.3 矩阵的相似变换与分解 ;2、正交分解 也称为“QR”分解，将一个矩阵表示成一个正交矩阵和一个上三角矩阵的乘积，A=Q*R 用qr()函数可得到分解后的两个矩阵，其调用形式为 [Q,R]=qr(A);3、奇异值分解 [U,S,V]=svd(A) 其中U和V是正交矩阵，而S是对角矩阵。即 A=U*S*V’ 奇异值分解是矩阵分析的有力工具。Matlab的一些函数也是基于此实现的，包括求伪逆函数pinv(A)，求秩函数rank(A)，求欧式范数norm(A,2)，求条件数cond(A);5.3 数值微积分问题 ; y=magic(6) y = 35 1 6 26 19 24 3 32 7 21 23 25 31 9 2 22 27 20 8 28 33 17 10 15 30 5 34 12 14 16 4 36 29 13 18 11 diff(y) ans = -32 31 1 -5 4 1 28 -23 -5 1 4 -5 -23 19 31 -5 -17 -5 22 -23 1 -5 4 1 -26 31 -5 1 4 -5;5.3.2 数值积分;[y,n]=sum(F,a,b); 等宽矩形法求定积分 [y,n]=trapz(F,a,b); 梯形法求定积分 [y,n]=quad(F,a,b,tol); 自适应辛普森积分 [y,n]=quadl(F,a,b,tol); 自适应Lobatto积分;例子：求无穷定积分;例子：;5.4 常微分方程的数值解法 ;求解常微分方程组的数值方法是多种多样的，如常用的Euler法、Runge-Kutta法、Adams线性多步法、Gear法等。;MATLAB求解常微分方程的函数如下：;ode23()和ode45()最常用，采用自适应变步长求解方法 [t,x]=ode23(方程函数名,[t0,tf],x0,选项) [t,x]=ode45(方程函数名,[t0,tf],x0,选项) 方程函数名为描述系统状态方程的M函数的名称，该函数名应该用引号括起来。 t0和tf分别为起始和终止计算时间，x0为系统的初始状态变量的值。 返回值为求解的时间变量和相应的状态列向量构成的矩阵转置;方程函数名的编写格式是固定的： function xdot=方程函数名(t,x) 其中t为时间变量，x为方程的状态变量，xdot代表状态变量的导数。 注意，及时微分方程是非时变的，也应该在函数输入变量列表中写上t占位。;在lorenzeq.m函数中：function xdot=lorenzeq(t,x)xdot=zeros(3,1);xdot=[-8/3*x(1)+x(2)*x(3);-10*x(2)+10*x(3); -x(1)*x(2)+28*x(2)-x(3)];谢 谢！