我不确定事件处理是否scipy.integrate.solve_ivp工作正常。在下面的例子中,我对一个导数进行积分,得到一个三次多项式,它的根在x=-6,x=-2和x=2。我设置了一个事件函数,返回y,在x值处为零。我希望在解决方案的t_events属性中看到三个条目,但是我只看到一个条目,尽管很明显解决方案跨越x轴三次。在

我是不是做错了什么?在

复制代码示例:def fprime(x, y):

return 3 * x**2 + 12 * x - 4

def event(x, y):

return y

import numpy as np

from scipy.integrate import solve_ivp

sol = solve_ivp(fprime, (-8, 4), np.array([-120]), t_eval=np.linspace(-8, 4, 10), events=[event])

上述代码导致:

^{pr2}$

```

问题是从sol.y数组中可以看到应该有三个0(有三个符号改变),但是只记录了一个事件。在

Scipy/Numpy/Python版本信息:1.0.0 1.13.3 sys.version_info(major=3, minor=6, micro=0, releaselevel='final', serial=0)

[UPDATE]:如果我使用可选的max_step参数来求解_ivp,并使其足够小,那么我会看到所有三个根。似乎事件函数并不是在t帴eval步骤中调用的,而是仅在内部解算器步骤上调用,后者远远少于t帴eval步骤,并且最终会跳过某些根。这似乎不是很有用,因为你必须知道如何设置最大步长以避免丢失根。在

python scipy库函数solve用法_如何在中使用事件scipy.integrate.solve_ivp相关推荐

  1. 下列变量名不符合python命名规范的是_下列选项中不符合Python语言变量命名规则的是...

    [单选题]下列Python保留字中,不用于表示分支结构的是 [单选题]以下程序的输出结果是: j = '' for i in "12345": j += i + ',' print ...

  2. react中绑定点击事件_在React中绑定事件处理程序的最佳方法

    react中绑定点击事件 by Charlee Li 通过李李 在React中绑定事件处理程序的最佳方法 (The best way to bind event handlers in React) ...

  3. python双样本t检验_两样本t检验的scipy实现

    以真实商业案例为数据基础,课程内容围绕scipy.stats和statsmodels包的相关功能展开,从统计分析实战的角度出发详细介绍了如何在Python中完成数据描述.t检验.单因素方差分析.卡方检 ...

  4. python模块调用的用法_如何使用Python语言中的random模块调用方法

    Python语言中的random模块,包含了很多方法属性,可以直接调用这些方法和设置属性.根据不同的方法展示不同的内容.下面利用几个实例查看random模块中的方法,操作如下: 工具/原料 Pytho ...

  5. python查看库函数和方法_查看python库函数和三方库函数文档

    原博文 2020-03-08 16:08 − 1.打开终端,输入python -m pydoc -p 0,(0是让系统自动分配端口,也可以自己指定端口)会返回Server ready at http: ...

  6. python 加载动图_在浏览器中使用TensorFlow.js和Python构建机器学习模型(附代码)...

    大数据文摘授权转载自数据派THU 作者:MOHD SANAD ZAKI RIZVI 本文主要介绍了: TensorFlow.js (deeplearn.js)使我们能够在浏览器中构建机器学习和深度学习 ...

  7. 下列数据类型中python不支持的是_ 下列选项中 ,Python 不支持的数据类型有 ( ) 。_学小易找答案...

    [单选题] 下列标识符中 , 合法的是 ( ) . [简答题]说明轴承代号7204AC表达的含义. [判断题]type() 函数可以查看变量的数据类型. ( ) [名词解释]限界 [单选题]体育教学 ...

  8. python input函数详解_对Python3中的input函数详解

    下面介绍python3中的input函数及其在python2及pyhton3中的不同. python3中的ininput函数,首先利用help(input)函数查看函数信息: 以上信息说明input函 ...

  9. python怎么横着输出_对python3中, print横向输出的方法详解

    对python3中, print横向输出的方法详解 Python 2 : print打印的时候,如果结尾有逗号,打出来时候不会换行.但是在python3里面就不行了. Python3: 3.0的pri ...

最新文章

  1. 调用API的SDK相关知识:实现回调函数.
  2. 2013多校训练赛第三场 总结
  3. Electron中与Vue集成流程
  4. TFboys:使用Tensorflow搭建深层网络分类器
  5. 【Transformer】CrossFormer:A versatile vision transformer based on cross-scale attention
  6. 系出名门Android(7) - 控件(View)之ZoomControls, Include, VideoView, WebView, RatingBar, Tab
  7. java已知一个二叉树_#二叉树复习#
  8. mysql 坚向变横向_(排班表一)使用SQL语句使数据从坚向排列转化成横向排列
  9. docker的php教程https,Docker搭建php环境教程详解
  10. Spring Cloud 分布式链路跟踪 Sleuth + Zipkin + Elasticsearch
  11. DEV-c++中sort的用法
  12. Windows系统安装Mysql前运行库依赖
  13. 【数据结构与算法】详解什么是哈希表,并用代码手动实现一个哈希表
  14. 原子结构示意图全部_原子结构示意图和元素及元素周期表
  15. KEIL MDK平台 S3C2440 编译链接、烧写调试
  16. MATLAB绘制地形图和等高线图
  17. nuvoton uboo2013引导流程 1 - 配置
  18. 大学生php实训总结_php实训报告.doc
  19. 【mysql报错】Data truncation: Data too long for column ‘XXX‘ at row 1
  20. 校招java实习生基础知识面试题

热门文章

  1. framebuffer驱动详解0——framebuffer介绍
  2. 字符串循环右移的一道题目
  3. Android 位置服务——BaiduLocation的使用
  4. MapReduce编程中常用的字符操作
  5. 什么时候需要好教材和好教材是怎样的
  6. 15.使用using和try/finally来做资源清理
  7. 软件项目管理的质量保证
  8. new 与 malloc区别
  9. 设计模式的理解的总结
  10. CPU加了缓存后,有人急了~