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

hduscimc

#include<iostream>
#include<cmath>
& \3 Q. ]" P0 s4 k0 U% v3 _% ^using namespace std;
int main()
{
        int i;
        float *x;
  Z% K: b7 z  P: f+ u        float c[3][4] ={0.101,2.304,3.555,1.183,
% M8 Y3 l& }! \: p2 X# j                                        -1.347,3.712,4.623,2.137,
, [' c# y; d: n# h) l                                        -2.835,1.072,5.643,3.035};
( P9 b4 m5 z' h& Q        float *ColPivot(float *,int);
* d! k2 e2 ?, c" D' N        x=ColPivot(c[0],3);
        for( i=0;i<=2;i++ )
                cout<<"x("<<i<<")="<<x[i]<<endl;
}

float *ColPivot( float *c,int n )
+ v8 n7 u- M" d" J- A{
        int i,j,t,k;
        float *x,p;
        x=new float[n*sizeof(float)];
        for( i=0;i<=n-2;i++)
        {
                k=i;
( B# V4 d1 \; r& h; k. \, ]. y, p                for(j=i+1;j<=n-1;j++)
& _# }. r6 D) l6 [                        if(fabs(*(c+j*(n+1)+i))>(fabs(*(c+k*(n+1)+i))))
                                k=j;
                        if(k!=i)
                                for( j=i;j<=n;j++ )
                                {
                                        p=*(c+i*(n+1)+j);
# ~8 n5 a8 N( D, ~                                        *(c+i*(n+1)+j)=*(c+k*(n+1)+j);
                                        *(c+k*(n+1)+j)=p;
: ~; f. y3 K3 j# I& B, }/ e# X                                }
1 U0 h3 _1 Y, G. f& B6 J7 M0 ~                for( j=i+1;j<=n-1;j++ )
                {
                        p=(*(c+j*(n+1)+i))/(*(c+i*(n+1)+i));
! ], j3 N% d8 E+ a                        for( t=i;t<=n;t++ )
                                *(c+j*(n+1)+t)-=p*(*(c+i*(n+1)+t));
9 g9 D1 D2 c: |5 W4 o; n                }
        }
0 s8 W( ?1 F, w  S& e        for( i=n-1;i>=0;i--)
        {
+ p* Z/ k9 e/ |9 n; U& I1 M                for( j=n-1;j>=i+1;j--)
, l, J: U4 s: y" v  K7 @+ D                        (*(c+i*(n+1)+n))-=x[j]*(*(c+i*(n+1)+j));
                x[i]=*(c+i*(n+1)+n)/(*(c+i*(n+1)+i));
4 H! e( e- h" Q$ f  M        }
        return x;
}

xiaozhangwh

谢谢分享，来试一下

deven1985

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

Rekcahpu

刚学高斯消元

0.9清1.8清2.7清

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

弘道

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