Gauss顺序消去法

木辰月生

2. %M函数gauss.m. ~; L% d. r: L. x$ o$ o% P4 D7 \
   
4. function x=gauss(a,b)( t% k; {1 w# L/ m
   
6. [n,m]=size(b);
   
7. 3 q4 }1 v9 _( J+ W9 H
8. if size(a)==[n,n]! X# F\\" |2 D' n/ l6 ?
   
10. for i=1:(n-1). p. l5 h; A( X0 b1 p) |5 B2 [; l
    
12. if a(i,i)==0,error('Divide by zero');end$ Y0 o3 v9 w, ?6 \. Y* o
    
14. b((i+1):n,:)=b((i+1):n,:)-a((i+1):n,i)*b(i,:)/a(i,i);
    
15. ( k1 N3 I- t+ H7 U, a+ R: X# }
16. a((i+1):n,i:n)=a((i+1):n,i:n)-...
    
17. , v0 Q' n4 a1 E
18. a((i+1):n,i)*a(i,i:n)/a(i,i);[a,b]
    
19. ) s8 l+ z7 [; Q4 @# p
20. end
    
21. ' b  b6 M9 l8 h# [
22. if a(n,n)==0,error('Divide by zero');end
    
23.   T& m9 g1 c5 k  N. Y
24. x(n,:)=b(n,:)/a(n,n);: @2 }/ c5 X$ e. v\\" T9 u5 ^3 d. L
    
26. for i=n-1:-1:10 p& \3 L$ i8 Y2 ~1 v- Y; p) O+ b. ~5 L
    
28. x(i,:)=(b(i,:)-a(i,(i+1):n)*x((i+1):n,:))/a(i,i);$ V% E& I. H1 {* e; E
    
30. end6 P5 Q% M/ Y- d2 u  V\\" n
    
32. else,error('Dimension of matrix must agree'),end

%M函数gauss.m<font class=\\"jammer\\">. ~; L% d. r: L. x$ o$ o% P4 D7 \</font><br /> function x=gauss(a,b)<font class=\\"jammer\\">( t% k; {1 w# L/ m</font><br /> [n,m]=size(b);<br /> <span style=\\"display:none\\">3 q4 }1 v9 _( J+ W9 H</span> if size(a)==[n,n]<font class=\\"jammer\\">! X# F\\" |2 D' n/ l6 ?</font><br /> for i=1:(n-1)<font class=\\"jammer\\">. p. l5 h; A( X0 b1 p) |5 B2 [; l</font><br /> if a(i,i)==0,error('Divide by zero');end<font class=\\"jammer\\">$ Y0 o3 v9 w, ?6 \. Y* o</font><br /> b((i+1):n,:)=b((i+1):n,:)-a((i+1):n,i)*b(i,:)/a(i,i);<br /> <span style=\\"display:none\\">( k1 N3 I- t+ H7 U, a+ R: X# }</span> a((i+1):n,i:n)=a((i+1):n,i:n)-... <br /> <span style=\\"display:none\\">, v0 Q' n4 a1 E</span> a((i+1):n,i)*a(i,i:n)/a(i,i);[a,b]<br /> <span style=\\"display:none\\">) s8 l+ z7 [; Q4 @# p</span> end<br /> <span style=\\"display:none\\">' b  b6 M9 l8 h# [</span> if a(n,n)==0,error('Divide by zero');end<br /> <span style=\\"display:none\\">  T& m9 g1 c5 k  N. Y</span> x(n,:)=b(n,:)/a(n,n);<font class=\\"jammer\\">: @2 }/ c5 X$ e. v\\" T9 u5 ^3 d. L</font><br /> for i=n-1:-1:1<font class=\\"jammer\\">0 p& \3 L$ i8 Y2 ~1 v- Y; p) O+ b. ~5 L</font><br /> x(i,:)=(b(i,:)-a(i,(i+1):n)*x((i+1):n,:))/a(i,i);<font class=\\"jammer\\">$ V% E& I. H1 {* e; E</font><br /> end<font class=\\"jammer\\">6 P5 Q% M/ Y- d2 u  V\\" n</font><br /> else,error('Dimension of matrix must agree'),end

9 a2 S4 C/ q2 m8 n计算方法中的有关Gauss消去法的M文件程序代码

阳台7777

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

Agnes2_3_4_5_6

不错，帮顶！！~

sunying52zs

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

wangyslyh

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

zjlsmrl

能不能给解释一下呢？看不懂耶。。
$ Q) ], `5 O9 J: A4 q& c

yjh2010

到底应用在哪方面啊

yamarychen

回复 木辰月生 的帖子
/ t1 b) M/ {( `$ ]; W9 h4 ~/ w
+ o: a7 K, V  N' D) Y顶啊顶~~

   

meseeme

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

pcm

留个脚印```````