简单迭代法(C++版)

本程序为本人参考一些资料后编写,在visual c++6.0下调试通过。

[@more@]

///
/简单迭代法:f(x)=0改写成x=g(x)
/原理:x_(k+1)=g(x_k)
/有条件限制,可能不收敛
/改进:输出每一次迭代的结果,将迭代模块做全面
///
#include
#include
#include
using namespace std;
//
double fun(double x);
void SimpleIterative(double(*pf)(double x));
//
//
int main()
{
SimpleIterative(fun);
return 0;
}
//
//
double fun(double x)
{
return pow(x+1,1/3.);///f(x)=x^3-x-1转化成x^3=x+1
}
/
void SimpleIterative(double(*pf)(double x))
{
double x,x0;
int maxflag,k=1;
double err;
cout<>maxflag;
cout<>err;
cout<>x;
cout< do
{
x0=x;
x=pf(x0);
cout< ++k;
}while(fabs(x-x0)>err&&k if(k>=maxflag)
cout< else
{
cout< cout< cout< }
}
/
/

来自 “ ITPUB博客 ” ,链接:http://blog.itpub.net/11355887/viewspace-974945/,如需转载,请注明出处,否则将追究法律责任。

转载于:http://blog.itpub.net/11355887/viewspace-974945/

  • 0
    点赞
  • 4
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
好的,以下是一个简单的高斯赛德尔迭代法C++代码实现: ```c++ #include <iostream> #include <cmath> using namespace std; double a[3][3] = {{10, 2, 1}, {1, 5, 1}, {2, 3, 10}}; // 系数矩阵 double b[3] = {7, -8, 6}; // 常数向量 int main() { double x[3] = {0}; // 初始解 double x_new[3] = {0}; // 新解 int max_iter = 10000; // 最大迭代次数 double tol = 1e-8; // 精度要求 int k = 0; while (k < max_iter) { for (int i = 0; i < 3; i++) { double sum = 0; for (int j = 0; j < 3; j++) { if (i != j) { sum += a[i][j] * x[j]; } } x_new[i] = (b[i] - sum) / a[i][i]; } bool flag = true; for (int i = 0; i < 3; i++) { if (fabs(x_new[i] - x[i]) > tol) { flag = false; break; } } if (flag) { cout << "迭代次数: " << k << endl; cout << "解: " << x_new[0] << " " << x_new[1] << " " << x_new[2] << endl; break; } for (int i = 0; i < 3; i++) { x[i] = x_new[i]; } k++; } if (k == max_iter) { cout << "迭代次数达到最大值" << endl; } return 0; } ``` 上述代码中,我们首先定义了系数矩阵 `a` 和常数向量 `b`,然后设定了初始解 `x` 和新解 `x_new` 的初始值为0。接着,我们设定了最大迭代次数 `max_iter` 和精度要求 `tol`,并且开始进行迭代求解。在每次迭代中,我们按照高斯赛德尔迭代法的公式,计算出新的解 `x_new`,然后判断新解和旧解之间的误差是否满足精度要求,如果满足,就输出解并结束迭代,否则将新解赋值给旧解,并继续迭代。如果迭代次数达到最大值,说明迭代方法没有收敛,输出错误信息。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值