matlab有限差分一维导热,一维导热方程-有限差分法-matlab实现

一维导热方程-有限差分法-matlab实现

p最新精品文档，知识共享 第五次作业(前三题写在作业纸上) 一、用有限差分方法求解 一维非定常热传导方程， 初始条件和边界条件见说明 .pdf 文件，热扩散系数 const， 22TTtx1. 用 Tylaor展开法推导出 FTCS格式的差分方程 2. 讨论该方程的相容性和稳定性，并说明稳定性要求对求解差分方程的影响。 3. 说明该方程的类型和定解条件，如何在程序中实现这些定解条件。 4. 编写 M 文件求解上述方程，并用 适当的文字 对程序做出 说明 。(部分由网络搜索得到，添加，修改后得到。) function rechuandaopde 以下所用数据，除了 t 的范围我根据题目要求取到了 20000，其余均从 pdf 中得来 a0.00001;a 的取值 xspan0 1;x 的取值范围 tspan0 20000;t 的取值范围 ngrid100 10;分割的份数，前面的是 t 轴的，后面的是 x 轴的 fx0;初值 g1t100;边界条件一 g2t100;边界条件二 T,x,tpdesolutiona,f,g1,g2,xspan,tspan,ngrid;计算所调用的函数 x,tmeshgridx,t; meshx,t,T;画图 ， 并且把坐标轴名称改为 x， t， T xlabel39;x39; ylabel39;t39; zlabel39;T39; T输出温度矩阵 dttspan2/ngrid1;t 步长 h30003000/dt; 最新精品文档，知识共享 h90009000/dt; h1500015000/dt;3000,9000,15000 下 ， 温度分别在 T 矩阵的哪些行 T3000Th3000, T9000Th9000, T15000Th15000,输出三个时间下的温度分布 不再对三个时间下的温度 -长度曲线画图，其图像就是三维图的截面 稳定性讨论 ,傅里叶级数法 dxxspan2/ngrid2;x 步长 sta4*a*dt/dx2*sinpi/22; if sta0,sta2 fprintf39;nsn39;,39;有稳定性 39; else fprintf39;nsn39;,39;没有稳定性 39; error end 真实值计算 xe,te,Tetruesolutiona,f,g1,g2,xspan,tspan,ngrid; xe,temeshgridxe,te; meshxe,te,Te;画图 ， 并且把坐标轴名称改为 xe， te， Te xlabel39;xe39; ylabel39;te39; zlabel39;Te39; Te输出温度矩阵 误差计算 jmax1/dx1;网格点数 rmswuchajisuanT,Te,jmax rms输出误差 最新精品文档，知识共享 function rmswuchajisuanT,Te,jmax for j1jmax rmsTj-Tej2/jmax1/2 end functionUe,xe,tetruesolutiona,f,g1,g2,xspan,tspan,ngrid nngrid1;t 份数 mngrid2;x 份数 Uezerosngrid; xelinspacexspan1,xspan2,m;画网格 nbsp;telinspacetspan1,tspan2,n;画网格 for j2n for i2m-1 for g1m-1 Uej,i100-400/2*g-1/pi*sin2*g-1*pi*xej*exp-a*2*g-12*pi2*tei end end end function U,x,tpdesolutiona,f,g1,g2,xspan,tspan,ngrid nngrid1;t 份数 mngrid2;x 份数 hrangexspan/m-1;x 网格长度 nbsp;xlinspacexspan1,xspan2,m;画网格 nbsp;krangetspan/n-1; t 网格长度 tlinspacetspan1,tspan2,n;画网格 Uzerosngrid; U,1g1t;边界条件 U,mg2t; 最新精品文档，知识共享 U1,fx;初值条件 nbsp;差分计算 nbsp;for j2n for i2m-1 Uj,i1-2*a*k/h2*Uj-1,ia*k/h2*Uj-1,i-1a*k/h2*Uj-1,i1; end end 5. 将温度随时间变化情况用曲线表示 00 . 20 . 40 . 60 . 8100 . 511 . 52x 1 04020406080100 xtT6. 给出 3000、 9000、 15000三个时刻的温度分布情况，对温度随时间变化规律做说明。 T3000100.0000 nbsp;63.4362 nbsp;34.2299 nbsp;15.8021 nbsp;7.4641 nbsp;7.4641 nbsp;15.8021 nbsp;34.2299 nbsp;63.4362 100.0000 T9000100.0000 nbsp;81.6930 nbsp;65.6076 nbsp;53.6839 nbsp;47.3466 nbsp;47.3466 nbsp;53.6839 nbsp;65.6076 nbsp;81.6930 100.0000 T15000100.0000 nbsp;89.9415 nbsp;81.0962 nbsp;74.5310 nbsp;71.0378 nbsp;71.0378 nbsp;74.5310 nbsp;81.0962 nbsp;89.9415 100.0000 最新精品文档，知识共享 根据数据分析，在同一个 x点上温度随时间的增加而增加，但增幅变小。 x-T图形仍为抛物线，但随着时间的增加，极值变小，图像变得平缓。 7. 用计算数据说明，并结合差分方程余项，空间、时间间隔对求解精度影响。 数据量较大，且原理相同，我取一个向量演示一下。 保持空间间隔不变，修改时间间隔，时间间隔加大，得到的误差加大。 保持时间间隔不变，修改空间间隔，空间间隔加大，得到的误差加大。 修改空间间隔的误差在增量比修改时间间隔的大。 从方差余项上来看， (没有公式编辑器。。。。只能从ppt里粘 贴了)这个余项里的 t， x都在分母上，所以与误差成正比，且 x的次数应该是比 t高，故影响较大。 8. 用计算数据说明，稳定性要求对求解精度的影响。 修改稳定性，即修改 x和 t分的份数( ngrid)，之后看误差。 稳定性越高，解的精度越高。即在满足稳定性要求( a* t/ x20.5)时， a* t/x2越接近 0，误差越小。 从概念上理解，稳定性越好，对引入时间层误差的抑制能力越强。所以误差越小。 二、调用 MATLAB 函数完成上述计算 1. 编写 M 文件求解上述方程，并用 适当的文字 对程序做出 说明 。 function pdepediaoyong m0; xlinspace0,1,11;x 的网格 tlinspace0,20000,101;t 的网格 sol pdepem,pdefun,icfun,bcfun,x,t;调用函数 Tsol,,1;解 figure;画图 surfx,t,T xlabel39;x39; ylabel39;t39; zlabel39;T39; 最新精品文档，知识共享 dt20000/100;t 步长 h30003000/dt; h90009000/dt; h1500015000/dt;3000,9000,15000 下，温度分别在 T 矩阵的哪些行 T3000Th3000, T9000Th9000, T15000Th15000,输出三个时间下的温度分布 不再对三个时间下的温度 -长度曲线画图，其图像就是三维图的截面 function c,f,spdefunx,t,T,DuDxPDE 方程函数 c100000; fDuDx; s0; function u0icfunx初始条件函数 u00; function pl,ql,pr,qrbcfunxl,Tl,xr,Tr,t边界条件函数 pl Tl-100; ql 0; pr Tr-100; qr 0; 2. 将温度随时间变化情况用曲线表示。 最新精品文档，知识共享 00 . 20 . 40 . 60 . 8100 . 511 . 52x 1 040204060801001203. 给出 3000、 9000、 15000三个时刻的温度分布情况。 T3000 100.0000 nbsp;67.1058 nbsp;39.8611 nbsp;21.1973 nbsp;10.9885 nbsp;7.8279 nbsp;10.9885 nbsp;21.1973 nbsp;39.8611 nbsp;67.1058 100.0000 T9000100.0000 nbsp;83.4839 nbsp;68.6032 nbsp;56.8191 nbsp;49.2705 nbsp;46.6732 nbsp;49.2705 nbsp;56.8191 nbsp;68.6032 nbsp;83.4839 100.0000 T15000100.0000 nbsp;90.8310 nbsp;82.5601 nbsp;75.9972 nbsp;71.7845 nbsp;70.3330 nbsp;71.7845 nbsp;75.9972 nbsp;82.5601 nbsp;90.8310 100.0000 根据数据分析，在同一个 x点上温度随时间的增加而增加，但增幅变 小。 x-T图形仍为抛物线，但随着时间的增加，极值变小，图像变得平缓。 4. 用计算数据说明，空间、时间间隔对求解精度影响，并与有限差分法的计算结果做比较。 调用前面做出来的真实值 ,跟 pdepe做出来的值计算误差，再与有限差分法的误差比较。用 pdepe 函数求的误差远小于有限差分法，所以 pdepe 函数法更精确。 5. 用计算数据说明，有无稳定性要求，为什么若有，如何对求解精度的影响。 不知道这个 pdepe 函数的稳定性要用什么检验。傅里叶级数检验不适用。/p