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

hduscimc

#include<iostream>
, Y% e7 _$ h# I1 |( ^#include<cmath>
$ [! K7 K$ e6 n: s3 l) w: Pusing namespace std;
: K4 X; b! r& J5 ]int main()
{
        int i;
/ z& c, z- M5 H% m; L        float *x;
        float c[3][4] ={0.101,2.304,3.555,1.183,
/ m: i) p( \/ r2 G3 M                                        -1.347,3.712,4.623,2.137,
7 E2 r6 @; ]3 R& R. t1 Y3 b                                        -2.835,1.072,5.643,3.035};
6 ~, r) V; Q  m  x4 m        float *ColPivot(float *,int);
        x=ColPivot(c[0],3);
        for( i=0;i<=2;i++ )
  z- \8 U# T  }5 S% m                cout<<"x("<<i<<")="<<x[i]<<endl;
" @/ {5 e4 Z" A8 B}
  l) T/ g' L; T& l3 o2 Z' c! k
float *ColPivot( float *c,int n )
{
        int i,j,t,k;
        float *x,p;
        x=new float[n*sizeof(float)];
        for( i=0;i<=n-2;i++)
        {
6 F1 e: T: Z7 Z) m5 G                k=i;
$ o& ?/ \  |; }8 @( `3 g% I                for(j=i+1;j<=n-1;j++)
                        if(fabs(*(c+j*(n+1)+i))>(fabs(*(c+k*(n+1)+i))))
                                k=j;
                        if(k!=i)
& }% @6 {+ g8 U1 e- ], Z* R                                for( j=i;j<=n;j++ )
: g2 p4 D  T! x) n/ d5 P4 U* c                                {
                                        p=*(c+i*(n+1)+j);
                                        *(c+i*(n+1)+j)=*(c+k*(n+1)+j);
0 v$ \/ l* W! n0 k5 b5 E- m                                        *(c+k*(n+1)+j)=p;
                                }
                for( j=i+1;j<=n-1;j++ )
                {
                        p=(*(c+j*(n+1)+i))/(*(c+i*(n+1)+i));
                        for( t=i;t<=n;t++ )
7 f* x/ j8 P: H) K3 }- V                                *(c+j*(n+1)+t)-=p*(*(c+i*(n+1)+t));
                }
; }. w2 v! J/ k% `        }
  v& G$ \0 i/ b3 ^6 x) I! M        for( i=n-1;i>=0;i--)
' R0 ?( u$ |" ~. i! w% ^4 u        {
                for( j=n-1;j>=i+1;j--)
0 b5 A4 b5 e/ ?1 y                        (*(c+i*(n+1)+n))-=x[j]*(*(c+i*(n+1)+j));
% {4 T; l$ I: h/ h# P* q& y                x[i]=*(c+i*(n+1)+n)/(*(c+i*(n+1)+i));
3 [  H* b) K" `. Z. Z  M& x/ o        }
1 {, P6 K+ O+ m. d* K6 U' [        return x;
6 f/ |1 i8 h" ~}

xiaozhangwh

谢谢分享，来试一下

deven1985

，不错，前几天的数值分析课上也是做的这个！

Rekcahpu

刚学高斯消元

0.9清1.8清2.7清

wuyuwuyul 自己程序拿出来看看呗！干嘛还要体力A！

弘道

谢谢楼主……辛苦啦！………………