%修改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箱已使用
????