在 numpy 中表示一阶微分方程

2024-01-08

我有一个等式dy/dx = x + y/5和一个初始值,y(0) = -3.

我想知道如何使用 pyplot 绘制该函数的精确图形。

我也有一个x = np.linspace(0, interval, steps+1)我想用它作为 x 轴。所以我只寻找 y 轴值。

提前致谢。


为了完整起见,这种方程可以很容易地进行数值积分,使用scipy.integrate.odeint.

import numpy as np
from scipy.integrate import odeint
import matplotlib.pyplot as plt

# function dy/dx = x + y/5.
func = lambda y,x : x + y/5.
# Initial condition
y0 = -3  # at x=0
# values at which to compute the solution (needs to start at x=0)
x = np.linspace(0, 4, 101)
# solution
y = odeint(func, y0, x)
# plot the solution, note that y is a column vector
plt.plot(x, y[:,0])
plt.xlabel('x')
plt.ylabel('y')
plt.show()
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

在 numpy 中表示一阶微分方程 的相关文章

  • 删除网格线,但保留框架(matplotlib 中的 ggplot2 样式)

    使用 Matplotlib 我想删除图中的网格线 同时保留框架 即轴线 我已经尝试过下面的代码和其他选项 但我无法让它工作 如何在删除网格线的同时保留框架 我这样做是为了在 matplotlib 中重现 ggplot2 图 我在下面创建了一
  • 如何将 python 点列表转换为 numpy 图像数组?

    我有一个 python 点列表 x y 坐标 200 245 344 248 125 34 它表示二维平面上的轮廓 我想使用一些 numpy scipy 算法进行平滑 插值等 它们通常需要 numpy 数组作为输入 例如scipy ndim
  • 从 MySQL 将数字数据加载到 python/pandas/numpy 数组的最快方法

    我想从 MySQL 表中读取一些数字 双精度 即 float64 数据 数据大小约为 200k 行 MATLAB 参考 tic feature accel off conn database c fetch exec conn select
  • 自定义seaborn jointplot中的轴标签

    我似乎陷入了一个相对简单的问题 但在经过一小时的搜索和大量的实验后无法解决它 我有两个 numpy 数组x and y我正在使用seaborn的联合图来绘制它们 sns jointplot x y 现在我想分别将 x 轴和 y 轴标记为 X
  • 在 python 中高效、快速地迭代元组列表中超过 3600 万个项目

    首先 在任何人将其标记为重复之前 请阅读以下内容 我不确定迭代的延迟是否是由于尺寸巨大或我的逻辑造成的 我有一个必须迭代的用例3600 万件商品在元组列表中 我的主要要求是速度和效率 样本清单 how are you I am fine h
  • 有一些 numpy.map 吗?

    我可能在这里遗漏了一些明显的东西 但我遗漏了一个功能numpy map 这与 Python 的相同map函数 但将输出收集在numpy大批 例如 我可以有一个图像生成器genImage i 生成 2D 图像 大小 m n 基于单个输入 我想
  • matplotlib 中无效的 rgba arg“#”

    我无法弄清楚在尝试在 matplotlib 中创建散点图时如何使用颜色 我正在尝试绘制具有不同颜色点的多个散点图来显示集群 colors 12efff eee111 eee00f e00fff 123456 abc222 000000 12
  • 内存高效的随机数迭代器,无需替换

    我觉得这应该很容易 但经过多次搜索和尝试后我找不到答案 基本上 我有大量的物品 我想以随机顺序进行采样 而不需要更换 在本例中 它们是二维数组中的单元 我用于较小数组的解决方案不会转换 因为它需要对内存数组进行改组 如果我必须采样的数量很小
  • Python 多元简单线性回归

    注意这是not关于多元回归的问题 这是一个关于在 Python NumPy 2 7 中多次进行简单 单变量 回归的问题 我有两个m x n arrays x and y 这些行彼此对应 每对都是用于测量的 x y 点的集合 那是 plt p
  • 使用自定义颜色渐变填充两条线之间的区域

    我正在做一项几乎已经完成的作业 但我想对其添加一些小改动 尝试使用基于温度的颜色图而不是简单的颜色来填充两条线之间的区域 绘制线条的方式本质上使它们成为独立的实体 所以我知道我可能需要两个彼此相遇或重叠的颜色图来完成此任务 但我不太确定如何
  • 在 x 轴上操作日期 Pandas Matplotlib

    我有一组非常简单的数据 如下所示 我正在寻找一种方法来绘制此堆积条形图并格式化 x 轴 日期 使其从 1996 年 31 月 12 日开始 到 2016 年 31 月 12 日结束 增量为 365 天 我编写的代码绘制了每个日期 因此 x
  • “扩展”numpy ndarray 的好方法?

    有没有 扩展 numpy ndarray 的好方法 假设我有一个像这样的 ndarray 1 2 3 4 我希望每行通过填充零来包含更多元素 1 2 0 0 0 3 4 0 0 0 我知道一定有一些蛮力的方法可以做到这一点 比如构造一个带有
  • 使用 Matplotlib、PyQt 和 Threading 进行实时绘图导致 python 崩溃

    我一直在努力研究我的 Python 应用程序 但找不到任何答案 我有 PyQT GUI 应用程序 它使用 Matplotlib 小部件 GUI 启动一个新线程来处理 mpl 小部件的绘图 恐怕我现在通过从另一个线程访问 matplotlib
  • 如何计算伽罗瓦域上的numpy数组?

    我想在伽罗华域 GF4 上使用 numpy 数组 所以 我将 GF4 类设置为数组元素 它适用于数组 整数计算 但不适用于数组 数组计算 import numpy class GF4 object class for galois fiel
  • 如何在刻度标签和轴之间添加空间

    我已成功增加刻度标签的字体 但现在它们距离轴太近了 我想在刻度标签和轴之间添加一点呼吸空间 如果您不想全局更改间距 通过编辑 rcParams 并且想要更简洁的方法 请尝试以下操作 ax tick params axis both whic
  • 将 github 上的包安装到 Spyder 中

    我一直在尝试安装并导入mpl finance来自 github 的包 在我的 Spyder 环境中没有成功 我努力了 pip install e git https github com matplotlib mpl finance git
  • 如何将交互式 matplotlib 图形插入 tkinter 画布

    我正在尝试将交互式 matplotlib 图形 具有滑块 重置按钮和单选按钮的图形 放入 tkinter Canvas 中 我已成功添加非交互式图表 但当它变为交互式时找不到问题 我尝试将所有内容更改为使用 matplotlib Figur
  • 如何将 numpy.matrix 提高到非整数幂?

    The 运算符为numpy matrix不支持非整数幂 gt gt gt m matrix 1 0 0 5 0 5 gt gt gt m 2 5 TypeError exponent must be an integer 我想要的是 oct
  • Numpy 优化

    我有一个根据条件分配值的函数 我的数据集大小通常在 30 50k 范围内 我不确定这是否是使用 numpy 的正确方法 但是当数字超过 5k 时 它会变得非常慢 有没有更好的方法让它更快 import numpy as np N 5000
  • Pandas 与 Numpy 数据帧

    看这几行代码 df2 df copy df2 1 df 1 df 1 values 1 df2 ix 0 0 我们的教练说我们需要使用 values属性来访问底层的 numpy 数组 否则我们的代码将无法工作 我知道 pandas Data

