使用Scipy生成单位抛物信号遇AttributeError错误求助

错误原因 #

scipy.special是Scipy的独立子模块，不属于scipy.signal的下属模块，直接通过signal.special调用会触发AttributeError。

修正步骤与代码 #

1. 正确导入子模块 #

需要单独导入scipy.special，而非通过scipy.signal访问。

2. 生成单位抛物信号 #

根据需求选择以下两种方案：

方案一：生成标准单位抛物信号

如果需要的是形如 $y = t^2$ 且振幅归一化到1的抛物信号，代码如下：

from scipy import signal
import numpy as np
import matplotlib.pyplot as plt
import scipy.special as sp

plt.title('Unit Parabolic Signal')
plt.xlabel('Sample') 
plt.ylabel('Amplitude')

# 生成-5到5之间的100个采样点
t = np.linspace(-5, 5, 100)
# 生成归一化的单位抛物信号
parabolic_signal = (t ** 2) / np.max(t ** 2)

plt.stem(t, parabolic_signal, basefmt=" ")
plt.show()

方案二：使用抛物柱面函数pbdv

如果需求是基于抛物柱面函数生成信号，修正导入后调用即可：

from scipy import signal
import numpy as np
import matplotlib.pyplot as plt
import scipy.special as sp

plt.title('Unit Parabolic Signal (Using pbdv)')
plt.xlabel('Sample') 
plt.ylabel('Amplitude')

# 生成-5到5之间的50个采样点
t = np.linspace(-5, 5, 50)
# 调用pbdv函数，取第一个返回值（对应抛物柱面函数的主解）
parabolic_signal = sp.pbdv(0, t)[0]
# 归一化到单位振幅
parabolic_signal = parabolic_signal / np.max(np.abs(parabolic_signal))

plt.stem(t, parabolic_signal, basefmt=" ")
plt.show()

关键说明 #

• scipy.special.pbdv返回包含两个数组的元组，对应抛物柱面函数的两个解，通常取第一个元素（索引0）即可。
• 归一化步骤用于将信号振幅调整为单位值，符合“单位抛物信号”的要求。

内容的提问来源于stack exchange，提问作者ITS_HJ 24