www.pudn.com > Correlation.Multiple.Cluster.rar > plot_laplacian.m

function plot_laplacian(number_clusters, power, AOA_deg, sigma_deg, ... 
		      delta_phi_deg, plot_number, colour) 
% plot_laplacian(number_clusters, power, AOA_deg, sigma_deg, 
%                delta_phi_deg, plot_number, colour) 
% 
% Plots the multi-cluster Laplacian PAS 
% 
% 
% STANDARD DISCLAIMER 
% 
% CSys is furnishing this item "as is". CSys does not provide any 
% warranty of the item whatsoever, whether express, implied, or 
% statutory, including, but not limited to, any warranty of 
% merchantability or fitness for a particular purpose or any 
% warranty that the contents of the item will be error-free. 
% 
% In no respect shall CSys incur any liability for any damages, 
% including, but limited to, direct, indirect, special, or 
% consequential damages arising out of, resulting from, or any way 
% connected to the use of the item, whether or not based upon 
% warranty, contract, tort, or otherwise; whether or not injury was 
% sustained by persons or property or otherwise; and whether or not 
% loss was sustained from, or arose out of, the results of, the 
% item, or any services that may be provided by CSys. 
% 
% (c) Laurent Schumacher, AAU-TKN/IES/KOM/CPK/CSys - September 2002 
 
phi_deg   = -180:1:180; 
maximum   = 0; 
 
phi_rad   = phi_deg*pi/180; 
sigma_rad = sigma_deg*pi/180; 
 
for (k = 1:1:number_clusters) 
  mask = .5.*(sign(delta_phi_deg(k)-phi_deg+1e-6)+... 
              sign(delta_phi_deg(k)+phi_deg+1e-6)); 
  PAS = (power(k)/(sqrt(2)*sigma_rad(k))).*exp(((-1).*sqrt(2).*abs(phi_rad))./sigma_rad(k)); 
  if (AOA_deg(k) > 0) 
    limit = find(phi_deg==(180-AOA_deg(k))); 
    mask  = [mask(limit+1:size(mask,2)),mask(1:limit)]; 
    PAS = [PAS(limit+1:size(PAS,2)),PAS(1:limit)]; 
  else 
    limit = find(phi_deg==(-180-AOA_deg(k))); 
    mask  = [mask(limit:size(mask,2)),mask(1:limit-1)]; 
    PAS = [PAS(limit:size(PAS,2)),PAS(1:limit-1)]; 
  end; 
  maximum = max(maximum,max(PAS)); 
  subplot(2,1,plot_number),plot(phi_deg, mask.*PAS, colour); 
  hold on; 
end; 
axis([-180 180 0 1.1*maximum]); 
grid; 
ylabel('Power (Linear)'); 
if (plot_number == 1) 
  title('Power Azimuth Spectrum'); 
else 
  xlabel('\phi [degrees]'); 
end;