线性方程组求解的数值实验报告课程.docVIP

线性方程组求解的数值实验 一、课题名称：双对角线线性方程组的数值实验 二、班级和姓名 张争（学号 0910190161） 三、实验问题摘要： 考虑一种特殊的对角线元素不为零的双对角线线性方程组（以n=7为例） = 写出解一般的n（奇数）阶方程组的程序（不要用消元过程，因为不用它可以十分方便的解出这个方程组）。 四、过程中的问题及处理方法： 解方程组时应当先解出x1和xn,r然后根据公式求出其他的解，由于矩阵的特殊性，应该从头尾开始解，最后解出。编程序时使用循环语句要分情况讨论。 五、计算公式： = (kn/2) = (kn/2) 六、结构程序设计 输入矩阵阶数输入各系数解出前（n-1）/2个解解出后（n-1）/2个解解出第（n+1）/2个解 七、程序运行结果： 这里令所有的系数等于1，实验结果如下： 八、主要程序段介绍： 解方程函数：由于矩阵具有好的对称性，解的形式也一样，用循环语句。但是第一个和最后一个解以及中间一个解的形式不同，需要单独解出，程序如下： void jiefangcheng() { int i; x[1]=b[1]/d[1]; for(i=2;i(i+1)/2;i++) x[i]=(b[i]-a[i-1]*x[i-1])/d[i]; x[n]=b[n]/d[n]; for(i=n-1;i(i+1)/2;i--) x[i]=(b[i]-a[i]*x[i+1])/d[i]; x[(n+1)/2]=(a[(n-1)/2]*x[(n-1)/2]-a[(n+1)/2]*x[(n+3)/2])/d[(n+1)/2]; } 九、源程序如下： //*****************张争******学号0910190161********** #includeiostream.h #define N 10000 int n; int a[N],b[N],d[N],x[N]; void input() { int i; cout请输入方程组的阶数：; cinn; cout请输入方程组的系数：endl; for(i=1;i=n;i++) { coutdi=; cind[i]; coutbi=; cinb[i]; } for(i=1;in;i++) { coutai=; cina[i]; } } void jiefangcheng() { int i; x[1]=b[1]/d[1]; for(i=2;i(i+1)/2;i++) x[i]=(b[i]-a[i-1]*x[i-1])/d[i]; x[n]=b[n]/d[n]; for(i=n-1;i(i+1)/2;i--) x[i]=(b[i]-a[i]*x[i+1])/d[i]; x[(n+1)/2]=(a[(n-1)/2]*x[(n-1)/2]-a[(n+1)/2]*x[(n+3)/2])/d[(n+1)/2]; } void output() { int i; cout方程组的解为：endl; for(i=1;i=n;i++) coutxi=x[i]endl; } void main() { input(); jiefangcheng(); output(); }Acknowledgments The authors would like to thank Johns Hopkins University for the TC-1 cells. This work was supported by a National Health Research Institutes intramural grant (IV-103-PP-22) and grants from the National Science Council, which were awarded to Y.C. Song (NSC 99-2321-B-400-004-MY3) and S.J. Liu (NSC 103-2321-B-400-008). Author Contributions Y.C.S. and S.J.L. designed the studies. Y.C.S. performed the research and analyzed the data. Y.C.S. and S.J.L. wrote the manuscript. Additional Information C57BL/6 mice were immunized subcutaneously (s.c.) once with