资源描述
细心整理一、 题目用向前差分格式计算如下热传导方程的初边值问题确定其精确解为二、 考虑的问题作为模型，考虑一维热传导方程：1.1其中是正常数，是给定的连续函数。此时此刻考虑其次类初边值问题的差分靠近： 初始条件：1.2边值条件：，1.3假设和在相应区域光滑，并且在满足相容条件，使上述问题有惟一充分光滑的解。三、 网格剖分取空间步长和时间步长，其中都是正整数。用两族平行直线和将矩形域分割成矩形网格，网格节点为。以表示网格内点集合，即位于开矩形的网点集合；表示全部位于闭矩形的网点集合；是网格界点集合。其次，用表示定义在网点的函数，四、 建立差分格式将方程在节点离散化， 1.4对充分光滑的解，由Taylor展式：1.51.61.71.5移项得：1.81.61.7相加得：1.9将1.81.9代入1.4得：1.10其中，舍去，得到靠近1.1的向前格式差分方程：， 1.11其中，记那么由1.4 由1.11 五、 截断误差3.边界条件在此题中，六、 稳定性分析用傅里叶方法对差分格式进展稳定性分析以表示网比，将1.11改写成便于计算的形式： 此题中以代入，得消去，那么知增长因子 由，得即 只需 解得 所以向前差分格式的稳定性条件是七、 结论 抛物型方程的有限差分法的步骤大致可以归纳如下： 1.对区域进展网格剖分 2.在离散结点建立相应的差分格式 3.处理初边值条件 4.进展稳定性分析由此题可以总结出，抛物型方程的有限差分法所得的数值解能够较好地靠近方程的精确 解，且区域剖分得越细，即步长越小，数值解与精确解的误差就越小，数值解越靠近精确解。附录MATLAB程序：取，那么，满足稳定性条件另取，那么，亦满足稳定性条件另取，那么，亦满足稳定性条件format longa=2;l=1;T=1;N=10;M=400;h=l/N;to=T/M;r=(a*to)/h2;for j=1:N+1 x(j)=(j-1)*h; for k=1:M+1 t(k)=(k-1)*to; u(j,k)=exp(x(j)+2*t(k); endendu %求解精确解for j=1:N+1 x(j)=(j-1)*h; us(j,1)=exp(x(j);endfor k=1:M+1 t(k)=(k-1)*to; us(1,k)=exp(2*t(k); us(N+1,k)=exp(1+2*t(k);endfor k=2:M+1 for j=2:N us(j,k)=r*us(j-1,k-1)+(1-2*r)*u(j,k-1)+r*us(j+1,k-1); endendus %求解数值解for k=1:M+1 for j=1:N+1 R(j,k)=abs(u(j,k)-us(j,k); endendR %计算误差Rmax=max(max(R) %求误差的最大值精确解与数值解的比拟：x=0:0.1:1;hold onplot(x,u(:,M+1),b);plot(x,us(:,M+1),y);title(t=1，h=1/10，=1/400时精确解和数值解的比拟)text(0.05,21,蓝：精确解);text(0.05,20,黄：数值解);hold off取不同步长时的误差比拟：x=0:1/10:1;y=0:1/20:1;z=0:1/40:1;hold onplot(x,R(:,M+1),b);hold offM分别取10，20，40结论 抛物型方程的有限差分法的步骤大致可以归纳如下： 1.对区域进展网格剖分 2.在离散结点建立相应的差分格式 3.处理初边值条件 4.进展稳定性分析由此题可以总结出，抛物型方程的有限差分法所得的数值解能够较好地靠近方程的精确 解，且区域剖分得越细，即步长越小，数值解与精确解的误差就越小，数值解越靠近精确解。班级：信计1 班名字：周 晓虹学号：20083710
展开阅读全文