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

hduscimc

#include<iostream>
#include<cmath>
) G) P& |" B0 V( u! I" z$ j* Tusing namespace std;
& k- r5 J4 s4 eint main()
{
        int i;
# L4 Y5 k) P) z- m; d4 I$ u: }        float *x;
7 P: G, u$ l- 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);
# {8 e% \. c5 y. d3 W9 S2 u3 R        x=ColPivot(c[0],3);
8 U0 X# C2 t" X! d. K        for( i=0;i<=2;i++ )
, G1 {4 s8 C# [                cout<<"x("<<i<<")="<<x[i]<<endl;
& B& Q, }2 s' S7 \& {6 \. d+ I}

8 N. T% ?- f& I4 Bfloat *ColPivot( float *c,int n )
7 r' k5 j2 C5 X" ^8 n3 s{
$ N9 k( P( r4 C' @0 i" Q9 A        int i,j,t,k;
        float *x,p;
$ ^8 e3 p  O. }         x=new float[n*sizeof(float)];
7 w! [) H; p$ ^% e, M- k. @( o! T        for( i=0;i<=n-2;i++)
! @9 @; ?6 o; B  R0 E- `: u! J. J$ t, q        {
                k=i;
1 _/ _' G/ l% o4 ?( w, e  d* P                for(j=i+1;j<=n-1;j++)
2 `9 ~1 ^( `. D0 _! k+ @                        if(fabs(*(c+j*(n+1)+i))>(fabs(*(c+k*(n+1)+i))))
                                k=j;
                        if(k!=i)
# o; T/ w" V8 ~* Y; h$ G, A- p0 U                                for( j=i;j<=n;j++ )
                                {
+ B, s8 _% u/ e/ ?                                        p=*(c+i*(n+1)+j);
                                        *(c+i*(n+1)+j)=*(c+k*(n+1)+j);
                                        *(c+k*(n+1)+j)=p;
% _4 d) `- j* S                                }
' j/ Q3 h9 h( X2 c: b" q# O6 l& h                for( j=i+1;j<=n-1;j++ )
+ v! r5 d% e' n- @% j7 @                {
8 ]( }. |. {, C( R                        p=(*(c+j*(n+1)+i))/(*(c+i*(n+1)+i));
8 h. Q" s  Z# b: `                        for( t=i;t<=n;t++ )
                                *(c+j*(n+1)+t)-=p*(*(c+i*(n+1)+t));
/ T- {( y  d# D                }
/ V! a! ^) O) F& i        }
  Z* x& s) M  t( t' o4 a        for( i=n-1;i>=0;i--)
; a1 Z$ E; h( B. u        {
                for( j=n-1;j>=i+1;j--)
( |- [* m) W; S7 w& W( b                        (*(c+i*(n+1)+n))-=x[j]*(*(c+i*(n+1)+j));
                x[i]=*(c+i*(n+1)+n)/(*(c+i*(n+1)+i));
( c2 A" h! x3 {! c2 k2 }2 i        }
, Y+ l: Z. ?: L, ?$ ]        return x;
" i7 X" a3 ~) m2 I5 Q* M}

xiaozhangwh

谢谢分享，来试一下

deven1985

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

Rekcahpu

刚学高斯消元

0.9清1.8清2.7清

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

弘道

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