fixed spectograms

This commit is contained in:
ti_mo
2026-05-20 10:16:46 +02:00
parent 1561016689
commit 2fa2bf7f9a
2 changed files with 11 additions and 11 deletions
+1 -1
View File
@@ -1,4 +1,4 @@
calculations_only = true;
calculations_only = false;
task1;
task2;
+10 -10
View File
@@ -17,22 +17,22 @@ duration4 = length(trace4.iq) / sample_rate;
% 2.
% get all four plot values
[~,F1,T1,P1] = spectrogram(trace1.iq, window, overlap, fft_precision, sample_rate, "centered", "power");
[~,F2,T2,P2] = spectrogram(trace2.iq, window, overlap, fft_precision, sample_rate, "centered", "power");
[~,F3,T3,P3] = spectrogram(trace3.iq, window, overlap, fft_precision, sample_rate, "centered", "power");
[~,F4,T4,P4] = spectrogram(trace4.iq, window, overlap, fft_precision, sample_rate, "centered", "power");
[S1,F1,T1] = spectrogram(trace1.iq, window, overlap, fft_precision, sample_rate, "centered");
[S2,F2,T2] = spectrogram(trace2.iq, window, overlap, fft_precision, sample_rate, "centered");
[S3,F3,T3] = spectrogram(trace3.iq, window, overlap, fft_precision, sample_rate, "centered");
[S4,F4,T4] = spectrogram(trace4.iq, window, overlap, fft_precision, sample_rate, "centered");
% begin calibrated color scale
P1_dB = 10 * log10(P1);
P2_dB = 10 * log10(P2);
P3_dB = 10 * log10(P3);
P4_dB = 10 * log10(P4);
P1_dB = 20 * log10(abs(S1).^2);
P2_dB = 20 * log10(abs(S2).^2);
P3_dB = 20 * log10(abs(S3).^2);
P4_dB = 20 * log10(abs(S4).^2);
all_P_dB = [P1_dB(:); P2_dB(:); P3_dB(:); P4_dB(:)];
% -50 to +50 so its like the example fig
p_dB_min = min(all_P_dB) - 50;
p_dB_max = max(all_P_dB) + 50;
p_dB_min = min(all_P_dB);
p_dB_max = max(all_P_dB);
% end power calibrated color scale
if ~calculations_only