資源簡介
Smart Antennas for Wireless Applications.This code creates Figure 7.2, plot of a Bartlett.MATLAB code

代碼片段和文件信息
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%?????????????????????????????????????????????????????????%%
%%?******************************************************??%%
%%?*?Smart?Antennas?for?Wireless?Applications?w/?Matlab?*??%%
%%?******************************************************??%%
%%?????????????????????????????????????????????????????????%%
%%?Chapter?7:?Ex?7.5???????????????????????????????????????%%
%%?????????????????????????????????????????????????????????%%
%%?Author:?Frank?Gross?????????????????????????????????????%%
%%?McGraw-Hill?2005???????????????????????????????????????%%
%%?Date:??1/26/2004????????????????????????????????????????%%
%%?????????????????????????????????????????????????????????%%
%%?This?code?creates?Figure?7.2?plot?of?a?Bartlett????????%%
%%?????Pseudospectrum?for?theta1=?-5?-10?&?theta2?=?5?10?%%?
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%??????
%%-------------------?Define?Variables:--------------------%%
%?M?-?#?of?elements?in?array????????????????????????????????%
%?sig2?-?noise?variance?????????????????????????????????????%
%?th1?th2?-?angles?of?arrival??????????????????????????????%
%?a1?a2?a?-?array?steering?vectors????????????????????????%
%?A?-?Steering?vector?matrix????????????????????????????????%
%?Rss?-?Source?correlation?matrix???????????????????????????%
%?Rxx?-?Array?correlation?matrix????????????????????????????%
%?P?-?Bartlett?Psuedospectra????????????????????????????????%
%%---------------------------------------------------------%%
%%-----?Given?Values?-----%%
M?=?6;??sig2?=?0.1;
disp(‘?‘)
theta?=?input(‘??What?are?the?AOA?pairs?for?this?plot?in?degrees??Enter?5?or?10:??‘);
disp(‘?‘)
if?theta?==?5
????th1?=?-5*pi/180;?th2?=?5*pi/180;
elseif?theta?==?10
????th1?=?-10*pi/180;?th2?=?10*pi/180;
end
%%-----?Create?array?steering?vectors?a1?&?a2?-----%%
a1?=?[];??a2?=?[];
i?=?1:M;
%?a1?=?exp(1j*(i-1)*pi*sin(th1));
%?a2?=?exp(1j*(i-1)*pi*sin(th2));
a1?=?exp(1j*(i-1)*4*pi/16);
a2?=?exp(-1j*(i-1)*4*pi/16);
%%-----?Create?A?Rss?&?Rxx?-----%%
A?=?[a1.‘?a2.‘];??
%?Rss?=?[1?0;0?1];???Rxx?=?A*Rss*A‘?+?sig2*eye(M);
Rss?=?[1?0;0?1];
Rxx?=?A*Rss*A‘?+?sig2*eye(M);
%%-----?Create?Bartlett?Psuedospectrum?-----%%
for?k?=?1:180;
???th(k)?=?-pi/6?+?pi*k/(3*180);
???clear?a
???a?=?[];
???for?jj?=?1:M
??????a?=?[a?exp(1j*(jj-1)*pi*sin(th(k)))];
???end
???Pb(k)?=?real(a*Rxx*a‘);
end
Pb_dB?=?10*log10(Pb/max(Pb));
%%-----?Plot?Results?-----%%
plot(th*180/piPb_dB‘k‘)
grid?on
title([‘\bfFigure?7.2?-?Bartlett?Psuedospectra?for?\theta?=?+/-‘num2str(theta)‘^o‘])
xlabel(‘Angle‘)
ylabel(‘|P(\theta)|?(dB)‘)
axis([-30?30?-30?10])
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????2657??2013-12-28?11:23??sa_Ex7_5.m
-----------?---------??----------?-----??----
?????????????????2657????????????????????1
- 上一篇:用matlab實現的圖像二值化的三種方法
- 下一篇:DOA估計MATLAB代碼
評論
共有 條評論