MATLAB求解偏微分方程扩散方程有限差分法

MATLAB求解偏微分方程（扩散方程）有限差分法，处理偏微分方程。

clear
clc
close all 

dx=0.05;    %x方向的步长
dt=0.001;   %t方向的步长
r=dt/（dx^2）;  %计算r的值
x=0:dx:1;     %得到x的序列
t=0:dt:0.2;     %得到t的序列

M=length（x）-1;
N=length（t）-1;

Phi=ones（N+1M+1）;
Phi（1:）=100;       %设置初值条件：Phi（x0）=100
Phi（2:N+11）=0;     %设置边界条件：Phi（0t）=0
Phi（2:N+1M+1）=0;   %设置边界条件：Phi（1t）=0

%根据差分方程，计算Phi的数值解
for k=1:N
    for i=2:M
        Phi（k+1i）=（1-2*r）*Phi（ki）+r*（Phi（ki-1）+Phi（ki+1））;
    end
end

[xt]=meshgrid（xt）;
mesh（xtPhi）     %绘制（xtPhi）的三维图
xlabel（‘x‘）
ylabel（‘t‘）
zlabel（‘\Phi（xt）‘）
title（‘扩散方程的数值模拟‘）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         639  2015-09-27 10:13  MATLAB求解偏微分方程（扩散方程）有限差分法 源程序代码\MATLAB_code.m
     文件          78  2015-09-27 10:07  MATLAB求解偏微分方程（扩散方程）有限差分法 源程序代码\欢迎交流.txt
     目录           0  2015-09-27 10:14  MATLAB求解偏微分方程（扩散方程）有限差分法 源程序代码\