一维热传导方程数值解法及matlab实现分离变量法和有限差分法

一维热传导方程的Matlab解法分离变量法和有限差分法问题描述实验原理分离变量法实验原理有限差分法实验目的利用分离变量法和有限差分法解热传导方程问题利用matlab进行建模构建图形研究不同的情况下采用何种方法从更深层次上理解热量分布与时间、空间分布关系。模拟与仿真作业(1)分离变量法(代码)x00.1*pipi;y00.041;x,tmeshgridx,y;s0;mlengthj;matlab可计算的最大数相当于无穷fori1mss200*1--1i/i*pi*sini*x.*exp-i2*t;end;surfx,t,s;xlabelx,ylabelt,zlabelT;title分离变量法(无穷);axis0pi010100;所得到的三维热传导图形为有限差分法uzeros10,25;t1xpi构造一个1025列的矩阵(初始化为0)用于存放时间t和变量x横坐标为x纵坐标为ts1/25/pi/102;fprintf稳定性系数S为n;disps;fori29ui,1100;end;forj125u1,j0;u10,j0;end;forj124fori29ui,j1s*ui1,j1-2*s*ui,js*ui-1,j;endenddispu;x,tmeshgrid125,110;surfx,t,u;xlabelt,ylabelx,zlabelT;title有限差分法解;所得到的热传导图形为(2)i分离变量法(取前100项和)x00.1*pipi;y00.041;x,tmeshgridx,y;s0;fori1100ss200*1--1i/i*pi*sini*x.*exp-i2*t;end;surfx,t,u;xlabelx,ylabelt,zlabelT;title分离变量法;axis0pi010100;所得到的热传导图形为Ii有限差分法根据(1)我们有如下图结论比较可得这两幅图基本相同，有限差分法和分离变量法对本题都适应(3)第一种情况(取无穷项)在原来程序代码的基础上加上disps,6;可得出第六列(即xpi/2)处温度随时间的变化情况x00.1*pipi;y00.041;x,tmeshgridx,y;s0;mlengthj;matlab可计算的最大数，相当于无穷fori1mss200*1--1i/i*pi*sini*x.*exp-i2*t;end;surfx,t,s;xlabelx,ylabelt,zlabelT;title分离变量法(无穷);axis0pi010100;disps,6;我们得到如下一组数据当温度低于50度是时间为t23.5*0.040.94第二种情况(取前100项和)在原来程序代码的基础上加上disps,6;可得出第六列(即xpi/2)处温度随时间的变化情况x00.1*pipi;y00.041;x,tmeshgridx,y;r0.04/0.1*pi2;fprintf稳定性系数S为dispr;s0;fori1100ss200*1--1i/i*pi*sini*x.*exp-i2*t;end;surfx,t,s;xlabelx,ylabelt,zlabelT;title分离变量法;axis0pi010100;disps,6;当温度低于50度是时间为t23.5*0.040.94第三种情况(有限差分法)在原来程序代码的基础上加上dispu5,;可得出第五行(即xpi/2)处温度随时间的变化情况uzeros10,25;t1xpi10行25列横坐标为x纵坐标为ts1/25/pi/102;fprintf稳定性系数S为n;disps;fori29ui,1100;end;forj125u1,j0;u10,j0;end;forj124fori29ui,j1s*ui1,j1-2*s*ui,js*ui-1,j;endenddispu;x,tmeshgrid125,110;surfx,t,u;xlabelt,ylabelx,zlabelT;title有限差分法解;dispu5,;得到如下结果我们知19列为50.350520列是数据为47.8902所以时间t为20*0.040.78结论比较一二三种情况，我们得到不同的时间，这是由于1、加和不同一种为100，一种为无穷；2、采用的方法不同一种为分离变量法，一种为有限差分法造成的。第一题完解根据课本知识我们知道稳定性系数S是衡量由差分格式所得的数据是否稳定的因数，如果S1/2，我们得到的数据就会不稳定，所画出的图形失真也就会很大。一下就S1/2，两种情况进行讨论(以第一题为例进行讨论)当S1/2,时取步长k1/25,hpi/20,S1.6211uzeros20,25;t1xpi20行25列横坐标为x纵坐标为ts1/25/pi/202;fprintf稳定性系数S为n;disps;fori219ui,1100;end;forj125u1,j0;u10,j0;end;forj124fori219ui,j1s*ui1,j1-2*s*ui,js*ui-1,j;endenddispu;x,tmeshgrid125,120;surfx,t,u;xlabelt,ylabelx,zlabelT;title有限差分法解s1/2时;所得数据为所得到的图形为显而易见，我们得到的图形与实际的热传导图形相差太大，这是由于S1/2，得到的数据是不稳定的，所以画出的图形也是不准确的。第三题完成根据课上所学知识，我们有如下方程为便于解释做题，我们令a1lpix;下面开始求解分离变量法根据课上所讲其中2000,0,00,0,0,0txxxxxltuauxltuutuxxl我们有如下代码x00.1*pipi;y00.410;x,tmeshgridx,y;u0;mlengthj;matlab可计算的最大数，相当于无穷fori0muu8*-1i/pi*2*i12*sin2*i1/2*x.*exp-2*i12/4*t;end;surfx,t,u;xlabelx,ylabelt,zlabelT;title分离变量法(无穷);dispu;得到如图所示的热传导现象有限差分法uzeros20,100;t1xpi20行100列横坐标为x纵坐标为ts1/100/pi/202;fprintf稳定性系数S为n;disps;fori120ui,1i/20*pi;;end;forj1100u1,j0;endforj199fori219ui,j1s*ui1,j1-2*s*ui,js*ui-1,j;endendforj1100u20,ju19,j;end;dispu;x,tmeshgrid1100,120;surfx,t,u;xlabelt,ylabelx,zlabelT;title有限差分法解;我们得到如图所示的热传导方程结论比较可得由以上两种方法作出的三维图形基本相同，符合热传导的热量分布随时间和空间的变化规律第四题完成