clear all;
close all;

max_speed_mph = 150;
constellation = 16;    % 16 = 16-QAM, 64 = 64-QAM
f = 120;				  % frequency in MHz
Tg = 300e-6;   % guard interval = 300 us
percent_subcarrier_spacing = 7;   % 9% based on DVB-H at 600 MHz for 100 mph

v = max_speed_mph * 1609 / 3600;  %v in m/s
current_modeA_Ts = 26.66e-3;
lambda = 3e8 / (f * 1e6)
fd_max = v / lambda
subcarrier_spacing = fd_max / (percent_subcarrier_spacing / 100)
Tu = 1/subcarrier_spacing
relative_symbol_rate = current_modeA_Ts / (Tg + Tu);
current_modeA_bit_rate = [38.2e3; 55.0e3];    % 16-QAM, 64-QAM for mode A, spectrum occupancy = 20kHz
current_Nsub = 461;

for i=1:2 
	for b=30:10:200
   	Nsub(b/10,1) = round(b * 1e3 / subcarrier_spacing);
   	bit_rate(b/10,i) = (Nsub(b/10,1) / current_Nsub) * current_modeA_bit_rate(i,1) * relative_symbol_rate;
   end
end

Nsub
bit_rate

plot(bit_rate)