有1条曲线,我想找到经过这个曲线过零点的所有点。python示例代码

2024-01-16 13:28

本文主要是介绍有1条曲线,我想找到经过这个曲线过零点的所有点。python示例代码,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

有1条曲线,我想找到经过这个曲线过零点的所有点。python示例代码

import numpy as np
import matplotlib.pyplot as plt# 假设这是你的曲线数据
x = np.linspace(-10, 10, 1000)  # 生成x轴的值
y = np.sin(x)  # 生成y轴的值,即曲线的值# 找到y值符号变化的位置
## np.sign() 如果数组值大于0则返回1,如果数组值小于0则返回-1,如果数组值为0则返回0
## np.diff() 差分
sign_changes = np.where(np.diff(np.sign(y)))[0]  # eg: array([ 28, 185, 342, 499, 656, 813, 970], dtype=int64)# 过零点的x坐标
zero_crossings_x = x[sign_changes] # array([-9.43943944, -6.2962963 , -3.15315315, -0.01001001,  3.13313313, 6.27627628,  9.41941942])# 由于过零点可能不会精确落在采样点上,我们可以通过线性插值来估计过零点的精确位置
zero_crossings_y = []
for i in sign_changes:# 线性插值zero_crossing = y[i] - (y[i+1] - y[i]) / (x[i+1] - x[i]) * y[i]zero_crossings_y.append(zero_crossing)# 打印结果
print("过零点的x坐标:", zero_crossings_x)
print("过零点的y坐标 (应接近0):", zero_crossings_y)# 绘制曲线和过零点
plt.plot(x, y, label='Curve')
plt.scatter(zero_crossings_x, zero_crossings_y, color='red', label='Zero Crossings')
plt.axhline(0, color='black', linewidth=0.5)
plt.legend()
plt.show()
## 输出结果:
过零点的x坐标: [-9.43943944 -6.2962963  -3.15315315 -0.01001001  3.13313313  6.276276289.41941942]
过零点的y坐标 (应接近0): [0.02932150337482438, -2.8198199947780667e-07, 0.02312027717766069, -1.6716403914568412e-07, 0.016918687676300323, -1.485976876530043e-07, 0.010716883956857083]

加粗样式

这篇关于有1条曲线,我想找到经过这个曲线过零点的所有点。python示例代码的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/612737

相关文章

python: 多模块(.py)中全局变量的导入

文章目录 global关键字可变类型和不可变类型数据的内存地址单模块(单个py文件)的全局变量示例总结 多模块(多个py文件)的全局变量from x import x导入全局变量示例 import x导入全局变量示例 总结 global关键字 global 的作用范围是模块(.py)级别: 当你在一个模块(文件)中使用 global 声明变量时,这个变量只在该模块的全局命名空

活用c4d官方开发文档查询代码

当你问AI助手比如豆包,如何用python禁止掉xpresso标签时候,它会提示到 这时候要用到两个东西。https://developers.maxon.net/论坛搜索和开发文档 比如这里我就在官方找到正确的id描述 然后我就把参数标签换过来

【Python编程】Linux创建虚拟环境并配置与notebook相连接

1.创建 使用 venv 创建虚拟环境。例如,在当前目录下创建一个名为 myenv 的虚拟环境: python3 -m venv myenv 2.激活 激活虚拟环境使其成为当前终端会话的活动环境。运行: source myenv/bin/activate 3.与notebook连接 在虚拟环境中,使用 pip 安装 Jupyter 和 ipykernel: pip instal

poj 1258 Agri-Net(最小生成树模板代码)

感觉用这题来当模板更适合。 题意就是给你邻接矩阵求最小生成树啦。~ prim代码:效率很高。172k...0ms。 #include<stdio.h>#include<algorithm>using namespace std;const int MaxN = 101;const int INF = 0x3f3f3f3f;int g[MaxN][MaxN];int n

【机器学习】高斯过程的基本概念和应用领域以及在python中的实例

引言 高斯过程(Gaussian Process,简称GP)是一种概率模型,用于描述一组随机变量的联合概率分布,其中任何一个有限维度的子集都具有高斯分布 文章目录 引言一、高斯过程1.1 基本定义1.1.1 随机过程1.1.2 高斯分布 1.2 高斯过程的特性1.2.1 联合高斯性1.2.2 均值函数1.2.3 协方差函数(或核函数) 1.3 核函数1.4 高斯过程回归(Gauss

【学习笔记】 陈强-机器学习-Python-Ch15 人工神经网络(1)sklearn

系列文章目录 监督学习:参数方法 【学习笔记】 陈强-机器学习-Python-Ch4 线性回归 【学习笔记】 陈强-机器学习-Python-Ch5 逻辑回归 【课后题练习】 陈强-机器学习-Python-Ch5 逻辑回归(SAheart.csv) 【学习笔记】 陈强-机器学习-Python-Ch6 多项逻辑回归 【学习笔记 及 课后题练习】 陈强-机器学习-Python-Ch7 判别分析 【学

计算机毕业设计 大学志愿填报系统 Java+SpringBoot+Vue 前后端分离 文档报告 代码讲解 安装调试

🍊作者:计算机编程-吉哥 🍊简介:专业从事JavaWeb程序开发,微信小程序开发,定制化项目、 源码、代码讲解、文档撰写、ppt制作。做自己喜欢的事,生活就是快乐的。 🍊心愿:点赞 👍 收藏 ⭐评论 📝 🍅 文末获取源码联系 👇🏻 精彩专栏推荐订阅 👇🏻 不然下次找不到哟~Java毕业设计项目~热门选题推荐《1000套》 目录 1.技术选型 2.开发工具 3.功能

nudepy,一个有趣的 Python 库!

更多资料获取 📚 个人网站:ipengtao.com 大家好,今天为大家分享一个有趣的 Python 库 - nudepy。 Github地址:https://github.com/hhatto/nude.py 在图像处理和计算机视觉应用中,检测图像中的不适当内容(例如裸露图像)是一个重要的任务。nudepy 是一个基于 Python 的库,专门用于检测图像中的不适当内容。该

代码随想录冲冲冲 Day39 动态规划Part7

198. 打家劫舍 dp数组的意义是在第i位的时候偷的最大钱数是多少 如果nums的size为0 总价值当然就是0 如果nums的size为1 总价值是nums[0] 遍历顺序就是从小到大遍历 之后是递推公式 对于dp[i]的最大价值来说有两种可能 1.偷第i个 那么最大价值就是dp[i-2]+nums[i] 2.不偷第i个 那么价值就是dp[i-1] 之后取这两个的最大值就是d

pip-tools:打造可重复、可控的 Python 开发环境,解决依赖关系,让代码更稳定

在 Python 开发中,管理依赖关系是一项繁琐且容易出错的任务。手动更新依赖版本、处理冲突、确保一致性等等,都可能让开发者感到头疼。而 pip-tools 为开发者提供了一套稳定可靠的解决方案。 什么是 pip-tools? pip-tools 是一组命令行工具,旨在简化 Python 依赖关系的管理,确保项目环境的稳定性和可重复性。它主要包含两个核心工具:pip-compile 和 pip