资源描述：
《用向前差分格式计算初边值问题.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、用向前差分格式计算初边值问题信计1班周晓虹一、题目用向前差分格式计算如下热传导方程的初边值问题已知其精确解为二、考虑的问题作为模型，考虑一维热传导方程：…………（1.1）其中是正常数，是给定的连续函数。现在考虑第二类初边值问题的差分逼近：初始条件：…………（1.2）边值条件：，，………（1.3）假设和在相应区域光滑，并且在满足相容条件，使上述问题有惟一充分光滑的解。三、网格剖分取空间步长和时间步长，其中都是正整数。用两族平行直线和将矩形域分割成矩形网格，网格节点为。以表示网格内点集合，即位于开矩形的网点集合；表示所有位于闭矩形的网点集合；是网格界点集合。其次，
2、用表示定义在网点的函数，用向前差分格式计算初边值问题信计1班周晓虹一、建立差分格式将方程在节点离散化，，…………（1.4）对充分光滑的解，由Taylor展式：…………（1.5）…………（1.6）…………（1.7）（1.5）移项得：…………（1.8）（1.6）（1.7）相加得：…………（1.9）将（1.8）（1.9）代入（1.4）得：…………（1.10）其中，舍去，得到逼近（1.1）的向前格式差分方程：，……（1.11）其中，，记用向前差分格式计算初边值问题信计1班周晓虹则由（1.4）由（1.11）一、截断误差（3）.边界条件在本题中，，，，，二、稳定性分析用傅
3、里叶方法对差分格式进行稳定性分析以表示网比，将（1.11）改写成便于计算的形式：（本题中）以代入，得消去，则知增长因子由，得用向前差分格式计算初边值问题信计1班周晓虹即只需解得所以向前差分格式的稳定性条件是一、结论抛物型方程的有限差分法的步骤大致可以归纳如下：1.对区域进行网格剖分2.在离散结点建立相应的差分格式3.处理初边值条件4.进行稳定性分析由本题可以总结出，抛物型方程的有限差分法所得的数值解能够较好地逼近方程的精确解，且区域剖分得越细，即步长越小，数值解与精确解的误差就越小，数值解越逼近精确解。附录MATLAB程序：取，，则，满足稳定性条件另取，，则，
4、亦满足稳定性条件另取，，则，亦满足稳定性条件formatlonga=2;l=1;T=1;N=10;M=400;h=l/N;to=T/M;r=(a*to)/h^2;forj=1:N+1用向前差分格式计算初边值问题信计1班周晓虹x(j)=(j-1)*h;fork=1:M+1t(k)=(k-1)*to;u(j,k)=exp(x(j)+2*t(k));endendu%求解精确解forj=1:N+1x(j)=(j-1)*h;us(j,1)=exp(x(j));endfork=1:M+1t(k)=(k-1)*to;us(1,k)=exp(2*t(k));us(N+1,k)
5、=exp(1+2*t(k));endfork=2:M+1forj=2:Nus(j,k)=r*us(j-1,k-1)+(1-2*r)*u(j,k-1)+r*us(j+1,k-1);endendus%求解数值解fork=1:M+1forj=1:N+1R(j,k)=abs(u(j,k)-us(j,k));endendR%计算误差Rmax=max(max(R))%求误差的最大值精确解与数值解的比较：x=0:0.1:1;holdonplot(x,u(:,M+1),'b');plot(x,us(:,M+1),'y');title('t=1，h=1/10，τ=1/400时精
6、确解和数值解的比较')text(0.05,21,'蓝：精确解');text(0.05,20,'黄：数值解');用向前差分格式计算初边值问题信计1班周晓虹holdoff取不同步长时的误差比较：x=0:1/10:1;y=0:1/20:1;z=0:1/40:1;holdonplot(x,R(:,M+1),'b');holdoffM分别取10，20，40结论抛物型方程的有限差分法的步骤大致可以归纳如下：1.对区域进行网格剖分2.在离散结点建立相应的差分格式3.处理初边值条件4.进行稳定性分析由本题可以总结出，抛物型方程的有限差分法所得的数值解能够较好地逼近方程的精确解
7、，且区域剖分得越细，即步长越小，数值解与精确解的误差就越小，数值解越逼近精确解。班级：信计1班名字：周晓虹学号：
