/*

超松弛迭代求解线性方程组

*/

#include

#include

#define MAX_N 20  /*方程最大维数*/

#define MAXREPT 100

#define EPSILON 0.00001  /*求解精度*/

int main()

{

int n,i,j,k;

double err,w;

double a[MAX_N][MAX_N],b[MAX_N][MAX_N],c[MAX_N],g[MAX_N],x[MAX_N],nx[MAX_N];

printf("input n value(dim ofAX=C):");  /*输入方程的维数*/

while(scanf("%d",&n)==1&&n>0&&n<=MAX_N){  /*以0结束输入*/

printf("now input the matrix a(i),i=0,...,%d:\n",n-1);

for(i=0;i

for(j=0;j

scanf("%lf",&a[i][j]);

printf("now input the matrix c(i),i=0,..,%d:\n",n-1);

for(i=0;i

scanf("%lf",&c[i]);

printf("now input the w value:");

while(scanf("%lf",&w)!=1||w<1||w>2)

printf("w must between 1 and 2,please enter again.\n");

for(i=0;i

for(j=0;j

b[i][j]=-a[i][j]/a[i][i];

g[i]=c[i]/a[i][i];

}

for(i=0;i

for(j=0;j

nx[j]=g[j];

for(j=0;j

for(k=0;k

nx[j]+=b[j][k]*nx[k];  /*迭代*/

for(k=j+1;k

nx[j]+=b[j][k]*x[k];

nx[j]=(1-w)*x[j]+w*nx[j];

}

err=0;

for(j=0;j

if(err

err=fabs(nx[j]-x[j]);  /*误差计算*/

for(j=0;j

x[j]=nx[j];

}

if(err

printf("solve problem:x_i=\n");

for(i=0;i

printf("%lf ",x[i]);

printf("\n");

}

continue;

printf("after %d repeat,no result...\n",MAXREPT);

}

return 0;

}

输入松弛因子后就停止了   下面的程序无法实现怎么办?

迭代法的c语言程序,松弛迭代法C程序相关推荐

  1. 用c语言编制牛顿法程序,求解试用newton法求函数,YTU 2405: C语言习题 牛顿迭代法求根...

    2405: C语言习题 牛顿迭代法求根 时间限制: 1 Sec  内存限制: 128 MB 提交: 562  解决: 317 题目描述 用牛顿迭代法求根.方程为ax3+bx2+cx+d=0.系数a,b ...

  2. 多重网格法-松弛迭代法-二维泊松方程-python实现

    这几天在家躲避疫情,闲来无事,写了这个多重网格法求解泊松方程的算法的代码. 多重网格法可能是目前为止解泊松方程最快的算法,n个格点需要n次计算就可以收敛,而快速傅里叶变换的收敛速度是n*logn, 共 ...

  3. jacobi matlab程序,jacobi迭代法实验MATLAB程序数值分析

    jacobi迭代法实验MATLAB程序数值分析 例1. 求线性方程组 得近似解.精确解为x*=[3,2,1]'. 解:对方程进行移项就得 记为Ax=b,或写为x=B0 x+f,其中 取初始值,代入原方 ...

  4. C++学习日记#1.1——四元一次方程组进行Sor松弛迭代法求解(现已推向n元一次方程组进行Sor松弛迭代法求解)

    在这篇文章中,不会有太多的讲解,主要是因为Sor迭代法和高斯-赛德尔(gauss-seidle)迭代法的程序只有迭代公式的一点点区别,当松弛因子为1的时候,Sor迭代法和高斯-赛德尔(gauss-se ...

  5. 基于交替迭代法的交直流混合系统潮流计算matlab程序iEEE9节点系统算例

    基于交替迭代法的交直流混合系统潮流计算matlab程序iEEE9节点系统算例 由于受端负荷量持续增大,直流输电受端系统的电压稳定性能自然引起了人们的关注.目前,交直流混合电网的电压稳定分析大致分为两类 ...

  6. c语言牛顿迭代法求解非线性方程组,利用牛顿迭代法求解非线性方程组

    最近一个哥们,是用牛顿迭代法求解一个四变量方程组的最优解问题,从网上找了代码去改进,但是总会有点不如意的地方,迭代的次数过多,但是却没有提高精度,真是 最近一个哥们,是用牛顿迭代法求解一个四变量方程组 ...

  7. c语言写输入汉字输出区位码程序,2017年计算机应用基础检测试题

    2017年计算机应用基础检测试题 计算机科学是一门包含各种各样与计算和信息处理相关主题的系统学科,从抽象的算法分析.形式化语法等等.下面是小编整理的关于计算机应用基础检测试题及答案,希望大家认真阅读! ...

  8. c语言编程 构建围墙,c语言程序设计朝盛 综合程序练习题.ppt

    c语言程序设计朝盛 综合程序练习题 综合程序练习题 综合程序练习题 大型程序设计的要求 综合程序练习题 综合程序练习题 大型程序设计题目分析: 1 打字练习 2 电子琴 3 简易画图板 4 贪吃蛇 5 ...

  9. python是c语言写的吗-先学C语言还是Python?资深程序员往往是这样建议的!

    编程是一门手艺, 学习编程首先要有目的,w3cschool先带大家了解一下自己需要学习什么吧! 0.学一门语言 其实学什么语言都可以,因为语言都是互通的,C语言也好, 对于初学者来说,Python 也 ...

最新文章

  1. window 添加环境变量
  2. SQL进阶教程 | 史上最易懂SQL教程 5小时零基础成长SQL大师
  3. python什么是数据结构_〖Python〗-- 数据结构
  4. sql limit 的用法
  5. vue从其它页面返回_vue页面按返回键等跳转重定向判断
  6. weblogic 8.1 安装并配置数据源
  7. SAP UI5 bindItem will cause OData refresh
  8. 深入理解Java的整型类型:如何实现2+2=5?
  9. igmp是哪个层协议_通俗易懂网络协议(IP)
  10. java开发微信公众平台备忘
  11. CentOS 5.4 制作 Python 2.6 RPM 包的方法
  12. 更适合智能家庭使用的新 Wi-Fi 技术问世了
  13. java微信支付超时_Java微信支付之关闭订单
  14. qq批量登录软件_QQ账号永久冻结
  15. 紫书刷题记录 UVa12100 打印队列
  16. 使用Arduino+L298N控制光驱两项四线步进电机
  17. 问题:During startup program exited with code 0xc0000135
  18. python列表的负数索引
  19. windows源文件名称大于文件系统支持的长度无法删除问题
  20. 公益中国系列活动进社区之 “健康进社区”系列活动

热门文章

  1. 浅谈软件设计的七大原则
  2. 0成本信息差项目:地图标注该如何做?
  3. 视频和PPT已公开 | 第八届中国互联网测试大会·2021 MTSC 上海站
  4. 永信至诚助“海南省首届网络安全大赛”决赛圆满收官
  5. MP4视频测试URL地址,亲测有效
  6. Python:好租网商务楼数据爬取
  7. 秒云与趋动科技联合发布容器云平台与GPU资源池化整体解决方案
  8. 4、SRS4.0源代码分析之RTMP推流处理
  9. mysql解决中文乱码问题
  10. 【大型电商项目开发】商品上架-es应用到商品上架-35