maker系列--最速下降法求线性方程

引用 MVH 2005-7-28 17:50

介绍一种求多元方程组的方法
come from shumo.com
%简单的GAUSS消去法
%a为系数方阵，b为右端向量

1. function x = gauss(a,b)
   
2. [n,m] = size(b);
   
3. if size(a) == [n,n]
   
4.     for i = 1 : (n - 1)
   
5.         if a(i,i) == 0 error('Divide by zero');end
   
6.         b((i + 1):n,:) = b((i + 1):n,:) - a((i + 1):n,i)*b(i,:)/a(i,i);
   
7.         a((i + 1):n,i:n) = a((i + 1):n,i:n) - ...
   
8.             a((i + 1):n,i)*a(i,i:n)/a(i,i);
   
9.         [a,b]
   
10.     end
    
11.     if a(n,n) == 0,error('Divide by zero');end
    
12.     x(n,:) = b(n,:)/a(n,n);
    
13.     for i = n-1:-1:1
    
14.         x(i,:) = b(i,:) - a(i,(i + 1):n)*x((i + 1):n,:)/a(i,i);
    
15.     end
    
16. else,error('Divede of matrix must agree');
    
17. end

复制代码

%请大家赐教其他的好方法

引用 xinyuan110 2006-5-23 11:11

不错，不错谢谢了！！！

查看全部评论(2)