昆明市住房和城乡建设局网站上看的,极速网站建设定制价格,英文网站建设的请示怎么写,wordpress批量改文章链接之前写过一篇博客#xff0c;介绍如何使用差分格式求解热传导方程
今天打开博客#xff0c;突然发现评论区被这篇文章霸屏了 询问实验室的小伙伴才知#xff0c;原来是被可爱的建模学子们攻占了
经过简单的了解#xff0c;发现今年建模的A题的核心就是求解一个热传导方程…之前写过一篇博客介绍如何使用差分格式求解热传导方程
今天打开博客突然发现评论区被这篇文章霸屏了 询问实验室的小伙伴才知原来是被可爱的建模学子们攻占了
经过简单的了解发现今年建模的A题的核心就是求解一个热传导方程因此之前所写文章的程序基本可以算是神助攻了完全可以帮助大家构建解题程序的基本框架。
但是
数学建模比赛考验的就是大家的学习能力以及解决问题的能力我只提供简要思路不要做伸手党直接找我要代码不要问我怎么写请对得起你将来获得的一等奖。
再有就是你要先确保你能看懂我已经提供的源代码否则下面我提供的思路你可能看了也白看。 这里统一对大家的问题做一个回答
本题适用差分解法吗
求解偏微分方程的方法中差分方法 和 有限元 是两类最主流的方法。
差分方法的优点是原理简单但是只能求解规则区域的数值解。
有限元背后的理论相对难很多但能够求解不规则区域问题。
本题适用哪种方法解答取决于你的模型假设
本题中涉及多种介质的热传导的求解我的建议是如果不想给自己找麻烦的话将每种介质层假设成规则的矩形。
既然带求解区域是矩形了那么本题使用差分方法来求解更加合适关于差分方法你可以随便找一本介绍偏微分方程数值解的书都有介绍。
当然你如果将模型假设定义为更符合实际的不规则问题能做出来当然是亮点但切记不要搬起石头砸自己的脚毕竟建模时间紧任务重。
是否适用于多层壁热传导
评论区有人问是否适用于多层壁热传导吗
其实就是问这个程序能否求解A题嘛O(∩_∩)O哈哈~
答案当然是能但显然不能直接拿来用给几点提示。
思路1
你单拿出其中一层来求解和我提供的算例已经没有本质区别了。
因此你可以一层一层的求解。先求第一层的数值解第一层的结果一有第二层的边界条件也就有了于是第二层也可以求了。
这样做的潜在问题是第一层中求解的误差必定会传递到第二层去数学上可能不太完美但是好理解代码改动也少。
思路2
我认为数学上更好的方式肯定是整体一起求解但这就有点困难了。
这样做时你需要对每一层边界在系数矩阵的相应位置处都按照边界处的对应关系进行相应处理。
这需要你对差分方法有着很好的理解如果我提供的代码你无法完全看懂建议就不要考虑了。
关于边界条件
构造的差分格式是保证解满足对应的方程但其实满足给定方程的解有无穷多种。
而边界条件的作用其实就是找出你想要的那个解。
之前文章中给出的算例包含的边界条件是
u(x,0)
u(0,t) 和 u(1,t)
在A题中右侧初始温度好像是没有的也就是u(1,t)没有
首先你要知道的是求解需要的边界条件并不一定非得是这几个
但是少了一个边界条件你就要想办法补上一个边界条件, 边界条件也不一定是已知函数的表达式导数的表达式也是可以的当然代码是一定需要相应修改的。
比如没有u(1,t)你可以想办法构造 du(0,t)/dx 或 du(0,t)/dt
这就看你如何理解原问题了建议查阅文献看看别人使用的是哪种边界条件相应的对代码进行修改。当然也可以通过模型假设将问题向你期待的边界条件上面靠。 思路提示目前就是这些了
最近在经营我的github觉得本文对你有帮助的话进来逛逛吧顺便给觉得好玩的项目点颗星~嘿嘿~
同时也欢迎大家关注我的团队(微信搜索公众号LSGO软件技术团队)我快毕业了希望团队能有一些新的血液~
最后祝大家都能获得好成绩~
