hc981.m

matlab数字信号处理代码

%《数字信号处理教程——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')                            % 置图形背景色为白