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

hduscimc

#include<iostream>
1 A0 @) t: V& Z#include<cmath>
using namespace std;
( }) ?* d$ A8 f. g& tint main()
8 p; y1 b& L# y4 v1 Z{
        int i;
5 d! X5 J& h5 _- O( |3 D        float *x;
( b1 {: H6 ?- q" L2 s, e        float c[3][4] ={0.101,2.304,3.555,1.183,
, S$ U% A5 v$ R  ?* }; V  p" w                                        -1.347,3.712,4.623,2.137,
                                        -2.835,1.072,5.643,3.035};
! O: f3 h3 u8 J# x: z2 V6 ]9 P        float *ColPivot(float *,int);
, A, ]  ~2 J/ Q7 }5 B4 O        x=ColPivot(c[0],3);
& F% M& z/ o# @1 D        for( i=0;i<=2;i++ )
6 X8 o. _0 |$ V* [8 B/ \                cout<<"x("<<i<<")="<<x[i]<<endl;
) N% L/ I6 O' Z2 U: z}
' K1 M  }& E. y1 o! J% A/ D$ b
float *ColPivot( float *c,int n )
+ m# M. i9 l0 J  G, b{
        int i,j,t,k;
  y: K7 Y/ B) a! U        float *x,p;
        x=new float[n*sizeof(float)];
+ s* V0 W! U( V# [        for( i=0;i<=n-2;i++)
        {
! E8 G( a5 N% }7 b0 f5 ?! A                k=i;
                for(j=i+1;j<=n-1;j++)
* T$ s2 P6 {. ]                        if(fabs(*(c+j*(n+1)+i))>(fabs(*(c+k*(n+1)+i))))
" k# U" c( B6 _# Y- B                                k=j;
                        if(k!=i)
4 u7 `' W& [/ A                                for( j=i;j<=n;j++ )
                                {
                                        p=*(c+i*(n+1)+j);
$ \5 {9 P% E1 S6 y4 v                                        *(c+i*(n+1)+j)=*(c+k*(n+1)+j);
$ n# X% }& ?' E                                        *(c+k*(n+1)+j)=p;
/ T* a! L# }7 ~  c. R                                }
                for( j=i+1;j<=n-1;j++ )
8 C/ V" [6 I! e                {
                        p=(*(c+j*(n+1)+i))/(*(c+i*(n+1)+i));
; c, b! c9 ~" j                        for( t=i;t<=n;t++ )
5 s+ C' m2 ^  j                                *(c+j*(n+1)+t)-=p*(*(c+i*(n+1)+t));
                }
8 a4 q. l# P; U* K9 u        }
        for( i=n-1;i>=0;i--)
, ~9 v0 T& r( B. N        {
                for( j=n-1;j>=i+1;j--)
                        (*(c+i*(n+1)+n))-=x[j]*(*(c+i*(n+1)+j));
                x[i]=*(c+i*(n+1)+n)/(*(c+i*(n+1)+i));
        }
4 q9 R7 ]/ K: G) B  {2 e        return x;
7 ^3 L- s' s! }* N* }4 |}

xiaozhangwh

谢谢分享，来试一下

deven1985

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

Rekcahpu

刚学高斯消元

0.9清1.8清2.7清

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

弘道

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