Gauss顺序消去法

木辰月生

2. %M函数gauss.m; D. M. m- s7 S# D( J
   
4. function x=gauss(a,b)
   
5. 9 Q  |2 ^! y/ {; r: _3 ?' v7 _\\" G
6. [n,m]=size(b);  M7 v! |8 T. M1 M) ?5 x
   
8. if size(a)==[n,n]
   
9. 2 e3 V% `\\" w3 |( o
10. for i=1:(n-1)
    
11. + t5 p& H+ x: _: T0 |  M  U3 U
12. if a(i,i)==0,error('Divide by zero');end9 {1 \. o2 ]. F. o, ]% H% H+ N. o
    
14. b((i+1):n,:)=b((i+1):n,:)-a((i+1):n,i)*b(i,:)/a(i,i);/ b3 }  m$ m/ R/ H3 [* g
    
16. a((i+1):n,i:n)=a((i+1):n,i:n)-...
    
17. 7 Q- G1 @5 j# Z' O& {7 @. l+ z/ f0 _
18. a((i+1):n,i)*a(i,i:n)/a(i,i);[a,b]
    
19. * x, A# E. ~! ]6 @: f4 o' I
20. end
    
21. ; Q$ ^: y$ t! h5 K$ v8 h
22. if a(n,n)==0,error('Divide by zero');end
    
23. 5 d+ j, Q- Z! |
24. x(n,:)=b(n,:)/a(n,n);. X6 j0 q2 u( F4 {; t: c, {
    
26. for i=n-1:-1:1
    
27. - G: }9 i7 v# c, Q
28. x(i,:)=(b(i,:)-a(i,(i+1):n)*x((i+1):n,:))/a(i,i);5 j  g# T* j: E% Z& f
    
30. end
    
31. * l) G1 R7 F% A: C6 g3 X+ i* S. e
32. else,error('Dimension of matrix must agree'),end

%M函数gauss.m<font class=\\"jammer\\">; D. M. m- s7 S# D( J</font><br /> function x=gauss(a,b)<br /> <span style=\\"display:none\\">9 Q  |2 ^! y/ {; r: _3 ?' v7 _\\" G</span> [n,m]=size(b);<font class=\\"jammer\\">  M7 v! |8 T. M1 M) ?5 x</font><br /> if size(a)==[n,n]<br /> <span style=\\"display:none\\">2 e3 V% `\\" w3 |( o</span> for i=1:(n-1)<br /> <span style=\\"display:none\\">+ t5 p& H+ x: _: T0 |  M  U3 U</span> if a(i,i)==0,error('Divide by zero');end<font class=\\"jammer\\">9 {1 \. o2 ]. F. o, ]% H% H+ N. o</font><br /> b((i+1):n,:)=b((i+1):n,:)-a((i+1):n,i)*b(i,:)/a(i,i);<font class=\\"jammer\\">/ b3 }  m$ m/ R/ H3 [* g</font><br /> a((i+1):n,i:n)=a((i+1):n,i:n)-... <br /> <span style=\\"display:none\\">7 Q- G1 @5 j# Z' O& {7 @. l+ z/ f0 _</span> a((i+1):n,i)*a(i,i:n)/a(i,i);[a,b]<br /> <span style=\\"display:none\\">* x, A# E. ~! ]6 @: f4 o' I</span> end<br /> <span style=\\"display:none\\">; Q$ ^: y$ t! h5 K$ v8 h</span> if a(n,n)==0,error('Divide by zero');end<br /> <span style=\\"display:none\\">5 d+ j, Q- Z! |</span> x(n,:)=b(n,:)/a(n,n);<font class=\\"jammer\\">. X6 j0 q2 u( F4 {; t: c, {</font><br /> for i=n-1:-1:1<br /> <span style=\\"display:none\\">- G: }9 i7 v# c, Q</span> x(i,:)=(b(i,:)-a(i,(i+1):n)*x((i+1):n,:))/a(i,i);<font class=\\"jammer\\">5 j  g# T* j: E% Z& f</font><br /> end<br /> <span style=\\"display:none\\">* l) G1 R7 F% A: C6 g3 X+ i* S. e</span> else,error('Dimension of matrix must agree'),end

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

阳台7777

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

Agnes2_3_4_5_6

不错，帮顶！！~

sunying52zs

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

wangyslyh

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

zjlsmrl

能不能给解释一下呢？看不懂耶。。
" P6 H. U  q- N. @- Z4 K' C

yjh2010

到底应用在哪方面啊

yamarychen

回复 木辰月生 的帖子
: ]! r9 Q5 I( V  [( m) o
顶啊顶~~

' a! x# \: ~4 r- X, j8 C1 f2 Q   

meseeme

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

pcm

留个脚印```````