资源描述:
《东南大学几代实验报告.docx》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、14-15-2《几何与代数》数学实验报告学号:__________姓名:_________得分:___________实验1求解线性方程组实验内容:用MATLAB求解如下线性方程组Ax=b,其中:b=[02114102]T.实验目的:1.了解MATLAB软件,学会MATLAB软件的一些基本操作.2.熟悉MATLAB软件的一些数值计算功能.3.练习编写简单的MATLAB程序.实验原理:1.对于齐次线性方程组Ax=b,根据gramer法则,其解为Xi=Di/D;2.当A可逆,即
2、A
3、≠0时,方程组Ax=b的解为X=inv(A)*B;3.当A可逆,即
4、A
5、≠0时,对增广矩阵[
6、A,b]进行初等行变换,把它化为最简形,则原矩阵b对应的地方即为A-1*b,取最简形的最后一列赋给X即得线性方程组的解。实验方案一:利用Cramer法则求解:1.在MATLAB命令窗口中输入如下命令:>>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=[0;2;1;1;4;1;0;2];%输入矩阵A
7、与i>>D=det([a,b,c,d,e,f,g,h]),%把行列式的值求出来赋给DD=10241>>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]),%把A的第i列替换成iD1=-26704D2=26704D
8、3=-6222D4=2203D5=7651D6=506D7=-446D8=2616>>X1=D1/D,X2=D2/D,X3=D3/D,X4=D4/D,X5=D5/D,X6=D6/D,X7=D7/D,X8=D8/D,实验方案二:利用矩阵方程求解:2.在MATLAB命令窗口中输入如下命令:>>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=[
9、0;0;0;0;0;0;1;8];>>A=[a,b,c,d,e,f,g,h],B=[0;2;1;1;4;1;0;2],%输入矩阵A和BA=1100000012100000013100000014100000015100000016100000017100000018B=B=02114102>>X=inv(A)*B,%把A逆与b相乘的结果赋给X,即为线性方程组的解实验方案三:利用Gauss消元法求解:3.在MATLAB命令窗口中输入如下命令:>>formatrat,%(采用有理数数据类型)>>A=[1,1,0,0,0,0,0,0,0;1,2,1,0,0,0,0,0,2;0
10、,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,4;0,0,0,0,1,6,1,0,1;0,0,0,0,0,1,7,1,0;0,0,0,0,0,0,1,8,2]A=110000000121000002013100001001410001000151004000016101000001710000000182>>rref(A)%将A的增广矩阵[A,b]化为简化阶梯形矩阵实验结果1:X1=-2.6076X2=2.6076X3=-0.6076X4=0.2151X5=0.7471X6=0.0494X7=-0.0436X8=0
11、.2554>>实验结果2:X=-2.60762.6076-0.60760.21510.74710.0494-0.04360.2554>>实验结果3:ans=10000000-897/34401000000897/34400100000-209/3440001000037/17200001000257/3440000010021/42500000010-26/5970000000147/184>>所以,该线性方程组的解分别为:X1=-897/344;X2=897/344;X3=-209/344;X4=37/172;X5=257/344;X6