Gauss顺序消去法

木辰月生

2. %M函数gauss.m/ n& D  @% q9 j3 [
   
4. function x=gauss(a,b)
   
5. + Q8 {* Y' V  f7 \
6. [n,m]=size(b);
   
7. ! b! Y& T! ?$ d, s, u
8. if size(a)==[n,n]
   
9. # R, P- C  b/ x
10. for i=1:(n-1)- L6 g2 R) d+ n0 f
    
12. if a(i,i)==0,error('Divide by zero');end2 S/ f: q+ {; f6 Q
    
14. b((i+1):n,:)=b((i+1):n,:)-a((i+1):n,i)*b(i,:)/a(i,i);4 g0 P6 r. y7 R' W, N& a
    
16. a((i+1):n,i:n)=a((i+1):n,i:n)-...
    
17. 3 P# @3 e+ t0 l, r2 x, K7 |, W$ z
18. a((i+1):n,i)*a(i,i:n)/a(i,i);[a,b]
    
19. 2 F' |7 v8 Y  Z& Y5 ?7 S8 w; o# s
20. end
    
21. : E- S  C6 F2 H2 c\\" M
22. if a(n,n)==0,error('Divide by zero');end
    
23. 1 [/ s8 o! C7 u
24. x(n,:)=b(n,:)/a(n,n);- d9 @% n2 Z8 A- m\\" c
    
26. for i=n-1:-1:1
    
27. % }0 Y+ w3 V/ ]2 d) m
28. x(i,:)=(b(i,:)-a(i,(i+1):n)*x((i+1):n,:))/a(i,i);2 x4 }0 p0 T( X$ ~
    
30. end
    
31. 3 @% b; h/ D\\" w) @6 n
32. else,error('Dimension of matrix must agree'),end

%M函数gauss.m<font class=\\"jammer\\">/ n& D  @% q9 j3 [</font><br /> function x=gauss(a,b)<br /> <span style=\\"display:none\\">+ Q8 {* Y' V  f7 \</span> [n,m]=size(b);<br /> <span style=\\"display:none\\">! b! Y& T! ?$ d, s, u</span> if size(a)==[n,n]<br /> <span style=\\"display:none\\"># R, P- C  b/ x</span> for i=1:(n-1)<font class=\\"jammer\\">- L6 g2 R) d+ n0 f</font><br /> if a(i,i)==0,error('Divide by zero');end<font class=\\"jammer\\">2 S/ f: q+ {; f6 Q</font><br /> b((i+1):n,:)=b((i+1):n,:)-a((i+1):n,i)*b(i,:)/a(i,i);<font class=\\"jammer\\">4 g0 P6 r. y7 R' W, N& a</font><br /> a((i+1):n,i:n)=a((i+1):n,i:n)-... <br /> <span style=\\"display:none\\">3 P# @3 e+ t0 l, r2 x, K7 |, W$ z</span> a((i+1):n,i)*a(i,i:n)/a(i,i);[a,b]<br /> <span style=\\"display:none\\">2 F' |7 v8 Y  Z& Y5 ?7 S8 w; o# s</span> end<br /> <span style=\\"display:none\\">: E- S  C6 F2 H2 c\\" M</span> if a(n,n)==0,error('Divide by zero');end<br /> <span style=\\"display:none\\">1 [/ s8 o! C7 u</span> x(n,:)=b(n,:)/a(n,n);<font class=\\"jammer\\">- d9 @% n2 Z8 A- m\\" c</font><br /> for i=n-1:-1:1<br /> <span style=\\"display:none\\">% }0 Y+ w3 V/ ]2 d) m</span> x(i,:)=(b(i,:)-a(i,(i+1):n)*x((i+1):n,:))/a(i,i);<font class=\\"jammer\\">2 x4 }0 p0 T( X$ ~</font><br /> end<br /> <span style=\\"display:none\\">3 @% b; h/ D\\" w) @6 n</span> else,error('Dimension of matrix must agree'),end

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

阳台7777

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

Agnes2_3_4_5_6

不错，帮顶！！~

sunying52zs

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

wangyslyh

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

zjlsmrl

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

yjh2010

到底应用在哪方面啊

yamarychen

回复 木辰月生 的帖子

: D- a5 P0 x. {7 O9 S: f/ m: E1 A顶啊顶~~

/ Z3 ?1 P7 O1 w   

meseeme

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

pcm

留个脚印```````