pri传统分选算法

此matlab 程序主要是针对3部重频固定的雷达信号进行分选，针对许多程序中与原理图不合理的地方进行了修改，符合了原理，经仿真可以准确分选出3部雷达信号。

%Original algorithm based on PRI transform applies to interleaved pulse  交错脉冲
%train with constant PRIs 1、sqrt（2）、sqrt（5）
%parameter  number of pulses N=1000;range of PRI
%[taumintaumax]=[010];numbers of PRI bins K=201   PRI箱数K 
%本程序是对雷达信号分选算法中的PRI变换所进行的仿真，主要针对3部重频固定的雷达进行信号分选
%!!!根据原理图更改的效果不如之前的好
clear all
clc
t1=0:333;
t2=0.1:sqrt（2）:（0.1+332*sqrt（2））;
t3=0.2:sqrt（5）:（0.2+332*sqrt（5））;
t=[t1 t2 t3];
clear t1 t2 t3
t=sort（t）;%升序排列，打乱次序方便分选
N=length（t）;
K=201;
taumin=0;
taumax=10;  
b=（taumax-taumin）/K;%b是箱宽，
C=zeros（1K）;
D=zeros（1K）;
for i=1:K
    tauk（i）=（i-1/2）*（taumax-taumin）/K+taumin;%tauk是第k个PRI箱的中心坐标=（k-0.5）*b+taumin
end
n=2;
while n<=N
    m=n-1;
    while m>=1
        tau=t（n）-t（m）;
        if （tau>taumin）&（tau<=taumax）
            for k=1:K
                if （tau>（tauk（k）-b/2））&（tau<=（tauk（k）+b/2））
                    D（k）=D（k）+exp（2*pi*t（n）*j/tau）;
                    C（k）=C（k）+1;%c（k）是每个小箱中脉冲个数
                end
            end
        elseif tau>taumax
            break; %源代码是m=m-1 根据原理应该是这个
        else
            ;%源代码是n=n+1 根据原理应该是这个             
        end
        m=m-1;
    end
    n=n+1;
end
plot（taukabs（D））
axis（[0 10 0 800]）
hold on         
X=225./tauk;                                                     %门限函数 观察时间原则
plot（taukX‘r-‘）
xlabel（‘tau‘）
ylabel（‘|D（k）|‘）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1507  2018-09-11 20:31  tradtionpritransform.m
     文件     3301091  2018-09-11 16:26  传统PRI变换法流程图.jpg
     文件         145  2018-09-11 20:35  新建文本文档.txt