%《數字信號處理教程——MATLAB釋義與實現》第九章例9.8.1程序hc981
% 稀疏收發天線陣列的方向圖生成
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
ht1=[1,1];hr1=[1,0,1,0,1,0,1,0,1,0,1,0,1,0,1];  % 給定稀疏收發陣元
he1=conv(ht1,hr1);                              % 有效陣列等于收發陣列的卷積
[Hr1,w]=freqz(hr1,1,'whole');                   % 求稀疏接收陣列的頻率特性（方向圖）
[He1,w]=freqz(he1,1,'whole');                   % 求有效收發陣列的頻率特性（方向圖）
plot(w-pi,fftshift(abs(Hr1))/16,'-.',w-pi,fftshift(abs(He1))/16) 
% 把頻率軸的中心移到原點畫圖，除以16，是為把零頻處的幅特性歸一化為1
legend('稀疏接收陣列的方向圖','收發雙向有效方向圖') %
set(gcf,'color','w')                            % 置圖形背景色為白