Gauss顺序消去法

木辰月生

2. %M函数gauss.m
   
3. % F2 e4 s) X4 y- K& a& m3 b
4. function x=gauss(a,b)7 V3 k/ p. X3 o1 i; I6 t2 t
   
6. [n,m]=size(b);
   
7. 5 U# q3 z\\" l: T8 Q, ~& i
8. if size(a)==[n,n]
   
9. * D7 N& n* Z3 c9 O
10. for i=1:(n-1); [. B1 ^+ e' R& y5 t9 k. B; u
    
12. if a(i,i)==0,error('Divide by zero');end
    
13. 2 L4 v; T' R' H8 S% a' Z0 \
14. b((i+1):n,:)=b((i+1):n,:)-a((i+1):n,i)*b(i,:)/a(i,i);
    
15. 6 O  b& M5 Q+ Q4 F5 _% P& ]2 L
16. a((i+1):n,i:n)=a((i+1):n,i:n)-... 4 L5 k+ [& Y+ N0 \+ V. _
    
18. a((i+1):n,i)*a(i,i:n)/a(i,i);[a,b]0 s, d6 |/ M0 f
    
20. end  B, v) k) `. y8 F' R0 r
    
22. if a(n,n)==0,error('Divide by zero');end) }. l/ m7 G0 M\\" X) y
    
24. x(n,:)=b(n,:)/a(n,n);
    
25. ) @; L: I( @! g) h5 i$ l. I
26. for i=n-1:-1:1
    
27. ( g7 h* n$ @# L* D/ V5 z7 j
28. x(i,:)=(b(i,:)-a(i,(i+1):n)*x((i+1):n,:))/a(i,i);
    
29. 3 |; m: p5 X2 V+ g0 A7 |
30. end
    
31. \\" l& P1 q* H/ S9 {: a\\" r$ H
32. else,error('Dimension of matrix must agree'),end

%M函数gauss.m<br /> <span style=\\"display:none\\">% F2 e4 s) X4 y- K& a& m3 b</span> function x=gauss(a,b)<font class=\\"jammer\\">7 V3 k/ p. X3 o1 i; I6 t2 t</font><br /> [n,m]=size(b);<br /> <span style=\\"display:none\\">5 U# q3 z\\" l: T8 Q, ~& i</span> if size(a)==[n,n]<br /> <span style=\\"display:none\\">* D7 N& n* Z3 c9 O</span> for i=1:(n-1)<font class=\\"jammer\\">; [. B1 ^+ e' R& y5 t9 k. B; u</font><br /> if a(i,i)==0,error('Divide by zero');end<br /> <span style=\\"display:none\\">2 L4 v; T' R' H8 S% a' Z0 \</span> b((i+1):n,:)=b((i+1):n,:)-a((i+1):n,i)*b(i,:)/a(i,i);<br /> <span style=\\"display:none\\">6 O  b& M5 Q+ Q4 F5 _% P& ]2 L</span> a((i+1):n,i:n)=a((i+1):n,i:n)-... <font class=\\"jammer\\">4 L5 k+ [& Y+ N0 \+ V. _</font><br /> a((i+1):n,i)*a(i,i:n)/a(i,i);[a,b]<font class=\\"jammer\\">0 s, d6 |/ M0 f</font><br /> end<font class=\\"jammer\\">  B, v) k) `. y8 F' R0 r</font><br /> if a(n,n)==0,error('Divide by zero');end<font class=\\"jammer\\">) }. l/ m7 G0 M\\" X) y</font><br /> x(n,:)=b(n,:)/a(n,n);<br /> <span style=\\"display:none\\">) @; L: I( @! g) h5 i$ l. I</span> for i=n-1:-1:1<br /> <span style=\\"display:none\\">( g7 h* n$ @# L* D/ V5 z7 j</span> x(i,:)=(b(i,:)-a(i,(i+1):n)*x((i+1):n,:))/a(i,i);<br /> <span style=\\"display:none\\">3 |; m: p5 X2 V+ g0 A7 |</span> end<br /> <span style=\\"display:none\\">\\" l& P1 q* H/ S9 {: a\\" r$ H</span> else,error('Dimension of matrix must agree'),end

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

阳台7777

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

Agnes2_3_4_5_6

不错，帮顶！！~

sunying52zs

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

wangyslyh

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

zjlsmrl

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

yjh2010

到底应用在哪方面啊

yamarychen

回复 木辰月生 的帖子
  S5 Q6 H& o" d- j4 A: D6 a# E" |# L
  f& x$ }" q1 t+ G) \0 c顶啊顶~~
0 m7 E: d1 ?2 y7 F! }+ w
   

meseeme

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

pcm

留个脚印```````