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

hduscimc

#include<iostream>
#include<cmath>
! [4 k3 G& w3 F) m$ jusing namespace std;
  f1 m9 ?) C1 X" y7 P- Fint main()
$ }3 M3 q! [& r& [' r{
8 X5 I5 n! K: g* s        int i;
        float *x;
        float c[3][4] ={0.101,2.304,3.555,1.183,
                                        -1.347,3.712,4.623,2.137,
* p  i, v0 O# M$ E5 I2 h                                        -2.835,1.072,5.643,3.035};
        float *ColPivot(float *,int);
6 j: K% {$ J  Y7 P6 s" J- X        x=ColPivot(c[0],3);
4 G! w3 R, v3 ^: X9 @        for( i=0;i<=2;i++ )
                cout<<"x("<<i<<")="<<x[i]<<endl;
% V% S4 x5 ]) u) n! i! l/ i: V8 ~+ |}
4 A' ?9 n# h; Y! W/ |5 W
float *ColPivot( float *c,int n )
8 G5 x' D$ A! x5 }  G+ C8 q% P  a{
        int i,j,t,k;
8 k, T( ^0 Q' D5 T: c        float *x,p;
; j1 p: [& o$ y3 m5 a         x=new float[n*sizeof(float)];
- @, k6 @7 T1 Y        for( i=0;i<=n-2;i++)
        {
9 o) u0 g# C7 `  W! a                k=i;
: K% n4 S1 X9 S, O5 c5 s                for(j=i+1;j<=n-1;j++)
( L/ c3 {7 S6 ]2 \% `                        if(fabs(*(c+j*(n+1)+i))>(fabs(*(c+k*(n+1)+i))))
1 P6 i8 y5 D; m" g; W5 t' R                                k=j;
                        if(k!=i)
5 I/ a6 q; Y$ L) y                                for( j=i;j<=n;j++ )
9 a9 B5 z5 f0 a: ~) g. u: o) ^* ^                                {
                                        p=*(c+i*(n+1)+j);
                                        *(c+i*(n+1)+j)=*(c+k*(n+1)+j);
                                        *(c+k*(n+1)+j)=p;
& H5 I6 j4 `: I                                }
                for( j=i+1;j<=n-1;j++ )
                {
+ o& B" a/ u9 j                        p=(*(c+j*(n+1)+i))/(*(c+i*(n+1)+i));
( q8 J: w% ^' j) H  S! |# Y& w. C+ P4 W                        for( t=i;t<=n;t++ )
! |3 n7 W" b7 ^, v* @6 E! ~; k                                *(c+j*(n+1)+t)-=p*(*(c+i*(n+1)+t));
                }
        }
        for( i=n-1;i>=0;i--)
        {
3 o! `* z# l6 E5 \& P8 C                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));
        }
        return x;
' q4 q  B. J0 d0 v  g9 S}

xiaozhangwh

谢谢分享，来试一下

deven1985

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

Rekcahpu

刚学高斯消元

0.9清1.8清2.7清

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

弘道

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