抛物线差分格式求解matlab程序

抛物线差分格式求解，包括一维古典显格式，DFF格式，CN格式，局部一维方法和预测校正格式的求解；文件中有具体的题目及解法说明并还有matlab程序供大家参考；易读易懂。

%运用Crank-Nicolson差分格式求解二维扩散方程的初值问题；
clear all
ti=0.00001;%时间步长
h=0.01;%空间步长
t=0.004;
k=t/ti;
m=1/h+1;
g=ti/h^2; %网格比;
%初值条件及边界条件
w=zeros（mm）;
u=zeros（mm）;
for i=（0.4/h+1）:（0.6/h+1）
    for j=（0.3/h+2）:（0.7/h）
        u（ij）=1;
    end
end
for i=（0.3/h+2）:（0.7/h）
    for j=（0.4/h+1）:（0.6/h+1）
        u（ij）=1;
    end
end
tic
for l=1:k
u=Gauss_Seidel（ug）;%Gauss_Seidel迭代
end
toc
t=toc
mesh（0:h:10:h:1u）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        479  2010-11-18 22:11  程序\CN.m

     文件        613  2010-11-18 17:51  程序\DFF.m

     文件        455  2010-11-18 16:40  程序\draw\CN.m

     文件        619  2010-11-14 19:24  程序\draw\DFF.m

     文件        223  2010-11-18 16:40  程序\draw\draw.m

     文件        425  2010-11-14 16:31  程序\draw\Gauss_Seidel.m

     文件        614  2010-11-18 16:44  程序\draw\gdxs.m

     文件       1192  2010-11-18 10:04  程序\draw\ycjz.m

     文件        929  2010-11-14 19:53  程序\draw\ywjb.m

     文件        310  2008-04-29 17:21  程序\draw\zuigan.m

     文件        425  2010-11-14 16:31  程序\Gauss_Seidel.m

     文件        608  2010-11-18 17:52  程序\gdxs.m

     文件       1261  2010-11-18 17:53  程序\ycjz.m

     文件        970  2010-11-18 17:53  程序\ywjb.m

     文件        310  2008-04-29 17:21  程序\zuigan.m

     文件     864256  2011-01-26 10:15  第二次作业.doc

     目录          0  2010-11-18 22:13  程序\draw

     目录          0  2010-11-18 22:14  程序

----------- ---------  ---------- -----  ----

               873689                    18