15-16-2几何与代数数学实验.docVIP

15-16-2《几何与代数》 数学实验报告 学号： 座机电话号码 姓名： 鲁罗林 得分： . 实验一：利用MATLAB用三种不同的方法求解线性方程组。其中，常数项列向量的分量是你的学号，系数矩阵为 方法一：利用Cramer法则求解； Trial a [1,1,0,0,0,0,0,0];b [1,2,1,0,0,0,0,0];c [0,1,3,1,0,0,0,0];d [0,0,1,4,1,0,0,0];e [0,0,0,1,5,1,0,0];f [0,0,0,0,1,6,1,0];g [0,0,0,0,0,1,7,1];h [0,0,0,0,0,0,1,8];i [7,1,1,1,5,3,3,4]; Trial D det [a,b,c,d,e,f,g,h] Trial D1 det [i,b,c,d,e,f,g,h] ,D2 det [a,i,c,d,e,f,g,h] ,D3 det [a,b,i,d,e,f,g,h] ,D4 det [a,b,c,i,e,f,g,h] ,D5 det [a,b,c,d,i,f,g,h] ,D6 det [a,b,c,d,e,i,g,h] ,D7 det [a,b,c,d,e,f,i,h] ,D8 det [a,b,c,d,e,f,g,i] , Trial X1 D1/D;X2 D2/D;X3 D3/D;X4 D4/D;X5 D5/D;X6 D6/D;X7 D7/D;X8 D8/D; Trial X1,X2,X3,X4,X5,X6,X7,X8 x1 17.0227 x2 -10.0227 x3 4.0227 x4 -1.0453 x5 1.1586 x6 0.2524 x7 0.3269 x8 0.4591 方法二：作为矩阵方程求解； Trial format rat Trial A [1,1,0,0,0,0,0,0;1,2,1,0,0,0,0,0;0,1,3,1,0,0,0,0;0,0,1,4,1,0,0,0;0,0,0,1,5,1,0,0;0,0,0,0,1,6,1,0;0,0,0,0,0,1,7,1;0,0,0,0,0,0,1,8] Trial B [7;1;1;1;5;3;3;4] Trial X inv A *B X 17.0227 -10.0227 4.0227 -1.0453 1.1586 0.2524 0.3269 0.4591 方法三：利用Gauss消元法求解。 A [1,1,0,0,0,0,0,0,7;1,2,1,0,0,0,0,0,1;0,1,3,1,0,0,0,0,1;0,0,1,4,1,0,0,0,1;0,0,0,1,5,1,0,0,5;0,0,0,0,1,6,1,0,3;0,0,0,0,0,1,7,1,3;0,0,0,0,0,0,1,8,4] Trial rref A ans 1 0 0 0 0 0 0 0 5260/309 0 1 0 0 0 0 0 0 -3097/309 0 0 1 0 0 0 0 0 1243/309 0 0 0 1 0 0 0 0 -323/309 0 0 0 0 1 0 0 0 95/82 0 0 0 0 0 1 0 0 26/103 0 0 0 0 0 0 1 0 17/52 0 0 0 0 0 0 0 1 118/257 x1 17.0227 x2 -10.0227 x3 4.0227 x4 -1.0453 x5 1.1586 x6 0.2524 x7 0.3269 x8 0.4591 实验二：“eigshow”是Matlab中平面线性变换的演示函数。对于矩阵，键入eigshow（），分别显示不同的单位向量及经变换后的向量。用鼠标拖动旋转，可以使产生一个单位圆，并显示所产生的轨迹。 分别对矩阵，考察单位向量变化时，变换后所得向量的轨迹，回答下列问题，并用代数方法解释。（1）问：和会不会在同一直线上？如果和在同一直线上，它们的长度之比是多少？（2）对哪些矩阵，和的轨迹有公共交点？（3）对哪些矩阵，和的转向相同，哪些相反？（4）你还发现什么有什么规律？（5）你能用代数知识解释这些现象吗？ 答： A,B,C,D中x,y会在一条直线上，长度之比是矩阵的行列式，M中x,y不会再同一条直线上。 A,B,D,M中转向相同，C中相反。 行列式为负的矩阵中转向相反。 不能！