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

hduscimc

#include<iostream>
#include<cmath>
using namespace std;
" Z6 G, L6 ~4 Y( h4 j/ t. u+ ]( Q/ vint main()
{
        int i;
        float *x;
* _9 V7 L# O2 K8 [8 c$ c( @        float c[3][4] ={0.101,2.304,3.555,1.183,
" U9 `0 v, `6 D& j& x: z                                        -1.347,3.712,4.623,2.137,
$ [- r' k, u8 |# I: R( j                                        -2.835,1.072,5.643,3.035};
% S; r. Y" ]4 V3 }% e        float *ColPivot(float *,int);
$ i" q2 Z/ y! E) g" v        x=ColPivot(c[0],3);
        for( i=0;i<=2;i++ )
% Y% i) }) e5 \5 s+ O                cout<<"x("<<i<<")="<<x[i]<<endl;
' x0 C" B: h" R}

float *ColPivot( float *c,int n )
{
9 L" y5 Q+ S8 o3 |6 b$ x6 s        int i,j,t,k;
! s7 T. d1 y; ]6 G        float *x,p;
3 ^! m: v% M8 S7 [         x=new float[n*sizeof(float)];
        for( i=0;i<=n-2;i++)
: ^1 n; B# N/ k. Z* X. _2 V        {
                k=i;
5 [6 s* @; R7 E                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)
3 k6 v6 W( A, H; H" t) P% B                                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;
  s/ V3 j4 \# n. I! Z5 I0 O4 _& O) K                                }
! k9 |; D  v* f2 {                for( j=i+1;j<=n-1;j++ )
                {
+ A6 q3 P8 L$ l* V( m3 s                        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));
                }
1 H, [, o) a" a6 k- ?9 H        }
        for( i=n-1;i>=0;i--)
# D4 t$ G; Z3 ^7 l1 o7 c- q        {
. `9 h- H1 q: d* k9 G5 P! U6 F                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));
        }
        return x;
  L) |. x6 I) i, a9 j6 _}

xiaozhangwh

谢谢分享，来试一下

deven1985

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

Rekcahpu

刚学高斯消元

0.9清1.8清2.7清

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

弘道

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