Gauss顺序消去法

木辰月生

2. %M函数gauss.m- S  R  ?6 f/ s9 c9 H3 _
   
4. function x=gauss(a,b)
   
5. 2 U2 m: o; _8 _# B9 P8 W& l0 m: Y
6. [n,m]=size(b);
   
7. 6 b1 }) l3 L3 _# g* I3 I
8. if size(a)==[n,n]) }4 p. s4 y6 E, _4 a
   
10. for i=1:(n-1)
    
11. ) K% {3 D! m! D& S- n
12. if a(i,i)==0,error('Divide by zero');end
    
13. % l6 e9 a% [2 H# t3 }4 c1 u$ r* J
14. b((i+1):n,:)=b((i+1):n,:)-a((i+1):n,i)*b(i,:)/a(i,i);
    
15. / L# W; l! _  @( z, q  p6 ]0 J! g
16. a((i+1):n,i:n)=a((i+1):n,i:n)-...
    
17. $ C3 ?6 z# y4 E+ {8 e' v
18. a((i+1):n,i)*a(i,i:n)/a(i,i);[a,b]0 c2 B/ X( f3 c
    
20. end. l) v- i+ c; {2 `* x
    
22. if a(n,n)==0,error('Divide by zero');end
    
23. - u& ~$ a4 x( y3 g+ f) T
24. x(n,:)=b(n,:)/a(n,n);% h' {  f) p  V+ L; \6 J. q
    
26. for i=n-1:-1:1
    
27. + B: p$ j6 l* S8 }4 f. G; q
28. x(i,:)=(b(i,:)-a(i,(i+1):n)*x((i+1):n,:))/a(i,i);0 k  Z# w/ D4 ^, R
    
30. end' g; B& c6 l( ?4 X
    
32. else,error('Dimension of matrix must agree'),end

%M函数gauss.m<font class=\\"jammer\\">- S  R  ?6 f/ s9 c9 H3 _</font><br /> function x=gauss(a,b)<br /> <span style=\\"display:none\\">2 U2 m: o; _8 _# B9 P8 W& l0 m: Y</span> [n,m]=size(b);<br /> <span style=\\"display:none\\">6 b1 }) l3 L3 _# g* I3 I</span> if size(a)==[n,n]<font class=\\"jammer\\">) }4 p. s4 y6 E, _4 a</font><br /> for i=1:(n-1)<br /> <span style=\\"display:none\\">) K% {3 D! m! D& S- n</span> if a(i,i)==0,error('Divide by zero');end<br /> <span style=\\"display:none\\">% l6 e9 a% [2 H# t3 }4 c1 u$ r* J</span> b((i+1):n,:)=b((i+1):n,:)-a((i+1):n,i)*b(i,:)/a(i,i);<br /> <span style=\\"display:none\\">/ L# W; l! _  @( z, q  p6 ]0 J! g</span> a((i+1):n,i:n)=a((i+1):n,i:n)-... <br /> <span style=\\"display:none\\">$ C3 ?6 z# y4 E+ {8 e' v</span> a((i+1):n,i)*a(i,i:n)/a(i,i);[a,b]<font class=\\"jammer\\">0 c2 B/ X( f3 c</font><br /> end<font class=\\"jammer\\">. l) v- i+ c; {2 `* x</font><br /> if a(n,n)==0,error('Divide by zero');end<br /> <span style=\\"display:none\\">- u& ~$ a4 x( y3 g+ f) T</span> x(n,:)=b(n,:)/a(n,n);<font class=\\"jammer\\">% h' {  f) p  V+ L; \6 J. q</font><br /> for i=n-1:-1:1<br /> <span style=\\"display:none\\">+ B: p$ j6 l* S8 }4 f. G; q</span> x(i,:)=(b(i,:)-a(i,(i+1):n)*x((i+1):n,:))/a(i,i);<font class=\\"jammer\\">0 k  Z# w/ D4 ^, R</font><br /> end<font class=\\"jammer\\">' g; B& c6 l( ?4 X</font><br /> else,error('Dimension of matrix must agree'),end

2 e, @# B5 l# e& u' x/ @; R7 f+ [计算方法中的有关Gauss消去法的M文件程序代码

阳台7777

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

Agnes2_3_4_5_6

不错，帮顶！！~

sunying52zs

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

wangyslyh

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

zjlsmrl

能不能给解释一下呢？看不懂耶。。
$ Y5 l  Q" U5 K" n/ s

yjh2010

到底应用在哪方面啊

yamarychen

回复 木辰月生 的帖子

/ [& D3 H  I8 B: y2 W  u* F8 t* C% ?7 m顶啊顶~~
8 w. _* H4 D1 X" f4 R/ |
: I2 C3 Z8 k, g! v7 P   

meseeme

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

pcm

留个脚印```````