计算方法实验报告.doc

计算方法 实验6 计算方法基本算法练习实验6 计算方法算法练习课程实验的训练重点在于根据数学模型提出求解的数值计算方法，编写程序上机计算结果。通过实验，掌握数学问题求解的数值方法及与此有关的理论，注意程序设计的可行性和复杂性，考虑算法时结合计算机的特点来实现并能根据结果程序结果进行算法分析。 实验目的(1) 根据数学模型提出求解的数值计算方法(2) 掌握控制语句与二维数组的使用 实验内容1（P59 例3）用雅克比迭代法法解方程组精度为0.000001 迭代初始值x1,x2,x3,x4分别为1，0，1，0 实验内容2高斯-塞德尔迭代法法解上面的方程组，迭代初始值x1,x2,x3,x4分别为1，0，1，0精度为0.000001 实验要求(1) 认真分析题目。(2) 进行算法设计。(3) 编写程序代码(4) 上机调试程序。(5) 保存和打印出程序的运行结果，并结合程序进行分析。计算方法 实验报告模板实验报告6： 计算方法基本算法练习时间：第7周 星期2第二大节 2009.3.31地点：实验楼311 座号： 48 姓名： 吴忠山 学号： 0640330218 班级：06计算2班实验内容第1题：1思路阐述： 解线性方程组的迭代法与解线性方程组的迭代法类似，是把非奇异线性方程组化为等价方程组，对任意给定的初始向量，有迭代公式可得X（k+1）=BX（k）+g，得到向量序列X（k），当X（k）收敛与a时，a即为所给线性方程组的解。2程序代码及其运行结果 #include stdio.h#include conio.h#include math.hmain() int n,i,j,h; float a1010,b10,x5050,p,r; printf(输入维数 n:); scanf(%d,&n); printf(输入:n); for(i=0;in;i+) for(j=0;jn;j+) printf(a%d%d,i,j); scanf(%f,&aij); printf(b%d,i); scanf(%f,&bi); printf(输入的增广矩阵为：n); for(i=0;in;i+) for(j=0;jn;j+) printf(%f ,aij); printf(%fn,bi); printf(输入精度 p：); scanf(%f,&p); printf(输入的精度为：%fn,p); for(i=0;in;i+) printf(输入 x%d0：,i); scanf(%f,&xi); printf(输入的初值为：); for(i=0;in;i+) printf(x%d0=%fn,i,xi0); h=0; dofor(i=0;in;i+) r=0; for(j=0;jp); printf(答案为：); for(i=0;in;i+) printf(%f ,xih); printf(n迭代次数为%d,h+1); getch();3结论分析或说明： 雅各比迭代法简单，但运算量大，迭代次数多，这是因为他的更新速度慢，要算出所有的X（k+1）才能算X（k+2）。实验内容第2题：1思路阐述： 高斯-赛德尔迭代法是在雅各比的基础上，进行了改进，每算出一个新分量便立即用它取代对应的旧分量进行迭代。2程序代码及其运行结果 #include stdio.h#include conio.h#include math.hmain() int n,i,j,h; float a1010,b10,x50,p,r,c,d; printf(输入维数 n:); scanf(%d,&n); printf(输入:n); for(i=0;in;i+) for(j=0;jn;j+) printf(a%d%d,i,j); scanf(%f,&aij); printf(b%d,i); scanf(%f,&bi); printf(输入的增广矩阵为：n); for(i=0;in;i+) for(j=0;jn;j+) printf(%f ,aij); printf(%fn,bi); printf(输入精度 p：); scanf(%f,&p); printf(输入的精度为：%fn,p); for(i=0;in;i+) printf(输入 x%d0：,i); scanf(%f,&xi); printf(输入的初值为：); for(i=0;in;i+) printf(x%d=%fn,i,xi); h=0; do for(i=0;in;i+) r=0; for(j=0;jp); printf(答案为：); for(i=0;in;i+)