信号分选SDIF序列差直方图算法附matlab代码.zip

clc; clear all; close all; % x=0.5; % T_pluse_width=2e-6; % T_detection_threshold=x*(1e-6/T_pluse_width);%单位 s_noise=0.05 * randn(1, 100000); x1=zeros(1,100000); x2=zeros(1,100000); x3=zeros(1,100000); %雷达1 % PRI1=197;%us % T1=2e-6;%脉冲宽度 % x1_start=0.0001; % x1_range=100000; % x1(x1_start:PRI1:x1_range)=1;%x(2:pri_3:x_range)+1; % x1=s_noise+x1; % 雷达2 % PRI2=332;%us % T2=2e-6;%脉冲宽度 % x2_start=0.00001; % x2_range=100000; % x2(x2_start:PRI2:x2_range)=1;%x(2:pri_3:x_range)+1; % x2=s_noise+x2; % 雷达3 PRI3=515;%us T3=2e-6;%脉冲宽度 x3_start=1; x3_range=100000; x3(x3_start:PRI3:x3_range)=1;%x(2:pri_3:x_range)+1; % x3=s_noise+x3; x=x3+s_noise; for i=1:100000%去噪 if x(1,i)<0.5 x(1,i)=0; end end %脉冲个数 s_number=0; for i=1:100000 if x(1,i)~=0 s_number=s_number+1; end end %到达时间 y_daodashijian=zeros(1, s_number); a=1;%此程序所有a均用来暂存数据 for i=1:100000 if x(1,i)~=0 y_daodashijian(1,a)=i; a=a+1; end end %到达时间排序(冒泡排序） a=0;%此程序所有a均用来暂存数据 for i=1:( s_number-1) for j=1: s_number-i if y_daodashijian(1,i)>y_daodashijian(1,i+1) a=y_daodashijian(1,i+1); y_daodashijian(1,i+1)=y_daodashijian(1,i); y_daodashijian(1,i)=a; end end end figure; % stem(y_daodashijian,'Marker','none');%验证是否排好序 for t=1:1 s_jiange=zeros(1,s_number-t); for i=1:(s_number-t) s_jiange(1,i)=y_daodashijian(1,i+t)-y_daodashijian(1,i); end %正常令5us为一个箱，10000us应该有2000箱子，但是观察到相邻脉冲间隔 % 最大不到200us，所以设置40个箱 s_xiang=zeros(1,10000) for i=1:10000 for j=1:(s_number-t) % if (s_jiange(1,j)>=5*(i-1))&&( s_jiange(1,j)<(i*5)) if s_jiange(1,j)==i s_xiang(1,i)=s_xiang(1,i)+1; end end end s_jieguo=zeros(1,10000); for i=1:10000 s_jieguo(1,i)=s_xiang(1,i); end figure(t); stem(s_jieguo,'r','Marker','none'); hold on; x_kebiancanshu=0.4; k=0.15; C=t; xpri=(1e-10:0.1:1200); y=x_kebiancanshu*(s_number-C)*exp(-xpri/k/2000);%coeff1*time*10^6 %N为直方图小区间个数（箱数） plot(xpri,y) %门限 title(['第',num2str(t),'阶SDIF直方图']) xlabel('脉冲到达时间（TOA)间隔差值/us');ylabel('统计数量/个'); axis([0 1200 0 500]); end