加权六点格式求解抛物型偏微分方程（Matlab实现）

最新推荐文章于 2023-04-21 22:11:57 发布

原创

最新推荐文章于 2023-04-21 22:11:57 发布 · 1.8k 阅读

24 ·

CC 4.0 BY-SA版权

文章标签：

#机器学习 #算法

本文介绍了一种使用加权六点格式求解抛物型偏微分方程的方法，并通过Matlab实现了该算法。具体包括定义时间及空间步长、初始化边界条件、构建迭代矩阵等步骤。

加权六点格式求解抛物型偏微分方程（Matlab实现）

clear
clc
at=[0,1];
ax=[0,1];
t1=0.01;
x1=0.01;
theta=0.5;
th=at(1):t1:at(2);
xh=ax(1):x1:ax(2);
n1=length(xh);
n2=length(th);
u=zeros(n2,n1);
syms x t;
bot=20;%时间边界
boxr=0;%空间左边界
boxl=100;%空间右边界
%下面求边界值
bt=zeros(1,n1);
for i=1:n1
   bt