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

1、线性方程组求解的数值实验 线性方程组求解的数值实验一、课题名称：双对角线线性方程组的数值实验 一、课题名称：双对角线线性方程组的数值实验二、班级和姓名： 二、班级和姓名：09101900 张争（学号 张争（学号 0910190161）三、实验问题摘要： 三、实验问题摘要：考虑一种特殊的对角线元素不为零的双对角线线性方程组（以 n=7 为例）=                76 65

2、54 4 33 22 11da da da d ad ad ad                7654321xxxxxxx                7654321bbbbbbb写出解一般的 n（奇数）阶方程组的程序（不要用消元过程，因为不用它可以十分方便的解出这个方程组） 。四、过程中的问题及处理方法： 四、过程中的问题及处理方法：解方程组时

3、应当先解出 x1 和 xn,r 然后根据公式求出其他的解，由于矩阵的特殊性，应该从头尾开始解，最后解出 。编程序时使用循环语句要分情况讨论。21  n x五、计算公式： 五、计算公式：= 1 x11 db(kn/2)kk k kk dx a b x 1   21212123212121     nn n n n nn dx a x a bx六、结构程序设计 六、结构程序设计输入矩阵阶数 输入矩阵阶数 输入各系

4、数 输入各系数 解出前（ 解出前（n-1）/2 个解 个解 解出后（ 解出后（n-1）/2 个解 个解 解出第 解出第    （n+1）/2 个解 个解七、程序运行结果： 七、程序运行结果：这里令所有的系数等于 这里令所有的系数等于 1，实验结果如下： ，实验结果如下：int a[N],b[N],d[N],x[N];void input(){int i;cout>n;cout>d[i];cout>b[i];}f

5、or(i=1;i>a[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]*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<<“方程组的解