clear?all
clc
%開始分選子程序

PRI1=150e-006;%單位us
a=0.05;
b=a*PRI1*2*（rand-0.5）;
t1=1:126;
%?z1=PRI1*t1+0.000223+b;
z1=PRI1*t1+2.23e-006;
%?z11=z1+z1.*（1-2*rand（1126））*0.01;
PRI2=240e-006;
t2=1:78;
z2=PRI2*t2+11e-006;
PRI3=370e-006;
t3=1:50;
z3=PRI3*t3+0.120e-006;
sig=[z1z2z3];

dd=1e-3;%dd=0;?dd=1e-3;?dd=1e-2??????????????????????????????????????????%設置抖動程度0和0.1%和1%
jitter=（1-2*rand（1254））*dd;
sig=sig+sig.*jitter;
sig=sort（sig）;
%以上程序用于生成交錯脈沖序列到達時間，交錯脈沖序列由三部雷達組成

sort_para.pnum=1;
sort_para.tnum=254;
sort_para.toa=sig;
sdif_para.pri_fixr=5e-006;

sdif_para.step=3;
sdif_para.pri_id=1;
sdif_para.gate=0.1;
%[sort_bank]=sdif_picture（sort_parasdif_para）

time=37e-005;%采樣時間最大重復周期10倍，
fs=250e006;%采樣頻率
toa_sdif=sort_para.toa;
tnum_sdif=sort_para.tnum;%調出要使用修改的參數
sig_num=1;
pri_num=0;
pri=0;
pri_total=0;
sig_total=0;
sort_id=0;
c=1;
tem_pri=0;
gate_num=0;
pri_oncenum=0;
pri_oncevalue=0;
pri_twicenum=0;
pri_twicevalue=0;
pri_search=0;
pri_searchnum=0;
search_num=1;
search_id=0;
tem_toa=0;
sortok_id=0;
sortok_id1=0;
%可能要調節的參數
coeff1=0.45;%0.8;
coeff2=0.3;
series=3;%最多做c級差
L=6;


ppri_num=0;
ppri_fnum=0;
ppri=zeros（18）;



%逐步提取脈寬對應TOA數據進行信號分選
for?i=1:sort_para.pnum??%%信號脈寬個數
????%變量初始化
????pri_num=0;
????pri=0;
????pri_total=0;
????sig_total=0;
????pri_ratio=0;
????
????clear?pri_save;
????pri_save.pri=0;????????????????????????????%保存一級差脈沖重復周期
????pri_save.count=0;??????????????????????????%保存各脈沖重復周期出現個數
????pri_save.pri_num=0;????????????????????????%保存一級差脈沖重復周期數
????pri_save.pri_total=0;
????
????%進行多級差計算
????for?step_i=1:sdif_para.step??%%?step:分選級數
????????%如果TOA個數小于門限，跳出此層循環，準備提取下一個脈寬內的數據
????????if?sort_para.tnum（i）<5+step_i???%%?tnum:信號到達時間個數
????????????break;
????????else
????????????pri_num=0;
????????????pri=0;
????????????pri_total=0;
????????????sig_total=0;
????????????pri_ratio=0;
????????????pri_gate_num=0;
????????????
????????
????????????child_id=0;
????????????child_num=0;
????????????child_fnum=0;
????????
????????????%TOA個數大于門限，則準備分選先進行pri分類統計
????????????for?j=1:sort_para.tnum（i）-step_i???%%%第step_i級一共有sort_para.tnum（i）-step_i個間隔
????????????????%先求step_i級TOA差
????????????????tem_pri=sort_para.toa（ij+step_i）-sort_para.toa（ij）;
????????????????????????????????
????????????????if?pri_num>0?%?pri數
????????????????????for?k=1:pri_num
????????????????????????%得到pri容限?
????????????????????????if?sdif_para.pri_id==1????%%分選重復周期容限?固定
????????????????????????????pri_range=sdif_para.pri_fixr;
????????????????????????else???%%分選重復周期容限?變化
????????????????????????????pri_range=floor（pri（k）*sdif_para.pri_ranr/100）;
????????????????????????end
????????????????????????%求pri下限
????????????????????????if?pri（k）????????????????????????????pri_low=0;
????????????????????????else
????????????????????????????pri_low=pri（k）-pri_range;
????????????????????????end
????????????????????????%求pri上限
????????????????????????p