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

hduscimc

#include<iostream>
! n3 I9 p# P0 j7 Z+ b#include<cmath>
using namespace std;
4 @  o8 p/ p& [* ?( ~8 oint main()
% @, u# \! K" w* t" C{
        int i;
; r% D, J4 ^3 M* L" x9 v        float *x;
        float c[3][4] ={0.101,2.304,3.555,1.183,
5 y6 _4 m) Y1 x* J; n$ Q                                        -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);
        for( i=0;i<=2;i++ )
3 U! T2 L3 }9 ?5 f. `                cout<<"x("<<i<<")="<<x[i]<<endl;
}
/ r/ X- A- w) L- a1 J' E+ Y
8 X  G) Z4 k0 ^( k3 Mfloat *ColPivot( float *c,int n )
+ D5 G0 I+ t: E' o# O  W{
) Q# i4 G; i! _5 \8 X9 E2 I        int i,j,t,k;
        float *x,p;
        x=new float[n*sizeof(float)];
        for( i=0;i<=n-2;i++)
, U( {; h1 o) l8 g' I) Q+ d        {
, @, {' q$ T+ M' f8 H& ?# r                k=i;
                for(j=i+1;j<=n-1;j++)
                        if(fabs(*(c+j*(n+1)+i))>(fabs(*(c+k*(n+1)+i))))
7 `" |1 W# L2 q; F: X3 r% K                                k=j;
                        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;
& R5 I7 I/ H; v! t! o% u                                }
+ L; g6 w  n+ f& v                for( j=i+1;j<=n-1;j++ )
( k; Z, u/ K* y( }, i4 o& K) H                {
8 d8 C7 p9 y. x+ x" R4 s                        p=(*(c+j*(n+1)+i))/(*(c+i*(n+1)+i));
/ G! k2 n8 p4 s6 N                        for( t=i;t<=n;t++ )
7 d* ~* \% e* }$ d& J                                *(c+j*(n+1)+t)-=p*(*(c+i*(n+1)+t));
- I7 X; Z4 V' j9 ]                }
  K  r/ H; b, E4 \9 A        }
        for( i=n-1;i>=0;i--)
        {
                for( j=n-1;j>=i+1;j--)
) Z" u) B- [; r- ~                        (*(c+i*(n+1)+n))-=x[j]*(*(c+i*(n+1)+j));
+ b+ |* O" T. V6 H: n: `; e3 C                x[i]=*(c+i*(n+1)+n)/(*(c+i*(n+1)+i));
" f8 a0 x7 a: p- T. r        }
        return x;
' e) P+ a: F/ R. x7 R6 U8 l5 {# K}

xiaozhangwh

谢谢分享，来试一下

deven1985

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

Rekcahpu

刚学高斯消元

0.9清1.8清2.7清

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

弘道

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