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

hduscimc

#include<iostream>
#include<cmath>
) ]0 a9 v0 V4 k/ a$ h) Zusing namespace std;
' ?9 z) g" D, Q) q' |int main()
{
        int i;
* L9 |, Z) X) f4 W% i  n: E        float *x;
        float c[3][4] ={0.101,2.304,3.555,1.183,
                                        -1.347,3.712,4.623,2.137,
                                        -2.835,1.072,5.643,3.035};
        float *ColPivot(float *,int);
        x=ColPivot(c[0],3);
2 ^* L( I5 t! X2 _% d' [9 d        for( i=0;i<=2;i++ )
                cout<<"x("<<i<<")="<<x[i]<<endl;
}
9 p+ f( m. M( ]* l  _9 a7 }
float *ColPivot( float *c,int n )
{
4 P7 i- V+ g! o. G        int i,j,t,k;
* B* U! b. S0 ?        float *x,p;
        x=new float[n*sizeof(float)];
        for( i=0;i<=n-2;i++)
        {
                k=i;
* \8 c5 H, n' W9 F5 P* s. g+ g3 J5 D                for(j=i+1;j<=n-1;j++)
                        if(fabs(*(c+j*(n+1)+i))>(fabs(*(c+k*(n+1)+i))))
& i, {+ F8 L& g8 x: f/ w( V* {* U% N                                k=j;
                        if(k!=i)
                                for( j=i;j<=n;j++ )
; P; `0 p: Q: a1 B$ e                                {
                                        p=*(c+i*(n+1)+j);
                                        *(c+i*(n+1)+j)=*(c+k*(n+1)+j);
                                        *(c+k*(n+1)+j)=p;
                                }
9 A# w% @2 L/ m( v  B                for( j=i+1;j<=n-1;j++ )
. `  l# y$ j+ O* C0 Z                {
: b4 N1 ^) N+ g6 G% a! x                        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));
/ t% m( y  X) ?% H4 f                }
        }
  p3 V4 _- @9 m- O2 C/ A        for( i=n-1;i>=0;i--)
        {
: Q8 u$ h2 k! l/ Q* l8 r* n* T                for( j=n-1;j>=i+1;j--)
                        (*(c+i*(n+1)+n))-=x[j]*(*(c+i*(n+1)+j));
+ `9 O& P* z. U5 t2 \0 v4 ^/ ~8 ~                x[i]=*(c+i*(n+1)+n)/(*(c+i*(n+1)+i));
        }
        return x;
; t7 m# Q, i0 P" ~6 }9 d, @}

xiaozhangwh

谢谢分享，来试一下

deven1985

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

Rekcahpu

刚学高斯消元

0.9清1.8清2.7清

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

弘道

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