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

hduscimc

#include<iostream>
#include<cmath>
' _0 C! M" w! |: Y5 r; nusing namespace std;
( ?& K- s3 ?# B- ^/ k' R) O5 ]8 Cint main()
5 q! m- B5 ~5 h# C5 r{
        int i;
        float *x;
        float c[3][4] ={0.101,2.304,3.555,1.183,
                                        -1.347,3.712,4.623,2.137,
3 I; D  e* Z. O; R! d& Q* h6 ]" n                                        -2.835,1.072,5.643,3.035};
        float *ColPivot(float *,int);
        x=ColPivot(c[0],3);
        for( i=0;i<=2;i++ )
- ~- M( T$ _3 P4 w                cout<<"x("<<i<<")="<<x[i]<<endl;
# C( X+ d7 |8 O6 {2 B- _* a}

float *ColPivot( float *c,int n )
{
        int i,j,t,k;
+ e: c9 k/ x3 t! V: R        float *x,p;
7 B/ s' v  {0 H3 I6 X2 }, p" U+ v         x=new float[n*sizeof(float)];
& R9 z: o7 e0 b7 P: r        for( i=0;i<=n-2;i++)
        {
                k=i;
* E* Q4 f" e& w( O# N/ V                for(j=i+1;j<=n-1;j++)
                        if(fabs(*(c+j*(n+1)+i))>(fabs(*(c+k*(n+1)+i))))
$ O; w5 F' B4 _4 i. c                                k=j;
                        if(k!=i)
" ~6 P3 G4 _8 ^% X" T/ ~                                for( j=i;j<=n;j++ )
                                {
                                        p=*(c+i*(n+1)+j);
                                        *(c+i*(n+1)+j)=*(c+k*(n+1)+j);
7 i/ S8 ?8 W7 m, A8 ~                                        *(c+k*(n+1)+j)=p;
                                }
5 d  _" ^, T2 t+ ~* _5 C' |9 q1 @1 v3 n                for( j=i+1;j<=n-1;j++ )
                {
                        p=(*(c+j*(n+1)+i))/(*(c+i*(n+1)+i));
1 i1 i9 F  i" L# `5 ^3 G$ o                        for( t=i;t<=n;t++ )
                                *(c+j*(n+1)+t)-=p*(*(c+i*(n+1)+t));
1 c! u; f6 e6 h: T# ]9 z$ P                }
6 w: V* _: e) x1 U2 y) w; I        }
. R- W& ^# \3 \$ R        for( i=n-1;i>=0;i--)
: d9 y' u6 J) Y        {
2 d3 H' q4 Q  V! F# z5 X8 w% g                for( j=n-1;j>=i+1;j--)
                        (*(c+i*(n+1)+n))-=x[j]*(*(c+i*(n+1)+j));
  {! _/ g1 |/ k: D  @6 W8 D                x[i]=*(c+i*(n+1)+n)/(*(c+i*(n+1)+i));
        }
" o2 ?% y2 \$ X5 s$ p9 [/ J        return x;
4 {2 Q% Y$ n4 `' V: [' M}

xiaozhangwh

谢谢分享，来试一下

deven1985

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

Rekcahpu

刚学高斯消元

0.9清1.8清2.7清

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

弘道

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