- 1、本文档共6页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
东南大学几何与代数matlab实验报告1
数学实验报告
学号, 姓名: 武亚文 , 得分:
实验1 求解线性方程组
实验内容: 用MATLAB求解如下线性方程组Ax = b, 其中
A =, b = T.
实验目的: 1.MATLAB软件, 学会MATLAB软件的一些基本操作;
2. 熟悉MATLAB软件的一些数值计算功能,如在矩阵方面的运算;
3. 练习编写简单的MATLAB程序。
实验原理: 1. 对于齐次线性方程组Ax = b, 根据gramer法则, 其解为Xi=Di/D;
2. 当A可逆,即|A|≠0时, 方程组Ax = b的解为X=inv(A)*B;
3. 当……时, 对增广矩阵[A, b]进行初等行变换, 把它化为最简形, 则原矩阵b对应的地方即为A-1*b,取最简形的最后一列赋给X即得线性方程组的解。
实验方案: 1. 在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=[0 4 2 1 1 7 0 2]; %输入矩阵A与b
X=[]; %声明一个名为X的空矩阵
D=det(A); %把A的行列式的值求出来赋给D
for i=1:8 %开始循环语句
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=;
A(:,i)=b; %把A的第i列替换成b
X=[X,det(A)/D]; %把求出的每个值填入空矩阵X对应的列中
i=i+1; %让i递增一控制循环的结束
end %结束for循环
X=X’%输出矩阵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=[0 4 2 1 1 7 0 2];%输入矩阵A和b
X=inv(A)*b%把A逆与b相乘的结果赋给X,即为线性方程组的解
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=[0 4 2 1 1 7 0 2]; %输入矩阵A和b
C=rref([A,b]); %将A的增广矩阵[A,b]化为最简形并赋给C
X=C(:,9) %取C的第九列并赋给X,即为线性方程组的解
实验结果: 1. 2. 3.
X =
-15187/187
7783/115
-4515/107
4238/175
-2116/163
7011/1009
-401/163
727/815
X =
-15187/187
7783/115
-4515/107
4238/175
-2116/163
7011/1009
-401/163
727/815
X =
. -15187/187
7783/115
-2152/51
557/23
-701/54
674/97
-123/50
33/37
5
对实验结果的分析:
在上述3种方案中,前两种的结果完全相同,而与第三种方案的结果
稍有不同,但差距极小,只要是由于不同的计算方法导致运算过程中近似程度不同,导致最终的精度不同。
实验2 研究三个平面的位置关系
内容: 用MATLAB研究下面的3个平面
(1: x + y + z = 1
(2: (x + y = 2
(3: x + t2z = t
当t取何值时交于一点? 当t取何值时交于一直线? 当t取何值
您可能关注的文档
- 三峡蓄水后典型河段分形维数的变化分析.pdf
- 三年级数学上册.千米的认识78页例3和例4666.ppt
- 三江并流区干热河谷成因新探析.pdf
- 三江源区生态移民社会适应问题的调查与思考.pdf
- 三线图解及其水文地质解释.pdf
- 三维地震勘探数据处理中的静校正方法及其应用.pdf
- 三维块体离散元可变形计算方法研究田振农.pdf
- 三维设计高二英语同步课件Module5SectionⅢOtherPartsoftheModule(外研版选修7).ppt
- 三草酸合铁(Ⅲ)酸钾的制备和组成分析探究.doc
- 上海市奉贤区2014届高三1月调研.doc
- 四川省德阳市罗江中学2025届高三考前热身化学试卷含解析.doc
- 山东省枣庄现代实验学校2025届高三下学期第五次调研考试化学试题含解析.doc
- 吉林省长春市十一高中等九校教育联盟2025届高三一诊考试生物试卷含解析.doc
- 2025届江苏省盐城市伍佑中学高考仿真模拟化学试卷含解析.doc
- 2025届广西贺州中学高考冲刺押题(最后一卷)生物试卷含解析.doc
- 安徽省池州市贵池区2025届高三第一次模拟考试生物试卷含解析.doc
- 宁夏银川一中2025届高三(最后冲刺)化学试卷含解析.doc
- 广东省广州市增城区四校联考2025届高考压轴卷化学试卷含解析.doc
- 2025届邯郸市第一中学高考生物必刷试卷含解析.doc
- 2025届安徽省安庆市石化第一中学高考仿真卷化学试卷含解析.doc
文档评论(0)