数学建模社区-数学中国

标题: [讨论]高斯消元法的程序对比 [打印本页]
作者: god    时间: 2005-1-19 17:01
标题: [讨论]高斯消元法的程序对比
下面我收集的高斯消元法的程序。还是请大家看看，并提点意见。% @/ y6 p+ [6 G7 U, Q9 a! W: `" X
[attach]966[/attach]

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

1.68 KB, 下载次数: 78, 下载积分: 体力 -2 点

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

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

1.09 KB, 下载次数: 67, 下载积分: 体力 -2 点

[讨论]高斯消元法的程序对比
作者: chzjao    时间: 2005-9-24 11:45
标题: 发下我编的,请各位帮忙修改下
[attach]2037[/attach]
, |/ ~- {& }5 Q2 `
[此贴子已经被作者于2005-9-24 11:46:03编辑过]

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

1.14 KB, 下载次数: 41, 下载积分: 体力 -2 点

[讨论]高斯消元法的程序对比
作者: jiangzsm    时间: 2006-4-17 14:42

还没看,感谢分享.

作者: lfer2008    时间: 2006-11-26 13:42

谢谢

作者: ylsi41    时间: 2006-12-4 23:50

谢谢

作者: ladeng321    时间: 2007-1-2 19:38
谢谢分享
作者: fansmiling    时间: 2007-6-18 12:17

8错8错

作者: LH611    时间: 2007-12-18 21:57
谢谢分享，要又AS写的就更好了！
作者: 水木年华zzu    时间: 2007-12-19 19:04
总的来说写的不错!继续努力
作者: shuainaruto    时间: 2009-4-11 11:05
感谢分享，下来研究一下
作者: hduscimc    时间: 2009-6-21 19:04
#include<iostream>
* C# `! }, e. u% G#include<cmath>7 w& W5 V. ?0 _2 Y: J' J- _$ p; ~
using namespace std;0 M# C  F. g2 ~7 p* U8 k9 s8 J
int main()
+ s; w" e9 u+ P5 r{
* d" O8 o" K( I! I0 L        int i;2 ?8 N5 `* m; V! y& Y* ~# z
        float *x;0 W  j0 ?' |  h9 l5 Z8 }
        float c[3][4] ={0.101,2.304,3.555,1.183,9 N! K4 V* p" l9 [! [. |9 l( f( U
                                        -1.347,3.712,4.623,2.137,
1 H9 J- Z- S% y7 L                                        -2.835,1.072,5.643,3.035};
; L/ \  L  G8 K1 v* _2 h        float *ColPivot(float *,int);/ ]6 @# G) G6 x. h/ w5 H: V0 n- T
        x=ColPivot(c[0],3);" [% p2 A# _& @+ e8 {& N
        for( i=0;i<=2;i++ )
) l1 w2 G8 ~3 [/ ^) G" V, `                cout<<"x("<<i<<")="<<x[i]<<endl;
  _. Y1 R& q; i}. l2 j4 Z3 |8 L) U: o
" s/ ?: m; k2 d
float *ColPivot( float *c,int n ). m0 Y# d. \% |9 a7 G( t
{
1 `) K* y. I$ Q8 l1 z        int i,j,t,k;
! k4 d6 |9 S- P( h        float *x,p;
& c0 e0 Q7 i- Z8 O5 o2 W% y- q         x=new float[n*sizeof(float)];
9 S; i3 f# g# h6 H! D        for( i=0;i<=n-2;i++)
: a  |+ `6 Y) J        {
& i8 Y7 k2 z. t, R                k=i;( q4 m0 d5 P0 F- m* H9 h& }1 K
                for(j=i+1;j<=n-1;j++)
6 R3 Q! w4 A8 K: `8 F! ]4 L                        if(fabs(*(c+j*(n+1)+i))>(fabs(*(c+k*(n+1)+i))))
; t+ G  b) N; A3 w4 V) _6 R                                k=j;
" I  r' }: Y$ m) U; e% y9 M                        if(k!=i)
3 {+ d: T( ]$ Y3 j1 i                                for( j=i;j<=n;j++ )" J+ ]4 J$ \# G3 ]/ H6 {
                                {' [. G2 ^0 w9 C' k5 c0 \
                                        p=*(c+i*(n+1)+j);( c5 k) |- h7 }6 X# X
                                        *(c+i*(n+1)+j)=*(c+k*(n+1)+j);1 H8 \) z! P. k" T( }2 _& [) P! A
                                        *(c+k*(n+1)+j)=p;
$ j! u1 u. s6 Z                                }
, Z' u# f. N! v+ b) C7 K7 h) u                for( j=i+1;j<=n-1;j++ )( @7 ~. B1 O+ R; ~
                {1 [' B* V- g- f+ ^3 @
                        p=(*(c+j*(n+1)+i))/(*(c+i*(n+1)+i));1 @( w  ?0 t% @4 Z) M
                        for( t=i;t<=n;t++ )' }* A8 ]7 ~+ ~0 _
                                *(c+j*(n+1)+t)-=p*(*(c+i*(n+1)+t));
3 }* ~8 b! j: K1 @8 L/ R                }) Y+ v, i: ~7 M' V2 h: k; G) X
        }' I) C4 q% [/ S3 x
        for( i=n-1;i>=0;i--)# n& Y1 B; X& J: H# \) m
        {$ U& ?4 ^  P3 |% ]; g
                for( j=n-1;j>=i+1;j--)
6 R9 @" H/ P' `- H3 ]% a) N+ Z* }* }                        (*(c+i*(n+1)+n))-=x[j]*(*(c+i*(n+1)+j));
& N& [  O$ Y8 A6 K& d  z$ r: Z                x[i]=*(c+i*(n+1)+n)/(*(c+i*(n+1)+i));+ ~" ^2 l* p( C, d& r/ I
        }
3 Q; e* `' ?  j% J' y% M0 \% @        return x;
; F* e& ~; R. P% ^+ y" S% k. I}
作者: xiaozhangwh    时间: 2009-7-1 21:56
谢谢分享，来试一下
作者: deven1985    时间: 2010-1-1 09:06
，不错，前几天的数值分析课上也是做的这个！
作者: Rekcahpu    时间: 2010-11-23 13:44
刚学高斯消元
作者: 0.9清1.8清2.7清    时间: 2012-8-31 16:04
wuyuwuyul 自己程序拿出来看看呗！干嘛还要体力A！
作者: 弘道    时间: 2014-7-29 09:12
谢谢楼主……辛苦啦！………………



欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5