rpp.m

经典书籍MIMO Signals And Systems的源代码

% Program rpp.m calculates the radiation power
% pattern (rpp) in dB (down of maximum) of
% arbitrary length thin wire dipole
clear
N = 2000;
dlam = input('Enter ratio antenna_length/lambda = ');
dBrange = input('Enter dB range of polar plot = ');

theta = [0:2*pi/N:(2*pi)];
rp = zeros(1,length(theta)); 
a = pi*dlam;
for k=1:length(theta)
      if abs(sin(theta(k))) < 10^(-dBrange/20)
      rp(k) = 10^(-dBrange/20);
   else
   rp(k)= abs((cos(a*cos(theta(k)))-cos(a))/sin(theta(k)));
      end   
end
rp = rp/max(rp);

%plot(theta, rp)
logpolar(theta, rp,'k', dBrange)