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

hduscimc

#include<iostream>
#include<cmath>
7 ?4 E7 m% f0 k6 |7 s8 [) C- H2 Husing namespace std;
int main()
7 s5 S5 A; W2 Q+ Z0 k  ^' S{
9 q  U$ M# j4 S& J8 ~9 C& N        int i;
. b3 W5 i* ]3 E$ S$ a        float *x;
+ W) [& P. r: f3 L/ W- l# f        float c[3][4] ={0.101,2.304,3.555,1.183,
                                        -1.347,3.712,4.623,2.137,
* x5 g! |: z! x: D, K# q3 ^                                        -2.835,1.072,5.643,3.035};
8 x+ U( c5 O5 V5 e$ e5 m        float *ColPivot(float *,int);
, u( i1 x6 C: J' D: v+ f        x=ColPivot(c[0],3);
        for( i=0;i<=2;i++ )
0 t. w- m, H+ m* |                cout<<"x("<<i<<")="<<x[i]<<endl;
2 \  K: X& i; R% l3 }}

" s: T* }' T) J7 S4 [! F; m; Cfloat *ColPivot( float *c,int n )
{
3 ]% S2 Y! g3 T9 F* r" E        int i,j,t,k;
        float *x,p;
        x=new float[n*sizeof(float)];
5 |, l0 ?. O' M: S+ u$ p, @        for( i=0;i<=n-2;i++)
        {
                k=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)
                                for( j=i;j<=n;j++ )
- j% A1 I. `: I$ M, Z) I: k5 \/ E                                {
0 Y8 E# Q8 ~7 c* `. W! @) I8 |                                        p=*(c+i*(n+1)+j);
                                        *(c+i*(n+1)+j)=*(c+k*(n+1)+j);
                                        *(c+k*(n+1)+j)=p;
- z' N& g) t2 I* |, `1 w                                }
7 l/ R- C9 o/ s3 m                for( j=i+1;j<=n-1;j++ )
                {
$ s( w' s3 D. T' a                        p=(*(c+j*(n+1)+i))/(*(c+i*(n+1)+i));
                        for( t=i;t<=n;t++ )
                                *(c+j*(n+1)+t)-=p*(*(c+i*(n+1)+t));
! E5 F; ~% A% M6 H0 E                }
. p+ B% _1 b/ L1 [. `  w! v        }
! o, P" K# A9 I; {% i  Y        for( i=n-1;i>=0;i--)
7 P' q% N9 P. Y2 C" N7 e0 e        {
                for( j=n-1;j>=i+1;j--)
                        (*(c+i*(n+1)+n))-=x[j]*(*(c+i*(n+1)+j));
                x[i]=*(c+i*(n+1)+n)/(*(c+i*(n+1)+i));
2 O. V* k% Y; {% {1 f7 B2 |. O1 N        }
1 L, i$ x8 t1 ?; q        return x;
}

xiaozhangwh

谢谢分享，来试一下

deven1985

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

Rekcahpu

刚学高斯消元

0.9清1.8清2.7清

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

弘道

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