function [x,N]= BJ(A,b,x0,d,eps,M)
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(1,1);
if(sum(d) ~= n)
disp('分块错误!');
return;
end
bnum = length(d);
bs = ones(bnum,1);
for i=1:(bnum-1)
bs(i+1,1)=sum(d(1:i))+1;
%获得对角线上每个分块矩阵元素索引的起始值
end
DB = zeros(n,n);
for i=1:bnum
endb = bs(i,1)+d(i,1)-1;
DB(bs(i,1):endb,bs(i,1):endb)=A(bs(i,1):endb,bs(i,1):endb);
%求A的对角分块矩阵
end
for i=1:bnum
endb = bs(i,1)+d(i,1)-1;
invDB(bs(i,1):endb,bs(i,1):endb)=inv(DB(bs(i,1):endb,bs(i,1):endb));
%求A的对角分块矩阵的逆矩阵
end
N = 0;
tol = 1;
while tol>=eps
x = invDB*(DB-A)*x0+invDB*b; %由于LB+DB=DB-A
N = N+1; %迭代步数
tol = norm(x-x0); %前后两步迭代结果的误差
x0 = x;
if(N>=M)
disp('Warning: 迭代次数太多,可能不收敛!');
return;
end
end
等天晴i
- 粉丝: 5935
- 资源: 10万+
最新资源
- 一对一MybatisProgram.zip
- matlab实现基于多源信息融合的水果虫灾智能检测系统研究-虫灾检测-智能检测系统-matlab-机器学习
- Docker从零走向实战视频(上).zip
- 五相电机邻近四矢量SVPWM模型-MATLAB-Simulink仿真模型包括: (1)原理说明文档(重要):包括扇区判断、矢量作用时间计算、矢量作用顺序及切时间计算、PWM波的生成; (2)输出部分仿
- 保险箱检测51-YOLO(v5至v11)、COCO、CreateML、Paligemma、TFRecord、VOC数据集合集.rar
- simpack软件与ansys,abqus联合仿真求解车桥耦合,地震波浪荷载联合仿真分析,全教程模型 1. abaqus-simpack车轨耦合振动分析 2. abaqus-simpack车轨桥耦合
- 单电阻采样 基于单电阻采样的相电流重构算法 keil完整工程 单电阻采样 f103的单电阻,完整工程,带文档,带硬件资料 f3平台的单电阻完整工程,代码详细注释 还有微芯的单电阻smo代码加文档
- IP102中分离出来的害虫数据集,使用yolov9标注
- 华为OD+考试真题+实现过程
- IP102中分离出来的害虫数据集,使用yolov8标注
- 充电器检测9-YOLO(v5至v11)、COCO、Create充电器检测9L、Paligemma、TFRecord、VOC数据集合集.rar
- 基于stm32的通信系统,sim800c与服务器通信,无线通信监测,远程定位,服务器通信系统,gps,sim800c,心率,温度,stm32 由STM32F103ZET6单片机核心板电路、DS18B2
- jQuery信息提示插件
- Docker从零走向实战视频(下).zip
- 华为OD+真题及解析+智能驾驶
- 程序名称:悬架设计计算程序 开发平台:基于matlab平台 计算内容:悬架偏频刚度挠度;螺旋弹簧,多片簧,少片簧,稳定杆,减震器的匹配计算;悬架垂向纵向侧向力学、纵倾、侧倾校核等;独立悬架杠杆比,等效
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