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

hduscimc

#include<iostream>
7 F( ?# T  v% `* G#include<cmath>
. C9 f/ E. E2 T" Musing namespace std;
int main()
- Q" T9 p& P, i" Y" S{
) @! O8 F- j$ l; {. e6 n        int i;
        float *x;
        float c[3][4] ={0.101,2.304,3.555,1.183,
  M% u2 z" Q7 X& V+ W, v                                        -1.347,3.712,4.623,2.137,
. Y* I; r. ~5 D7 c                                        -2.835,1.072,5.643,3.035};
$ E) D& T7 S4 t! t# ~/ n        float *ColPivot(float *,int);
9 Q: r  f0 t# ^6 n) W        x=ColPivot(c[0],3);
) t4 M0 b( g/ R        for( i=0;i<=2;i++ )
                cout<<"x("<<i<<")="<<x[i]<<endl;
' p6 V7 G' E8 H: V, M}

5 _1 M$ f( k# l5 F5 ?! P7 @# Hfloat *ColPivot( float *c,int n )
' c. E/ d& y# g) Y6 Y3 s5 J+ I{
; J3 x$ [1 K8 a. A4 Q        int i,j,t,k;
$ w$ z8 I0 B9 G8 I2 Y$ j        float *x,p;
8 S  Z+ E8 T' Z6 ~! f/ N: F         x=new float[n*sizeof(float)];
  F# ], u  _, ?8 b0 O& R        for( i=0;i<=n-2;i++)
        {
  H! c! k5 N8 K/ |  m" j                k=i;
1 y+ g1 t& \# r' N6 r" w& Q                for(j=i+1;j<=n-1;j++)
( q* o" n1 a% \  H% s                        if(fabs(*(c+j*(n+1)+i))>(fabs(*(c+k*(n+1)+i))))
3 s6 s4 y5 m" i. e! s  w                                k=j;
" K1 E; K; n6 w7 ^7 Y% N' z                        if(k!=i)
                                for( j=i;j<=n;j++ )
0 t+ J1 }6 i/ W3 l9 q) e, l                                {
                                        p=*(c+i*(n+1)+j);
                                        *(c+i*(n+1)+j)=*(c+k*(n+1)+j);
9 N; r2 l7 ]3 E4 o                                        *(c+k*(n+1)+j)=p;
                                }
                for( j=i+1;j<=n-1;j++ )
5 t7 V, V( m, z5 {- X                {
3 C& s0 {, B/ g0 X7 t' ]( M                        p=(*(c+j*(n+1)+i))/(*(c+i*(n+1)+i));
! L. O+ j6 T' D                        for( t=i;t<=n;t++ )
9 I/ c+ `" T9 I( J7 ~/ t; A! Q0 U" K5 y                                *(c+j*(n+1)+t)-=p*(*(c+i*(n+1)+t));
$ u# J1 w1 B: S, b( ^# R                }
, S& Q& _, a8 w        }
* }$ S" @; [1 f9 \. W) J* S        for( i=n-1;i>=0;i--)
" u/ ]9 ]9 _$ F0 L; Q        {
3 e/ G( q7 y6 L8 \; f                for( j=n-1;j>=i+1;j--)
- n3 b' l# _5 \7 R                        (*(c+i*(n+1)+n))-=x[j]*(*(c+i*(n+1)+j));
, l) D+ p+ U& v7 }+ U* P7 ^% ^' E4 t$ }                x[i]=*(c+i*(n+1)+n)/(*(c+i*(n+1)+i));
# J; C4 E4 J+ j) J        }
4 \% ~  L; U. r% {4 `        return x;
# {4 c: K! _' N; c; P. T}

xiaozhangwh

谢谢分享，来试一下

deven1985

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

Rekcahpu

刚学高斯消元

0.9清1.8清2.7清

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

弘道

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