| 查看: 1021 | 回复: 3 | |||
| 【悬赏金币】回答本帖问题，作者lzwangs将赠送您 15 个金币 | |||
[求助]
Matlab求解带参数的偏微分方程时出现 out of memory已有1人参与
|
|||
|
求带参数反应扩散方程 \frac{\partial u}{\partial t}=D\frac{\partial^2 u}{\partial t^2}+\mu u(1-u) , 边界条件为诺依曼边界条件，初值条件见程序附件部分。\mu>0 是参数，用的是后向差分法（本程序解不带参数的方程经过测试完全没有问题）。 程序中将 \mu 设为 syms 格式，最后的数值解中也应该有参数 \mu. 问题1： 时间方面迭代不到10次，会出现 out of memory。 原因应该就是引入参数的符号运算和数值运算结合在一起造成的，求问题原因和解决方案? 问题2： 用了 vpa转换，迭代6次（和目标次数相差太远）的时候不会 out of memory，但求出来的是关于 \mu 的多项式，且没有合并同类项，太长导致计算机无法显示，求解决方案！ 附件中是m代码，请高手留意￥￥￥￥￥￥￥￥￥部分！ 诚求高手指点交流！ |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：libolin3@tal.com - 附件 1 : test1.m
2018-05-16 23:40:41, 2.7 K
» 猜你喜欢
博士论文被抄袭
已经有43人回复
-
"颜宁:基础研究应顶天立地"能做到基础研究同时顶天立地的才是牛人
已经有6人回复
-
大学直属学院卸任的副院长退休后还享受副处级的养老待遇吗？
已经有4人回复
-
投稿RSC旗下杂志突然间看不到投稿状态，里面投稿的文章信息不见了，有遇到过的吗？
已经有3人回复
-
求论文投稿期刊推荐
已经有4人回复
-
AC审稿
已经有6人回复
-
聊天
已经有13人回复
-
南京大学-广州大学联合招聘博士后 欢迎广大优秀人才！！！
已经有13人回复
-
博士白读了
已经有52人回复
-
with editor 两个月了，什么原因？
已经有16人回复
somomo91
专家顾问 (职业作家)
-
专家经验: +3 - 应助: 144 (高中生)
- 金币: 19352.2
- 红花: 115
- 帖子: 3854
- 在线: 1053.2小时
- 虫号: 2399301
- 注册: 2013-04-03
- 专业: 动力学与控制
- 管辖: 计算模拟
2楼2018-05-23 01:49:24
送红花一朵|
我让显示 u（N+1,:）（即关于时间变量最后一节点处的所有值） 的时候，其中的元素（除了边界上的 u（N+1,1）和u（N+1,J+1） 为 0，这是诺依曼边界条件决定的，不需计算），其余项是关于 \mu 的多项式，很多项是二次的（屏幕上能显示出来的部分全是二次），并没有合并同类项，或者舍去系数特别小的项，结果太长（好像超过两万多字符了），光 u（N+1,2）一个值太长而导致屏幕无法全部输出，这应该是是导致 out of memory 的原因（因为涉及多项式的存储），请问有没有好的方法解决这一问题，我知道有人编过这样的程序，人家解的是好几个方程组成的方程组，但是无奈搞不到人家的代码，自己也一直想不到好的解决方案。不知您能否想想有什么好的办法没有？不胜感激！  |
3楼2018-05-23 09:11:51
qqinly
专家顾问 (著名写手)
至人无己，神人无功，圣人无名
-
专家经验: +68 - 应助: 119 (高中生)
- 金币: 15820.8
- 红花: 29
- 帖子: 2612
- 在线: 588.8小时
- 虫号: 2666806
- 注册: 2013-09-20
- 性别: GG
- 专业: 信号理论与信号处理
- 管辖: 程序语言
4楼2018-06-16 09:41:21
5