scipy.solve_ivp:在第二次发生事件时终止轨迹

如何解决scipy.solve_ivp:在第二次发生事件时终止轨迹

我想解决返回初始状态所需的时间。作为一个玩具示例,让我们看一个简单的谐波振荡器d / dt(x,y)=(y,-x),初始条件为(x,y)=(1,0)。

我希望代码返回T = 6.28 ...(=2π),因为这是返回y为正的初始条件所花费的时间。当达到这一点时,它应该停止集成。

import numpy as np
from scipy.integrate import solve_ivp

def f(t,x):
    return [x[1],-x[0]]

def event(t,x):
    return x[0]

event.direction = 1
event.terminal = True # Try false,this gives the result but integrates beyond the event I need.

initial_condition = [0,1]

sol = solve_ivp(f,[0,10],initial_condition,method='RK45',dense_output=True,\
                events = event,rtol=1e-6,atol=3e-6)

print(sol.t_events)

麻烦是将event.terminal = True设置为立即启动积分,因为事件在t = 0时得到满足。是否有办法让积分仅在两次事件都得到满足后才停止?

解决方法

与您的解决方法类似,我将 1e-100 添加到事件而不是初始状态。这样做的好处是根本不影响轨迹,甚至理论上也不影响。

当然,这确实意味着您未来的事件会提前一小部分发生,但如果这对您很重要,则取决于您的情况。

对于您编写的示例代码,这将使您的事件如下:

def event(t,x):
    return x[0] + 1e-100

我尝试过的一件事不起作用是在集成过程中修改终止参数。我想如果我从 event.terminal = False 开始,然后在事件定义集 event.terminal = True 中,那么它将在第一个事件期间继续(在 t=0 时),并在第二个事件处终止。但是,如果在集成过程中更改了更改的参数,它似乎会忽略。

(这个问题与 Solve_ivp integration gets stuck if initial condition triggers event with terminal = True 相关,但我没有标记任何东西的声誉。)

,

为以后参考,这里是我一直在使用的解决方法:

将初始条件设置为事件发生前的一小部分,例如zoo::na.locf(),而不是library(data.table) dtA <- data.table(Time = seq(as.ITime("00:15:00"),by = 900L,length.out = 96L),Values1 = seq(1L,by = 2L,Values2 = seq(2L,Values3 = seq(1.5,by = 1.0,Category = rep("A",96L)) dtB <- data.table(Time = seq(as.ITime("00:15:00"),Values1 = seq(to = 192L,Values2 = seq(to = 191L,Values3 = seq(to = 97.5,Category = rep("B",96L)) dt <- rbind(dtA,dtB) 。即使对于更高自由度的混沌系统,它似乎也不会影响轨迹,但是它不会在t = 0时触发事件。

另一个可行的解决方法是,在一段时间内不进行任何事件集成,然后将结果用作进一步集成的初始条件,这次是在事件已就绪的情况下进行。我怀疑这个答案对任何人都有用,但我只是为了完整性而发布了它。

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

相关推荐


依赖报错 idea导入项目后依赖报错,解决方案:https://blog.csdn.net/weixin_42420249/article/details/81191861 依赖版本报错:更换其他版本 无法下载依赖可参考:https://blog.csdn.net/weixin_42628809/a
错误1:代码生成器依赖和mybatis依赖冲突 启动项目时报错如下 2021-12-03 13:33:33.927 ERROR 7228 [ main] o.s.b.d.LoggingFailureAnalysisReporter : *************************** APPL
错误1:gradle项目控制台输出为乱码 # 解决方案:https://blog.csdn.net/weixin_43501566/article/details/112482302 # 在gradle-wrapper.properties 添加以下内容 org.gradle.jvmargs=-Df
错误还原:在查询的过程中,传入的workType为0时,该条件不起作用 &lt;select id=&quot;xxx&quot;&gt; SELECT di.id, di.name, di.work_type, di.updated... &lt;where&gt; &lt;if test=&qu
报错如下,gcc版本太低 ^ server.c:5346:31: 错误:‘struct redisServer’没有名为‘server_cpulist’的成员 redisSetCpuAffinity(server.server_cpulist); ^ server.c: 在函数‘hasActiveC
解决方案1 1、改项目中.idea/workspace.xml配置文件,增加dynamic.classpath参数 2、搜索PropertiesComponent,添加如下 &lt;property name=&quot;dynamic.classpath&quot; value=&quot;tru
删除根组件app.vue中的默认代码后报错:Module Error (from ./node_modules/eslint-loader/index.js): 解决方案:关闭ESlint代码检测,在项目根目录创建vue.config.js,在文件中添加 module.exports = { lin
查看spark默认的python版本 [root@master day27]# pyspark /home/software/spark-2.3.4-bin-hadoop2.7/conf/spark-env.sh: line 2: /usr/local/hadoop/bin/hadoop: No s
使用本地python环境可以成功执行 import pandas as pd import matplotlib.pyplot as plt # 设置字体 plt.rcParams[&#39;font.sans-serif&#39;] = [&#39;SimHei&#39;] # 能正确显示负号 p
错误1:Request method ‘DELETE‘ not supported 错误还原:controller层有一个接口,访问该接口时报错:Request method ‘DELETE‘ not supported 错误原因:没有接收到前端传入的参数,修改为如下 参考 错误2:cannot r
错误1:启动docker镜像时报错:Error response from daemon: driver failed programming external connectivity on endpoint quirky_allen 解决方法:重启docker -&gt; systemctl r
错误1:private field ‘xxx‘ is never assigned 按Altʾnter快捷键,选择第2项 参考:https://blog.csdn.net/shi_hong_fei_hei/article/details/88814070 错误2:启动时报错,不能找到主启动类 #
报错如下,通过源不能下载,最后警告pip需升级版本 Requirement already satisfied: pip in c:\users\ychen\appdata\local\programs\python\python310\lib\site-packages (22.0.4) Coll
错误1:maven打包报错 错误还原:使用maven打包项目时报错如下 [ERROR] Failed to execute goal org.apache.maven.plugins:maven-resources-plugin:3.2.0:resources (default-resources)
错误1:服务调用时报错 服务消费者模块assess通过openFeign调用服务提供者模块hires 如下为服务提供者模块hires的控制层接口 @RestController @RequestMapping(&quot;/hires&quot;) public class FeignControl
错误1:运行项目后报如下错误 解决方案 报错2:Failed to execute goal org.apache.maven.plugins:maven-compiler-plugin:3.8.1:compile (default-compile) on project sb 解决方案:在pom.
参考 错误原因 过滤器或拦截器在生效时,redisTemplate还没有注入 解决方案:在注入容器时就生效 @Component //项目运行时就注入Spring容器 public class RedisBean { @Resource private RedisTemplate&lt;String
使用vite构建项目报错 C:\Users\ychen\work&gt;npm init @vitejs/app @vitejs/create-app is deprecated, use npm init vite instead C:\Users\ychen\AppData\Local\npm-