python 光谱曲线的二阶微分

光谱曲线的二阶微分可以通过Python中的SciPy库来实现。具体步骤如下：

1. 导入需要的库和数据

import numpy as np
import matplotlib.pyplot as plt
from scipy.signal import savgol_filter
from scipy.misc import derivative

# 假设光谱数据存储在x和y中
x = np.array([1, 2, 3, 4, 5])
y = np.array([2, 4, 5, 4, 2])

2. 平滑处理

# 使用Savitzky-Golay滤波器平滑数据
y_smooth = savgol_filter(y, window_length=5, polyorder=2)

3. 计算二阶导数

# 计算二阶导数
y_second_derivative = derivative(savgol_filter, x, dx=1, n=2, args=(y,))

4. 绘制结果

# 绘制原始数据和平滑后的数据
plt.plot(x, y, 'o', label='Original data')
plt.plot(x, y_smooth, label='Smoothed data')

# 绘制二阶导数
plt.plot(x, y_second_derivative, label='Second derivative')

plt.legend()
plt.show()