线性代数实报告1.docVIP

线性代数实验报告 学号******, 姓名: *******, 得分: 实验1 求解线性方程组 实验内容: 用MATLAB求解如下线性方程组Ax = b, 其中 A =, b = T. 实验目的: 1.了解matlab软件，学会matlab软件的基本操作； 2.熟悉matlab软件的数值计算功能； 3.练习编写简单的matlab 程序； 4.掌握matlab解线性方程组的原理，学会使用。 实验原理: 1. 对于方程组Ax = b，由克拉默法则解得； 2. 对于方程组Ax = b，得x=a\b，用矩阵除法解得； 3. 利用矩阵的初等变换, 对增广矩阵[A, b]进行行初等变换解得。 实验方案: 1. 在MATLAB命令窗口中输入如下命令: a1=[5;1;0;0;0;0;0;0]; a2=[6;5;1;0;0;0;0;0]; a3=[0;6;5;1;0;0;0;0]; a4=[0;0;6;5;1;0;0;0]; a5=[0;0;0;6;5;1;0;0]; a6=[0;0;0;0;6;5;1;0]; a7=[0;0;0;0;0;6;5;1]; a8=[0;0;0;0;0;0;6;5]; b=[2;1;1;0;9;1;1;0]; A=[a1,a2,a3,a4,a5,a6,a7,a8]; D=det(A);X=[]; for i=1:8 A=[a1,a2,a3,a4,a5,a6,a7,a8]; A(:,i)=b;X=[X,det(A)/D]; i=i+1; end format rat ,X 2. 在MATLAB命令窗口中输入如下命令: A=[5,6,0,0,0,0,0,0; 1,5,6,0,0,0,0,0; 0,1,5,6,0,0,0,0; 0,0,1,5,6,0,0,0; 0,0,0,1,5,6,0,0; 0,0,0,0,1,5,6,0; 0,0,0,0,0,1,5,6; 0,0,0,0,0,0,1,5]; b=[2;1;1;0;9;1;1;0];format rat,X=A\b 3. 在MATLAB命令窗口中输入如下命令: A=[5,6,0,0,0,0,0,0; 1,5,6,0,0,0,0,0; 0,1,5,6,0,0,0,0; 0,0,1,5,6,0,0,0; 0,0,0,1,5,6,0,0; 0,0,0,0,1,5,6,0; 0,0,0,0,0,1,5,6; 0,0,0,0,0,0,1,5]; b=[2;1;1;0;9;1;1;0];B=[A,b]; format rat,C=rref(B); X=C(:,9) 实验结果: 1. X = Columns 1 through 6 7515/166 -5048/135 3163/133 -2335/174 7111/985 -1411/619 1405/1628 -281/1628 Columns 7 through 8 -123/301 424/1505 2. X = 7515/166 -5048/135 3163/133 -2335/174 7111/985 -2300/1009 1405/1628 -281/1628 3. X = 2671/59 -5048/135 1855/78 -2335/174 823/114 -155/68 63/73 -29/168 对实验结果的分析: 上述3种方案所得的结果不一致, 这是因为方案一和方案二都是直接求得，而方案三则是先将B化为行最简矩阵后求得，中间的过程不同，导致答案也不同。 实验2 研究三个平面的位置关系 实验内容: 用MATLAB研究下面的3个平面 p1: x + y + z = 1 p2: ?x + y = 2 p3: 2x + t2z = t 当t取何值时交于一点? 当t取何值时交于一直线? 当t取何值时没有公共的交点? 并在每一种情形下, 用MATLAB在同一个坐标系内绘制出这3个平面的图形(其中, 没有公共的交点的情况, 只要给t取一个适当的值即可)。 实验目的: 1.掌握练习编写简单的mat