tct-doa估计

利用TCT算法实现的对多个信号进行doa估计。

clear;
clc;
Sensor=8;                                     %天线阵元数
bw1=1e7;                                      % 信号带宽
bw2=1e6;
bw3=51e6;
f1=1e8;
f2=1e7;
f3=1e7;
T1=1e-2;                                      %信号脉冲宽度  
L=1024;                                       %采样点数
dt1=T1/L;                                     %采样间隔
gc=3e8;                                       %光速
snr=10;                                       %信躁比
degrad=pi/180;
p=3;                                          %信号源个数 
angle=10;           
angle2=20;                                    %信号源角度
angle3=-40;  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%构造采样序列
t1=（0:L-1）*dt1;                               %时域采样点
sig1=exp（j*pi*bw1*t1.*t1/T1+j*2*pi*f1*t1）;    %产生宽带信号
sig2=exp（j*pi*bw2*t1.*t1/T1+j*2*pi*f2*t1）;  
sig3=exp（j*pi*bw3*t1.*t1/T1+j*2*pi*f3*t1）;  
sig1=awgn（sig1snr）;
sig2=awgn（sig2snr）;                          %加白噪音
sig3=awgn（sig3snr）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%宽带信号频率
F=f1+bw1/（T1）*t1;                           
f1=F（128）;
f2=F（2*128）;
f3=F（3*128）;
f4=F（4*128）;
f5=F（5*128）;
f6=F（6*128）;
f7=F（7*128）;
f8=F（8*128）;
F=[f1f2f3f4f5f6f7f8];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%构造采样数据矩阵
Lambda=gc/f8;                                 %波长
d=0.5*Lambda;                                 %阵元间距离
tmp2=d*[0:Sensor-1 ]‘;            

v11=sig1（1:128）;
v12=sig1（129:256）;
v13=sig1（257:384）;
v14=sig1（385:512）;
v15=sig1（513:640）;
v16=sig1（641:768）;
v17=sig1（769:896）;
v18=sig1（897:1024）;
v21=sig2（1:128）;
v22=sig2（129:256）;
v23=sig2（257:384）;
v24=sig2（385:512）;
v25=sig2（513:640）;
v26=sig2（641:768）;
v27=sig2（769:896）;
v28=sig2（897:1024）;
v31=sig3（1:128）;
v32=sig3（129:256）;
v33=sig3（257:384）;
v34=sig3（385:512）;
v35=sig3（513:640）;
v36=sig3（641:768）;
v37=sig3（769:896）;
v38=sig3（897:1024）;

s1=[（v11）;（v21）;v31];
s2=[（v12）;（v22）;v32];
s3=[（v13）;（v23）;v33];
s4=[（v14）;（v24）;v34];
s5=[（v15）;（v25）;v35];
s6=[（v16）;（v26）;v36];
s7=[（v17）;（v27）;v37];
s8=[（v18）;（v28）;v38];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%构造方向矩阵
A1=[exp（-j*2*pi*f1/gc*tmp2*sin（angle*degrad））exp（-j*2*pi*f1/gc*tmp2*sin（angle2*degrad））exp（-j*2*pi*f1/gc*tmp2*sin（angle3*degrad））]; 
A2=[exp（-j*2*pi*f2/gc*tmp2*sin（angle*degrad））exp（-j*2*pi*f2/gc*tmp2*sin（angle2*degrad））exp（-j*2*pi*f2/gc*tmp2*sin（angle3*degrad））]; 
A3=[exp（-j*2*pi*f3/gc*tmp2*sin（angle*degrad））exp（-j*2*pi*f3/gc*tmp2*sin（angle2*degrad））exp（-j*2*pi*f3/gc*tmp2*sin（angle3*degrad））]; 
A4=[exp（-j*2*pi*f4/gc*tmp2*sin（angle*degrad））exp（-j*2*pi*f4/gc*tmp2*sin（angle2*degrad））exp（-j*2*pi*f4/gc*tmp2*sin（angle3*degrad））]; 
A5=[exp（-j*2*pi*f5/gc*tmp2*sin（angle*degrad））exp（-j*2*pi*f5/gc*tmp2*sin（angle2*degrad））exp（-j*2*pi*f5/gc*tmp2*sin（angle3*degrad））]; 
A6=[exp（-j*2*pi*f6/gc*tmp2*sin（angle*degrad））exp（-j*2*pi*f6/gc*tmp2*sin（angle2*degrad））exp（-j*2*pi*f6/gc*tmp2*sin（angle3*degrad））]; 
A7=[exp（-j*2*pi*f7/gc*

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       6738  2015-05-04 11:11  tctchange.m

----------- ---------  ---------- -----  ----

                 6738                    1