線性代数实验报告3.docVIP

数学实验报告 学号: , 得分: 实验2 研究三个平面的位置关系 用MATLAB研究下面的3个平面 p1: x + y + z = 1 p2: ?x + y = 2 p3: x + t2z = t 当t取何值时交于一点? 当t取何值时交于一直线? 当t取何值时没有公共的交点? 并在每一种情形下, 用MATLAB在同一个坐标系内绘制出这3个平面的图形公共的交点. 实验目的：1. 练习编写简单的MATLAB程序. 2. 掌握用MATLAB软件绘制简单图形的方法. 实验原理: 联立这3个平面x + y + z = 1 ?x + y = 2 2x + t2z = t , 令A = (1,-1,2)T, b =(1,2,t)T, B = [A, b], 则原问题转化为线性方程组 实验方案: 1. 在MATLAB命令窗口中输入如下命令: y=-20:1:20;z=y; [Y,Z]=meshgrid(y,z); X1=ones(size(Y))-Y-Z; X2=-2*ones(size(Y))+Y; X3=(-ones(size(Y))-Z)/2; surf(X1,Y,Z),hold on, mesh(X2,Y,Z),mesh(X3,Y,Z) 2. 在MATLAB命令窗口中输入如下命令: y=-20:1:20;z=y; [Y,Z]=meshgrid(y,z); X1=ones(size(Y))-Y-Z; X2=-2*ones(size(Y))+Y; X3=(ones(size(Y))-Z)/2; surf(X1,Y,Z),hold on, mesh(X2,Y,Z),mesh(X3,Y,Z) 3. y=-20:1:20;z=y; [Y,Z]=meshgrid(y,z); X1=ones(size(Y))-Y-Z; X2=-2*ones(size(Y))+Y; X3=(-3*ones(size(Y))-9*Z)/2; surf(X1,Y,Z),hold on, mesh(X2,Y,Z),mesh(X3,Y,Z) 实验结果: 1、 2、 3、 对实验结果的分析: 实验二： “eigshow”是Matlab中平面线性变换的演示函数。如果是一个矩阵，键入eigshow（），就可分别显示不同的单位向量列及经变换后的向量。用鼠标拖动旋转，可以使产生一个单位圆，并显示的轨迹。分别对矩阵，考察单位向量变化时，变换后所得向量的轨迹，回答下列问题，并给出代数解释：（1）问：和会不会在同一直线上？（2）如果和在同一直线上，它们的长度之比是多少？（3）用其它的矩阵试一下，你还发现有什么规律？ A A=[1,2;2,1];eigshow(A)（1）x和y在同一条直线；（2）此时x和y的长度比为1或1/3 B B=[3,1;2,3];eigshow(B) (1)x和y在同一条直线；（2）x和y的长度比为 C C=[3,-1;2,3];eigshow(C) 此时x和y不在同一直线上。 （3）用其它的矩阵试一下，你还发现有什么规律？ 设A= = +=1 =a+b =c+d 化简得: 会在一直线上 不会在一直线上  附录 实验报告参考模板 1