快速PRI信号分选

%修改PRI变换法分选程序
%作用：对已有脉冲序列进行修改PRI分选，绘制PRI谱图，以及门限
%需要先生成雷达信号
%作者：巫书航
%时间：2007年11月13日


tmin=0.5;tmax=10;              %PRI取值范围设定 
K=200;                        %PRI箱个数，K越大，PRI变换图越精确
doum=0.05;                    %脉冲抖动上限
b=（tmax-tmin）/K;              %每个PRi箱宽度，此处也代表求得PRI的精度
k=[1:K];
tk=（k-0.5）*b+tmin;            %每个PRI箱的中心坐标
for tk1=1:length（tk）
bk（tk1）=max（[b2*doum*tk（tk1）]）;                 %每个PRI箱宽度
end
nn=length（C）;
Zi0=0.03;



Dk=zeros（size（tk））;           %PRI变换函数Dk
Ck=zeros（size（tk））;           %标记脉冲数函数
Ok=zeros（size（tk））;           %时间起点函数
JI=zeros（size（tk））;           %标记函数

    for nn=2:NN
            for mm=1:nn-1
                    tao=C（nn）-C（mm）;                     %对脉冲到达时间进行重复相减计算差值 
                                     
                    if tao>（1-doum）*tmin&tao<（1+doum）*tmax               %防止超限制的处理    
                        
                        for l=1:K
                            
                        if tao>tk（l）-bk（l）/2&tao                            Ck（l）=Ck（l）+1;
                            if JI（l）==0Ok（l）=C（nn）;JI（l）=1;
                            end                                  %若PRI箱第一次使用，则选时间起点为C（nn）并标记PRI箱已使用