%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