追赶法求解块三对角矩阵以及matlab代码实现

数值分析课的作业

单从题目来看，矩阵A和矩阵B都是三对角矩阵

参考资料：块三对角矩阵方程的追赶法及其应用 - 豆丁网Docin

推导过程如下：

求解matlab程序：

function x=tridiagsolver_block(n)
c=(1:n^2);
for i=1:n,a{i}=c(i:i+n-1);b{i}=a{i}';
end
B=diag(repmat([4], 1, n))+diag(repmat([-1], 1, n-1), 1)+diag(repmat([-1], 1, n-1), -1);%矩阵B
G{1}=B;
D{1}=-inv(G{1});
y{1}=inv(B)*b{1};
for i=2:n-1,G{i}=B+D{i-1};D{i}=-inv(G{i});y{i}=inv(G{i})*(b{i}+y{i-1});
end
G{n}=B+D{n-1};
D{n}=-inv(G{n});
y{n}=inv(G{n})*(b{n}+y{n-1});
x{n}=y{n};
for i=n-1:-1:1,x{i}=y{i}-D{i}*x{i+1};
end

注意，此处要使用cell数据类型将矩阵存储到A中，所以均采用“{}”，最后调用函数输入n就行啦^_^

追赶法求解块三对角矩阵以及matlab代码实现相关推荐

三阶魔方自动求解及动态可视化matlab代码
三阶魔方自动求解及动态可视化matlab代码思路与步骤总结思考参考链接源代码第一次写博客,想总结分享下以前做过的一些有趣的东西,目的是为了回望过去与展望未来,同时为了提高自己的写作表达能力 ...
货物配送问题的matlab,免疫算法求解配送中心选址问题matlab代码
免疫算法求解配送中心选址问题matlab代码所属分类:matlab例程开发工具:matlab 文件大小:31KB 下载次数:3 上传日期:2020-12-01 16:11:00 上传者:代码大 ...
追赶法求解三对角线性方程组的MATLAB程序
function [M] = after_method(a,b,c,g,tolerance) %% 追赶法求解三对角矩阵 % a为三对角矩阵左下对角线上的值 % b为三对角矩阵中间对角线上的值 % c ...
matlab 读取图片后分区域编号_你的第一个有限元求解器——仅十行MATLAB代码
有限元分析话题中有不少讨论有限元求解器的问题,但大都停留在概念层面,未见实际代码.望本文能略起抛砖引玉之作用. 以下代码是基于MATLAB编写. 问题描述考虑一平面有界区域 ,设其边界为 .我们求解 ...
模拟退火法和蚂蚁优化算法求解TSP问题（Matlab代码实现）
目录 1 概述 2 蚂蚁优化算法运行结果 3 模拟退火法求解结果 4 Matlab代码 1 概述旅行商问题被定义为一个推销员在所有城市的旅行﹐以最低的成本回到最初的城市.TSP问题的计算规模随着 ...
用追赶法求方程组c语言程序,数值计算——追赶法求解三对角方程组（附代码）...
目录追赶法基础理论追赶法c++程序代码程序运行结果源码文件下载地址追赶法基础理论在数值计算中,对三次样条曲线插值和用差分方法求解常微分方程边值问题时,通常会遇到Ax=d三对角形式的方程组: ...
【路径优化】基于帝企鹅算法求解TSP问题（Matlab代码实现）
目录 1 帝企鹅算法 2 旅行商问题(TSP) 3 运行结果 4 参考文献 5 Matlab代码实现 1 帝企鹅算法帝企鹅优化算法(emperor penguin optimizer,EPO)是Ga ...
【优化求解】基于灰狼算法GWO求解最优目标matlab代码
1 简介 Mirjalili 等人提出了一种新的群体智能算法---灰狼优化算法(GWO),并通过多个基准测试函数进行测试,从结果上验证了该算法的可行性,通过对比,GWO 算法已被证明在算法对函数求解精 ...
基于改进禁忌搜索算法求解TSP问题（Matlab代码实现）
目录 1 概述 2 改进禁忌搜索算法 3 运行结果 4 参考文献 5 Matlab代码实现 1 概述当城市数量较少时,理论上可以通过穷举法来列举出最优方案,然而当城市数量较多时,所有路线之和将呈指数 ...

追赶法求解块三对角矩阵以及matlab代码实现

追赶法求解块三对角矩阵以及matlab代码实现相关推荐

最新文章

热门文章