MATLAB与线性代数基本运算数学建模.pptxVIP

MATLAB与线性代数基本运算数学建模会计学第1页/共24页一、矩阵的基本输入在MATLAB命令窗口输入： A＝[1，2，3；2，3，4]或 A＝[ 1 2 3 2 3 4 ]第2页/共24页二、产生特殊矩阵的函数ones 创建全1矩阵zeros 创建零矩阵 eye 创建单位矩阵rand（randn） 创建随机矩阵round四舍五入运算length(A) 矩阵的长度size(A) 矩阵的尺寸第3页/共24页三、矩阵的函数输入A=rand(2,3)B=randn(2,3)C=round(10*randn(2,3))D=eye(5)第4页/共24页 矩阵的剪裁从一个矩阵中取出若干行（列）构成新矩阵称为剪裁，“：”是非常重要的剪裁工具。例如，键入：A=[1 2 3;4 5 6;7 8 9]; A(3, :) %A的第三行输出：ans= 7 8 9第5页/共24页 矩阵的剪裁键入：B=A(2:3, :) %A的2，3行输出：B= 4 5 6 7 8 9键入：A(:, 1) %A的第一列输出：ans= 1 4 7第6页/共24页 矩阵的剪裁键入： A=[1 2 3;4 5 6;7 8 9]; C=A(1:2, [1 3]) 输出：C= 1 3 4 6还有A(1:2:3, 3:-1:1), 第7页/共24页矩阵的拼接将几个矩阵接在一起称为拼接，左右拼接行数要相同，上下拼接列数要相同。键入：D=[C, zeros(2,1)]输出：D= 1 3 0 4 6 0 第8页/共24页矩阵的拼接键入：E=[D;eye (2),ones(2,1)]输出：E = 1 3 0 4 6 0 1 0 1 0 1 1第9页/共24页四、矩阵的基本运算 1、矩阵的加、减与数乘：+， -， *2、矩阵的乘法：*3、矩阵的转置：4、方阵的幂运算：^5、方阵的逆：inv6、方阵的行列式：det7、矩阵的秩：rank第10页/共24页五、求线性方程组的唯一解 问题：设A为n阶可逆矩阵，求方程组Ax＝b的解。方法一：x=inv(A)* b 或：x=A^-1* b方法二：U＝rref（[A，b]）U为矩阵[A，b]的行最简形第11页/共24页例：求非齐次线性方程组的唯一解。 第12页/共24页解：在MATLAB命令窗口输入：A=[2,1,2,4;-14,17,-12,7;7,7,6,6;-2,-9,21,-7]b=[5;8;5;10]x=inv(A)*bx=A^-1*bU=rref([A,b])第13页/共24页六、求线性方程组的通解 问题：求方程组Ax＝b的通解方法一：U＝rref（[A,b])方法二：Ax=b的特解：x0=A\b Ax=0的通解：x=null(A,’r’)第14页/共24页例：求非齐次线性方程组的通解。 第15页/共24页解：在MATLAB命令窗口输入：A=[2,4,-1,4,16;-3,-6,2,-6,-23;3,6,-4,6,19;1,2,5,2,19]b=[-2;7;-23;43]U=rref([A,b]) x0=A\bx=null(A,r) 第16页/共24页七、分析向量组的线性相关性把向量以列的形式放入矩阵A中：A＝[a1,a2,a3,…，am][R,s]＝rref（A)R 为矩阵 A 的行最简形s 为矩阵R的基准元素所在列数所构成的行向量第17页/共24页例： 已知向量组，，求出它的最大无关组，并用该最大无关组来线性表示其它向量。，，第18页/共24页解：在MATLAB命令窗口输入：a1=[1;1;0;2;2]; % 输入5个列向量a2=[3;4;0;8;3];a3=[2;3;0;6;1];a4=[9;3;2;1;2];a5=[6;-2;2;-9;2];A=[a1,a2,a3,a4,a5]; % 构造矩阵A[R,s]=rref(A) 第19页/共24页八、求方阵的特征值和特征向量r=eig(A)r为矩阵A的所有特征值所构成的列向量[P,D]=eig(A)D为对角矩阵，对角线上元素为A的所有特征值；P的列向量是A的属于对应特征值的单位特征向量。第20页/共24页九、化二次型为标准形例 用正交变换法将以下二次型化为标准形。第21页/共24页解：在MATLAB命令窗口输入：A=[1,0,0;0,2,2;0,2,2]; % 输入二次型的矩阵[P,D]=eig(A) 第22页/共24页十、小结矩阵的输入= [ ] , ; round rand randn矩阵的基本运算 + - * ‘ ^ inv rank det 求解方程组 rref([A,b]) null(A,’r’)