pritransform.m

PRI变换法

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


tmin=0.5;tmax=10;              %PRI取值范围设定 
K=201;                        %PRI箱个数，K越大，PRI变换图越精确
b=(tmax-tmin)/K;              %每个PRi箱宽度，此处也代表求得PRI的精度
kk=[1:K];
tk=(kk-0.5)*b+tmin;            %每个PRI箱的中心坐标
NN=length(C);
                      %用来标记当前PRI箱


Dk=zeros(size(tk));           %PRI变换函数Dk
Ck=zeros(size(tk));           %标记脉冲数函数
  
    for nn=2:NN
            for mm=1:nn-1
                    tao=C(nn)-C(mm);                     %对脉冲到达时间进行重复相减计算差值 
                                     
                    if tao>=tmin&tao<=tmax,               %防止超限制的处理    
                        
                        for l=1:K
                            
                        if tao<=tk(l)+b/2&tao>=tk(l)-b/2,
                             Ck(l)=Ck(l)+1;
                            Dk(l)=Dk(l)+exp(2*pi*i*C(nn)/tao);   %对在当前PRI箱内的到达时间差值，累计Dk值                                  
                        end  
                        
                        end
                        
                    end
            end
    end
   

for ll=1:K
 Dk(ll)=abs(Dk(ll));                %对Dk取模
end

Dk1=zeros(size(tk));
%设定门限
aa=0.3;
bb=0.3;
cc=6;
for v=1:length(tk) 
    
       Men1(v)=aa*TT/tk(v);                     %代表以tt为PRI的脉冲个数  150是观察时间
       Men2(v)=bb*Ck(v);                     %为了消除子谐波
       Men3=cc*(TT*(length(C)/TT)^2*b)^0.5;    %为消除噪声
       Dk1(v)=max([Men1(v),Men2(v),Men3]);
   end
plot(tk,Dk,'k',tk,Dk1,'k');                      %做出PRI变换图