- 1、本文档共6页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
数学实验报告
学号, 姓名: 武亚文 , 得分:
实验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取何值
您可能关注的文档
- “汽车”空调常用维修工具.doc
- 自-动搬运分拣机器人比赛方案2011-v2.1.doc
- 逸科住宅设备冲压线“自-动化”方案-.pptx
- “汽车”美容业现状与发展.doc
- 自-动点胶机说明书.doc
- 自-动分拣垃圾桶相关报告材料.doc
- 自-动分拣系统实验报告.doc
- 自-动生产线报告.doc
- 自-动输送分拣系统招标书.doc
- 自-动分拣机毕业论文.docx
- 人音版七年级音乐上册(简谱)第一单元《 彩色的中国》教学设计教案1000字.docx
- (2024年秋季版)七年级道德与法治下册 第一单元 人与人之间 第三课 交往有艺术 第2框 表达与倾听教案 教科版.docx
- 口语交际:商量 教学设计-2023-2024学年语文二年级上册统编版.docx
- 八年级生物上册 第2章 人体生命活动的调节 第3节 基因与人体性状同步教案 沪教版.docx
- 初中主题班会活动策划实施方案内容.docx
- 第2课 我是生活小主人(教案)湘美版 美术三年级上册.docx
- 6《景阳冈》教学设计-2023-2024学年语文五年级下册统编版.docx
- 北师大版六年级数学上册第四单元第5课时这月我当家(1)(同步教学设计).docx
- Module 4 Unit 3 Language practice (2) & Around the world & Module task-教学设计 (1) 2023-2024学年外研版九年级英语上册.docx
- Unit 4 Section B 2a-2e教学设计 2023-2024学年人教版八年级英语下册.docx
文档评论(0)