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

hduscimc

#include<iostream>
8 T9 N( t- U6 S1 ?) t#include<cmath>
/ `  f  Y& B1 r: f$ K  zusing namespace std;
% k& ?: `9 F1 o. r: cint main()
! {* ~6 E  y. v5 p( X{
/ H3 M! F8 H% x& c2 Y6 [/ k        int i;
9 j! P! i4 e6 z6 b* i& c        float *x;
        float c[3][4] ={0.101,2.304,3.555,1.183,
                                        -1.347,3.712,4.623,2.137,
% G6 X6 d6 }& R6 R                                        -2.835,1.072,5.643,3.035};
        float *ColPivot(float *,int);
        x=ColPivot(c[0],3);
) t- `. ~% _( }! t; ]* R/ w  N! t. O% m        for( i=0;i<=2;i++ )
# R4 y* M! c0 a# o2 R                cout<<"x("<<i<<")="<<x[i]<<endl;
* G; C* A; b: T! P% C" C0 X$ W}

& p" D5 R  I. A( r( s8 Q0 t' t* w( @float *ColPivot( float *c,int n )
0 H0 z* M7 I! f{
- j; w& e2 ^6 d, w        int i,j,t,k;
        float *x,p;
        x=new float[n*sizeof(float)];
        for( i=0;i<=n-2;i++)
/ N1 ]9 Q  M2 K        {
  r7 d2 k& a) N( f7 V' J                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;
. ]1 V/ P) Y# I                        if(k!=i)
# T$ Z! }. Y; J1 l( O2 }                                for( j=i;j<=n;j++ )
  [( y9 ?! |4 S8 l0 s: ]7 L                                {
. X8 |6 Q# A. F+ g/ g, U: N                                        p=*(c+i*(n+1)+j);
                                        *(c+i*(n+1)+j)=*(c+k*(n+1)+j);
+ ]1 D6 T4 _7 J! P5 @. R# c$ |) T                                        *(c+k*(n+1)+j)=p;
                                }
  O* B" J- D9 c# j- j# ^$ ]                for( j=i+1;j<=n-1;j++ )
                {
* n! t. y0 j' U  ~* F- d4 n8 T                        p=(*(c+j*(n+1)+i))/(*(c+i*(n+1)+i));
                        for( t=i;t<=n;t++ )
# C7 X2 x) G  `6 f9 ?3 g. Y                                *(c+j*(n+1)+t)-=p*(*(c+i*(n+1)+t));
, X- V" M/ [, D( A8 a, v                }
9 H8 P+ k& g6 f) g        }
        for( i=n-1;i>=0;i--)
        {
* z( ^' }! J6 T- q  n$ \                for( j=n-1;j>=i+1;j--)
$ U5 T8 _; j+ l' O" |) p. a7 l$ T                        (*(c+i*(n+1)+n))-=x[j]*(*(c+i*(n+1)+j));
/ ~: f" v7 x9 t9 @; ?. C( B+ t8 g                x[i]=*(c+i*(n+1)+n)/(*(c+i*(n+1)+i));
        }
+ ]' r0 b9 M& d# O1 B        return x;
) n4 S9 }6 O: f2 P3 ^}

xiaozhangwh

谢谢分享，来试一下

deven1985

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

Rekcahpu

刚学高斯消元

0.9清1.8清2.7清

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

弘道

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