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
