数学建模社区-数学中国

标题: [讨论]高斯消元法的程序对比 [打印本页]
作者: god    时间: 2005-1-19 17:01
标题: [讨论]高斯消元法的程序对比
下面我收集的高斯消元法的程序。还是请大家看看，并提点意见。
1 j. U$ s' d/ N# M+ o: M1 e[attach]966[/attach]

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

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

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

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

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

[讨论]高斯消元法的程序对比
作者: chzjao    时间: 2005-9-24 11:45
标题: 发下我编的,请各位帮忙修改下
[attach]2037[/attach]
6 P0 C* H. \. q: I! o
[此贴子已经被作者于2005-9-24 11:46:03编辑过]

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

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>
* P) D1 p) e8 C#include<cmath>" |1 \6 s7 r1 t$ K$ {! S& _
using namespace std;
7 ]9 b6 x! t# V' c( Fint main()
0 I0 b; u5 ^4 M. C{
# r5 g  v9 l( [. s, U: c        int i;, J1 w2 V- }+ k3 v
        float *x;
* D! }6 s1 C0 n5 _3 ~        float c[3][4] ={0.101,2.304,3.555,1.183,) x' d/ T% h' o- z0 S$ W
                                        -1.347,3.712,4.623,2.137,
* P" L$ D# R& m; C5 O  W$ a                                        -2.835,1.072,5.643,3.035};. M3 D' L& E% \8 R5 ^+ s) p& q
        float *ColPivot(float *,int);
; x7 @3 O+ C0 J& X        x=ColPivot(c[0],3);2 c, U1 R8 p& C. k- r' D3 y
        for( i=0;i<=2;i++ )
5 m; k1 L( ], b) s0 h! N; [- h                cout<<"x("<<i<<")="<<x[i]<<endl;
# B& N) c8 X2 d$ E. g8 v! i}
7 h+ m8 T1 ^  k$ q. [, g2 x: I" D. O9 e. C
float *ColPivot( float *c,int n )/ V0 z- h, i6 X3 F
{4 X, V) W/ }2 u# R+ Z
        int i,j,t,k;( v+ C$ g: C+ }
        float *x,p;2 w5 D6 q% N3 v  v9 f
        x=new float[n*sizeof(float)];
+ B# ^! b9 J$ z( D% g8 K; D        for( i=0;i<=n-2;i++)
* h- i& G2 U$ H  [) I0 n3 A7 b6 S        {- s$ }9 }4 F; T' w. R- d% \& V
                k=i;
) y! O6 A6 T/ ^                for(j=i+1;j<=n-1;j++)
9 _2 b# f2 `5 T& I! u/ ~( D" k1 M                        if(fabs(*(c+j*(n+1)+i))>(fabs(*(c+k*(n+1)+i))))
: f% v6 p# K" ?8 r                                k=j;
7 q9 T! c3 H& L* i                        if(k!=i)7 E( X' Q. ~/ j3 q
                                for( j=i;j<=n;j++ )
& W" I0 g& q7 A, ]% l7 W9 I( [                                {
4 R; l7 j+ {3 y7 F: E9 E7 v$ p) H5 i                                        p=*(c+i*(n+1)+j);
/ {# k8 Z6 o5 c& X  z" V% u3 C+ f) V4 v                                        *(c+i*(n+1)+j)=*(c+k*(n+1)+j);* l7 t% F' A2 ~! o) m8 G! c1 |
                                        *(c+k*(n+1)+j)=p;7 t% d" U& y! Z* @( v4 ]- Q, C
                                }
5 H1 T! b7 v$ \& E3 E9 }7 s; W  D                for( j=i+1;j<=n-1;j++ )
0 n' P) ?" [5 w1 [                {
' v& [9 y6 c$ _0 q0 h                        p=(*(c+j*(n+1)+i))/(*(c+i*(n+1)+i));' \- \5 w2 Y$ f
                        for( t=i;t<=n;t++ )' P& X/ W" Y- D- [6 J8 h  o
                                *(c+j*(n+1)+t)-=p*(*(c+i*(n+1)+t));
) @/ |/ p' g! f+ k, s$ ?( s                }
1 Z! v% _: f8 ]9 G* N        }1 ?7 P, @. D3 \1 q9 _
        for( i=n-1;i>=0;i--)5 U+ f0 D4 U# J8 M0 n5 H
        {
2 ~7 S2 [) \% C; Z2 m3 a! W                for( j=n-1;j>=i+1;j--)
5 b4 H3 u8 v2 y& B) z1 L" i                        (*(c+i*(n+1)+n))-=x[j]*(*(c+i*(n+1)+j));
8 A; k, `% ?2 M6 [                x[i]=*(c+i*(n+1)+n)/(*(c+i*(n+1)+i));
. P, w1 c* u: x; A5 Z        }# A) [; f5 @. r! ]6 V7 y
        return x;0 [9 T$ D: p- V  u% ?4 u
}
作者: 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