MATLAB Listing

%Figure 5.2;‘Mobile satellite communications – Principles and Trends’ by M.Richharia (MR);

%Polarisation loss versus Mobile and satellite antenna axial ratios; 

%Equation (4.17) Ohmori, Wakana and Kawase (1998)-See MR’s book for reference;

AmdB=linspace(0,10,10);%beamwidth

figure(‘color’,’white’);

hold;

Am = 10.^ (AmdB/20);

AsdB=0;

As = 10^ (AsdB/20);

ploss0 = (Am.^2 + As^2 + 2*Am*As)./((1+Am.^2)*(1+As ^2));

AsdB=1;

As = 10^ (AsdB/20);

ploss1 = (Am.^2 + As^2 + 2*Am*As)./((1+Am.^2)*(1+As ^2));

AsdB=2;

As = 10^ (AsdB/20);

ploss2 = (Am.^2 + As^2 + 2*Am*As)./((1+Am.^2)*(1+As ^2));

AsdB=3;

As = 10^ (AsdB/20);

ploss3 = (Am.^2 + As^2 + 2*Am*As)./((1+Am.^2)*(1+As ^2));

AsdB=1;

As = 10^ (AsdB/20);

%ploss5 = ((1+Am*As).^2)./((1+Am.^2)*(1+As ^2));

%ploss5 = (Am.^2 + As^2 – 2*Am*As)./((1+Am.^2)*(1+As ^2));

AsdB=2;

As = 10^ (AsdB/20);

%ploss6 = ((1+Am*As).^2)./((1+Am.^2)*(1+As ^2));

%ploss6 = (Am.^2 + As^2 – 2*Am*As)./((1+Am.^2)*(1+As ^2));

AsdB=3;

As = 10^ (AsdB/20);

%ploss7 = ((1+Am*As).^2)./((1+Am.^2)*(1+As ^2));

%ploss7 = (Am.^2 + As^2 – 2*Am*As)./((1+Am.^2)*(1+As ^2));

plot(AmdB,10*log10(ploss0),’k’,AmdB,10*log10(ploss1),’k’,AmdB,10*log10(ploss2),’k’,AmdB,10*log10(ploss3),’k’);

%plot (AmdB,10*log10(ploss5),AmdB,10*log10(ploss6),AmdB,10*log10(ploss7));

%legend (‘AsdB=0′,’AsdBmax=1′,’AsdBmax=2′,’AsdBmax=3′,’AsdBmin=1′,’AsdBmin=2′,’AsdBmin=3’);

%legend (‘Axial ratio – satellite (dB)’);

axis([0 10 -inf inf]);

xlabel(‘Axial Ratio – Mobile (dB)’);

ylabel(‘Polarization Loss (dB)’);

grid on;

%zoom;