三个关于matlab机算的matlab机算.ppt

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 实验3 向量组的相关性及方程组的通解 * 实验3 向量组的相关性及方程组的通解 * * 在MATLAB命令窗口中输入： la07 运行结果为： 最大线性无关组为：a1 a2 a4 A0 = 1 3 9 1 4 3 0 0 2 2 8 1 2 3 2 a3= -1*a1+ 1*a2+ 0*a4 a5= 3*a1+ -2*a2+ 1*a4 实验3 向量组的相关性及方程组的通解 * 实验3 向量组的相关性及方程组的通解 * % 计算带符号变量的齐次线性方程组的解 clear syms k % 定义符号变量k A=[1-2*k,3,3,3;3,2-k,3,3;3,3,2-k,3;3,3,3,11-k]; % 给系数矩阵赋值 D=det(A); % 算出系数矩阵的行列式D kk=solve(D); % 解方程“D＝0”，得到解kk，即k值 for i=1:4 AA=subs(A,k,kk(i)); % 分别把k值代入系数矩阵A中 fprintf(当k=); disp(kk(i)); % 显示k的取值 fprintf(基础解系为：\n); disp(null(AA)) % 计算齐次线性方程组“Ax=0”的基础解系 end 实验3 向量组的相关性及方程组的通解 * 在MATLAB命令窗口中输入： la08 运算结果为： 当k=7/2 基础解系为： [ 1] [ 2] [ 2] [ -2] 当k=14 基础解系为： [ 1] [ 2] [ 2] [ 5] 当k=-1 基础解系为： [ -1, -1] [ 1, 0] [ 0, 1] [ 0, 0] 当k=-1 基础解系为： [ -1, -1] [ 1, 0] [ 0, 1] [ 0, 0] 实验3 向量组的相关性及方程组的通解 * 实验1 矩阵的基本运算 实验2 行列式与方程组的求解 实验3 向量组的相关性及方程组的通解 实验4 特征向量与二次型 实验5 线性代数的几何概念与Matlab作图 基础篇 * 实验4 特征向量与二次型 实验目的: 掌握Matlab软件对向量组正交化的方法； 掌握Matlab软件求方阵特征值和特征向量的方法； 掌握Matlab软件分析方阵是否可对角化的方法； 掌握Matlab软件化二次型为标准型的方法； 掌握Matlab软件分析对称阵是否正定的方法； 了解MATLAB软件关于矩阵分解的命令。 * 实验4 特征向量与二次型 4.2实验指导: 本实验利用MATLAB软件来求解: 矩阵的特征值和特征向量； 分析矩阵的对角化问题； 给出二次型标准化的方法； 分析二次型的正定性； 了解矩阵的分解命令。 表4.1给出了与本实验相关的MATLAB命令或函数。 * * 实验4 特征向量与二次型 解：在MATLAB命令窗口输入： a1=[1;2;3]; a2=[-1;1;2]; a3=[5;1;0]; A=[a1,a2,a3] P=orth(A) % 将矩阵A的列向量组正交规范化， % P的列构成了空间V的一个标准正交基 % P的列数反应了空间V的维数 运算结果为： P = -0.9266 0.3359 -0.3116 -0.4343 -0.2105 -0.8358 从矩阵P的列数可以看出，原向量组是线性相关的，它生成的空间是二维的。 * 解：在MATLAB的M文件编辑器中，编写la09.m文件，它给出了三种求矩阵特征值的方法： % 矩阵特征值的求解方法 clear A=[2,-2,-20,-19;-2,16,-9,11;-8,4,-6,1;0,-8,-4,-7]; %方法一： syms k