首页 > 后端开发 > Python教程 > 如何在 Python 中使用线性插值在绘图上找到精确的零和非零截距?

如何在 Python 中使用线性插值在绘图上找到精确的零和非零截距?

Susan Sarandon
发布: 2024-10-21 07:12:29
原创
782 人浏览过

How to Find Exact Zero and Non-Zero Intercepts on Plots Using Linear Interpolation in Python?

查找曲线与零的交点

在 Python 中,当值不是一个值时,从图中获取精确的 y 轴值可能会很困难。整数。本文解决了这个问题,并提出了一种基于线性插值的解决方案。

给定两个数组(vertical_data 和gradient(Temperature_data)),使用 plt.plot 生成绘图。然而,该图显示的 y 值接近但不完全为零。

用于根估计的线性插值

估计 numpy 数组的精确根,可以使用简单的线性插值方法。以下代码演示了如何查找任意曲线的零值:

<code class="python">import numpy as np

def find_roots(x, y):
    s = np.abs(np.diff(np.sign(y))).astype(bool)
    return x[:-1][s] + np.diff(x)[s]/(np.abs(y[1:][s]/y[:-1][s])+1)

x = .4+np.sort(np.random.rand(750))*3.5
y = (x-4)*np.cos(x*9.)*np.cos(x*6+0.05)+0.1

z = find_roots(x,y)

plt.plot(x,y)
plt.plot(z, np.zeros(len(z)), marker="o", ls="", ms=4)</code>
登录后复制

此代码标识曲线的根,并将它们在精确的 y 值为零处将它们绘制为圆。

非零截距

通过修改求根的线,可以使用相同的方法来查找曲线与任何非零 y 值 (y0) 的交点:

<code class="python">z = find_roots(x,y-y0)</code>
登录后复制

两条曲线的交点

线性插值法也可以用来求两条曲线的交点。通过找到两条曲线之间差异的根,我们可以估计它们的交点:

<code class="python">y2 = (x - 2) * np.cos(x * 8.) * np.cos(x * 5 + 0.03) + 0.3

z = find_roots(x,y2-y1)

plt.plot(x,y1)
plt.plot(x,y2, color="C2")
plt.plot(z, np.interp(z, x, y1), marker="o", ls="", ms=4, color="C1")</code>
登录后复制

以上是如何在 Python 中使用线性插值在绘图上找到精确的零和非零截距?的详细内容。更多信息请关注PHP中文网其他相关文章!

来源:php
本站声明
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系admin@php.cn
作者最新文章
热门教程
更多>
最新下载
更多>
网站特效
网站源码
网站素材
前端模板