Matlab实现有限差分法解二维热传导问题

利用Matlab解二维热传导问题，主要采用了有限差分法，同时使用追赶法解对角矩阵，包括相应函数、例程及图像

function x=chase（Af）
%x = chase（Af），A为对角矩阵，f为右侧常量x为方程的解
%该函数正常工作前提为A的行列式不为零
n=rank（A）;
[LU]=LU（An）;
sum=0;
y=zeros（n1）;
y（1）=f（1）;
for k=2:n
    for j=1:k-1
        sum=sum+L（kj）*y（j）;
    end
    y（k）=f（k）-sum;
    sum=0;
end
x=zeros（n1）;
x（n）=y（n）/U（nn）;
k=n-1;
while（k>=1）
    for j=k+1:n
        sum=sum+U（kj）*x（j）;
    end
    x（k）=（y（k）-sum）/U（kk）;
    sum = 0;k=k-1;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        454  2019-07-30 09:50  Code\chase.m

     文件      29019  2019-07-30 09:41  Code\example.jpg

     文件       1771  2019-07-30 09:49  Code\example.m

     文件        332  2019-07-30 09:40  Code\LU.m

     目录          0  2019-07-30 09:49  Code

----------- ---------  ---------- -----  ----

                31576                    5