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求解偏微分方程（擴散方程）有限差分法?源程序代碼\