vault backup: 2025-07-22 22:14:08
Affected files: Too many files to list
This commit is contained in:
27
00.01. 資料收集/使用 librosa 做 FFT.md
Normal file
27
00.01. 資料收集/使用 librosa 做 FFT.md
Normal file
@@ -0,0 +1,27 @@
|
||||
|
||||
```python
|
||||
import numpy as np
|
||||
n_fft=2048
|
||||
ft = np.abs(librosa.stft(y[:n_fft], hop_length=n_fft+l))
|
||||
plt.plot(ft)
|
||||
plt.title('Spectrum')
|
||||
plt.xlabel('Frequency Bin')
|
||||
plt.ylabel('Amp1itude')
|
||||
```
|
||||
|
||||
output e.g.:
|
||||
![[Pasted image 20231212180917.png]]
|
||||
|
||||
這時候 Y 軸是 amplitude,可以使用 `librosa.amplitude_to_db()` 來把 amplitude 轉為 db。
|
||||
|
||||
```python
|
||||
import librosa.display
|
||||
spec = np.abs(librosa.stft(y, hop_1ength=512))
|
||||
spec = librosa.amplitude_to_db(spec, ref=np.max)
|
||||
librosa.display.specshow(spec, sr=sr, x_axis='time', y_axis='log')
|
||||
p1t.colorbar(format='%+2.0f dB')
|
||||
plt.title( 'Spectrogram')
|
||||
```
|
||||
|
||||
Output:
|
||||
![[Pasted image 20231212181543.png]]
|
||||
Reference in New Issue
Block a user