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

hduscimc

#include<iostream>
+ |$ B) a2 K/ T' Y) |6 A) w, `3 }#include<cmath>
4 \9 ?# \2 C# eusing namespace std;
( ^4 t+ R+ f3 X2 i5 b1 U- p) S# c! Nint main()
1 c& F* v2 E* i  L, I0 w8 g* S{
0 v# e& ?# x3 `. d4 h6 A        int i;
        float *x;
  g" C- l1 b5 M; O% X. M8 K6 x        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};
/ o$ _( }* l3 g, A        float *ColPivot(float *,int);
9 q; R- J  M4 v: s- U- D        x=ColPivot(c[0],3);
7 [% h2 W- q3 ~; Y: g; \2 Z( j, l        for( i=0;i<=2;i++ )
5 Z# G/ ?( m& p+ ]# ?                cout<<"x("<<i<<")="<<x[i]<<endl;
% n9 M* E/ t' A5 ^% d# x}
1 n, s+ a6 ^( j- P# z% h$ F
float *ColPivot( float *c,int n )
{
        int i,j,t,k;
7 p5 k3 s9 G2 f( B        float *x,p;
        x=new float[n*sizeof(float)];
        for( i=0;i<=n-2;i++)
        {
                k=i;
1 n, }9 B. b! ^  e% r" p9 k                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)
                                for( j=i;j<=n;j++ )
% e0 a0 @1 G; A' U8 \" ^- X                                {
                                        p=*(c+i*(n+1)+j);
& ?; l, [: @/ b4 w( _, y/ B' |. H( K) ]                                        *(c+i*(n+1)+j)=*(c+k*(n+1)+j);
                                        *(c+k*(n+1)+j)=p;
                                }
                for( j=i+1;j<=n-1;j++ )
4 s, @7 c/ I, W; G                {
8 S7 m5 G5 v% ~, i2 f, }0 R                        p=(*(c+j*(n+1)+i))/(*(c+i*(n+1)+i));
                        for( t=i;t<=n;t++ )
                                *(c+j*(n+1)+t)-=p*(*(c+i*(n+1)+t));
                }
. ]: j4 X7 U" P; L% \  n( f        }
        for( i=n-1;i>=0;i--)
        {
' \" a' R3 a3 l& c* C                for( j=n-1;j>=i+1;j--)
! V2 y4 C: F7 F7 P                        (*(c+i*(n+1)+n))-=x[j]*(*(c+i*(n+1)+j));
                x[i]=*(c+i*(n+1)+n)/(*(c+i*(n+1)+i));
        }
2 ]1 @( k, Y  u& ^        return x;
}

xiaozhangwh

谢谢分享，来试一下

deven1985

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

Rekcahpu

刚学高斯消元

0.9清1.8清2.7清

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

弘道

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