怎样编制能自动选取主元又能手动选取主元的求解线性方程组的Gauss消去过程

jan52

怎样编制一个能自动选取主元，又能手动选取主元的求解线性
方程组的Gauss消去过程?
麻烦写个程序!

[ 本帖最后由 lxq 于 2006-9-28 21:23 编辑 ]

AaronSpark

1. function x=gauss(A,b)
   
2. %x=gauss(A,b)
   
3. n=length(A);
   
4. a=[A,b];
   
5. for k=1:n-1
   
6.     maxa=max(abs(a(k:n,k)));
   
7.     if maxa==0
   
8.         return;
   
9.     end
   
10.     for i=k:n
    
11.         if abs(a(i,k))==maxa
    
12.             y=a(i,k:n+1);a(i,k:n+1)=a(k,k:n+1);a(k,k:n+1)=y;
    
13.             break;
    
14.         end
    
15.     end
    
16.     for i=k+1:n
    
17.         l(i,k)=a(i,k)/a(k,k);
    
18.         a(i,k+1:n+1)=a(i,k+1:n+1)-l(i,k).*a(k,k+1:n+1);
    
19.     end
    
20. end
    
21. 
    
22. %回代
    
23. if a(n,n)==0
    
24.     return
    
25. end
    
26. x(n)=a(n,n+1)/a(n,n);
    
27. for i=n-1:-1:1
    
28.     x(i)=(a(i,n+1)-sum(a(i,i+1:n).*x(i+1:n)))/a(i,i);
    
29. end

复制代码