Gauss顺序消去法

木辰月生

2. %M函数gauss.m
   
3. ) V. E2 A! L* M
4. function x=gauss(a,b)
   
5. : |! Y9 `' V, s0 i
6. [n,m]=size(b);
   
7. 9 ~& a! j! O3 I) {$ [* U5 s2 ?
8. if size(a)==[n,n]
   
9. 6 b6 i+ b( p& @8 H; n) G) x6 y
10. for i=1:(n-1)
    
11. # o1 U9 X1 A6 m% Q8 L- R
12. if a(i,i)==0,error('Divide by zero');end\\" @! k6 o# z8 O3 \
    
14. b((i+1):n,:)=b((i+1):n,:)-a((i+1):n,i)*b(i,:)/a(i,i);
    
15. * d& G- k7 K* l: s' E2 l4 Y
16. a((i+1):n,i:n)=a((i+1):n,i:n)-... / ^( V* [, P; O8 C& `/ @
    
18. a((i+1):n,i)*a(i,i:n)/a(i,i);[a,b]
    
19. / v* w, z' [( p0 _, M1 N
20. end  S) Y% b* h/ a
    
22. if a(n,n)==0,error('Divide by zero');end# M9 q- v% y8 j) \% e1 j# `
    
24. x(n,:)=b(n,:)/a(n,n);! I# ]0 \# e$ Q% |\\" J
    
26. for i=n-1:-1:1
    
27. ) ]) e, x% N% @+ D, ]. s6 K
28. x(i,:)=(b(i,:)-a(i,(i+1):n)*x((i+1):n,:))/a(i,i);7 {0 n2 K; S0 `) L
    
30. end
    
31. 5 a: K: _9 i3 B3 C+ ^' a9 [2 K
32. else,error('Dimension of matrix must agree'),end

%M函数gauss.m<br /> <span style=\\"display:none\\">) V. E2 A! L* M</span> function x=gauss(a,b)<br /> <span style=\\"display:none\\">: |! Y9 `' V, s0 i</span> [n,m]=size(b);<br /> <span style=\\"display:none\\">9 ~& a! j! O3 I) {$ [* U5 s2 ?</span> if size(a)==[n,n]<br /> <span style=\\"display:none\\">6 b6 i+ b( p& @8 H; n) G) x6 y</span> for i=1:(n-1)<br /> <span style=\\"display:none\\"># o1 U9 X1 A6 m% Q8 L- R</span> if a(i,i)==0,error('Divide by zero');end<font class=\\"jammer\\">\\" @! k6 o# z8 O3 \</font><br /> b((i+1):n,:)=b((i+1):n,:)-a((i+1):n,i)*b(i,:)/a(i,i);<br /> <span style=\\"display:none\\">* d& G- k7 K* l: s' E2 l4 Y</span> a((i+1):n,i:n)=a((i+1):n,i:n)-... <font class=\\"jammer\\">/ ^( V* [, P; O8 C& `/ @</font><br /> a((i+1):n,i)*a(i,i:n)/a(i,i);[a,b]<br /> <span style=\\"display:none\\">/ v* w, z' [( p0 _, M1 N</span> end<font class=\\"jammer\\">  S) Y% b* h/ a</font><br /> if a(n,n)==0,error('Divide by zero');end<font class=\\"jammer\\"># M9 q- v% y8 j) \% e1 j# `</font><br /> x(n,:)=b(n,:)/a(n,n);<font class=\\"jammer\\">! I# ]0 \# e$ Q% |\\" J</font><br /> for i=n-1:-1:1<br /> <span style=\\"display:none\\">) ]) e, x% N% @+ D, ]. s6 K</span> x(i,:)=(b(i,:)-a(i,(i+1):n)*x((i+1):n,:))/a(i,i);<font class=\\"jammer\\">7 {0 n2 K; S0 `) L</font><br /> end<br /> <span style=\\"display:none\\">5 a: K: _9 i3 B3 C+ ^' a9 [2 K</span> else,error('Dimension of matrix must agree'),end

7 g: X* ~" o: ]计算方法中的有关Gauss消去法的M文件程序代码

柏厷ふ

好啊好啊，非常好

木辰月生

回复 mjw3918 的帖子
& }. X2 \, C6 ]2 E  i
( U1 s0 ?" @- ]& P. `看不懂难道就说别人有病啊··那毕加索的画你看得懂吗
   

SDHZTAO

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

不学无术

哦～～

qian7981

鉴定完毕！  

mjw3918

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

vicky9269

哦～～

wzzh

鉴定完毕！  

dtlxblxb2

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