MATLAB Listing

%Figure 3.44 (a) and (b); ‘Mobile satellite communication; Principles and trends’ by M.Richharia (MR);

%Data source:ITU-R report; see MR’s book for reference;

%Estimated topospheric attenuation due to Oxygen and water for an elevation angle of 10 deg, single transit, at various altitudes;

%Matlab source code written by AR/MR:Year 2000;

altitude=[0 1 2 3 4 5 6];%in km

altitude1=[0 1 2 3 ]

altitude2=[0 1 2 3 ]

altitude3=[0 1 2 3 4]

oxygenAttenuation15=[0.19 0.16 0.14 0.12 0.1 0.08 0.07];%Attenuation (dB) for 1.5GHz;

oxygenAttenuation20=[0.36 0.3 0.26 0.22 0.18 0.16 0.13];%Attenuation (dB) for 20GHz;

oxygenAttenuation40=[1.52 1.29 1.09 0.92 0.78 0.66 0.56];%Attenuation (dB) for 40GHz;

waterAttenuation15=[0.0016 0.001 0.0006 0.0004];%Attenuation (dB) for 1.5GHz;

waterAttenuation20=[1.35 0.92 0.62 0.42];%Attenuation (dB) for 20GHz;

waterAttenuation40=[1.17 0.75 0.47 0.3];%Attenuation (dB) for 40GHz;

cloudAttenuation15=[0.02 0.02 0.01 0];%Attenuation (dB) for 1.5GHz;

cloudAttenuation20=[3.48 3.48 1.74 0];%Attenuation (dB) for 20GHz;

cloudAttenuation40=[13.9 13.9 7 0];%Attenuation (dB) for 40GHz;

rainAttenuation15=[0.009 0.008 0.006 0.004 0];%Attenuation (dB) for 1.5GHz;

rainAttenuation20=[16.3 14.2 11.3 6.9 0];%Attenuation (dB) for 20GHz;

rainAttenuation40=[43 37.3 29.6 18.3 0];%Attenuation (dB) for 40GHz;

% Plot figure;

figure(‘color’,’white’);

hold;

plot(altitude,oxygenAttenuation15,’k’,altitude1, waterAttenuation15,’k*–‘,altitude,oxygenAttenuation20,’k’,altitude,oxygenAttenuation40,’k’,altitude1,waterAttenuation20,’k*–‘,altitude1,waterAttenuation40,’k*–‘);

legend(‘Oxygen attenuation’, ‘Water attenuation’)

xlabel(‘Altitude (km)’);

ylabel(‘Attenuation (dB)’);

grid on;

figure(‘color’,’white’);

hold;

plot(altitude2,cloudAttenuation15,’k’,altitude3, rainAttenuation15,’k*–‘,altitude2,cloudAttenuation20,’k’,altitude2,cloudAttenuation40,’k’,altitude3,rainAttenuation20,’k*–‘,altitude3,rainAttenuation40,’k*–‘);

legend(‘Cloud attenuation’, ‘Rain attenuation’)

xlabel(‘Altitude (km)’);

ylabel(‘Attenuation (dB)’);

grid on;