数学建模社区-数学中国

标题: [讨论]高斯消元法的程序对比 [打印本页]

---

作者: god    时间: 2005-1-19 17:01
标题: [讨论]高斯消元法的程序对比
下面我收集的高斯消元法的程序。还是请大家看看，并提点意见。
[attach]966[/attach]

[讨论]高斯消元法的程序对比.rar

2005-1-19 17:00 上传

点击文件名下载附件

下载积分: 体力 -2 点

1.68 KB, 下载次数: 78, 下载积分: 体力 -2 点

[讨论]高斯消元法的程序对比

[讨论]高斯消元法的程序对比.rar

2005-1-19 17:01 上传

点击文件名下载附件

下载积分: 体力 -2 点

1.09 KB, 下载次数: 67, 下载积分: 体力 -2 点

[讨论]高斯消元法的程序对比

---

作者: chzjao    时间: 2005-9-24 11:45
标题: 发下我编的,请各位帮忙修改下
[attach]2037[/attach]
! D2 L$ {$ B2 [; U. P

[此贴子已经被作者于2005-9-24 11:46:03编辑过]

[讨论]高斯消元法的程序对比.txt

2005-9-24 11:44 上传

点击文件名下载附件

下载积分: 体力 -2 点

1.14 KB, 下载次数: 41, 下载积分: 体力 -2 点

[讨论]高斯消元法的程序对比

---

作者: jiangzsm    时间: 2006-4-17 14:42

还没看,感谢分享.

---

作者: lfer2008    时间: 2006-11-26 13:42

谢谢

---

作者: ylsi41    时间: 2006-12-4 23:50

谢谢

---

作者: ladeng321    时间: 2007-1-2 19:38
谢谢分享

---

作者: fansmiling    时间: 2007-6-18 12:17

8错8错

---

作者: LH611    时间: 2007-12-18 21:57
谢谢分享，要又AS写的就更好了！

---

作者: 水木年华zzu    时间: 2007-12-19 19:04
总的来说写的不错!继续努力

---

作者: shuainaruto    时间: 2009-4-11 11:05
感谢分享，下来研究一下

---

作者: hduscimc    时间: 2009-6-21 19:04
#include<iostream>
: T5 q9 Y3 ]3 R- g- ^6 ~$ Q#include<cmath>
- A7 R. d5 ^8 q* L: tusing namespace std;
! m& t: @5 X! V* T/ s' kint main()
& h0 J5 A7 Q( P) h2 R+ g{
  v7 q9 l$ G5 H        int i;
+ ?1 j; Z2 s+ ^1 V  r        float *x;
        float c[3][4] ={0.101,2.304,3.555,1.183,
                                        -1.347,3.712,4.623,2.137,
                                        -2.835,1.072,5.643,3.035};
$ S0 d) u& l3 _# @        float *ColPivot(float *,int);
* y, }) s, Q' V! D, b9 m        x=ColPivot(c[0],3);
        for( i=0;i<=2;i++ )
2 a, d5 D9 \  \                cout<<"x("<<i<<")="<<x[i]<<endl;
}
7 A" S% _6 b( u/ R9 m) k
- z6 t. E* _) ]+ ^$ Ofloat *ColPivot( float *c,int n )
{
" l* h3 Q% d! X4 l* v. [        int i,j,t,k;
        float *x,p;
6 E( d0 H- l) A) w9 [         x=new float[n*sizeof(float)];
        for( i=0;i<=n-2;i++)
        {
$ B4 M$ k: Z5 g, c# A* Z, \1 k                k=i;
                for(j=i+1;j<=n-1;j++)
                        if(fabs(*(c+j*(n+1)+i))>(fabs(*(c+k*(n+1)+i))))
% l! H- x/ c1 ^5 J% ?# f2 G: E                                k=j;
                        if(k!=i)
9 X& _+ @7 ^! j                                for( j=i;j<=n;j++ )
9 _4 i$ y+ x1 p8 m  c                                {
1 s) a% ]1 l  m$ q0 ^                                        p=*(c+i*(n+1)+j);
- y2 j# {8 a7 {' e3 l( I9 ]                                        *(c+i*(n+1)+j)=*(c+k*(n+1)+j);
" N& I; `6 V$ [2 f8 w                                        *(c+k*(n+1)+j)=p;
+ Z, r: L6 n% T" w, b/ `5 H                                }
                for( j=i+1;j<=n-1;j++ )
                {
                        p=(*(c+j*(n+1)+i))/(*(c+i*(n+1)+i));
  m/ J# b) l1 U3 [$ o) d                        for( t=i;t<=n;t++ )
2 A4 u1 G( n- R3 T3 q: x                                *(c+j*(n+1)+t)-=p*(*(c+i*(n+1)+t));
                }
        }
+ W. [3 D8 p1 T# P8 w5 y; S( a6 w        for( i=n-1;i>=0;i--)
        {
, _% Z1 L5 ?" Y1 g                for( j=n-1;j>=i+1;j--)
                        (*(c+i*(n+1)+n))-=x[j]*(*(c+i*(n+1)+j));
$ y4 k: W0 I  Y( g. W                x[i]=*(c+i*(n+1)+n)/(*(c+i*(n+1)+i));
        }
" J* P: x$ x  l# e4 ~        return x;
7 W) }) o2 P' W. r$ E! v}

---

作者: xiaozhangwh    时间: 2009-7-1 21:56
谢谢分享，来试一下

---

作者: deven1985    时间: 2010-1-1 09:06
，不错，前几天的数值分析课上也是做的这个！

---

作者: Rekcahpu    时间: 2010-11-23 13:44
刚学高斯消元

---

作者: 0.9清1.8清2.7清    时间: 2012-8-31 16:04
wuyuwuyul 自己程序拿出来看看呗！干嘛还要体力A！

---

作者: 弘道    时间: 2014-7-29 09:12
谢谢楼主……辛苦啦！………………

---

欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)

Powered by Discuz! X2.5