简单迭代法的代码实现.doc

/*简单迭代法的代码实现*/ #include #include #include using namespace std; double e 2.座机电话号码8284; 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 cout a endl; a b; b f a ; k++; if k 100 cout 迭代失败！（可能是函数不收敛） endl; return ; cout b endl; return; int main cout 请输入初始值x0和要求得结果的精度：; double x,d; cin x d; SimpleDiedai x,d ; return 0; /*牛顿迭代法的代码实现*/ #include #include #include using namespace std; double e 2.座机电话号码8284; 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 cout a endl; a b; b f a ; k++; if k 100 cout 迭代失败！（可能是函数不收敛） endl; return ; cout b endl; return; int main cout 请输入初始值x0和要求得结果的精度：; double x,d; cin x d; NewtonDiedai x,d ; return 0; /*雅可比算法的代码实现*/ #include #include #include #include using namespace std; //函数求数组中的最大值 double MaxOfList vector x double max x[0]; int n x.size ; for int i 0;i n;i++ if x[i] max max x[i]; return max; //雅可比迭代公式 void Jacobi vector vector A,vector B,int n vector X n,0 ; vector Y n,0 ; vector D n,0 ; int k 0; //记录循环次数 do X Y; for int i 0;i n;i++ double tem 0; for int j 0;j n;j++ if i! j tem + A[i][j]*X[j]; Y[i] B[i]-tem /A[i][i]; cout left setw 8 Y[i] ; cout endl; k++; if k 100 cout 迭代失败！（可能是函数不收敛） endl; return ; for int a 0;a n;a++ D[a] X[a]-Y[a]; while MaxOfList D 0.00001 || MaxOfList D -0.00001 ; return ; int main int n; cout 请输入方程组未知数的个数n：; cin n; cout endl; vector vector A n,vector n,0 ; vector B n,0 ; cout 请输入方程组的系数矩阵： endl; for int i 0;i n;i++ for int j 0;j n;j++ cin A[i][j]; cout endl; cout 请输入方程组的值向量： endl; for int k 0;k n;k++ cin B[k]; cout endl; cout 您输入的方程组为： endl; for int a 0;a n;a++ for int b 0;b n;b++ cout A[a][b] ; cout B[a] endl; cout endl; cout 由雅可比迭代公式求的方程组的解为： endl; Jaco