Gauss顺序消去法

木辰月生

2. %M函数gauss.m
   
3. 6 E3 Y  n& f) o
4. function x=gauss(a,b)
   
5. 0 `! n7 x. G' E  \4 N( {
6. [n,m]=size(b);
   
7. 5 y% R% [! ?6 i8 S9 M. g' g
8. if size(a)==[n,n]/ n6 o0 t/ d+ H4 Z! ^
   
10. for i=1:(n-1)
    
11. * R2 X/ b! }* r6 E( x- }% b
12. if a(i,i)==0,error('Divide by zero');end. D0 I6 R; u/ l1 h& x
    
14. b((i+1):n,:)=b((i+1):n,:)-a((i+1):n,i)*b(i,:)/a(i,i);\\" x  U/ V3 T( O& u
    
16. a((i+1):n,i:n)=a((i+1):n,i:n)-... - ^( P\\" t! c\\" c2 E% L$ o4 `
    
18. a((i+1):n,i)*a(i,i:n)/a(i,i);[a,b]: N7 o) i8 R, k* i
    
20. end
    
21. / S: r9 X* t4 M( Y0 I$ u6 c0 V( }
22. if a(n,n)==0,error('Divide by zero');end; y0 [+ |1 H$ {
    
24. x(n,:)=b(n,:)/a(n,n);
    
25. 6 s) l+ H2 {7 I# {* b/ C
26. for i=n-1:-1:1
    
27. ! O2 }; n\\" F7 @7 u# F
28. x(i,:)=(b(i,:)-a(i,(i+1):n)*x((i+1):n,:))/a(i,i);
    
29. / H7 S& J7 f5 i+ B: I
30. end& `) |; \0 ^- }+ K2 h& }% b
    
32. else,error('Dimension of matrix must agree'),end

%M函数gauss.m<br /> <span style=\\"display:none\\">6 E3 Y  n& f) o</span> function x=gauss(a,b)<br /> <span style=\\"display:none\\">0 `! n7 x. G' E  \4 N( {</span> [n,m]=size(b);<br /> <span style=\\"display:none\\">5 y% R% [! ?6 i8 S9 M. g' g</span> if size(a)==[n,n]<font class=\\"jammer\\">/ n6 o0 t/ d+ H4 Z! ^</font><br /> for i=1:(n-1)<br /> <span style=\\"display:none\\">* R2 X/ b! }* r6 E( x- }% b</span> if a(i,i)==0,error('Divide by zero');end<font class=\\"jammer\\">. D0 I6 R; u/ l1 h& x</font><br /> b((i+1):n,:)=b((i+1):n,:)-a((i+1):n,i)*b(i,:)/a(i,i);<font class=\\"jammer\\">\\" x  U/ V3 T( O& u</font><br /> a((i+1):n,i:n)=a((i+1):n,i:n)-... <font class=\\"jammer\\">- ^( P\\" t! c\\" c2 E% L$ o4 `</font><br /> a((i+1):n,i)*a(i,i:n)/a(i,i);[a,b]<font class=\\"jammer\\">: N7 o) i8 R, k* i</font><br /> end<br /> <span style=\\"display:none\\">/ S: r9 X* t4 M( Y0 I$ u6 c0 V( }</span> if a(n,n)==0,error('Divide by zero');end<font class=\\"jammer\\">; y0 [+ |1 H$ {</font><br /> x(n,:)=b(n,:)/a(n,n);<br /> <span style=\\"display:none\\">6 s) l+ H2 {7 I# {* b/ C</span> for i=n-1:-1:1<br /> <span style=\\"display:none\\">! O2 }; n\\" F7 @7 u# F</span> x(i,:)=(b(i,:)-a(i,(i+1):n)*x((i+1):n,:))/a(i,i);<br /> <span style=\\"display:none\\">/ H7 S& J7 f5 i+ B: I</span> end<font class=\\"jammer\\">& `) |; \0 ^- }+ K2 h& }% b</font><br /> else,error('Dimension of matrix must agree'),end

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

阳台7777

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

Agnes2_3_4_5_6

不错，帮顶！！~

sunying52zs

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

wangyslyh

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

zjlsmrl

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

yjh2010

到底应用在哪方面啊

yamarychen

回复 木辰月生 的帖子
' H6 O. |4 L/ ?& l# y% X
顶啊顶~~

- U6 w6 q- [$ Q$ T/ D/ J: [   

meseeme

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

pcm

留个脚印```````