db10 = @(x) 10 * log10(x);
db20 = @(x) 20 * log10(x);
inv_db10 = @(x) 10.^(x/10);
inv_db20 = @(x) 10.^(x/20);
Lnorm = -220; % [dBc/Hz], Normalized Phase Noise Floor
L1_f = -114; % [dBc/Hz], Normalized 1/f Noise, Measured at 10kHz offset, Normalized to 1GHz
frf_meas = 1000e6; %[Hz]
foffset_meas = 10e3; %[Hz]
frf = 2.1e9; %[Hz]
fpfd = 0.2e6; %[Hz]
% plot the results (except at DC)
semilogx(f_Hz(2:end), g_dBc1Hz(2:end), 'r', 'Linewidth', 2)
Lout_1 = ( Lnorm + db10(fpfd) + db20(frf/fpfd) ) * ones( size( f_Hz(2:end) ) );
semilogx(f_Hz(2:end), Lout_1, 'g', 'Linewidth', 2)
Lout_2 = L1_f + db10(foffset_meas./f_Hz(2:end)) + db20(frf/frf_meas);
semilogx(f_Hz(2:end), Lout_2, 'c', 'Linewidth', 2)
Lout_3 = db10( inv_db10(Lout_1) + inv_db10(Lout_2) );
semilogx(f_Hz(2:end), Lout_3, 'y', 'Linewidth', 2)