array.m

用matlab编写的电磁学方面的函数和图像

% array.m - gain computation for 1D equally-spaced isotropic array
%
% Usage: [g, phi] = array(d, a, Nph);
%
% d   = array spacing in units of lambda
% a   = row vector of array weights 
% Nph = number of azimuthal angles in [0,pi] (actually, Nph+1 angles)
%
% g   = row vector of gain values evaluated at phi
% phi = row vector of (Nph + 1) equally-spaced angles over [0, pi] (in radians)
%
% notes: computes g(phi) = |A(psi)|^2, where A(psi) = sum_n a(n)z^n, 
%        with z = e^(j*psi) and psi = 2*pi*d*cos(phi),
%
%        normalizes g to unity maximum,
%        [g,phi] can be passed into gain plotting functions DBZ,ABZ,
%        e.g., dbz(phi,g), abz(phi,g)
%
%        to compute the gain of a scanned array use SCAN or STEER first, e.g.
%        [g, phi] = array(d, scan(a,psi0), Nph);
%        [g, phi] = array(d, steer(d,a,phi0), Nph);
%
%        uses the I2SP function DTFT
%        
% see also UNIFORM, BINOMIAL, TAYLOR, DOLPH, SECTOR, MULTIBEAM

% S. J. Orfanidis - 1997 - www.ece.rutgers.edu/~orfanidi/ewa

function [g, phi] = array(d, a, Nph)

if nargin==0, help array; return; end        

phi = (0 : Nph) * pi / Nph;     % equally-spaced over [0,pi]
        
psi = 2 * pi * d * cos(phi);

A = dtft(a, -psi);              % array factor, note dsft(a,psi)=dtft(a,-psi)
g = abs(A).^2;                  % power gain                                     
g = g/max(g);                   % normalized to unity maximum