传真(Fax)编程

开发平台：

Matlab

fig10_5.m：源码内容

							% Use this code to produce figure 10.5a and 10.5b based on equation 10.34
clear all
close all
eps = 0.00001;
k = 2*pi;
theta = -pi : pi / 10791 : pi;
var = sin(theta);
nelements = 8;
d = 1;         %  d = 1;
num = sin((nelements * k * d * 0.5) .* var);
if(abs(num) <= eps)
   num = eps;
end
den = sin((k* d * 0.5) .* var);
if(abs(den) <= eps)
   den = eps;
end
pattern = abs(num ./ den);
maxval = max(pattern);
pattern = pattern ./ maxval;
figure(1)
plot(var,pattern)
xlabel('sine angle - dimensionless')
ylabel('Array pattern')
grid
figure(2)
plot(var,20*log10(pattern))
axis ([-1 1 -60 0])
xlabel('sine angle - dimensionless')
ylabel('Power pattern [dB]')
grid;
figure(3)
theta = theta +pi/2;
polar(theta,pattern)
title ('Array pattern')

						