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（‘扩散方程的数值模拟‘）
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