| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 897 人关注过本帖
标题:[求助]高斯消元法的c程序
只看楼主 加入收藏
pp786195
Rank: 1
等 级:新手上路
帖 子:16
专家分:0
注 册:2007-9-20
收藏
 问题点数:0 回复次数:4 
[求助]高斯消元法的c程序
vc++6.0编辑器怎么运行c程序啊,不会用,郁闷。
高斯消元法的思路是什么,现在只看书,一个程序也不会编,真难。
希望高手指点一下
搜索更多相关主题的帖子: 高斯 
2007-10-31 14:17
cdutli
Rank: 1
等 级:新手上路
帖 子:4
专家分:0
注 册:2007-10-28
收藏
得分:0 
数值分析书上有撒


2007-10-31 15:06
cdutli
Rank: 1
等 级:新手上路
帖 子:4
专家分:0
注 册:2007-10-28
收藏
得分:0 
void Gauss_cal(Matrix * pAug_matrix)
{
for(int k = 0; k < ((*pAug_matrix).size()-1); ++k) { //!!max(k)
for(int i = k+1; i < (*pAug_matrix).size(); ++i) {
float l = (*pAug_matrix)[i][k] / (*pAug_matrix)[k][k] ;//!!
(*pAug_matrix)[i][k] = 0;
for(int j = k+1; j < (*pAug_matrix)[i].size(); ++j) {
(*pAug_matrix)[i][j] -= l * (*pAug_matrix)[k][j];
}
}
}
}

//Matrix为一个二维数组

2007-10-31 15:13
冰烨
Rank: 1
来 自:不咔界
等 级:新手上路
威 望:1
帖 子:135
专家分:0
注 册:2007-10-30
收藏
得分:0 

具体思路看数值分析课本吧....
下面是我学数分时候写过的一个"选主元素法"解方程的程序,可以参考一下
#include<stdio.h>
#include<math.h>
#define M 3
#define N 3
main(){
double a[M][M]={10,-1,-2,-1,10,-2,-1,-1,5};//系数
double b[M]={7.2,8.3,4.2};
int k,i,j,l;
double d,t;
k=0;
for(k=0;k<N;k++){
d=a[k][k];
l=k;
i=k+1;
for(i=k+1;i<N;i++){
if(fabs(a[i][k])>fabs(d)){
d=a[i][k];
l=i;
}
}

if(d!=0){
if(l!=k){
for(j=k;j<N;j++){
t=a[l][j];
a[l][j]=a[k][j];
a[k][j]=t;
}
t=b[k];
b[k]=b[l];
b[l]=t;
}
}
for(j=k+1;j<N;j++)
a[k][j]=a[k][j]/a[k][k];
b[k]/=a[k][k];
for(i=k+1;i<N;i++)
for(j=k+1;j<N;j++)a[i][j]=a[i][j]-a[i][k]*a[k][j];
for(i=k+1;i<N;i++)b[i]=b[i]-a[i][k]*b[k];
}
for(i=N-1;i>=0;i--){
for(j=i+1;j<N;j++)
b[i]-=a[i][j]*b[j];
}
for(i=0;i<N;i++)
printf("%lf ",b[i]);
getchar();getchar();
return 0;
}


要么死去,要么奋斗!
2007-10-31 15:21
pp786195
Rank: 1
等 级:新手上路
帖 子:16
专家分:0
注 册:2007-9-20
收藏
得分:0 
楼上的果真是高手

我要学好c语言
2007-11-03 20:18
快速回复:[求助]高斯消元法的c程序
数据加载中...
 
   



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

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