实验1 MATLAB初步.pptVIP

函数 eye %生成单位矩阵 Y = eye(n) %生成n×n单位阵 Y = eye(m,n) %生成m×n单位阵 Y = eye(4,3) Y = 1 0 0 0 1 0 0 0 1 0 0 0 Y = eye(3,4) Y = 1 0 0 0 0 1 0 0 0 0 1 0 MATLAB初步 函数 ones %生成全1矩阵 Y = ones(2,3) Y = 1 1 1 1 1 1 函数 rand %生成均匀分布随机矩阵 Y = rand(n) %生成n×n随机矩阵，其元素在 %(0，1）内 MATLAB初步 函数 blkdiag %产生以a,b,c,d,…为对角线元素的矩阵 out = blkdiag(1,2,3,4) MATLAB初步 4. 矩阵的基本数学运算 矩阵的加减法使用“+”、“-”运算符，格式与数字运算完全相同，但要求加减的两矩阵同阶。 (1)矩阵的加和减 a=[1 2 3;2 3 4;3 4 5]; b=[1 1 1;2 2 2;3 3 3]; c=a+b d=a-b c = 2 3 4 4 5 6 6 7 8 d = 0 1 2 0 1 2 0 1 2 MATLAB初步 (2)矩阵的乘法 矩阵的乘法使用运算符“*”，要求若A为i × j阶， 则B必须为j×k阶时，A和B才可以相乘。 “.*” 则是每个元素对应相乘。 e=[b,[5 5 5]] e = 1 1 1 5 2 2 2 5 3 3 3 5 f=a*c f = 28 34 40 40 49 58 52 64 76 f=e*a=? f=a.*b=? MATLAB初步 （4）矩阵的逆运算 (3)．矩阵与常数间运算 常数与矩阵的运算即是同此矩阵的各元素之间 进行运算，如 数加是指每个元素都加上此常数， 数乘即是每个元素都与此数相乘。需要注意的是， 当进行数除时，常数通常只能做除数。 命令 inv inv(a) 例2.4 求矩阵A的逆 a=[2 1 -3 -1;3 1 0 7;-1 2 4 -2;1 0 -1 5]; format rat; ％以有理数的格式输出数据 MATLAB初步 a=[1,3;1,3]; inv(a) Warning: Matrix is singular to working precision. (Type warning off MATLAB:singularMatrix to suppress this warning.) ans = 1/0 1/0 1/0 1/0 a=[1,3;1,3;1,3]; inv(a) ??? Error using == inv Matrix must be square. MATLAB初步 (5)．矩阵的行列式运算 矩阵的行列式的值可由det函数计算得出。 a1=det(a); a2=det(inv(a)); a1*a2 ans = 1.0000 .例1.3.8 求上例中的A矩阵及其逆的行列式之积 MATLAB初步 （6）矩阵的变维 函数：reshape （7）矩阵的变向 旋转： rot90 左右翻转： flipud fliplr 上下翻转: A=1:9; B=reshape(A,3,3) MATLAB初步 注意：矩阵元素赋值的先后顺序 是“列优先”原则。 即先赋值给第一列， 然后是第二列,… (8) 除法运算 Matlab提供了两种除法运算：左除（\）和右除 （/）。一般情况下，x=a\b是方程a*x =b的解， 而x=b/a是方程x*a=b的解。 a=[1 2 3; 4 2 6; 7 4 9] b=[4; 1; 2]; x=a\b 则显示：x= -1.5000 2.0000 0.5000 如果a为非