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

hduscimc

#include<iostream>
#include<cmath>
- x- B6 v0 w; e7 p  Cusing namespace std;
int main()
{
; X) S% y, |% ^4 s5 M; s1 j        int i;
        float *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};
5 `8 a. ^6 u2 @" K2 M1 k  X$ z' u        float *ColPivot(float *,int);
        x=ColPivot(c[0],3);
        for( i=0;i<=2;i++ )
+ n/ U* X% c" J! Z0 a1 `7 E                cout<<"x("<<i<<")="<<x[i]<<endl;
+ `3 D: V$ j/ R6 ]0 m}

float *ColPivot( float *c,int n )
& ?8 ^1 w9 D' z2 ~6 i{
' }, \5 }, V7 m  D8 \8 r: S        int i,j,t,k;
0 a9 l9 g) `/ @3 X$ }        float *x,p;
        x=new float[n*sizeof(float)];
; j9 x8 R- F) A- h        for( i=0;i<=n-2;i++)
        {
                k=i;
% k8 r3 i) D7 T, M0 E$ m                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)
7 D% y: A& G% X$ j' h8 n# ~                                for( j=i;j<=n;j++ )
                                {
                                        p=*(c+i*(n+1)+j);
+ a) K/ {- m! q# Y. a                                        *(c+i*(n+1)+j)=*(c+k*(n+1)+j);
4 l6 y, m0 {: e. v7 t6 w                                        *(c+k*(n+1)+j)=p;
                                }
                for( j=i+1;j<=n-1;j++ )
1 I% s1 M  s3 R; V                {
4 q+ o# q  z; p, g                        p=(*(c+j*(n+1)+i))/(*(c+i*(n+1)+i));
* w. I. L, I# V( \4 s; j7 _                        for( t=i;t<=n;t++ )
/ _$ b7 ~3 O- O( O                                *(c+j*(n+1)+t)-=p*(*(c+i*(n+1)+t));
                }
7 Q' u. l5 d; _8 e0 B6 s/ o        }
        for( i=n-1;i>=0;i--)
        {
                for( j=n-1;j>=i+1;j--)
* ~! m+ C* Q! t  c# `$ d                        (*(c+i*(n+1)+n))-=x[j]*(*(c+i*(n+1)+j));
8 l8 m0 m+ k( L, z) I" Z                x[i]=*(c+i*(n+1)+n)/(*(c+i*(n+1)+i));
$ m9 |6 S4 k2 f, \9 m        }
        return x;
, j$ x% ?* [# o, A}

xiaozhangwh

谢谢分享，来试一下

deven1985

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

Rekcahpu

刚学高斯消元

0.9清1.8清2.7清

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

弘道

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