10112几何与代数数学实验1.docVIP

10-11-2《几何与代数》 数学实验报告 学号： 姓名： 得分： . 要求：报告中应包含实验中你所输入的所有命令及运算结果，在第16周之前。 实验一：平板的稳态温度分布问题（线性方程组应用） 在热传导的研究中，一个重要的问题是确定一块平板的稳态温度分布。假定下图中的平板代表一条金属梁的截面，并忽略垂直于该截面方向上的热传导。 已知平板内部有9个节点，每个节点的温度近似等于与它相邻的四个节点温度的平均值，例如,；为避免出现分数，可写成。设4条边界上的温度分别等于每位同学学号的后四位的4倍，例如学号同学计算时，选择、、、。 求：（1）建立可以确定平板内节点温度的线性方程组； （2）用MATLAB软件求解该线性方程组 ； 方法一：利用Cramer法则求解； 方法二：作为矩阵求解； 方法三：利用Gauss消元法求解。（1）建立可以确定平板内节点温度的线性方程组；4T1-T2-T4=12 -T1+4T2-T5-T3=0 -T2+4T3-T6=16 -T1+4T4-T5-T7=12 -T2-T4+4T5-T6-T8=0 -T3-T5+4T6-T9=16 -T4+4T7-T8=16 -T5-T7+4T8-T9=4 -T6-T8+4T9=20 （2）用MATLAB软件求解该线性方程组 ； 方法一：利用Cramer法则求解； a=[4,-1,0,-1,0,0,0,0,0]; b=[-1,4,-1,0,-1,0,0,0,0]; c=[0,-1,4,0,0,-1,0,0,0]; d=[-1,0,0,4,-1,0,-1,0,0]; e=[0,-1,0,-1,4,-1,0,-1,0]; f=[0,0,-1,0,-1,4,0,0,-1]; g=[0,0,0,-1,0,0,4,-1,0]; h=[0,0,0,0,-1,0,-1,4,-1]; i=[0,0,0,0,0,-1,0,-1,4]; j=[12,0,16,12,0,16,16,4,20]; D=det([a,b,c,d,e,f,g,h,i]), D = 100352 D1=det([j,b,c,d,e,f,g,h,i]),D2=det([a,j,c,d,e,f,g,h,i]),D3=det([a,b,j,d,e,f,g,h,i]), D4=det([a,b,c,j,e,f,g,h,i]),D5=det([a,b,c,d,j,f,g,h,i]),D6=det([a,b,c,d,e,j,g,h,i]), D7=det([a,b,c,d,e,f,j,h,i]),D8=det([a,b,c,d,e,f,g,j,i]),D9=det([a,b,c,d,e,f,g,h,j]), D1 = 659456 D2 = 566272 D3 = 802816 D4 = 867328 D5 = 802816 D6 = 1039360 D7 = 802816 D8 = 738304 D9 = 946176 format rat 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 = 46/7 T2 = 79/14 T3 = 8 T4 = 121/14 T5 = 8 T6 = 145/14 T7 = 8 T8 = 103/14 T9 = 66/7 方法二：作为矩阵求解； format rat A=[4,-1,0,-1,0,0,0,0,0;-1,4,-1,0,-1,0,0,0,0;0,-1,4,0,0,-1,0,0,0;-1,0,0,4,-1,0,-1,0,0; 0,-1,0,-1,4,-1,0,-1,0;0,0,-1,0,-1,4,0,0,-1;0,0,0,-1,0,0,4,-1,0;0,0,0,0,-1,0,-1,4,-1; 0,0,0,0,0,-1,0,-1,4],B=[12;0;16;12;0;16;16;4;20], A = Columns 1 through 5 4 -1 0 -1 0 -1