Gauss顺序消去法

木辰月生

2. %M函数gauss.m
   
3. 9 r# M$ i; S& l5 s
4. function x=gauss(a,b)
   
5. : z  r; R; I$ ~; t2 G
6. [n,m]=size(b);
   
7. 3 @- x# j0 H: W, E: P  u
8. if size(a)==[n,n]
   
9. ( J- J! X2 z/ U3 g/ e6 F1 o
10. for i=1:(n-1)
    
11. % o' `\\" k' `$ e
12. if a(i,i)==0,error('Divide by zero');end: G: q* m/ ^+ f% x
    
14. b((i+1):n,:)=b((i+1):n,:)-a((i+1):n,i)*b(i,:)/a(i,i);
    
15. ' s0 ~% k' a! ~* ?/ N7 f
16. a((i+1):n,i:n)=a((i+1):n,i:n)-...
    
17.   ]7 w4 T- z\\" m5 O+ c1 x2 C
18. a((i+1):n,i)*a(i,i:n)/a(i,i);[a,b]$ ?4 m! Y! A- B' P% ~# M
    
20. end
    
21. 0 g# f+ v5 r8 m
22. if a(n,n)==0,error('Divide by zero');end
    
23. - [1 S1 f3 p: G
24. x(n,:)=b(n,:)/a(n,n);0 L6 a: [2 p! ?8 c! l9 |5 _
    
26. for i=n-1:-1:1, W0 x, f\\" t* k, N/ _! V. s
    
28. x(i,:)=(b(i,:)-a(i,(i+1):n)*x((i+1):n,:))/a(i,i);
    
29. 8 K! V4 p: W\\" A2 N
30. end
    
31. % F1 ^3 G- u2 q0 [6 a, r( \8 _
32. else,error('Dimension of matrix must agree'),end

%M函数gauss.m<br /> <span style=\\"display:none\\">9 r# M$ i; S& l5 s</span> function x=gauss(a,b)<br /> <span style=\\"display:none\\">: z  r; R; I$ ~; t2 G</span> [n,m]=size(b);<br /> <span style=\\"display:none\\">3 @- x# j0 H: W, E: P  u</span> if size(a)==[n,n]<br /> <span style=\\"display:none\\">( J- J! X2 z/ U3 g/ e6 F1 o</span> for i=1:(n-1)<br /> <span style=\\"display:none\\">% o' `\\" k' `$ e</span> if a(i,i)==0,error('Divide by zero');end<font class=\\"jammer\\">: G: q* m/ ^+ f% x</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\\">' s0 ~% k' a! ~* ?/ N7 f</span> a((i+1):n,i:n)=a((i+1):n,i:n)-... <br /> <span style=\\"display:none\\">  ]7 w4 T- z\\" m5 O+ c1 x2 C</span> a((i+1):n,i)*a(i,i:n)/a(i,i);[a,b]<font class=\\"jammer\\">$ ?4 m! Y! A- B' P% ~# M</font><br /> end<br /> <span style=\\"display:none\\">0 g# f+ v5 r8 m</span> if a(n,n)==0,error('Divide by zero');end<br /> <span style=\\"display:none\\">- [1 S1 f3 p: G</span> x(n,:)=b(n,:)/a(n,n);<font class=\\"jammer\\">0 L6 a: [2 p! ?8 c! l9 |5 _</font><br /> for i=n-1:-1:1<font class=\\"jammer\\">, W0 x, f\\" t* k, N/ _! V. s</font><br /> x(i,:)=(b(i,:)-a(i,(i+1):n)*x((i+1):n,:))/a(i,i);<br /> <span style=\\"display:none\\">8 K! V4 p: W\\" A2 N</span> end<br /> <span style=\\"display:none\\">% F1 ^3 G- u2 q0 [6 a, r( \8 _</span> else,error('Dimension of matrix must agree'),end

5 t7 i4 ]' i! L" m& e2 b- t计算方法中的有关Gauss消去法的M文件程序代码

阳台7777

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

Agnes2_3_4_5_6

不错，帮顶！！~

sunying52zs

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

wangyslyh

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

zjlsmrl

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

yjh2010

到底应用在哪方面啊

yamarychen

回复 木辰月生 的帖子

顶啊顶~~

   

meseeme

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

pcm

留个脚印```````