数值分析迭代法实验报告.docxVIP

数值分析实验报告第一题，雅克比迭代法 题目：略 解决方案：1，截图为：2，程序为 #includeiostreamusingnamespace std;int main(){int n,i,j,k,count;double A[100][100],b[100],x[100]={0},B[100][100],f[100],t[100]={0};cout***叶伟见 计1307 133001***endl; cout***这是用雅克比迭代法求解线性方程组的过程***endl;cout***本程序只能解决方程个数少于100的方程组***endl;cout***系数矩阵是一个n阶方阵,n为方程的个数***endl;cout请输入方程的个数n=;cinn;cout 请依次输入系数矩阵：endl;for(i=0;in;i++){for(j=0;jn;j++)cinA[i][j];}cout请输入方程组的常数矩阵：endl;for(i=0;in;i++){cinb[i];}cout请输入迭代次k的值k=;cink;for(i=0;in;i++){for(j=0;jn;j++){if(i==j)B[i][j]=0;elseB[i][j]=-A[i][j]/A[i][i];}}for(i=0;in;i++){f[i]=b[i]/A[i][i];}for(count=1;count=k;count++){for(i=0;in;i++){for(j=0;jn;j++){t[i]=t[i]+B[i][j]*x[j];}t[i]=t[i]+f[i];}for(i=0;in;i++){x[i]=t[i];t[i]=0;}}cout方程组的解集是：endl;for(i=0;in;i++){coutx[i] ;}return 0;}第二题，高斯赛德尔迭代法 题目，略 解决方案：1，截图2，程序#includeiostreamusingnamespace std;int main(){int n,i,j,k,count;double A[100][100],b[100],x[100]={0},B[100][100],f[100],t[100]={0}; cout***这是用高斯赛德尔迭代法求解线性方程组的过程***endl;cout***本程序只能解决方程个数少于100的方程组***endl;cout***系数矩阵是一个n阶方阵,n为方程的个数***endl;cout请输入方程的个数n=;cinn;cout 请依次输入系数矩阵：endl;for(i=0;in;i++){for(j=0;jn;j++)cinA[i][j];}cout请输入方程组的常数矩阵：endl;for(i=0;in;i++){cinb[i];}cout请输入迭代次k的值k=;cink;for(i=0;in;i++){for(j=0;jn;j++){if(i==j)B[i][j]=0;elseB[i][j]=-A[i][j]/A[i][i];}}for(i=0;in;i++){f[i]=b[i]/A[i][i];}for(count=1;count=k;count++){for(i=0;in;i++){for(j=0;jn;j++){ t[i]=t[i]+B[i][j]*x[j];}x[i]=t[i]+f[i]; t[i]=0;}}cout方程组的解集是：endl;for(i=0;in;i++){coutx[i] ;}return 0;}第三题，用二分法求解方程的零点 题目：用二分法求方程f(x)=4x*x*x+x*x-10 在区间[1，1.5]上的根，要求求出具有3位有效数的近似根 解决方法：1程序:#includeiostream#includecmathusing namespace std;double f(double x){double f;f=x*x*x+4*x*x-10;return f;}int main(){double a=1,b=1.5;int m;cout请输入零点值的有效位数:;cinm;while((b-a)=0.5*pow(10,-m)){if(f((a+b)/2)==0)break;elseif(f(a)*f((a+b)/2)0)b=(a+b)/2;elsea=(a+b)/2;}coutf零点的值是：(a+b)/2;return 0;} 2,截图：