原创力文档- 1、本文档共10页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
10-11-2《几何与代数》
数学实验报告
学号姓名：周志浩 得分： .
实验一：平板的稳态温度分布问题（线性方程组应用）
在热传导的研究中，一个重要的问题是确定一块平板的稳态温度分布。假定下图中的平板代表一条金属梁的截面，并忽略垂直于该截面方向上的热传导。
已知平板内部有9个节点，每个节点的温度近似等于与它相邻的四个节点温度的平均值，例如,；为避免出现分数，可写成。设4条边界上的温度分别等于每位同学学号的后四位的4倍，例如学号同学计算时，选择、、、。
求：（1）建立可以确定平板内节点温度的线性方程组；
（2）用MATLAB软件的三种方法求解该线性方程组 (请输出精确解（分数形式）) ；
方法一：利用Cramer法则求解；
方法二：作为逆矩阵的方法求解；
方法三：利用Gauss消元法即通过初等行变换求解。
实验部分
构造的线性方程：
方法一：Cramer法则
format rat
a1=[4,-1,0,-1,0,0,0,0,0,];a2=[-1,4,-1,0,-1,0,0,0,0];a3=[0,-1,4,0,0,-1,0,0,0];a4=[-1,0,0,4,-1,0,-1,0,0];a5=[0,-1,0,-1,4,-1,0,-1,0];a6=[0,0,-1,0,-1,4,0,0,-1];a7=[0,0,0,-1,0,0,4,-1,0];a8=[0,0,0,0,-1,0,-1,4,-1];a9=[0,0,0,0,0,-1,0,-1,4];b=[12,0,4,12,0,4,24,12,16];
D=det([a1,a2,a3,a4,a5,a6,a7,a8,a9]),
D =
100352
D1=det([b,a2,a3,a4,a5,a6,a7,a8,a9]),D2=det([a1,b,a3,a4,a5,a6,a7,a8,a9]),D3=det([a1,a2,b,a4,a5,a6,a7,a8,a9]),D4=det([a1,a2,a3,b,a5,a6,a7,a8,a9]),D5=det([a1,a2,a3,a4,b,a6,a7,a8,a9]),D6=det([a1,a2,a3,a4,a5,b,a7,a8,a9]),D7=det([a1,a2,a3,a4,a5,a6,b,a8,a9]),D8=det([a1,a2,a3,a4,a5,a6,a7,b,a9]),D9=det([a1,a2,a3,a4,a5,a6,a7,a8,b]),
D1 =
630784
D2 =
419328
D3 =
344064
D4 =
899584
D5 =
702464
D6 =
555520
D7 =
1060864
D8 =
935424
D9 =
774144
T1=D1/D,T2=D2/D,T3=D3/D,T4=D4/D,T5=D5/D,T6=D6/D,T7=D7/D,T8=D8/D,T9=D9/D
T1 =
44/7
T2 =
117/28
T3 =
24/7
T4 =
251/28
T5 =
7
T6 =
155/28
T7 =
74/7
T8 =
261/28
T9 =
54/7
方法二：逆矩阵a1=[4,-1,0,-1,0,0,0,0,0,];a2=[-1,4,-1,0,-1,0,0,0,0];a3=[0,-1,4,0,0,-1,0,0,0];a4=[-1,0,0,4,-1,0,-1,0,0];a5=[0,-1,0,-1,4,-1,0,-1,0];a6=[0,0,-1,0,-1,4,0,0,-1];a7=[0,0,0,-1,0,0,4,-1,0];a8=[0,0,0,0,-1,0,-1,4,-1];a9=[0,0,0,0,0,-1,0,-1,4]; b=[12,0,4,12,0,4,24,12,16];
A=[a1,a2,a3,a4,a5,a6,a7,a8,a9],B=[12,0,4,12,0,4,24,12,16],
A =
Columns 1 through 5
4 -1 0 -1 0
-1 4 -1 0 -1
文档评论（0）