有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的Darts库实现时间序列预测

《Python的Darts库实现时间序列预测》Darts一个集统计、机器学习与深度学习模型于一体的Python时间序列预测库,本文主要介绍了Python的Darts库实现时间序列预测,感兴趣的可以了解... 目录目录一、什么是 Darts?二、安装与基本配置安装 Darts导入基础模块三、时间序列数据结构与

Python正则表达式匹配和替换的操作指南

《Python正则表达式匹配和替换的操作指南》正则表达式是处理文本的强大工具,Python通过re模块提供了完整的正则表达式功能,本文将通过代码示例详细介绍Python中的正则匹配和替换操作,需要的朋... 目录基础语法导入re模块基本元字符常用匹配方法1. re.match() - 从字符串开头匹配2.

Python使用FastAPI实现大文件分片上传与断点续传功能

《Python使用FastAPI实现大文件分片上传与断点续传功能》大文件直传常遇到超时、网络抖动失败、失败后只能重传的问题,分片上传+断点续传可以把大文件拆成若干小块逐个上传,并在中断后从已完成分片继... 目录一、接口设计二、服务端实现(FastAPI)2.1 运行环境2.2 目录结构建议2.3 serv

C#实现千万数据秒级导入的代码

《C#实现千万数据秒级导入的代码》在实际开发中excel导入很常见,现代社会中很容易遇到大数据处理业务,所以本文我就给大家分享一下千万数据秒级导入怎么实现,文中有详细的代码示例供大家参考,需要的朋友可... 目录前言一、数据存储二、处理逻辑优化前代码处理逻辑优化后的代码总结前言在实际开发中excel导入很

通过Docker容器部署Python环境的全流程

《通过Docker容器部署Python环境的全流程》在现代化开发流程中,Docker因其轻量化、环境隔离和跨平台一致性的特性,已成为部署Python应用的标准工具,本文将详细演示如何通过Docker容... 目录引言一、docker与python的协同优势二、核心步骤详解三、进阶配置技巧四、生产环境最佳实践

Python一次性将指定版本所有包上传PyPI镜像解决方案

《Python一次性将指定版本所有包上传PyPI镜像解决方案》本文主要介绍了一个安全、完整、可离线部署的解决方案,用于一次性准备指定Python版本的所有包,然后导出到内网环境,感兴趣的小伙伴可以跟随... 目录为什么需要这个方案完整解决方案1. 项目目录结构2. 创建智能下载脚本3. 创建包清单生成脚本4

SpringBoot+RustFS 实现文件切片极速上传的实例代码

《SpringBoot+RustFS实现文件切片极速上传的实例代码》本文介绍利用SpringBoot和RustFS构建高性能文件切片上传系统,实现大文件秒传、断点续传和分片上传等功能,具有一定的参考... 目录一、为什么选择 RustFS + SpringBoot?二、环境准备与部署2.1 安装 RustF

Python实现Excel批量样式修改器(附完整代码)

《Python实现Excel批量样式修改器(附完整代码)》这篇文章主要为大家详细介绍了如何使用Python实现一个Excel批量样式修改器,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一... 目录前言功能特性核心功能界面特性系统要求安装说明使用指南基本操作流程高级功能技术实现核心技术栈关键函

python获取指定名字的程序的文件路径的两种方法

《python获取指定名字的程序的文件路径的两种方法》本文主要介绍了python获取指定名字的程序的文件路径的两种方法,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要... 最近在做项目,需要用到给定一个程序名字就可以自动获取到这个程序在Windows系统下的绝对路径,以下

使用Python批量将.ncm格式的音频文件转换为.mp3格式的实战详解

《使用Python批量将.ncm格式的音频文件转换为.mp3格式的实战详解》本文详细介绍了如何使用Python通过ncmdump工具批量将.ncm音频转换为.mp3的步骤,包括安装、配置ffmpeg环... 目录1. 前言2. 安装 ncmdump3. 实现 .ncm 转 .mp34. 执行过程5. 执行结