### Python 中频谱分析与可视化库推荐 以下是一些适用于频谱分析和可视化的 Python 库，涵盖从基础计算到高级交互式可视化的功能： #### 1. NumPy 和 SciPy NumPy 提供了数组操作的基础支持，而 SciPy 则扩展了这些功能以支持科学计算和信号处理[^1]。SciPy 的 `scipy.fftpack` 模块可以用于快速傅里叶变换（FFT），从而实现频谱分析。 ```python import numpy as np from scipy.fftpack import fft # 示例：生成正弦波并进行 FFT 分析 fs = 1000 # 采样频率 t = np.linspace(0, 1, fs, endpoint=False) # 时间向量 f = 5 # 正弦波频率 x = np.sin(2 * np.pi * f * t) # 正弦波 X = fft(x) # 快速傅里叶变换 freqs = np.fft.fftfreq(len(t), d=1/fs) # 频率轴 ``` #### 2. Matplotlib Matplotlib 是一个强大的绘图库，能够绘制静态、动态和交互式的频谱图。结合 NumPy 和 SciPy，可以轻松实现频谱的可视化[^2]。 ```python import matplotlib.pyplot as plt # 可视化频谱 plt.plot(freqs[:len(freqs)//2], np.abs(X[:len(X)//2])) plt.xlabel('Frequency (Hz)') plt.ylabel('Amplitude') plt.title('Spectrum Analysis') plt.grid() plt.show() ``` #### 3. librosa librosa 是一个专注于音频和音乐信号处理的库，提供了丰富的频谱分析工具，包括短时傅里叶变换（STFT）、梅尔频谱图等[^3]。 ```python import librosa import librosa.display # 加载音频文件 y, sr = librosa.load(librosa.ex('trumpet')) # 计算短时傅里叶变换（STFT） D = librosa.stft(y) # 可视化频谱图 librosa.display.specshow(librosa.amplitude_to_db(np.abs(D), ref=np.max), y_axis='log', x_axis='time', cmap='viridis') plt.colorbar(format='%+2.0f dB') plt.title('Spectrogram') plt.show() ``` #### 4. Bokeh 和 Plotly Bokeh 和 Plotly 是两个交互式可视化库，适合需要动态交互的频谱分析场景。Plotly 提供了简单易用的 API 来创建交互式图表[^4]。 ```python import plotly.graph_objects as go # 创建交互式频谱图 fig = go.Figure(data=[go.Scatter(x=freqs[:len(freqs)//2], y=np.abs(X[:len(X)//2]))]) fig.update_layout(title='Interactive Spectrum Analysis', xaxis_title='Frequency (Hz)', yaxis_title='Amplitude') fig.show() ``` #### 5. Scikit-learn Scikit-learn 是一个机器学习库，虽然其主要用途不是频谱分析，但可以通过特征提取和降维技术对频谱数据进行进一步分析和分类[^5]。 ```python from sklearn.decomposition import PCA # 对频谱数据进行主成分分析 pca = PCA(n_components=2) principalComponents = pca.fit_transform(np.abs(D).T) # 可视化降维后的结果 plt.scatter(principalComponents[:, 0], principalComponents[:, 1]) plt.xlabel('Principal Component 1') plt.ylabel('Principal Component 2') plt.title('PCA of Spectrogram Features') plt.show() ``` ---