Gauss顺序消去法

木辰月生

2. %M函数gauss.m; f% x; y- a! j. B$ q6 ]
   
4. function x=gauss(a,b)
   
5. % C2 H; ^, o) g* D4 u: N* b3 R
6. [n,m]=size(b);) R+ }* |- x( e4 O0 @
   
8. if size(a)==[n,n]: M/ }\\" D* X& |7 O
   
10. for i=1:(n-1)
    
11. + G5 q& a5 \\\" l* w  s/ {
12. if a(i,i)==0,error('Divide by zero');end
    
13. + d4 ?; N- a: B% W\\" N
14. b((i+1):n,:)=b((i+1):n,:)-a((i+1):n,i)*b(i,:)/a(i,i);
    
15. / P7 C7 r+ p5 T% q, j6 C1 c
16. a((i+1):n,i:n)=a((i+1):n,i:n)-...
    
17. , `1 U9 \6 m5 @1 h
18. a((i+1):n,i)*a(i,i:n)/a(i,i);[a,b]
    
19. + W3 K7 v' T9 `: }
20. end
    
21. $ I, l! m( o. c, W
22. if a(n,n)==0,error('Divide by zero');end
    
23. 9 b% }7 a, L6 e! E
24. x(n,:)=b(n,:)/a(n,n);
    
25. $ U# o\\" }5 y! b7 ^% b4 P# d
26. for i=n-1:-1:1* w7 Z. A: f8 S\\" C/ t
    
28. x(i,:)=(b(i,:)-a(i,(i+1):n)*x((i+1):n,:))/a(i,i);
    
29. 9 |% h7 _9 D# n. x7 N
30. end
    
31. 7 E* g\\" r) q- V9 a
32. else,error('Dimension of matrix must agree'),end

%M函数gauss.m<font class=\\"jammer\\">; f% x; y- a! j. B$ q6 ]</font><br /> function x=gauss(a,b)<br /> <span style=\\"display:none\\">% C2 H; ^, o) g* D4 u: N* b3 R</span> [n,m]=size(b);<font class=\\"jammer\\">) R+ }* |- x( e4 O0 @</font><br /> if size(a)==[n,n]<font class=\\"jammer\\">: M/ }\\" D* X& |7 O</font><br /> for i=1:(n-1)<br /> <span style=\\"display:none\\">+ G5 q& a5 \\\" l* w  s/ {</span> if a(i,i)==0,error('Divide by zero');end<br /> <span style=\\"display:none\\">+ d4 ?; N- a: B% W\\" N</span> b((i+1):n,:)=b((i+1):n,:)-a((i+1):n,i)*b(i,:)/a(i,i);<br /> <span style=\\"display:none\\">/ P7 C7 r+ p5 T% q, j6 C1 c</span> a((i+1):n,i:n)=a((i+1):n,i:n)-... <br /> <span style=\\"display:none\\">, `1 U9 \6 m5 @1 h</span> a((i+1):n,i)*a(i,i:n)/a(i,i);[a,b]<br /> <span style=\\"display:none\\">+ W3 K7 v' T9 `: }</span> end<br /> <span style=\\"display:none\\">$ I, l! m( o. c, W</span> if a(n,n)==0,error('Divide by zero');end<br /> <span style=\\"display:none\\">9 b% }7 a, L6 e! E</span> x(n,:)=b(n,:)/a(n,n);<br /> <span style=\\"display:none\\">$ U# o\\" }5 y! b7 ^% b4 P# d</span> for i=n-1:-1:1<font class=\\"jammer\\">* w7 Z. A: f8 S\\" C/ t</font><br /> x(i,:)=(b(i,:)-a(i,(i+1):n)*x((i+1):n,:))/a(i,i);<br /> <span style=\\"display:none\\">9 |% h7 _9 D# n. x7 N</span> end<br /> <span style=\\"display:none\\">7 E* g\\" r) q- V9 a</span> else,error('Dimension of matrix must agree'),end

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

柏厷ふ

好啊好啊，非常好

木辰月生

回复 mjw3918 的帖子

看不懂难道就说别人有病啊··那毕加索的画你看得懂吗
   

SDHZTAO

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

不学无术

哦～～

qian7981

鉴定完毕！  

mjw3918

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

vicky9269

哦～～

wzzh

鉴定完毕！  

dtlxblxb2

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