信号分选CDIF累计差值直方图算法matlab仿真程序

clc; clear all; close all; %对应第一部分输入信号 s_noise=0.05 * randn(1, 100000);%噪声 x1=zeros(1,100000); x2=zeros(1,100000); x3=zeros(1,100000); %三个输入信号 %雷达1 PRI1=197;%us x1_start=1;%脉冲初始时间 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=1; 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=x1+x2+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 %对脉冲到达时间排序(冒泡排序） 这里可以用sort直接进行排序 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 % stem(y_daodashijian,'Marker','none');验证是否排好序 s_jieguo=zeros(1,10000); for t=1:4%CDIF阶数 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 %令1us为一个箱，10000us有10000箱子，脉冲统计中点为1、2、3··· s_xiang=zeros(1,10000)%进行脉冲累计 for i=1:10000 for j=1:(s_number-t) if s_jiange(1,j)==i s_xiang(1,i)=s_xiang(1,i)+1; end end end s_jieguo_guocheng=zeros(1,10000);%存放中间值，因为n阶CDIF需要小于n的每一阶累计值和第n阶的值进行叠加 for i=1:10000 s_jieguo_guocheng(1,i)=s_xiang(1,i); end s_jieguo=s_jieguo+s_jieguo_guocheng; figure(t); %绘图 stem(s_jieguo,'r','Marker','none'); hold on; x_kebiancanshu=0.6; xpri=(1e-10:0.1:1200); y=x_kebiancanshu*1e-1./(xpri*1e-6);%coeff1*time*10^6%门限绘制 plot(xpri,y) %门限 title([num2str(t),'阶CDIF直方图']) xlabel('PRI估计值/us'); ylabel('统计值/V'); axis([0 1200 0 600]); end