• Create BookmarkCreate Bookmark
  • Create Note or TagCreate Note or Tag
  • PrintPrint
Share this Page URL
Help

Chapter 5. Radio Wave > MATLAB/Simulink Models

5.7. MATLAB/Simulink Models

5.7.1. Reflection Coefficient

%clear
m1=1;
m2=1;
e1=1;
e2=9;
n1=sqrt(m1/e1);
n2=sqrt(m2/e2);
thtaid=[0:90];
thtai=thtaid/180*pi;
thtar=thtai;
thtat=acos(cos(thtai).*sqrt(m1*e1/m2/e2));
ii=1;
for i=1:91
   if imag(thtat(i))~=0
      ii=ii+1;
   end
end
whitebg('white')
figure(1)
plot(thtaid,thtat/pi*180,'k')
grid
xlabel('Incident Angle (Degree)')
ylabel('Transmit Angle (Degree)')
R1=ones(91,1);
R2=ones(91,1);
T1=zeros(91,1);
T2=zeros(91,1);
R1(ii:91)=((-n1.*sin(thtai(ii:91))+n2.*sin(thtat(ii:91)))./(n1.*sin(thtai(ii:91))+n2.*sin(thtat(ii:91)))).^2;
R2(ii:91)=((n2.*sin(thtai(ii:91))-n1.*sin(thtat(ii:91)))./(n2.*sin(thtai(ii:91))+n1.*sin(thtat(ii:91)))).^2;
T1(ii:91)=(4*n1*n2.*sin(thtai(ii:91)).*sin(thtat(ii:91)))./(n1.*sin(thtai(ii:91))+n2.*sin(thtat(ii:91))).^2;
T2(ii:91)=(4*n1*n2.*sin(thtai(ii:91)).*sin(thtat(ii:91)))./(n2.*sin(thtai(ii:91))+n1.*sin(thtat(ii:91))).^2;
whitebg('white')
figure(2)
plot(thtaid,R1,'k',thtaid,T1,'k');
grid
xlabel('Angle (Degree)');
ylabel('Coefficient');
gtext('Parallel Reflection Coefficient')
gtext('Parallel Transmit Coefficient')
whitebg('white')
figure(3)
plot(thtaid,R2,'k',thtaid,T2,'k');
grid
xlabel('Angle (Degree)');
ylabel('Coefficient');
gtext('Perpendicular Reflection Coefficient')
gtext('Perpendicular Transmit Coefficient')


					  

5.7.2. Signal Strength

d=[3:50];
l1=0.32*39.37/12;
l2=0.12*39.37/12;
l3=0.052*39.37/12;
lu=0.5;
ll=3;
data1=[3.5 6; 3.5 14; 3.9 12; 5 11; 5 19; 5 24; 7 15; 7.3 13; 7.3 17; 
7.3 24; 7.3 27; 7.5 20; 7.5 18; 8.2 21; 8.2 27; 9.0 18; 9.0 28; 9.0 30; 
10 20; 10 26; 12 32; 12 35];
data2=[2.3 6; 2.3 12; 3.1 7; 3.1 11; 3.5 13; 4 12; 4 19; 4 24; 5 8;
5 12; 5 15; 5 21; 6 20; 7 16; 7 29; 8 22; 8 32; 9 25; 9 28; 10 24;
10 30; 13 23; 13 26; 13 27; 13 33; 13 39];
whitebg('white')
figure(1)
provptu=10*log10(l1^2/16/pi^2/lu./d.^2);
provpt=10*log10(l1^2/16/pi^2./d.^2);
provptl=10*log10(l1^2/16/pi^2/ll./d.^2);
plot(d,provptu,'k',d,provpt,'k',d,provptl,'k');
grid
xlabel('Distance (ft)')
ylabel('Attenuation (dB)')
whitebg('white')
figure(2)
provptu=10*log10(l2^2/16/pi^2/lu./d.^2);
provpt=10*log10(l2^2/16/pi^2./d.^2);
provptl=10*log10(l2^2/16/pi^2/ll./d.^2);
plot(d,provptu,'k',d,provpt,'k',d,provptl,'k',data1(:,1)*39.37/12,-data1(:,2)-40,'kx');
grid
xlabel('Distance (ft)')
ylabel('Attenuation (dB)')
whitebg('white')
figure(3)
provptu=10*log10(l3^2/16/pi^2/lu./d.^2);
provpt=10*log10(l3^2/16/pi^2./d.^2);
provptl=10*log10(l3^2/16/pi^2/ll./d.^2);
plot(d,provptu,'k',d,provpt,'k',d,provptl,'k',data2(:,1)*39.37/12,-data2(:,2)-46,'kx');
grid
xlabel('Distance (ft)')
ylabel('Attenuation (dB)')


					  


PREVIEW

                                                                          

Not a subscriber?

Start A Free Trial


  
  • Creative Edge
  • Create BookmarkCreate Bookmark
  • Create Note or TagCreate Note or Tag
  • PrintPrint