随机推荐

  • 在没有焦点的情况下监视 tkinter 中的按键

    我正在使用 tkinter 在 Python 中编写一个小型计时模块 在这方面 我想全局监控何时按下退出键来停止计时 不幸的是 tkinters bind 和 bind all 函数仅在窗口处于焦点时才拾取击键 我已经研究了其他几个用于记录
  • 尝试在android中以编程方式模拟滑动

    我希望我的应用程序在单击按钮时模拟滑动触摸事件 向上 向下 向左 向右 然后 TextView 将向下 向上滚动 我尝试过使用 Motion Event 但是在我分别调度 ACTION DOWN ACTION MOVE 和 ACTION U
  • 下载数据库文件到app目录

    我想在使用 DownloadManager 单击按钮时更新应用程序的 sqlite 数据库 但它说 java lang IllegalArgumentException 不是文件URI data user 0 com example lau
  • 使用 asyncio 逐行读取文件

    我希望在写入多个日志文件时读取它们并使用 asyncio 处理它们的输入 该代码必须在 Windows 上运行 根据我在 stackoverflow 和网络上的搜索了解到 异步文件 I O 在大多数操作系统上都很棘手 select例如 将无
  • 使用自定义字体时如何垂直纠正导航栏的 titleView 文本位置?

    我们在导航栏中的 titleView 使用自定义字体 不知何故 苹果总是把这种字体画得太高 当您在导航栏中使用自定义字体时 如何纠正这种奇怪的偏移 I used setTitleVerticalPositionAdjustment forB
  • R:Rvest - 隐藏了我不想要的文字

    我正在对这个网络进行网络抓取 http www falabella com pe falabella pe category cat40536 Climatizacion navAction push 我只需要产品的信息 品牌 产品名称 价
  • 如何使用 Python 在 Windows 中创建符号链接?

    我正在尝试创建symlinks在 Windows 8 上使用 Python 我发现这个帖子 https stackoverflow com questions 1447575 symlinks on windows这是我脚本的一部分 imp
  • 如何将现有的kafka主题分区分散到更多目录中?

    默认情况下 kafka使用一个目录来保存日志 为了提高性能 建议将更多磁盘挂载到代理 并将每个磁盘分配给一个目录 然后在server properties输入log dirs 作为逗号分隔的目录列表 文档称 分区将以循环方式分布在目录中 据
  • 经典 ASP 有 OR/M 吗?

    是否有可在经典 ASP 中使用的 OR M 对象关系映射器 即使是简化的类对象也会对处理简单的 CRUD 任务有很大帮助 是的 我知道 ASP NET 有很多 我将其中一些用于 ASP NET 站点 然而 这是一个使用 ASP VBScri
  • Firebase 的云功能 - 在Write 上发送电子邮件

    当有任何内容写入时 我正在尝试向我的电子邮件发送一封测试电子邮件 emails但电子邮件永远不会发送 并且功能日志为空 exports sendTestEmail functions database ref emails onWrite
  • 如何指定绘图应转到哪个图形?

    我打开了多个图形 我想在运行时独立更新它们 下面的玩具示例应该可以阐明我的意图 clf figure name a and b a and b should be plotted to this window hold on ylim 10
  • 缩小浏览器时,响应式页脚显示在部分内容上

    我对 CSS 和 HTML 还很陌生 据我所知 这个问题已被问过多次 但我似乎无法解决当前的问题 我一直在从事一些小项目 比如尝试重新创建具有额外功能的其他网站 我当前的任务是重新创建一个首页搜索区域 例如 google 实施与发布 我一直
  • SharePoint:看不到任何筛选器 Web 部件

    在 Office 2007 官方网站中 有大量对筛选 Web 部件的引用 当我尝试添加其中之一时 Sharepoint 中的 Web 部件列表不显示任何筛选器 Web 部件 如果有人经历过相同的问题和任何解决方案 请发帖 Thanks 在
  • Perl 中美元符号和数字符号如何一起工作?

    今天遇到一个问题 要求我在perl中确定数组的最大索引 我曾经这样做过 my array 1 2 3 print array array 1 但今天我偶然发现了这段代码 my array 1 2 3 print array array 我在
  • 使用 xattr 设置 Mac OSX 隔离属性

    StackOverflow 和其他地方有很多关于如何清除 Mac 隔离属性的信息 就我而言 我想设置它 这是为了测试我的应用程序是否已正确签名 以便用户在下载后会收到 不受信任的开发人员 警告 我的应用程序特别大 我们从大文件下载网站分发
  • JTable 行标题文本

    有人可以告诉我是否有一种方法可以添加带有文本的行标题 或者也许有某种类型的技巧 这将允许我们做这样的事情 行标题是 one way https stackoverflow com a 8187799 714968 更好的方法 http ti
  • 如何抑制“void*”到“foo*”转换的警告(通过 -fpermissive 减少错误)

    我正在尝试用 g 编译一些 c 代码 是的 故意的 我收到如下错误 例如 error invalid conversion from void to unsigned char fpermissive adding fpermissive到
  • Material-UI 和 gatsby 官方示例

    我根据官方示例使用 gatsby 和 Material UI 我想知道是否必须添加withRoot包装到每个使用 MUI 的组件 或者我是否可以包装我想要使用 MUI 的最高级别组件 我在渲染时遇到问题 不知道如何找出原因 在开发中一切都很
  • Kafka 消费者:如果未提交先前的消息偏移并且禁用自动提交,则希望再次读取相同的消息

    我已经关闭了自动提交 并且在阅读后也没有提交来自消费者的偏移量 检查的消费者滞后也保持不变 它确保偏移量不会被提交 但问题是 它正在消耗下一条消息 而不是再次使用相同的消息 我怎样才能一遍又一遍地阅读同一条消息 仅当先前的偏移量已提交时 我
  • 在 numpy 中表示一阶微分方程

    我有一个等式dy dx x y 5和一个初始值 y 0 3 我想知道如何使用 pyplot 绘制该函数的精确图形 我也有一个x np linspace 0 interval steps 1 我想用它作为 x 轴 所以我只寻找 y 轴值 提前