如何使用可用的曲率數據集找到曲線的切線？

Question

我有一個曲率數據集，我需要找到曲線的切線。 我的代碼如下，但不幸的是，我沒有得到所需的結果：

chData = efficient.get('Car.Road.y')  
fittedParameters = (np.gradient(chData_m_5[:],1)) # 999 values
plt.plot(chData[1:]) # orginally 1000 values 
plt.plot(fittedParameters)
plt.show()

output 是：

編輯1：

我對代碼進行了以下更改以獲得曲率的切線，但不幸的是，這離曲線有點遠。 請指導我與問題相關的問題解決方案。 謝謝！

fig, ax1 = plt.subplots()
chData_m = efficient.get('Car.Road.y')  

x_fit = chData_m.timestamps
y_fit = chData_m.samples

fittedParameters = np.polyfit(x_fit[:],y_fit[:],1)

f = plt.figure(figsize=(800/100.0, 600/100.0), dpi=100)
axes = f.add_subplot(111)

    # first the raw data as a scatter plot
axes.plot(x_fit, y_fit,  'D')

    # create data for the fitted equation plot
xModel = np.linspace(min(x_fit), max(x_fit))
yModel = np.polyval(fittedParameters, xModel)

    # now the model as a line plot
axes.plot(xModel, yModel)

axes.set_xlabel('X Data') # X axis data label
axes.set_ylabel('Y Data') # Y axis data label

    # polynomial derivative from numpy
deriv = np.polyder(fittedParameters)

    # for plotting
minX = min(x_fit)
maxX = max(x_fit)

    # value of derivative (slope) at a specific X value, so
    # that a straight line tangent can be plotted at the point
    # you might place this code in a loop to animate
pointVal = 10.0 # example X value
y_value_at_point = np.polyval(fittedParameters, pointVal)
slope_at_point = np.polyval(deriv, pointVal)

ylow = (minX - pointVal) * slope_at_point + y_value_at_point
yhigh = (maxX - pointVal) * slope_at_point + y_value_at_point

    # now the tangent as a line plot
axes.plot([minX, maxX], [ylow, yhigh])

plt.show()
plt.close('all') # clean up after using pyplot

output 是：

Answer 1

很可能只是一個縮放問題，我們可以通過為獨立於原始數據縮放的梯度創建雙軸來解決。 為了安全起見，我們還向np.gradient提供 x 值，以防它們的間距不均勻。

import matplotlib.pyplot as plt
import numpy as np

fig, ax1 = plt.subplots()

def func(x, a=0, b=100, c=1, n=3.5):
    return a + (b/(1+(c/x)**n))

x_fit = np.linspace(0.1, 70, 100)
y_fit = func(x_fit, 1, 2, 15, 2.4)

tang = np.gradient(y_fit, x_fit)

ax1.plot(x_fit, y_fit, c="blue", label="data")
ax1.legend()
ax1.set_ylabel("data")

ax2 = ax1.twinx()
ax2.plot(x_fit, tang, c="red", label="gradient")
ax2.legend()
ax2.set_ylabel("gradient")

plt.show()

樣品 output：

如果我們將其繪制在同一張圖中，該圖：