d******t 发帖数: 1114 | 1 【 以下文字转载自 Mathematics 讨论区 】
发信人: davidibt (david), 信区: Mathematics
标 题: 微分方程用MATLAB求解
发信站: BBS 未名空间站 (Wed Jul 11 18:42:00 2012, 美东)
如附件所示,用MATLAB求解怎么弄都求不出来。所以还望高手指点下,双黄包感谢 |
d*a 发帖数: 1863 | 2 考
你也太图省事了
这么复杂的传热方程要是能有解析解,这世界就太美好了 |
C********g 发帖数: 274 | 3 k是温度的函数,方程是非线性的,恐怕没有解析解。
如果k是常数,大概有解析解:根据叠加原理构造辅助解,然后分离变量;或者用格林
函数
【在 d******t 的大作中提到】 : 【 以下文字转载自 Mathematics 讨论区 】 : 发信人: davidibt (david), 信区: Mathematics : 标 题: 微分方程用MATLAB求解 : 发信站: BBS 未名空间站 (Wed Jul 11 18:42:00 2012, 美东) : 如附件所示,用MATLAB求解怎么弄都求不出来。所以还望高手指点下,双黄包感谢
|
d******t 发帖数: 1114 | 4 多谢楼上两位的指点。但是怎么构造辅助解啊?自己做的时候就知道没有解析解,所以
只有找数值解了,但是就是困在这步不知道怎么进行下一步 |
f*****e 发帖数: 2992 | 5 远离PDE。
【在 d******t 的大作中提到】 : 多谢楼上两位的指点。但是怎么构造辅助解啊?自己做的时候就知道没有解析解,所以 : 只有找数值解了,但是就是困在这步不知道怎么进行下一步
|
C********g 发帖数: 274 | 6 其次方程+非其次边界,非其次方程+其次边界。
不过原方程是非线性的,这些就用不上了。最快最省事的方式肯定是用软件做数值解,
随便一个CFD类的软件应该都行。软件用的熟的话,10分钟左右大概能完成。
如果用matlab,就用里面的pde toolbox,应该能解决
【在 d******t 的大作中提到】 : 多谢楼上两位的指点。但是怎么构造辅助解啊?自己做的时候就知道没有解析解,所以 : 只有找数值解了,但是就是困在这步不知道怎么进行下一步
|
k*******2 发帖数: 1429 | 7 Kirchhoff transform ?
【在 d******t 的大作中提到】 : 多谢楼上两位的指点。但是怎么构造辅助解啊?自己做的时候就知道没有解析解,所以 : 只有找数值解了,但是就是困在这步不知道怎么进行下一步
|
C********g 发帖数: 274 | 8 这个确实值得试一下,这个积分变换还是挺简单的,新的变量如果能解出来,温度就能
显示的得到。
不过LZ的边界条件有问题哦,r方向上有三个边界条件,z方向上一个都没有。。。
另外r方向上还有一个负的边界r=-0.002,这个坐标系很奇怪,我怀疑是弄错了,
除非已经做了坐标变换。
LZ去检查一下边界条件吧
【在 k*******2 的大作中提到】 : Kirchhoff transform ?
|