垂直切换的仿真（终）(2)

text(0,1,'\\leftarrow Wimax','FontSize',10)

title('the state of mobile station in LTE and Wimax network')

%############################END THE PROGRAM ##############################

2. is_in_LTE_coverage.m

function is_LTE=is_in_LTE_coverage(MS_coordinate) Tcide whether the mobile station is in the LTE network global LTE_BS_coordinate;

d=norm(MS_coordinate-LTE_BS_coordinate); if d<=1000&d>=35 is_LTE=true; else

is_LTE=false; end

3. is_in_Wimax_coverage.m

function is_Wimax=is_in_Wimax_coverage(MS_coordinate) Tcide whether the mobile station is in the Wimax network global Wimax_BS_coordinate;

d=norm(MS_coordinate-Wimax_BS_coordinate); if d<=300&d>=36 is_Wimax=true; else

is_Wimax=false; end

4. calculate_RSS.m

function RSS=calculate_RSS(MS_coordinate,MS_speed,s,record_time,Tc) êlculate the received signal strength switch s

case 'LTE'

slow_fading=calculate_slow_fading(MS_coordinate,'LTE',... record_time,Tc);

fast_fading=calculate_fast_fading(MS_speed,record_time,'LTE',Tc); RSS=33-slow_fading+fast_fading; otherwise

slow_fading=calculate_slow_fading(MS_coordinate,'Wimax',... record_time,Tc);

fast_fading=calculate_fast_fading(MS_speed,record_time,'Wimax',Tc); RSS=23-slow_fading+fast_fading; end

5. calculate_fast_fading.m

function fast_fading=calculate_fast_fading(MS_speed,t,s,Tc)

c=3.0*10^8; %the speed of electromagnetic wave fc1=2000*10^6; %the frequency of the carrier for LTE fc2=2500*10^6; %the frequency of the carrier for Wimax fm1=MS_speed*fc1/c; %the maximum of doppler frequency shift %in LTE network

fm2=MS_speed*fc2/c; %the maximum of doppler frequency shift N0=2; switch s

case 'LTE' fc=fc1;

fm=fm1; otherwise fc=fc2; fm=fm2; end N=50; for i=1:N

R_t(i)=jakes(N0,fm,t+i/N*Tc,fc); end

fast_fading=10*log10(mean(R_t.^2));

6. calculate_slow_fading.m

function slow_fading= calculate_slow_fading(MS_coordinate,s,record_time,Tc) êlcualte the large scale fading

global fc_LTE fc_Wimax LTE_BS_coordinate Wimax_BS_coordinate shadow_fading=calculate_shadow_fading(s,record_time,Tc); switch s

case 'LTE' f=fc_LTE;

d=norm(MS_coordinate-LTE_BS_coordinate); otherwise

f=fc_Wimax;

d=norm(MS_coordinate-Wimax_BS_coordinate); end

propagation_loss=32.44+20*log10(d/1000)+20*log10(f/10^6)+12+3+24; slow_fading=propagation_loss+shadow_fading;

7. calculate_shadow_fading.m

function shadow_fading=calculate_shadow_fading(s,record_time,Tc) global large_scale_time_interval large_scale_LTE large_scale_wimax ... large_scale_time_start large_scale_time_end N if record_time>large_scale_time_end

large_scale_time_start=large_scale_time_start+... large_scale_time_interval/2;

large_scale_time_end=large_scale_time_end+large_scale_time_interval/2; large_scale_LTE(1:5)=large_scale_LTE(6:10)

large_scale_LTE(6:10)=ShadowGenerate(3,0.5,N/2); large_scale_wimax(6:10)=ShadowGenerate(3,0.5,N/2); end switch s

case 'LTE'

shadow_fading=interp1(linspace(large_scale_time_start,... large_scale_time_end,N),large_scale_LTE,record_time,'cubic'); otherwise

shadow_fading=interp1(linspace(large_scale_time_start,...

large_scale_time_end,N),large_scale_wimax,record_time,'cubic'); end

8. ShadowGenerate.m

function corr = ShadowGenerate(standard_d,relation_coeff,N) % e.g. SFstd = 8dB lambda = 0.5;

forwcoeff = standard_d * sqrt(1-relation_coeff.^2); backwcoeff = [ones(1,N) -relation_coeff]; uncorr = randn(1,N);

corr = filter(forwcoeff,backwcoeff,uncorr);