test.m

宽带信号测向的几种直接处理算法

function test
fs=360; fl=80; fh=120;  %频率选择
% N=10;
L=3;
TargetNum=2;            %声源数目
ang=[5,12];            %目标方位
theta0=10.5;
dcss=2;
D=5;
Nf=33;
Td=0.00625;
x=16;                    %ULA阵元数
Snaps=512;              %快拍数
angles=[0:.1:20];        %宽带波束扫描范围(度)
f0=(fl+fh)/2;
% d=0.067;d1=0.2;d2=0.256; %计算阵元坐标
% px1=[0:4]'*d;
% px2=[0:4]'*d+px1(5)+d1;
% px3=[0:4]'*d+px2(5)+d2;
% px4=[0:4]'*d+px3(5)+d1;
% x=[px1;px2;px3;px4];
% for ang=-20:2:20
    SNR=5;
%     ang
% load('PSD715')
% data = SenSG(fs,fl,fh,Snaps,ang,SNR,x,PSDnoise,PSDsig);
data = SenSG(fs,fl,fh,Snaps,ang,SNR,x);
PCSS=CSS(data,angles,fl,fh,fs,dcss,theta0,x);
P=BASS(data,angles,fl,fh,fs,Td,L,Nf,D,x);
PFDM=FDM(data,angles,fl,fh,fs,L,Nf,D,x);
FP=FBASS(data,angles,fl,fh,fs,Td,L,Nf,D,x,theta0);
plot(angles,10*log10(PCSS/max(PCSS)),angles,10*log10(P/max(P)),angles,10*log10(PFDM/max(PFDM)),angles,10*log10(FP/max(FP)));
grid on
legend('CSS','BASS','FDM','FBASS','Location','Best')
line(ang(1),-5:.2:0);
line(ang(2),-5:.2:0);