计算传热作业1.docVIP

计算传热学课程大作业报告 PAGE \* MERGEFORMAT PAGE \* MERGEFORMAT 1 储运与建筑工程学院能源与动力工程系 计算传热学课程大作业报告 作业题目：代数方程组的求解 学生姓名：田 学 号： 专业班级：能动1 专业写全称，仿宋GB-2312 专业写全称，仿宋GB-2312，三号字体。 专业写全称，仿宋GB-2312 专业写全称，仿宋GB-2312，三号字体。 2017年9月23日  目 录 TOC \o 1-3 \h \z \u HYPERLINK \l _Toc345138880 一、计算题目 3 HYPERLINK \l _Toc345138881 二、离散方程 3 HYPERLINK \l _Toc345138887 三、程序设计 4 HYPERLINK \l _Toc345138888 3.1 高斯赛德尔迭代法 4 HYPERLINK \l _Toc345138889 3.2 TDMA法 5 四、 HYPERLINK \l _Toc345138893 程序及计算结果验证 6 HYPERLINK \l _Toc345138894 五、网格独立性考核 PAGEREF _Toc345138894 \h 7 HYPERLINK \l _Toc345138888 3.1 高斯赛德尔迭代法 7 HYPERLINK \l _Toc345138889 3.2 TDMA法8 HYPERLINK \l _Toc345138894 六、结果分析与结论 8 HYPERLINK \l _Toc345138888 3.1 高斯赛德尔迭代法 9 HYPERLINK \l _Toc345138889 3.2 TDMA法10 一、计算题目 分别用高斯赛德尔迭代和TDMA方法求解方程 （1） 在=-5，-1,0,1,5情况下的解，并表示在图中。 其中，=0，=0；=1，=1. 二、离散方程 采用控制容积法： 即 （2） （3） 假设均分网格，则有 上式则变为： （4） 即 （5） （6） 程序设计 3.1 高斯赛德尔迭代法 由已知公式 可设计高斯赛德尔迭代C语言程序如下： #include stdio.h #include stdlib.h #include math.h int main() { double e=0,x; int i,j,b,k,d; double a[100]; scanf(%lf%d,x,d); for(j=0;j1/x;j++) { a[j]=0; } b=1/x; a[b]=1; while(1){ for(i=0;i1/x-1;i++) { a[i+1]=((2-x*d)*a[i+2]+(2+x*d)*a[i])/4; printf(i = %d\n,i); if(i==1/x-3) e=a[i+2]; } if(fabs((a[i]-e))/a[i]0.00001) break; } for(k=0;k=1/x;k++) { printf(%lf ,a[k]); } system(pause); return 0; } 其中，输入项x和d分别代表步长和系数的值，即和的值，输出a[k]即为的值。 3.2 TDMA法 由已知公式 （7） 目标方程 （8） 联立（7）、（8）方程可得： （9） 则有 ， 利用端点条件 =0，=0；=1，=1 可得： , 到此，可设计C语言计算程序如下： #include stdio.h #include stdlib.h #include math.h int main() { double e = 0, x; int i, j, b, k, d; double a[100], m[100]; scanf(%lf%d, x, d