| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 908 人关注过本帖
标题:高手请进
只看楼主 加入收藏
SB
Rank: 1
等 级:新手上路
帖 子:3
专家分:0
注 册:2004-10-15
收藏
 问题点数:0 回复次数:5 
高手请进

这是我写的一个用行主元的高斯消去法解方程组的程序,编译已通过了,但不能得到正确答案或者有float point error:domain的错误提示,请各位大虾指导一下.

#define n 4 #include <math.h> void exchrow(double a[n+1][n+1],double b[n+1])/*选主元*/ { double temp,max; int i,j,k,m; for(k=1;k<=n-1;k++) { max=fabs(a[k][k]); for(i=k+1;i<=n;i++) if(fabs(a[i][k]>max)) {max=fabs(a[i][k]); m=i;} if(k!=m) {for(j=1;j<=n;j++) { temp=a[k][j]; a[k][j]=a[m][j]; a[m][j]=temp;} temp=b[m];b[m]=b[k];b[k]=temp;} } } void eliminate(double a[n+1][n+1],double b[n+1])/*消元*/ { int i,j,k; for(k=1;k<=n-1;k++) {for(i=k+1;i<=n;i++) {for(j=k;j<=n;j++) a[i][j]=a[i][j]-a[k][j]*a[i][k]/a[k][k]; b[i]=b[i]-b[k]*a[i][k]/a[k][k]; } } } void backsub(double a[n+1][n+1],double b[n+1])/*回代*/ { int i,j; double sum=0.0; b[n]=b[n]/a[n][n];/*解得的x[i]存放在b[i]中*/ for(i=n-1;i>=1;i--) {for(j=i+1;j<=n;j++) sum+=a[i][j]*b[j]; b[i]=(b[i]-sum)/a[i][i];} } main() {double a[n+1][n+1],b[n+1]; int i,j; printf("input a\n"); for(i=1;i<=n;i++) for(j=1;j<=n;j++) { printf("a[%d][%d]=",i,j); scanf("%f",&a[i][j]);} printf("input b\n"); for(i=1;i<=n;i++) {printf("b[%d]=",i); scanf("%f",&b[i]);} exchrow(a,b); eliminate(a,b); backsub(a,b); printf("root is\n"); for(i=1;i<=n;i++) printf("x[%d]=%f\n",i,b[i]); }

搜索更多相关主题的帖子: 高斯 max fabs double 
2004-10-15 19:48
SB
Rank: 1
等 级:新手上路
帖 子:3
专家分:0
注 册:2004-10-15
收藏
得分:0 

天哪,这里人气这么旺,怎么没人给我指点一下啊

2004-10-15 21:33
mingyue
Rank: 1
等 级:新手上路
帖 子:60
专家分:0
注 册:2004-5-13
收藏
得分:0 

没有仔细研究你的算法,只是大概看了一下语法方面的问题

main() {double a[n+1][n+1],b[n+1]; /*不能用变量来定义数组下标*/ int i,j; printf("input a\n"); for(i=1;i<=n;i++) for(j=1;j<=n;j++) { printf("a[%d][%d]=",i,j); scanf("%f",&a[i][j]);} printf("input b\n"); for(i=1;i<=n;i++) {printf("b[%d]=",i); scanf("%f",&b[i]);} exchrow(a,b); /*还要加上数组个数的参数*/ eliminate(a,b); /*同上*/ backsub(a,b); /*同上*/ printf("root is\n"); for(i=1;i<=n;i++) printf("x[%d]=%f\n",i,b[i]); }


2004-10-15 22:54
live41
Rank: 10Rank: 10Rank: 10
等 级:贵宾
威 望:67
帖 子:12442
专家分:0
注 册:2004-7-22
收藏
得分:0 
不是不想帮,我没学过数值分析。
2004-10-15 23:06
Knocker
Rank: 8Rank: 8
等 级:贵宾
威 望:47
帖 子:10454
专家分:603
注 册:2004-6-1
收藏
得分:0 

简单的语法错误!

scanf("%f",&a[i][j]);} scanf("%f",&b[i]);}

是%lf


九洲方除百尺冰,映秀又遭蛮牛耕。汽笛嘶鸣国旗半,哀伤尽处是重生。     -老K
治国就是治吏。礼义廉耻,国之四维。四维不张,国之不国。   -毛泽东
2004-10-15 23:23
SB
Rank: 1
等 级:新手上路
帖 子:3
专家分:0
注 册:2004-10-15
收藏
得分:0 

3楼大概没看见我开始的宏定义 5楼的意见很对.但还是有错误,谁能帮我看一下算法有什么错误

2004-10-16 20:23
快速回复:高手请进
数据加载中...
 
   



关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.062686 second(s), 7 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved