数学建模社区-数学中国

标题: [讨论]高斯消元法的程序对比 [打印本页]
作者: god    时间: 2005-1-19 17:01
标题: [讨论]高斯消元法的程序对比
下面我收集的高斯消元法的程序。还是请大家看看，并提点意见。
8 w6 T! }& @: M* j$ L[attach]966[/attach]

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

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

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

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

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

[讨论]高斯消元法的程序对比
作者: chzjao    时间: 2005-9-24 11:45
标题: 发下我编的,请各位帮忙修改下
[attach]2037[/attach]
' f5 n. m6 \! Q
[此贴子已经被作者于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>& t& L3 X: k4 u. _% H
#include<cmath>
$ H7 }& [  O7 q; M. }! E- Ausing namespace std;8 d0 B8 L! Y0 T
int main()
; M3 D1 {1 L4 J{+ D, {6 {+ s  i: h/ J  K
        int i;% I% ?" i- b% v. w+ H) |2 m
        float *x;
: T: B  o& S+ ~/ B4 X  C) L        float c[3][4] ={0.101,2.304,3.555,1.183,6 G6 L1 R7 B+ d" _
                                        -1.347,3.712,4.623,2.137," t9 j  l$ x" x
                                        -2.835,1.072,5.643,3.035};5 ?( ?) G% c& E* \8 \
        float *ColPivot(float *,int);5 m: V* w0 d& C1 u0 _
        x=ColPivot(c[0],3);
! _- ~7 h8 x8 P/ d4 s, n0 i0 F% d        for( i=0;i<=2;i++ )
) t* f, v5 H3 z5 g7 d/ r+ c                cout<<"x("<<i<<")="<<x[i]<<endl;
% P% W/ C2 c% q1 H7 _1 D}' t# s! {- p, K

, J9 B' A& e+ G& f9 Vfloat *ColPivot( float *c,int n )3 o' i1 T" K8 y) V/ U% {) F1 R
{
+ L, I$ Z4 ^0 J) ], m$ s& H        int i,j,t,k;
- l! ?$ e2 z9 R/ {) |        float *x,p;
5 Y7 V0 Y  S( R# T& s6 E, e) N         x=new float[n*sizeof(float)];3 a3 o% Z2 p; s. f; h- p3 P: x
        for( i=0;i<=n-2;i++)
' q( n7 e+ t" G( K9 x' S0 U        {
& w2 v  b# u% C                k=i;% M& L3 z: X) m
                for(j=i+1;j<=n-1;j++)
1 D; C& V% M& s$ A# ^2 e                        if(fabs(*(c+j*(n+1)+i))>(fabs(*(c+k*(n+1)+i))))' T+ D, l9 ~' L: T+ e( t
                                k=j;0 q1 G- x# h5 Q; D" B; d
                        if(k!=i)6 H( |6 e( H9 j0 x4 u
                                for( j=i;j<=n;j++ )
1 e/ M7 z) F* [9 f) V                                {
  n. q$ z* v& g                                        p=*(c+i*(n+1)+j);, H. }3 C8 t5 E; C
                                        *(c+i*(n+1)+j)=*(c+k*(n+1)+j);. Y) Y! {6 @' F  G
                                        *(c+k*(n+1)+j)=p;
  [3 C1 e! ^6 g) E0 l' k                                }
( e5 W. r' `9 s                for( j=i+1;j<=n-1;j++ )3 f7 {) T. m( j+ l$ e. _3 t
                {; Q: C. b, U; H7 z( K8 H  X3 W
                        p=(*(c+j*(n+1)+i))/(*(c+i*(n+1)+i));
1 w+ w; c5 j) y                        for( t=i;t<=n;t++ )
& R  S, J) c4 h- @+ d                                *(c+j*(n+1)+t)-=p*(*(c+i*(n+1)+t));
0 P; t  w$ I: _5 M9 [                }6 O: K/ P2 J  K: q1 |  ]
        }  O1 u3 W$ \/ l. o2 z
        for( i=n-1;i>=0;i--)
% O2 j. V5 f  |$ }7 D7 m        {
1 j1 G* z5 X* ?& q" b5 b+ U5 y                for( j=n-1;j>=i+1;j--)
+ c9 z. Q5 \& W, b4 ^* {2 b                        (*(c+i*(n+1)+n))-=x[j]*(*(c+i*(n+1)+j));
' M. f- O; `8 [. }: z1 _                x[i]=*(c+i*(n+1)+n)/(*(c+i*(n+1)+i));" I/ W1 N6 v6 }+ {  W
        }5 z4 W* Z4 o! }1 V% m
        return x;" @: j7 C" A3 e3 v- e) L; c  _
}
作者: 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