• 大小: 2KB
    文件类型: .rar
    金币: 1
    下载: 0 次
    发布日期: 2021-06-01
  • 语言: 其他
  • 标签: TCT  

资源简介

利用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


评论

共有 条评论