Gauss顺序消去法

木辰月生

2. %M函数gauss.m
   
3. 8 V1 v: B. s' |+ o/ V/ m3 U. }* |
4. function x=gauss(a,b)
   
5.   v! O. M, B6 c' V7 G
6. [n,m]=size(b);
   
7. + X9 G+ C! j1 W; X3 l
8. if size(a)==[n,n]
   
9. $ H7 i2 r* q2 w; l1 v+ u
10. for i=1:(n-1)) N$ `) ?, n! O7 g
    
12. if a(i,i)==0,error('Divide by zero');end3 N3 R  I\\" q5 `- a
    
14. b((i+1):n,:)=b((i+1):n,:)-a((i+1):n,i)*b(i,:)/a(i,i);0 p1 M5 g3 h8 B( Y; W
    
16. a((i+1):n,i:n)=a((i+1):n,i:n)-... ) O8 ^\\" K1 n* K( K, H( {3 j+ U
    
18. a((i+1):n,i)*a(i,i:n)/a(i,i);[a,b]
    
19. / s$ V7 k6 v+ H! f$ g0 V
20. end
    
21. 5 F8 d; `8 w/ U& z! q. r# G
22. if a(n,n)==0,error('Divide by zero');end* ]5 l8 l/ v+ P# U/ H) h
    
24. x(n,:)=b(n,:)/a(n,n);4 j9 q- A; F# M+ F6 A- c( l, I
    
26. for i=n-1:-1:19 I/ ^/ r) i4 \$ f( z
    
28. x(i,:)=(b(i,:)-a(i,(i+1):n)*x((i+1):n,:))/a(i,i);
    
29. 5 V* }, L8 \: Y$ B9 H- W
30. end
    
31. 5 {; V; a7 l9 o3 U0 g
32. else,error('Dimension of matrix must agree'),end

%M函数gauss.m<br /> <span style=\\"display:none\\">8 V1 v: B. s' |+ o/ V/ m3 U. }* |</span> function x=gauss(a,b)<br /> <span style=\\"display:none\\">  v! O. M, B6 c' V7 G</span> [n,m]=size(b);<br /> <span style=\\"display:none\\">+ X9 G+ C! j1 W; X3 l</span> if size(a)==[n,n]<br /> <span style=\\"display:none\\">$ H7 i2 r* q2 w; l1 v+ u</span> for i=1:(n-1)<font class=\\"jammer\\">) N$ `) ?, n! O7 g</font><br /> if a(i,i)==0,error('Divide by zero');end<font class=\\"jammer\\">3 N3 R  I\\" q5 `- a</font><br /> b((i+1):n,:)=b((i+1):n,:)-a((i+1):n,i)*b(i,:)/a(i,i);<font class=\\"jammer\\">0 p1 M5 g3 h8 B( Y; W</font><br /> a((i+1):n,i:n)=a((i+1):n,i:n)-... <font class=\\"jammer\\">) O8 ^\\" K1 n* K( K, H( {3 j+ U</font><br /> a((i+1):n,i)*a(i,i:n)/a(i,i);[a,b]<br /> <span style=\\"display:none\\">/ s$ V7 k6 v+ H! f$ g0 V</span> end<br /> <span style=\\"display:none\\">5 F8 d; `8 w/ U& z! q. r# G</span> if a(n,n)==0,error('Divide by zero');end<font class=\\"jammer\\">* ]5 l8 l/ v+ P# U/ H) h</font><br /> x(n,:)=b(n,:)/a(n,n);<font class=\\"jammer\\">4 j9 q- A; F# M+ F6 A- c( l, I</font><br /> for i=n-1:-1:1<font class=\\"jammer\\">9 I/ ^/ r) i4 \$ f( z</font><br /> x(i,:)=(b(i,:)-a(i,(i+1):n)*x((i+1):n,:))/a(i,i);<br /> <span style=\\"display:none\\">5 V* }, L8 \: Y$ B9 H- W</span> end<br /> <span style=\\"display:none\\">5 {; V; a7 l9 o3 U0 g</span> else,error('Dimension of matrix must agree'),end

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

阳台7777

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

Agnes2_3_4_5_6

不错，帮顶！！~

sunying52zs

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

wangyslyh

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

zjlsmrl

能不能给解释一下呢？看不懂耶。。
2 T/ b' o$ S& }# v& ]

yjh2010

到底应用在哪方面啊

yamarychen

回复 木辰月生 的帖子
" S. s/ p& Y; L8 W2 D7 a
: W" b$ R5 \7 z* X9 |# F$ e顶啊顶~~

" O; F. l+ `) E; \) M* b% b% o   

meseeme

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

pcm

留个脚印```````