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

hduscimc

#include<iostream>
8 Q! V9 \7 T+ J: Y% T#include<cmath>
using namespace std;
int main()
{
' D% K5 b$ y5 X0 x# W0 a9 ?        int i;
) k3 h$ Y6 Z; J; e( O) C3 D$ F# a        float *x;
7 b! H6 V/ A7 b' d* ?3 s        float c[3][4] ={0.101,2.304,3.555,1.183,
! X: I9 M/ {  z9 H                                        -1.347,3.712,4.623,2.137,
4 o# q1 E# W. C- h                                        -2.835,1.072,5.643,3.035};
; z! ~8 s3 Y/ e1 j5 ^        float *ColPivot(float *,int);
        x=ColPivot(c[0],3);
! a; X2 Z( R, h$ C3 k! C/ E        for( i=0;i<=2;i++ )
+ i0 Z3 p) b% @9 p( l3 c                cout<<"x("<<i<<")="<<x[i]<<endl;
}

/ j, u" K% }2 V, Y5 H0 F$ k& ]; [float *ColPivot( float *c,int n )
{
" t) U- }( \. p2 R9 O        int i,j,t,k;
        float *x,p;
8 Q1 g* ~  r, x9 Y* o, q         x=new float[n*sizeof(float)];
        for( i=0;i<=n-2;i++)
        {
                k=i;
                for(j=i+1;j<=n-1;j++)
+ L& x% `9 v& q. \  f$ Q                        if(fabs(*(c+j*(n+1)+i))>(fabs(*(c+k*(n+1)+i))))
                                k=j;
                        if(k!=i)
                                for( j=i;j<=n;j++ )
                                {
                                        p=*(c+i*(n+1)+j);
. Q( }& h6 P8 b5 Y" a                                        *(c+i*(n+1)+j)=*(c+k*(n+1)+j);
" z3 E9 b+ y! t                                        *(c+k*(n+1)+j)=p;
- `: h; o: @' `                                }
                for( j=i+1;j<=n-1;j++ )
                {
                        p=(*(c+j*(n+1)+i))/(*(c+i*(n+1)+i));
. t& k8 ~$ {: ~                        for( t=i;t<=n;t++ )
                                *(c+j*(n+1)+t)-=p*(*(c+i*(n+1)+t));
                }
        }
        for( i=n-1;i>=0;i--)
        {
( S! [$ q" B) a' h& X2 Q1 Q2 f                for( j=n-1;j>=i+1;j--)
' f/ u) x4 `, S6 O( \9 E' d% }                        (*(c+i*(n+1)+n))-=x[j]*(*(c+i*(n+1)+j));
( N" _, q( w% K1 ~# m; F+ I  Y& b                x[i]=*(c+i*(n+1)+n)/(*(c+i*(n+1)+i));
. b$ W& Q: ~& M# l* i        }
, t, e+ x+ J- H& q- i+ [) j5 b        return x;
# x: A! X+ v* K& C% K}

xiaozhangwh

谢谢分享，来试一下

deven1985

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

Rekcahpu

刚学高斯消元

0.9清1.8清2.7清

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

弘道

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