matlab中求解线性方程组的源代码

模糊数学在工程技术、管理科学、金融工程等领域应用中的很多问题都可以用模糊方程和模糊线性系统来描述。 但是,实现模糊方程和模糊线性系统的求解十分困难,对求解方法的研究一直以来都是重点,也是难点。 无论从理论研究还是从实际应用的角度来说,对模糊方程和模糊线性系统的求解研究都具有重要意义。 本文针对传统方法求解模糊方程和模糊线性系统在模糊数运算、隶属函数解析表示、模糊解判定等方面存在的困难,借助模糊结构元理论,相应地提出了一套模糊方程和模糊线性系统的求解方法。首先,利用两个单调函数的自反单调变换构造了等式限定算子,推广了等式限定运算,处理了存在负模糊情况下关于乘法运算的不可逆问题。 并将等式限定运算思想应用到求解模糊线性方程中,给出了模糊解的结构元表示方法和解存在的充要条件。同时,推广了模糊线性方程,研究了更一般的双重模糊线性方程。此外,还研究了关于矩形复模糊数和圆楔形复模糊数线性方程的求解问题。 其次,定义了幂模糊数和幂模糊数方程,基于结构元方法研究了幂模糊数运算和幂模糊数方程的求解。同时,实现了一元二次模糊方程的求解,利用区间[-1,1]上的单调函数将一元二次模糊方程的求解问题转化为二元二次参数方程组的求解问题,给出了二次模糊方程解存在的充要条件,并辅以数值例子。 最后,利用结构元技术提出了模糊线性系统的求解方法,给出了模糊解存在的充要条件,并辅以实例计算。由于该求解方法是借助[-1,1]上关于y轴对称的单调函数实现的,结果表明在解存在的判定上优于Embedding法。 同时,管理毕业论文www.yifanglunwen.com [-1,1]还研究了一类由模糊结构元线性生成的模糊线性系统,其求解特点是可转为经典线性系统,避免了参数的讨论。本文提出的模糊方程和模糊线性系统的结构元求解方法,极大地简化了模糊数运算的困难,实现了模糊解的判定和解析表达,为模糊数学基础理论问题的研究以及实际问题中的应用与推广奠定了基础。

function [xN]= BGS（Abx0depsM）
if nargin==4
    eps= 1.0e-6;
    M  = 10000;
elseif nargin<4
    error
    return
elseif nargin ==5
    M  = 10000;
end  

NS = size（A）;
n = NS（11）;
bnum = length（d）;
bs = ones（bnum1）;
for i=1:（bnum-1）
    bs（i+11）=sum（d（1:i））+1;
    %获得对角线上每个分块矩阵元素索引的起始值
end

DB = zeros（nn）;       
for i=1:bnum
    endb = bs（i1）+d（i1）-1;
    DB（bs（i1）:endbbs（i1）:endb）=A（bs（i1）:endbbs（i1）:endb）;
    %求A的对角分块矩阵
end
LB = -tril（A-DB）;      %求A的下三角分块阵
UB = -triu（A-DB）;      %求A的上三角分块阵

N = 0;
tol = 1;
while tol>=eps 
    invDL = inv（DB-LB）;
    x = invDL*UB*x0+invDL*b;    %块迭代公式
    N = N+1;
    tol = norm（x-x0）;
    x0 = x;
    if（N>=M）
        disp（‘Warning: 迭代次数太多，可能不收敛！‘）;
        return;
    end
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2012-11-08 19:59  解线性方程组的迭代法\
     文件         846  2008-04-10 00:14  解线性方程组的迭代法\BGS.m
     文件        1092  2008-04-10 00:13  解线性方程组的迭代法\BJ.m
     文件         884  2008-04-10 00:14  解线性方程组的迭代法\BSOR.m
     文件         521  2008-04-10 00:11  解线性方程组的迭代法\JOR.m
     文件         579  2006-11-04 17:53  解线性方程组的迭代法\SOR.m
     文件         679  2006-11-04 18:10  解线性方程组的迭代法\SSOR.m
     文件         447  2008-04-10 00:12  解线性方程组的迭代法\conjgrad.m
     文件         537  2008-04-10 00:08  解线性方程组的迭代法\crs.m
     文件         289  2008-04-10 00:12  解线性方程组的迭代法\fastdown.m
     文件         510  2006-11-04 14:10  解线性方程组的迭代法\gauseidel.m
     文件         624  2008-04-10 20:52  解线性方程组的迭代法\grs.m
     文件         521  2006-11-04 14:00  解线性方程组的迭代法\jacobi.m
     文件         458  2008-04-10 00:13  解线性方程组的迭代法\preconjgrad.m
     文件         348  2006-11-04 14:00  解线性方程组的迭代法\richason.m
     文件         535  2008-04-10 00:08  解线性方程组的迭代法\rs.m
     文件         565  2006-11-04 18:24  解线性方程组的迭代法\twostep.m