数学实验与Matlab-共205页.ppt

数学实验与Matlab 实 验 一 矩阵运算与Matlab命令 将表格写成矩阵形式 计 算 输入下面Matlab指令 A=[4 2 3;1 3 2;1 3 3;3 2 2], B=[35 20 60 45;10 15 50 40;20 12 45 20] C=A*B 请自行计算观看结果 Matlab基本指令 向量的创建和运算 1. 直接输入向量 》x1=[1 2 4],x2=[1,2,1],x3=x1’ 运行结果 x1 = 1 2 4 x2 = 1 2 1 x3 = 1 2 4 2.冒号创建向量 x1=3.4:6.7, x2=3.4:2:6.7 x3=2.6:-0.8:0 运算结果 x1 = 3.4000 4.4000 5.4000 6.4000 x2 = 3.4000 5.4000 x3 = 2.6000 1.8000 1.0000 0.2000 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窗口创建向量后并运行后，向量就存在于工作空间，可以被调用。 向量的运算 设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] 数的乘幂： 如 a^2 元素群运算(四则运算） 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] 观察运行结果 创建矩阵（符号矩阵的创建） 用指令“syms”说明符号变量。 syms a11 a12 a13 a14 a21 a22 a23 a24 a31 a32 a33 a34 b11 b12 b13 b14 b21 b22 b23 b24 b31 b32 b33 b34 A1=[a11 a12 a13 a14 ;a21 a22 a23 a24; a31 a32 a33 a34], B1=[b11 b12 b13 b14 ;b21 b22 b23 b24; b31 b32 b33 b34] 运行 矩阵的运算(矩阵的加减、数乘、乘积) C=A1+B1 D=A1-B1 syms c, cA=c*A1 A2=A1(:,1:3), B1 G=A2*B1 矩阵的运算(矩阵的加减、数乘、乘积) 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) 矩阵的运算(左除和右除) 求矩阵