Gauss顺序消去法

木辰月生

2. %M函数gauss.m+ v! R. J, k$ H0 O3 U& Y1 y
   
4. function x=gauss(a,b)
   
5. / j/ P9 n' D( u4 S/ E
6. [n,m]=size(b);
   
7. ' S9 A* j; [; A0 a$ x
8. if size(a)==[n,n]7 R% S0 J4 M# S) ?9 |
   
10. for i=1:(n-1)
    
11. 4 F\\" @9 L8 p0 Z' q4 Q( @6 [+ `0 t
12. if a(i,i)==0,error('Divide by zero');end
    
13. 0 z  e, M\\" a\\" a# Y/ X) `5 q8 r- g, w
14. b((i+1):n,:)=b((i+1):n,:)-a((i+1):n,i)*b(i,:)/a(i,i);
    
15. , h' U) j\\" U\\" @9 t+ M3 D
16. a((i+1):n,i:n)=a((i+1):n,i:n)-...
    
17. ( J& [, f- o+ G) o
18. a((i+1):n,i)*a(i,i:n)/a(i,i);[a,b]
    
19. + o0 y) ?1 a1 N  {) p! g
20. end9 m9 h2 ~. Q. p; b8 y
    
22. if a(n,n)==0,error('Divide by zero');end+ _( M( B4 y- v
    
24. x(n,:)=b(n,:)/a(n,n);9 L8 B0 \\\" b5 ]% t* Q* h. g
    
26. for i=n-1:-1:1
    
27. ; W0 ~0 N6 g7 U\\" n9 O( `
28. x(i,:)=(b(i,:)-a(i,(i+1):n)*x((i+1):n,:))/a(i,i);
    
29. . A/ E. t. S3 P
30. end
    
31. 4 ]! p7 T) ]1 G, l
32. else,error('Dimension of matrix must agree'),end

%M函数gauss.m<font class=\\"jammer\\">+ v! R. J, k$ H0 O3 U& Y1 y</font><br /> function x=gauss(a,b)<br /> <span style=\\"display:none\\">/ j/ P9 n' D( u4 S/ E</span> [n,m]=size(b);<br /> <span style=\\"display:none\\">' S9 A* j; [; A0 a$ x</span> if size(a)==[n,n]<font class=\\"jammer\\">7 R% S0 J4 M# S) ?9 |</font><br /> for i=1:(n-1)<br /> <span style=\\"display:none\\">4 F\\" @9 L8 p0 Z' q4 Q( @6 [+ `0 t</span> if a(i,i)==0,error('Divide by zero');end<br /> <span style=\\"display:none\\">0 z  e, M\\" a\\" a# Y/ X) `5 q8 r- g, w</span> b((i+1):n,:)=b((i+1):n,:)-a((i+1):n,i)*b(i,:)/a(i,i);<br /> <span style=\\"display:none\\">, h' U) j\\" U\\" @9 t+ M3 D</span> a((i+1):n,i:n)=a((i+1):n,i:n)-... <br /> <span style=\\"display:none\\">( J& [, f- o+ G) o</span> a((i+1):n,i)*a(i,i:n)/a(i,i);[a,b]<br /> <span style=\\"display:none\\">+ o0 y) ?1 a1 N  {) p! g</span> end<font class=\\"jammer\\">9 m9 h2 ~. Q. p; b8 y</font><br /> if a(n,n)==0,error('Divide by zero');end<font class=\\"jammer\\">+ _( M( B4 y- v</font><br /> x(n,:)=b(n,:)/a(n,n);<font class=\\"jammer\\">9 L8 B0 \\\" b5 ]% t* Q* h. g</font><br /> for i=n-1:-1:1<br /> <span style=\\"display:none\\">; W0 ~0 N6 g7 U\\" n9 O( `</span> x(i,:)=(b(i,:)-a(i,(i+1):n)*x((i+1):n,:))/a(i,i);<br /> <span style=\\"display:none\\">. A/ E. t. S3 P</span> end<br /> <span style=\\"display:none\\">4 ]! p7 T) ]1 G, l</span> else,error('Dimension of matrix must agree'),end

: a8 f- _2 p- r: y% X计算方法中的有关Gauss消去法的M文件程序代码

柏厷ふ

好啊好啊，非常好

木辰月生

回复 mjw3918 的帖子
( h+ x% q7 Z" L& f# i
8 V  v6 E* t. T) m& e8 C看不懂难道就说别人有病啊··那毕加索的画你看得懂吗
   

SDHZTAO

楼主，你写得实在是太好了。我惟一能做的，就只有把这个帖子顶上去这件事了

不学无术

哦～～

qian7981

鉴定完毕！  

mjw3918

我基本上是采用看英语文章的办法，先泛读，再精读，再一句一句看，最后再提纲挈领，总算是明白一点了，当然，也可能还是领悟错了。最后要说的一句话是：楼主，你很牛叉，希望你不是真的有病。   

vicky9269

哦～～

wzzh

鉴定完毕！  

dtlxblxb2

我回不抢呢 考虑再三 还是不抢了吧 ^_^