12 3 4 5 / 5 页下一页

查看: 23362|回复: 48

[代码资源] Gauss顺序消去法

[复制链接]

字体大小: 正常放大

木辰月生

4 主题	3 听众	25 积分

升级 21.05%

该用户从未签到

电梯直达

1^#

发表于 2010-9-7 21:24 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

%M函数gauss.m
( Z. P% d/ @4 M! I7 n3 f4 \: l2 j
function x=gauss(a,b)' k+ w( o$ a+ V4 w( f) Y; w) v
[n,m]=size(b);5 a$ K& ]# b6 @' B: E! C& N) f$ i
if size(a)==[n,n]
. f' ~, `( ?+ B. O& V* Q4 u
for i=1:(n-1)$ y& ~1 }2 t% p u0 @2 d
if a(i,i)==0,error('Divide by zero');end* a( z. M# z/ ^# g& J3 n
b((i+1):n,:)=b((i+1):n,:)-a((i+1):n,i)*b(i,:)/a(i,i);\\" p F3 y! I2 i0 X
a((i+1):n,i:n)=a((i+1):n,i:n)-... / u8 ^; w& n7 ~0 c! Y; P' E' a
a((i+1):n,i)*a(i,i:n)/a(i,i);[a,b]
7 k! c1 @\\" J9 R- p3 ?
end1 [+ p, o5 _\\" G$ O5 ?3 w' E% |8 [\\" l9 k
if a(n,n)==0,error('Divide by zero');end\\" h% U1 _/ c+ v* X\\" c
x(n,:)=b(n,:)/a(n,n);
; R% e% W8 ]1 ?
for i=n-1:-1:1
, {! ^' F. _( g6 [
x(i,:)=(b(i,:)-a(i,(i+1):n)*x((i+1):n,:))/a(i,i);
* B d9 \7 R1 Q! c3 o
end; E( {& F6 z2 f
else,error('Dimension of matrix must agree'),end

%M函数gauss.m 
( Z. P% d/ @4 M! I7 n3 f4 \: l2 j
function x=gauss(a,b)' k+ w( o$ a+ V4 w( f) Y; w) v

[n,m]=size(b);5 a$ K& ]# b6 @' B: E! C& N) f$ i

if size(a)==[n,n] 
. f' ~, `( ?+ B. O& V* Q4 u
for i=1:(n-1)$ y& ~1 }2 t% p  u0 @2 d

if a(i,i)==0,error('Divide by zero');end* a( z. M# z/ ^# g& J3 n

b((i+1):n,:)=b((i+1):n,:)-a((i+1):n,i)*b(i,:)/a(i,i);\\" p  F3 y! I2 i0 X

a((i+1):n,i:n)=a((i+1):n,i:n)-... / u8 ^; w& n7 ~0 c! Y; P' E' a

a((i+1):n,i)*a(i,i:n)/a(i,i);[a,b] 
7 k! c1 @\\" J9 R- p3 ?
end1 [+ p, o5 _\\" G$ O5 ?3 w' E% |8 [\\" l9 k

if a(n,n)==0,error('Divide by zero');end\\" h% U1 _/ c+ v* X\\" c

x(n,:)=b(n,:)/a(n,n); 
; R% e% W8 ]1 ?
for i=n-1:-1:1 
, {! ^' F. _( g6 [
x(i,:)=(b(i,:)-a(i,(i+1):n)*x((i+1):n,:))/a(i,i); 
* B  d9 \7 R1 Q! c3 o
end; E( {& F6 z2 f

else,error('Dimension of matrix must agree'),end

计算方法中的有关Gauss消去法的M文件程序代码

zan

转播0 淘帖0 分享0 收藏0 支持0 反对0 微信

使用道具举报

阳台7777

31 主题	3 听众	1770 积分

升级 77%

TA的每日心情

	无聊 2012-10-25 07:54

签到天数: 39 天

[LV.5]常住居民I

群组: 2012第三期美赛培训

2^#

发表于 2010-9-7 21:42 |只看该作者 |招呼Ta 关注Ta

帮你顶顶，暂时看不懂···

我也说一句

发表

使用道具举报

Agnes2_3_4_5_6

1 主题	3 听众	137 积分

升级 18.5%

该用户从未签到

群组: Matlab讨论组

3^#

发表于 2010-9-7 21:49 |只看该作者 |招呼Ta 关注Ta

不错，帮顶！！~

我也说一句

发表

使用道具举报

sunying52zs

5 主题	3 听众	109 积分

升级 4.5%

该用户从未签到

群组: 数学趣味、游戏、IQ等

群组: Matlab讨论组

群组: Linux推广

群组: 数学建摸协会

群组: LINGO

4^#

发表于 2010-9-7 21:57 |只看该作者 |招呼Ta 关注Ta

是不是很有用的我看不懂啊~~

我也说一句

发表

使用道具举报

wangyslyh

9 主题	3 听众	145 积分

升级 22.5%

TA的每日心情

	郁闷 2012-9-22 18:28

签到天数: 7 天

[LV.3]偶尔看看II

群组: 学术交流A

群组: 西安交大数学建模

群组: B题讨论群

5^#

发表于 2010-9-7 22:04 |只看该作者 |招呼Ta 关注Ta

帮你顶顶，暂时看不懂···

我也说一句

发表

使用道具举报

zjlsmrl

1 主题	3 听众	43 积分

升级 40%

该用户从未签到

自我介绍: 200 字节以内

不支持自定义 Discuz! 代码

群组: 数学建模保研联盟

6^#

发表于 2010-9-7 22:39 |只看该作者 |招呼Ta 关注Ta

能不能给解释一下呢？看不懂耶。。

我也说一句

发表

使用道具举报

yjh2010

6 主题	3 听众	52 积分

升级 49.47%

该用户从未签到

群组: Matlab讨论组

7^#

发表于 2010-9-7 23:21 |只看该作者 |招呼Ta 关注Ta

到底应用在哪方面啊

我也说一句

发表

使用道具举报

yamarychen

1 主题	3 听众	48 积分

升级 45.26%

该用户从未签到

8^#

发表于 2010-9-7 23:39 |只看该作者 |招呼Ta 关注Ta

回复 木辰月生 的帖子

顶啊顶~~

我也说一句

发表

使用道具举报

meseeme

0 主题	2 听众	64 积分

升级 62.11%

该用户从未签到

9^#

发表于 2010-9-8 08:00 |只看该作者 |招呼Ta 关注Ta

楼主那种裂纸欲出的大手笔，竟使我忍不住一次次的翻开楼主的帖子……

我也说一句

发表

使用道具举报

pcm

0 主题	2 听众	34 积分

升级 30.53%

该用户从未签到

国际赛参赛者

10^#

发表于 2010-9-8 12:00 |只看该作者 |招呼Ta 关注Ta

留个脚印```````

我也说一句

发表

使用道具举报

12 3 4 5 / 5 页下一页

返回列表

帐号		密码		只需要一步，快速开始		注册地址	找回密码

[代码资源] Gauss顺序消去法

浏览过的版块

新人进步奖

QQ

电话咨询

关于我们| 联系我们| 诚征英才| 对外合作| 产品服务|