Python中的 matplotlib.pyplot.angle_spectrum()
Matplotlib是Python中的一个库,它是 NumPy 库的数值数学扩展。 Pyplot是Matplotlib模块的基于状态的接口,它提供了一个类似 MATLAB 的接口。
matplotlib.pyplot.acorr()函数:
matplotlib 库的 pyplot 模块中的angle_spectrum()函数用于绘制角度谱。通常,它计算序列的角谱并绘制。
Syntax: angle_spectrum(x, Fs=2, Fc=0, window=mlab.window_hanning, pad_to=None, sides=’default’, **kwargs)
Parameters: This method accept the following parameters that are described below:
- x: This parameter is a sequence of data.
- Fs : This parameter is a scalar. Its default value is 2.
- window: This parameter take a data segment as an argument and return the windowed version of the segment. Its default value is window_hanning()
- sides: This parameter specifies which sides of the spectrum to return. This can have following values : ‘default’, ‘onesided’ and ‘twosided’.
- pad_to : This parameter contains the integer value to which the data segment is padded.
- Fc: This parameter is also contains the integer value to offsets the x extents of the plot to reflect the frequency range. Its default value is 0
Returns: This returns the following:
- spectrum :This returns the angle spectrum in radians.
- freqs :This returns the frequencies corresponding to the elements in spectrum.
- line : This returns the line created by this function.
The resultant is (spectrum, freqs, line)
下面的示例说明了 matplotlib.pyplot 中的 matplotlib.pyplot.angle_spectrum()函数:
示例 #1:
# Implementation of matplotlib.pyplot.angle_spectrum()
# function
import matplotlib.pyplot as plt
import numpy as np
np.random.seed(10**5)
dt = 0.0001
Fs = 1 / dt
geeks = np.array([24.40, 110.25, 20.05,
22.00, 61.90, 7.80,
15.00, 22.80, 34.90,
57.30])
nse = np.random.randn(len(geeks))
r = np.exp(-geeks / 0.05)
s = 0.1 * np.sin(2 * np.pi * geeks) + nse
# plot angle_spectrum
plt.angle_spectrum(s, Fs = Fs)
# Display the simple spectrum and
# angle_spectrum plot
plt.show()
输出:
示例 #2:
# Implementation of matplotlib.pyplot.angle_spectrum()
# function
import matplotlib.pyplot as plt
import numpy as np
np.random.seed(0)
dt = 0.01
Fs = 1 / dt
t = np.arange(0, 10, dt)
nse = np.random.randn(len(t))
r = np.exp(-t / 0.05)
cnse = np.convolve(nse, r)*dt
cnse = cnse[:len(t)]
s = 0.1 * np.sin(2 * np.pi * t) + cnse
# plot simple spectrum
plt.subplot(2, 1, 1)
plt.plot(t, s)
# plot angle_spectrum
plt.subplot(2, 1, 2)
plt.angle_spectrum(s, Fs = Fs)
# Display the simple spectrum and
# angle_spectrum plot
plt.show()
输出: