Gauss顺序消去法

木辰月生

2. %M函数gauss.m
   
3. * q\\" d: ?% u1 g' r0 c9 N
4. function x=gauss(a,b)
   
5. # c; `8 l8 ]\\" _3 K3 p6 O( @
6. [n,m]=size(b);
   
7. 8 W# K. u. |) o! Z6 W5 X8 f0 p
8. if size(a)==[n,n]: O, I! _6 M\\" N* i
   
10. for i=1:(n-1)
    
11. ( M4 J$ R/ ?& o& I; S
12. if a(i,i)==0,error('Divide by zero');end: d$ }\\" W\\" b; ]7 g
    
14. b((i+1):n,:)=b((i+1):n,:)-a((i+1):n,i)*b(i,:)/a(i,i);$ g\\" `8 _; |% x\\" y, g' m& ]9 w
    
16. a((i+1):n,i:n)=a((i+1):n,i:n)-... / K6 K- {2 Z0 x: {. I1 g
    
18. a((i+1):n,i)*a(i,i:n)/a(i,i);[a,b]* w# M( ?6 q+ x$ R, j3 Z
    
20. end
    
21. . t' |- Q5 j3 Y3 I% X' S
22. if a(n,n)==0,error('Divide by zero');end
    
23. 4 C; f0 o  W6 R4 X& w6 b\\" k0 A
24. x(n,:)=b(n,:)/a(n,n);& E/ x9 R$ |+ V. T  W
    
26. for i=n-1:-1:1# I; T' E3 ~$ b+ P( P
    
28. x(i,:)=(b(i,:)-a(i,(i+1):n)*x((i+1):n,:))/a(i,i);5 @) b0 _* J: {1 e
    
30. end/ F1 S4 j+ H+ h6 I* H
    
32. else,error('Dimension of matrix must agree'),end

%M函数gauss.m<br /> <span style=\\"display:none\\">* q\\" d: ?% u1 g' r0 c9 N</span> function x=gauss(a,b)<br /> <span style=\\"display:none\\"># c; `8 l8 ]\\" _3 K3 p6 O( @</span> [n,m]=size(b);<br /> <span style=\\"display:none\\">8 W# K. u. |) o! Z6 W5 X8 f0 p</span> if size(a)==[n,n]<font class=\\"jammer\\">: O, I! _6 M\\" N* i</font><br /> for i=1:(n-1)<br /> <span style=\\"display:none\\">( M4 J$ R/ ?& o& I; S</span> if a(i,i)==0,error('Divide by zero');end<font class=\\"jammer\\">: d$ }\\" W\\" b; ]7 g</font><br /> b((i+1):n,:)=b((i+1):n,:)-a((i+1):n,i)*b(i,:)/a(i,i);<font class=\\"jammer\\">$ g\\" `8 _; |% x\\" y, g' m& ]9 w</font><br /> a((i+1):n,i:n)=a((i+1):n,i:n)-... <font class=\\"jammer\\">/ K6 K- {2 Z0 x: {. I1 g</font><br /> a((i+1):n,i)*a(i,i:n)/a(i,i);[a,b]<font class=\\"jammer\\">* w# M( ?6 q+ x$ R, j3 Z</font><br /> end<br /> <span style=\\"display:none\\">. t' |- Q5 j3 Y3 I% X' S</span> if a(n,n)==0,error('Divide by zero');end<br /> <span style=\\"display:none\\">4 C; f0 o  W6 R4 X& w6 b\\" k0 A</span> x(n,:)=b(n,:)/a(n,n);<font class=\\"jammer\\">& E/ x9 R$ |+ V. T  W</font><br /> for i=n-1:-1:1<font class=\\"jammer\\"># I; T' E3 ~$ b+ P( P</font><br /> x(i,:)=(b(i,:)-a(i,(i+1):n)*x((i+1):n,:))/a(i,i);<font class=\\"jammer\\">5 @) b0 _* J: {1 e</font><br /> end<font class=\\"jammer\\">/ F1 S4 j+ H+ h6 I* H</font><br /> else,error('Dimension of matrix must agree'),end

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

阳台7777

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

Agnes2_3_4_5_6

不错，帮顶！！~

sunying52zs

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

wangyslyh

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

zjlsmrl

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

yjh2010

到底应用在哪方面啊

yamarychen

回复 木辰月生 的帖子
3 R* r3 C! P2 b: h
- ~! x8 ]0 ^$ r. N% ?顶啊顶~~

   

meseeme

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

pcm

留个脚印```````