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

hduscimc

#include<iostream>
! x7 j8 J5 S! V4 |5 w  v#include<cmath>
using namespace std;
* B% U! b: f8 fint main()
  \! `) J" X: a* f. @3 I{
% ]6 G% G6 i2 ^        int i;
! l( o( _- Z% m4 d        float *x;
        float c[3][4] ={0.101,2.304,3.555,1.183,
. A  H! g3 j5 a  ~3 p+ z3 X; U8 @                                        -1.347,3.712,4.623,2.137,
* e" ^1 Y4 W  w' U, r/ Y$ v" ^                                        -2.835,1.072,5.643,3.035};
        float *ColPivot(float *,int);
        x=ColPivot(c[0],3);
. E4 ]+ F+ h2 \, Z' N6 G5 v& ^        for( i=0;i<=2;i++ )
                cout<<"x("<<i<<")="<<x[i]<<endl;
; D+ Z+ x3 R! A( }" `% g}

float *ColPivot( float *c,int n )
{
4 y3 y5 L: R+ w: K2 \; R        int i,j,t,k;
        float *x,p;
1 H  [' I9 _; s' s, h( @         x=new float[n*sizeof(float)];
        for( i=0;i<=n-2;i++)
        {
                k=i;
7 w3 @) O- B: X' E6 _9 ^                for(j=i+1;j<=n-1;j++)
                        if(fabs(*(c+j*(n+1)+i))>(fabs(*(c+k*(n+1)+i))))
                                k=j;
9 f6 }6 P% `* w( h. I                        if(k!=i)
                                for( j=i;j<=n;j++ )
                                {
                                        p=*(c+i*(n+1)+j);
                                        *(c+i*(n+1)+j)=*(c+k*(n+1)+j);
                                        *(c+k*(n+1)+j)=p;
                                }
                for( j=i+1;j<=n-1;j++ )
! I4 y# H3 V( V. a/ ^0 z: L) g                {
                        p=(*(c+j*(n+1)+i))/(*(c+i*(n+1)+i));
/ \. I) H: j8 A9 ~' m2 T. V) x                        for( t=i;t<=n;t++ )
                                *(c+j*(n+1)+t)-=p*(*(c+i*(n+1)+t));
                }
        }
, j0 w! s: w! ?1 R. x        for( i=n-1;i>=0;i--)
; _) t  R* o/ _! Z/ m6 w7 X0 t        {
                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));
        }
3 d: c% u; F" R; m8 t- D8 n        return x;
}

xiaozhangwh

谢谢分享，来试一下

deven1985

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

Rekcahpu

刚学高斯消元

0.9清1.8清2.7清

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

弘道

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