- 1、本文档共6页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
东南大学几何与代数matlab实验报告(大一专用).
数学实验报告
学号: , 姓名: 吴雪松 , 得分:
实验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 9 0 0 9 3 2 0]; %输入矩阵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=[0 9 0 0 9 3 2 0];
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 9 0 0 9 3 2 0];%输入矩阵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 9 0 0 9 3 2 0]; %输入矩阵A和b
C=rref([A,b]); %将A的增广矩阵[A,b]化为最简形并赋给C
X=C(:,9) %取C的第九列并赋给X,即为线性方程组的解
实验结果: 1. 2. 3.
X =
34.9644
-29.1370
19.9534
-11.7717
6.4842
-1.9415
1.0372
-0.2074
X =
34.9644
-29.1370
19.9534
-11.7717
6.4842
-1.9415
1.0372
-0.2074
X =
34.9643
-29.1370
19.9535
-11.7717
6.4842
-1.9415
1.0372
-0.2075
对实验结果的分析:
在上述3种方案中,前两种的结果完全相同,而与第三种方案的结果
稍有不同,但差距极小,只要是由于不同的计算方法导致运算过程中近似程度不同,导致最终的精度不同。
实验2 研究三个平面的位置关系
实验内容: 用MATLAB研究下面的3个平面
p1: x + y + z = 1
p2: ?x + y = 2
p3: 2x + t2z = t
当t取何值时交于一点? 当t取何值时交于一直线? 当t取何值时没有公共的交点? 并在每一种情形下, 用MATLAB
您可能关注的文档
最近下载
- 《改造我们的学习》课件+2024-2025学年统编版高中语文选择性必修中册.pptx VIP
- T∕CECS 10028-2019 绿色建材评价 钢结构房屋用钢构件(可复制版).pdf
- 毕业典礼的活动策划方案.pptx VIP
- A4线缆标签打印模板.xls VIP
- 2022年高中学业水平考试试题及答案9科.docx VIP
- 临潭县住房和城乡建设局社会稳定风险评估工作实施方案.pdf VIP
- 新北区苏科版二年级上册劳动《6、树叶书签》教案.docx VIP
- 试验设计与数据处理(第三版)李云雁-第4章-回归分析.pptx VIP
- 西子奥的斯OH-CON电气原理图纸XAA21310S 一体化LMCB图纸.pdf
- A4标签打印模板.xls VIP
文档评论(0)