数理方程matlab实验_第1页
数理方程matlab实验_第2页
数理方程matlab实验_第3页
全文预览已结束

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

有限长细杆在第一种边界条件下的热传导定解与可视化刘强(2010001040005)摘要:本文主要研究了有限长细杆在第一种边界条件下的热传导问题以及用matlab编程将其可视化,通过图形变化进一步分析这种条件下热传导的特点。关键词: 热传导 边界条件 matlab一、 问题描述设有一根长度为L均匀的,内部无热源的热传导细杆,侧面绝热,其左右端均保持零度,初始温度分布已知。该定解问题应为:设,于是可以得到于是得到固有值问题很容易得到由此易解得由叠加原理得并且由初始条件可得假定我们取定,取且所得的解为:二、 程序设计用MATLAB求解,所用程序如下5:function = main() N=50;t=0:0.0001:1;x=0:0.21:20; w=Sum(N,t(1); h=plot(x,w,linewidth,5); axis(0,20,0,1.5); for n=2:length(t) %对不同的时间t绘制图形 w=Sum(N,t(n); set(h,ydata,w); drawnow; pause(0.1) end function u=Sum(N,t) %该函数用于计算时间固定时u的波形图 x=0:0.21:20;u=0; for k=1:2*N %累加 u=u+exp(-(k2*pi2*102/400*t)*2/k/pi*(cos(k*pi*10/20). -cos(k*pi*11/20)*sin(k*pi*x./20);End有限长细杆热传导解析解的等时距截图当t继续增大时的分布图三、 问题分析从matlab程序运行来看,热量一开始分布在(10,11)这个区间,这个与题设给定的初态是吻合的;而由后面的等时距截图(当然也可直接从matlab运行上来看)可一观察得到,热量不断地向着两边传播,因为中间部分在不断减少,而两边的在不断的增高;但是到了一定程度以后,可以看出,整个细杆上的热量都开始衰减,这个当然是由e指数来控制的。最后可以看出,在t趋于无穷大时,整条杆上的热量都趋于0,并且我们还可以发现,整个变化过程是先快后慢,究其根本原因就在于exp这一项的负指数变化规律。通过这次对数学物理方程的课程内容进行仿真可视化,让我对课本内容有了更加深刻的认识,但是对于matlab的掌握还是不够熟练,其中也参考了数学物理方程的matlab的解法与可视化的一些指导。参考文献 1 数学物理方程 李明奇,田太心.电子科技大学出版社,2010;2 数学物理方法 梁昆淼 北京: 高等教育版社,1998 P223-226; 3 MATLAB实用教程 徐金明(MA

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论