Gauss顺序消去法

木辰月生

2. %M函数gauss.m
   
3. & \9 [  u) R  B: t0 @
4. function x=gauss(a,b)- J) u5 }, ~1 x6 w
   
6. [n,m]=size(b);
   
7. 2 C. y; Z1 F6 o) W* f\\" I6 n$ b$ S
8. if size(a)==[n,n]
   
9. 4 y3 h) t, b. D& _
10. for i=1:(n-1)
    
11.   x( K4 G8 p) r/ a$ D0 B# {
12. if a(i,i)==0,error('Divide by zero');end4 h3 @3 t' g% ]! g0 i
    
14. b((i+1):n,:)=b((i+1):n,:)-a((i+1):n,i)*b(i,:)/a(i,i);
    
15. : d3 ]  b7 T! J& j( ?
16. a((i+1):n,i:n)=a((i+1):n,i:n)-... ( h& T. S) Y( H3 I3 S+ F
    
18. a((i+1):n,i)*a(i,i:n)/a(i,i);[a,b]( t; u& A; ^$ ~* b9 K\\" e& Y' y8 _
    
20. end
    
21. 3 ^7 A- I7 m8 ^  O4 v' y- c7 w
22. if a(n,n)==0,error('Divide by zero');end; y# p- L8 L; y' f: W0 Z/ ?
    
24. x(n,:)=b(n,:)/a(n,n);
    
25. , B! Q, R( Z$ ?\\" S7 N
26. for i=n-1:-1:1
    
27. 7 ^# J$ k) P6 e1 ^3 H* G0 F
28. x(i,:)=(b(i,:)-a(i,(i+1):n)*x((i+1):n,:))/a(i,i);3 S8 A9 W9 E3 t1 C  D( B( w$ Y
    
30. end: x% m& z/ G$ G
    
32. else,error('Dimension of matrix must agree'),end

%M函数gauss.m<br /> <span style=\\"display:none\\">& \9 [  u) R  B: t0 @</span> function x=gauss(a,b)<font class=\\"jammer\\">- J) u5 }, ~1 x6 w</font><br /> [n,m]=size(b);<br /> <span style=\\"display:none\\">2 C. y; Z1 F6 o) W* f\\" I6 n$ b$ S</span> if size(a)==[n,n]<br /> <span style=\\"display:none\\">4 y3 h) t, b. D& _</span> for i=1:(n-1)<br /> <span style=\\"display:none\\">  x( K4 G8 p) r/ a$ D0 B# {</span> if a(i,i)==0,error('Divide by zero');end<font class=\\"jammer\\">4 h3 @3 t' g% ]! g0 i</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\\">: d3 ]  b7 T! J& j( ?</span> a((i+1):n,i:n)=a((i+1):n,i:n)-... <font class=\\"jammer\\">( h& T. S) Y( H3 I3 S+ F</font><br /> a((i+1):n,i)*a(i,i:n)/a(i,i);[a,b]<font class=\\"jammer\\">( t; u& A; ^$ ~* b9 K\\" e& Y' y8 _</font><br /> end<br /> <span style=\\"display:none\\">3 ^7 A- I7 m8 ^  O4 v' y- c7 w</span> if a(n,n)==0,error('Divide by zero');end<font class=\\"jammer\\">; y# p- L8 L; y' f: W0 Z/ ?</font><br /> x(n,:)=b(n,:)/a(n,n);<br /> <span style=\\"display:none\\">, B! Q, R( Z$ ?\\" S7 N</span> for i=n-1:-1:1<br /> <span style=\\"display:none\\">7 ^# J$ k) P6 e1 ^3 H* G0 F</span> x(i,:)=(b(i,:)-a(i,(i+1):n)*x((i+1):n,:))/a(i,i);<font class=\\"jammer\\">3 S8 A9 W9 E3 t1 C  D( B( w$ Y</font><br /> end<font class=\\"jammer\\">: x% m& z/ G$ G</font><br /> else,error('Dimension of matrix must agree'),end

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

阳台7777

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

Agnes2_3_4_5_6

不错，帮顶！！~

sunying52zs

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

wangyslyh

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

zjlsmrl

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

yjh2010

到底应用在哪方面啊

yamarychen

回复 木辰月生 的帖子

1 `/ i) Y9 n3 ^3 a$ i$ W0 R! k顶啊顶~~

   

meseeme

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

pcm

留个脚印```````