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（‘擴散方程的數值模擬‘）
msgbox（‘MATLAB編程答疑，請加QQ:?1530497909‘‘MATLAB答疑‘‘help‘）
web?http://url.cn/NSFcAs?-browser

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2014-06-21?11:45??MATLAB求解偏微分方程（擴散方程）有限差分法?源程序代碼\
?????文件?????????740??2014-06-21?11:45??MATLAB求解偏微分方程（擴散方程）有限差分法?源程序代碼\MATLAB_code.m
?????文件?????????333??2014-05-17?20:48??MATLAB求解偏微分方程（擴散方程）有限差分法?源程序代碼\Read?Me.url
?????文件?????????294??2014-02-06?11:28??MATLAB求解偏微分方程（擴散方程）有限差分法?源程序代碼\聯系我們.txt