Matlab-杜睿.pptVIP

数学实验与Matlab 目 录 实验一：矩阵运算与Matlab命令; 实验二：函数可视化与Matlab作图; 实验三：插值和拟和; 实验四：微分、积分和常微分方程; 实验五：最优化方法; 如何撰写数学建模论文.  实 验 一 矩阵运算与Matlab命令 Matlab基本指令 向量的创建和运算 1. 直接输入向量 》x1=[1 2 4], x2=[1,2,1], x3=x1’ 3.生成线性等分向量 指令x=linspace(a,b,n) 在[a,b]区间产生 n 个等分点(包括端点) 例如：x=linspace(0,1,5) 结果 x = 0 0.2500 0.5000 0.7500 1.0000 工作空间 在Matlab窗口创建向量后并运行后，向量就存在于工作空间（Workspace），可以被调用。 向量的运算 设三维向量x=[x1 x2 x3]; y=[y1 y2 y3]; ，a, b为标量。 向量的数乘：a*x=[a*x1 a*x2 a*x3] 向量的平移： x+b=[x1+b x2+b x3+b] 向量和： x+y=[x1+y1 x2+y2 x3+y3] 向量差： x-y=[x1-y1 x2-y2 x3-y3] 向量对应元素的运算 x.*y=[x1*y1 x2*y2 x3*y3] (乘积) x./y=[x1/y1 x2/y2 x3/y3] (右除，右边的y做分母) x.\y=[y1/x1 y2/x2 y3/x3] (左除，左边的x做分母) x.^5=[x1^5 x2^5 x3^5] (乘幂) 2.^x=[2^x1 2^x2 2^x3] x.^y=[x1^y1 x2^y2 x3^y3] 函数计算 Matlab有许多内部函数，可直接作用于向量产生一个同维的函数向量。 如：x=linspace(0,4*pi,100);（产生100维向量x） y=sin(x); (y也自动为100维向量) y1=sin(x).^2; y2=exp(-x).*sin(x); 观察结果 创建矩阵（数值矩阵的创建） 直接输入法创建简单矩阵。 A=[1 2 3 4; 5 6 7 8; 9 10 11 12] B=[-1.3,sqrt(3);(1+2)*4/5,sin(5);exp(2),6] 观察运行结果 A = 1 2 3 4 5 6 7 8 9 10 11 12 B = -1.3000 1.7321 2.4000 -0.9589 7.3891 6.0000 矩阵的运算(矩阵的加减、数乘、乘积等) A, A_trans=A H=[1 2 3 ; 2 1 0 ; 1 2 3 ], K=[1 2 3 ; 2 1 0 ; 2 3 1] h_det=det(H), k_det=det(K), H_inv=inv(H), K_inv=K^-1 矩阵的运算(左除和右除) 左除“ \ ”： 求矩阵方程AX=B的解；（ A 、B的行要保持一致） 解为 X=A\B； 当A为方阵且可逆时有X=A\B=inv(A)*B； 右除“ / ”： 求矩阵方程XA=B的解 （A 、B的列要保持一致） 解为 X=B/A ， 当A为方阵且可逆时有X=B/A=B*inv(A) 矩阵的运算(左除和右除) 例1：求矩阵方程： 设A、B满足关系式：AB＝2B+A,求B。 其中A=[3 0 1; 1 1 0; 0 1 4]。 解：有(A-2I)B＝A 程序 ： A=[3 0 1; 1 1 0;0 1 4]; B=inv(A-2*eye(3))*A, B=(A-2*eye(3))\A 观察结果： 分块矩阵（矩阵的标识） 例2 取出A的1、3行和1、3列的交叉处元素构成新矩阵A1。 解：程序 A=[1 0 1 1 2;0 1 -1 2 3; 3 0 5 1 0;2 3 1 2 1], vr=[1, 3]; vc=[1, 3]; A1=A(vr, vc) 观察结果 分块矩阵（矩阵的标识） 例3 将A分为四块，并把它们赋值到矩阵B中，观察运行后的结果。 解：程序 A11=A(1:2,1:2),A12=A(1:2,3:5), A21=A(3:4,1:2),A22=A(3:4,3:5) B=[A11 A12; A21 A22]