clear
clc
close?all?

dx=0.05;????%x方向的步長(zhǎng)
dt=0.001;???%t方向的步長(zhǎng)
r=dt/（dx^2）;??%計(jì)算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;???????%設(shè)置初值條件：Phi（x0）=100
Phi（2:N+11）=0;?????%設(shè)置邊界條件：Phi（0t）=0
Phi（2:N+1M+1）=0;???%設(shè)置邊界條件：Phi（1t）=0

%根據(jù)差分方程，計(jì)算Phi的數(shù)值解
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（‘?dāng)U散方程的數(shù)值模擬‘）
msgbox（‘MATLAB編程答疑，請(qǐng)加QQ:?1530497909‘‘MATLAB答疑‘‘help‘）
web?http://url.cn/NSFcAs?-browser

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----

?????文件????????740??2014-06-21?11:45??MATLAB求解偏微分方程（擴(kuò)散方程）有限差分法?源程序代碼\MATLAB_code.m

?????文件????????333??2014-05-17?20:48??MATLAB求解偏微分方程（擴(kuò)散方程）有限差分法?源程序代碼\Read?Me.url

?????文件????????294??2014-02-06?11:28??MATLAB求解偏微分方程（擴(kuò)散方程）有限差分法?源程序代碼\聯(lián)系我們.txt

?????目錄??????????0??2014-06-21?11:45??MATLAB求解偏微分方程（擴(kuò)散方程）有限差分法?源程序代碼

-----------?---------??----------?-----??----

?????????????????1367????????????????????4