python – NumPy odeint输出额外的变量

内容导读

互联网集市收集整理的这篇技术教程文章主要介绍了python – NumPy odeint输出额外的变量，小编现在分享给大家，供广大互联网技能从业者学习和参考。文章包含1187字，纯文字阅读大概需要2分钟。

内容图文

在Numpy中使用odeint进行模拟期间保存中间变量的最简单方法是什么？

例如：

def dy(y,t)
    x = np.rand(3,1)
    return y + x.sum()

sim = odeint(dy,0,np.arange(0,1,0.1))

在模拟过程中保存x中存储数据的最简单方法是什么？理想情况下,在传递给odeint的t参数中指定的点.

解决方法:

通过一些注意事项来解决odeint的一种方便方法是在类中使用dy将调用包含在odeint中,使用dy作为另一种方法,并将self作为参数传递给dy函数.例如,

class WrapODE(object):
    def __init__(self):
        self.y_0 = 0.
        self.L_x = []
        self.timestep = 0
        self.times = np.arange(0., 1., 0.1)

    def run(self):
        self.L_y = odeint(
            self.dy,
            self.y_0, self.times,
            args=(self,))

    @staticmethod
    def dy(y, t, self):
        """"
        Discretized application of dudt

        Watch out! Because this is a staticmethod, as required by odeint, self
        is the third argument
        """
        x = np.random.rand(3,1)
        if t >= self.times[self.timestep]:
            self.timestep += 1
            self.L_x.append(x)
        else:
            self.L_x[-1] = x
        return y + x.sum()

要清楚,这是一个容易陷入困境的黑客.例如,除非odeint执行Euler步进,否则dy将被调用的次数多于您指定的时间步数.为了确保每个y得到一个x,if t> = self.times [self.timestep]：block中的monkey business选择一个数组中的一个点来存储来自时间向量的每个时间值的数据.您的特定应用可能会导致其他疯狂的问题.请务必为您的应用程序彻底验证此方法.

内容总结

以上是互联网集市为您收集整理的python – NumPy odeint输出额外的变量全部内容，希望文章能够帮你解决python – NumPy odeint输出额外的变量所遇到的程序开发问题。 如果觉得互联网集市技术教程内容还不错，欢迎将互联网集市网站推荐给程序员好友。

内容备注

版权声明：本文内容由互联网用户自发贡献，该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容， 请发送邮件至 gblab@vip.qq.com 举报，一经查实，本站将立刻删除。

内容手机端