@@ -81,6 +81,7 @@ audio_stft = zaf.stft(audio_signal,window_function,step_length);
81
81
audio_spectrogram = abs(audio_stft(2:window_length/2+1,:));
82
82
83
83
% Display the spectrogram in dB, seconds, and Hz
84
+ number_samples = length(audio_signal);
84
85
xtick_step = 1;
85
86
ytick_step = 1000;
86
87
figure
@@ -237,9 +238,10 @@ mel_filterbank = zaf.melfilterbank(sampling_frequency,window_length,number_mels)
237
238
mel_spectrogram = zaf.melspectrogram(audio_signal,window_function,step_length,mel_filterbank);
238
239
239
240
% Display the mel spectrogram in in dB, seconds, and Hz
241
+ number_samples = length(audio_signal);
240
242
xtick_step = 1;
241
243
figure
242
- zaf.melspecshow(mel_spectrogram, length(audio_signal) , sampling_frequency, window_length, xtick_step)
244
+ zaf.melspecshow(mel_spectrogram, number_samples , sampling_frequency, window_length, xtick_step)
243
245
title('Mel spectrogram (dB)')
244
246
```
245
247
@@ -339,8 +341,8 @@ imagesc(abs(cqt_kernel))
339
341
axis xy
340
342
colormap(jet)
341
343
title('Magnitude CQT kernel')
342
- xlabel('FFT length ')
343
- ylabel('CQT frequency ')
344
+ xlabel('FFT index ')
345
+ ylabel('CQT index ')
344
346
```
345
347
346
348
<img src =" images/cqtkernel.png " width =" 1000 " >
@@ -375,7 +377,7 @@ minimum_frequency = 55;
375
377
maximum_frequency = 3520;
376
378
cqt_kernel = zaf.cqtkernel(sampling_frequency,octave_resolution,minimum_frequency,maximum_frequency);
377
379
378
- % Compute the (magnitude) CQT spectrogram using the kernel
380
+ % Compute the CQT spectrogram using the kernel
379
381
time_resolution = 25;
380
382
cqt_spectrogram = zaf.cqtspectrogram(audio_signal,sampling_frequency,time_resolution,cqt_kernel);
381
383
@@ -577,10 +579,11 @@ window_function = sqrt([window_function2; window_function2(window_length/2:-1:1)
577
579
audio_mdct = zaf.mdct(audio_signal,window_function);
578
580
579
581
% Display the MDCT in dB, seconds, and Hz
582
+ number_samples = length(audio_signal);
580
583
xtick_step = 1;
581
584
ytick_step = 1000;
582
585
figure
583
- zaf.specshow(abs(audio_mdct),length(audio_signal) ,sampling_frequency,xtick_step,ytick_step)
586
+ zaf.specshow(abs(audio_mdct),number_samples ,sampling_frequency,xtick_step,ytick_step)
584
587
title('MDCT (dB)')
585
588
```
586
589
0 commit comments