简单迭代法及代码实现.docVIP

/*简单迭代法的代码实现*/ #includeiostream #includestring #includecmath using namespace std; double e=2.718281818284; double f(double x){ return pow(e,-1*x); } void SimpleDiedai(double x,double d){ double a=x; double b=f(a); int k=0;//记录循环的次数 while(((a-b)d) || ((a-b)-1*d)){ coutaendl; a=b; b=f(a); k++; if(k100){ cout迭代失败！（可能是函数不收敛）endl; return ; } } coutbendl; return; } int main(){ cout请输入初始值x0和要求得结果的精度：; double x,d; cinxd; SimpleDiedai (x,d); return 0; } /*牛顿迭代法的代码实现*/ #includeiostream #includestring #includecmath using namespace std; double e=2.718281818284; double f(double x){ double a=pow(e,-1*x); return x-(x-a)/(1+a); } void NewtonDiedai(double x,double d){ double a=x; double b=f(a); int k=0; //记录循环的次数 while(((a-b)d) || ((a-b)-1*d)){ coutaendl; a=b; b=f(a); k++; if(k100){ cout迭代失败！（可能是函数不收敛）endl; return ; } } coutbendl; return; } int main(){ cout请输入初始值x0和要求得结果的精度：; double x,d; cinxd; NewtonDiedai(x,d); return 0; } /*雅可比算法的代码实现*/ #includeiostream #includeiomanip #includestring #includevector using namespace std; //函数求数组中的最大值 double MaxOfList(vectordoublex){ double max=x[0]; int n=x.size(); for(int i=0;in;i++) if(x[i]max) max=x[i]; return max; } //雅可比迭代公式 void Jacobi(vectorvectordouble A,vectordouble B,int n){ vectordouble X(n,0); vectordouble Y(n,0); vectordouble D(n,0); int k=0; //记录循环次数 do{ X=Y; for(int i=0;in;i++){ double tem=0; for(int j=0;jn;j++){ if(i!=j) tem += A[i][j]*X[j]; } Y[i]=(B[i]-tem)/A[i][i]; coutleftsetw(8)Y[i] ; } coutendl; k++; if(k100){ cout迭代失败！（可能是函数不收敛）endl; return ; } for(int a=0;an;a++){ D[a]=X[a]-Y[a]; } }while( MaxOfList(D)0.00001 || MaxOfList(D)-0.00001); return ; } int main(){ int n; cout请输入方程组未知数的个数n：; cinn; coutendl; vectorvectordouble A(n,vectordouble(n,0)); vectordoubleB(n,0); cout请输入方程组的系数矩阵：endl; for(int i=0;in;i++){ for(int j=0;jn;j++){ cinA[i][j]; } } coutendl; cout请输入方程组的值向量：endl; fo