MATLAB Listing

%Figure 3.52; ‘Mobile satellite communication; Principles and trends’ by M.Richharia;

%Plot of supportable bit rate against satellite EIRP for various types of terminals;

%Matlab source code written by AR/MR;

rb=[2.4 24 64 128 256 512 1200 2400];%kbits/sec

rb=0:1:2400;

rb=rb*1000;%bits/sec

%C/No=Eb/No+10log(rb) : Eb/No=4.5 …(1)

EbNo=4.5;

CNo=EbNo+10*log10(rb);

rb=rb/1000;%convert rb to kbit/sec

%Pathloss=20(log(4pi*41680*1000/0.15) …(2)

Pathloss=20*log10(4*pi*40000*1000/0.2);

%EIRP=C/No+Path loss-G/T-228.6 …(3)

GT=0;

EIRP1=CNo+Pathloss-GT-228.6;

GT=-5;

EIRP2=CNo+Pathloss-GT-228.6;

GT=-15;

EIRP3=CNo+Pathloss-GT-228.6;

GT=-25;

EIRP4=CNo+Pathloss-GT-228.6;

%Next with varying fade margin

%Add fademargin to EIRP estimated in(eq 3)

GT=0;

fademargin=2;

EIRP5=CNo+Pathloss-GT-228.6 + fademargin;

GT=-5;

fademargin=3;

EIRP6=CNo+Pathloss-GT-228.6+ fademargin;

GT=-15;

fademargin=6;

EIRP7=CNo+Pathloss-GT-228.6+ fademargin;

GT=-25;

fademargin=10;

EIRP8=CNo+Pathloss-GT-228.6+ fademargin;

% Plot: figure(‘color’,’white’)

hold;

plot(rb,EIRP1,’k’,rb,EIRP2,’k’,rb,EIRP3,’k’,rb,EIRP4,’k’,rb,EIRP5,’k:’,rb,EIRP6,’k:’,rb,EIRP7,’k:’,rb,EIRP8,’k:’);

xlabel(‘Bit rate (Kbit/sec)’);

ylabel(‘Satellite EIRP (dBW)’);

axis([0 70 -10 50]);

figure(‘color’,’white’)

hold;

plot(rb,EIRP1,’k’,rb,EIRP2,’k’,rb,EIRP3,’k’,rb,EIRP4,’k’,rb,EIRP5,’k:’,rb,EIRP6,’k:’,rb,EIRP7,’k:’,rb,EIRP8,’k:’);

%legend (‘G/T=0′,’G/T=-5′,’G/T=-15′,’G/T=-25’);

xlabel(‘Bit rate (Kbit/sec)’);

ylabel(‘Satellite EIRP (dBW)’);

axis([100 2000 0 70]);

%grid on;

%zoom;