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

hduscimc

#include<iostream>
#include<cmath>
using namespace std;
4 _! \+ ]* \0 O4 i2 lint main()
  v0 P- I) p$ h! ?6 ~6 u6 v{
7 K* z. b# ^. F+ F7 X        int i;
3 Q8 R* G0 n! F        float *x;
% c% i2 }" l! F; {$ U1 g: t- @3 ?; D        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);
        x=ColPivot(c[0],3);
+ Z( Z7 Z+ l' I. W; h        for( i=0;i<=2;i++ )
                cout<<"x("<<i<<")="<<x[i]<<endl;
}

float *ColPivot( float *c,int n )
{
/ N4 n: F# S* S4 S. C& L8 M        int i,j,t,k;
& o& _4 Q  t9 K) r- n  |        float *x,p;
4 e5 v8 p  i  {7 S         x=new float[n*sizeof(float)];
2 N  j+ \4 U2 ]+ H  Q. @" ~$ C        for( i=0;i<=n-2;i++)
8 F; N& q1 `* r4 y4 U        {
; X# v- ?( s5 \! [; a& 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;
                        if(k!=i)
& A9 {8 `) w2 }( S' H0 b/ b                                for( j=i;j<=n;j++ )
7 b7 o( ~4 f2 j9 u6 u* s- r" L                                {
  V3 v& _3 X! Q  C% o3 Z* E                                        p=*(c+i*(n+1)+j);
                                        *(c+i*(n+1)+j)=*(c+k*(n+1)+j);
4 R* J3 U  }$ B4 f& R                                        *(c+k*(n+1)+j)=p;
% h% c2 Y' Z3 b7 P" |) {                                }
$ i- O6 X, J5 F! d3 e: s5 k                for( j=i+1;j<=n-1;j++ )
                {
. X; ]  E& q9 ^+ d6 j                        p=(*(c+j*(n+1)+i))/(*(c+i*(n+1)+i));
# V/ }; s0 w6 Q& h9 W0 X) e                        for( t=i;t<=n;t++ )
                                *(c+j*(n+1)+t)-=p*(*(c+i*(n+1)+t));
- O5 |, Q" _6 Y8 O: P                }
) z4 M$ I. r* z; O* x        }
: Z& R( k8 b' s7 [& f        for( i=n-1;i>=0;i--)
9 {+ e* {* G% @        {
                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));
* o* Q2 b( ]/ d4 N8 j7 j! D! I        }
        return x;
}

xiaozhangwh

谢谢分享，来试一下

deven1985

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

Rekcahpu

刚学高斯消元

0.9清1.8清2.7清

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

弘道

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