| 查看: 1311 | 回复: 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
» 猜你喜欢
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 (高中生)
- 金币: 15507.8
- 散金: 186
- 红花: 29
- 帖子: 2632
- 在线: 589.7小时
- 虫号: 2666806
- 注册: 2013-09-20
- 性别: GG
- 专业: 信号理论与信号处理
- 管辖: 程序语言
4楼2018-06-16 09:41:21
