matlab机算讲义.ppt

实验1 矩阵的基本运算 实验2 行列式与方程组的求解 实验3 向量组的相关性及方程组的通解 实验4 特征向量与二次型 实验5 线性代数的几何概念与Matlab作图;实验1 矩阵的基本运算;;;实验1 矩阵的基本运算;在本实验中用到MATLAB的运算符号及命令或函数列举如下： 1、运算符号: 2、命令或函数 表1.2给出了与本实验相关的MATLAB命令或函数。若要进一步了解和学习某个命令或函数的详细功能和用法时，MATLAB提供了一个help命令。 ;7;实验1 矩阵的基本运算;实验1 矩阵的基本运算;实验1 矩阵的基本运算;实验1 矩阵的基本运算;实验1 矩阵的基本运算;实验1 矩阵的基本运算;实验1 矩阵的基本运算;实验1 矩阵的基本运算;实验1 矩阵的基本运算;实验1 矩阵的基本运算;实验1 矩阵的基本运算;实验1 矩阵的基本运算;实验1 矩阵的基本运算;实验1 矩阵的基本运算;实验1 矩阵的基本运算 实验2 行列式与方程组的求解 实验3 向量组的相关性及方程组的通解 实验4 特征向量与二次型 实验5 线性代数的几何概念与Matlab作图;实验2 行列式与方程组的求解;实验2 行列式与方程组的求解;实验2 行列式与方程组的求解;实验2 行列式与方程组的求解;27;实验2 行列式与方程组的求解;实验2 行列式与方程组的求解;实验2 行列式与方程组的求解;实验2 行列式与方程组的求解;实验2 行列式与方程组的求解;;在MATLAB命令窗口中输入： la01 得到以下人机对话结果： 方程个数n＝5 n = 5 系数矩阵A=[6,2,3,4,5;2,-3,7,10,13;3,5,11,-16,21;2,-7,7,7,2;7,3,-5,3,10] A = 6 2 3 4 5 2 -3 7 10 13 3 5 11 -16 21 2 -7 7 7 2 7 3 -5 3 10 常数列向量b=[80;59;90;22;85] ; 显然，当方程组的系数矩阵不是方阵，或系数行列式等于零时，逆矩阵法和克莱姆法都不能实现方程组的求解，而初等行变换的方法适合各种线性方程组的求解，在下一个实验中将继续讨论rref命令的详细应用。关于矩阵的左除运算，有着多种的含义，在实验3和实验5中将逐步讨论它的其它功能。;实验2 行列式与方程组的求??;实验2 行列式与方程组的求解;实验2 行列式与方程组的求解;实验2 行列式与方程组的求解;实验2 行列式与方程组的求解;% 6.伴随矩阵求逆法：;实验2 行列式与方程组的求解;ans = 1.0e+015 * -4.5036 -4.5036 4.5036 9.0072 9.0072 -9.0072 -4.5036 -4.5036 4.5036 显然此结果是不正确的，因为A不可逆。;实验2 行列式与方程组的求解;45;实验2 行列式与方程组的求解;在MATLAB的M文件编辑器中，编写程序la04.m： % 验证行列式按行（列）展开公式 clear A=round(10*randn(5)); % 构造5阶随机数方阵 D=det(A); % 计算矩阵A的行列式 % 矩阵A按第一行元素展开：s=a11*A11+a12*A12+…+a15*A15 s=0; for i=1:5 T=A; T(1,:)=[]; % 删去阵矩第1行 T(:,i)=[]; % 删去矩阵第i列 % 此时，|T| 为矩阵A元素a1i的余子式 s=s+A(1,i)*(-1)^(1+i)*det(T); end e=D-s % 验算D与s是否相等;在MATLAB的命令窗口中输入： la04 计算结果为： e = 0;在MATLAB的M文件编辑器中，编写程序la05.m： % 计算5阶方阵A的第一行