Python中文
首页
教程
问答
标签
搜索
登录
注册
带两个耦合节点的scipy.integrate.ode?
回答此问题可获得
20
贡献值,回答如果被采纳可获得
50
分。
<p>我目前正在尝试使用SciPy的integrate.ode包来解决一对耦合的一阶ode:比如说<a href="http://en.wikipedia.org/wiki/Lotka%E2%80%93Volterra_equation" rel="noreferrer">Lotka-Volterra predator-prey equation</a>。然而,这意味着在集成循环期间,我必须在每次迭代中更新发送给方法的参数,并且仅仅跟踪上一个值并在每次迭代中调用<code>set_f_params()</code>似乎并没有起到作用。</p> <pre><code>hprev = Ho pprev = Po yh = np.zeros(0) yp = np.zeros(0) while dh.successful() and dp.successful() and dp.t < endtime and dh.t < endtime: hparams = [alpha, beta, pprev] pparams = [delta, gamma, hprev] dh.set_f_params(hparams) dp.set_f_params(pparams) dh.integrate(dh.t + stepsize) dp.integrate(dp.t + stepsize) yh = np.<a href="https://www.cnpython.com/list/append" class="inner-link">append</a>(yh, dh.y) yp = np.append(yp, dp.y) hprev = dh.y pprev = dp.y </code></pre> <p>我在通过<code>set_f_params</code>的每次迭代中设置的值似乎不会传播到回调方法,这并不奇怪,因为web上的任何示例似乎都不涉及传递给回调的“live”变量,但这是唯一一个可以将这些值获取到回调方法中的方法。</p> <p>有人对如何使用SciPy对这些ODEs进行数值集成有什么建议吗?</p>
0 条评论
分类:
Python问答
请先
登录
后评论
默认排序
时间排序
1 个回答
匿名
1天前
擅长:python、mysql、java
<p>我也有类似的问题。结果是,积分器并不是在每次调用integrate()时重新计算微分方程函数,而是在它自己的内部时间进行计算。我将积分器的max_step选项更改为与stepsize相同,这对我有效。</p>
请先
登录
后评论
针对此问题:
更多的回答
关注
89
关注
收藏
1
收藏,
216
浏览
网友 提问于 2天前
相关Python问题
尽管Python中的所有内容都是引用,为什么Python导师在没有指针的列表中绘制字符串和整数?
9 回答
尽管python中的表达式为false,但循环仍在运行
6 回答
尽管python代码正确,但从nifi ExecuteScript处理器获取语法错误
6 回答
尽管Python在Neovim中工作得很好,但插件不能识别Neovim中的Python主机
4 回答
尽管python字典包含了大量的条目,但它并没有增长
4 回答
尽管python说模块存在,为什么我会得到这个消息?
2 回答
尽管setuptools和控制盘是最新的,但无法识别singleversionexternallymanaged
1 回答
尽管stdout和stderr重定向,但未捕获错误消息
1 回答
尽管Tensorboard的事件太大,但Tensorboard的步骤太少了
9 回答
尽管tkinter上的变量已更改,但显示未更改
6 回答
尽管try/except使用Python进行单元测试时出现断言错误
1 回答
尽管URL是sam,但仍会抛出“达到最大重定向”
10 回答
尽管url有效,Pandas仍读取url的\u csv错误
2 回答
尽管while中存在时间延迟,但LINUX线程的CPU利用率为100%(1)
4 回答
尽管x0在范围内,Scipy优化仍会引发ValueError
4 回答
尽管xpath正确,但使用selenium单击链接仍不起作用
3 回答
尽管下载了ffmpeg并设置了路径变量python,但没有后端错误
6 回答
尽管下载了i,但找不到型号“fr”
9 回答
尽管下载了plotnine包,但未获取名为“plotnine”的模块时出错
9 回答
尽管为所有行指定了权重,网格(0)仍不起作用
8 回答