| 查看: 1243 | 回复: 3 | |||
[求助]
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代码，请高手留意￥￥￥￥￥￥￥￥￥部分！ 诚求高手指点交流！ |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : test1.m
2018-05-16 23:40:41, 2.7 K
» 猜你喜欢
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有4人回复
-
磺酰氟产物，毕不了业了！
已经有6人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有9人回复
-
Cas 72-43-5需要30g，定制合成，能接单的留言
已经有8人回复
-
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？
已经有8人回复
-
论文终于录用啦！满足毕业条件了
已经有25人回复
-
自荐读博
已经有3人回复
-
不自信的我
已经有5人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
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 (高中生)
- 金币: 15500.8
- 散金: 186
- 红花: 29
- 帖子: 2631
- 在线: 589.7小时
- 虫号: 2666806
- 注册: 2013-09-20
- 性别: GG
- 专业: 信号理论与信号处理
- 管辖: 程序语言
4楼2018-06-16 09:41:21