clear?all;clc;close?all;
%%%%??仿真數(shù)據(jù)??%%%%
T=0.25;M=100;tc=T/M;fc=1/tc;%%%?一個掃頻周期T?為250ms?；?M為每個周期采樣點100?；?fc采樣頻率
N=256;n=1:N;%%?一個相關(guān)積累時間64
%%%
fs=1;%%??信號頻率
c_db=-35;c=10^（c_db/20）;%%?噪聲基底
s_db=-30;s=10^（s_db/20）;%%?信號幅度
signal=s.*exp（-j*2*pi*fs*T*n）;
nosic=c.*random（‘Normal‘011N）;
x=nosic+signal;
%%%
Fx=abs（fftshift（fft（x）））;
freq=（-1/（2*T）:1/（T*N）:（1/（2*T）-1/（T*N）））;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%?CA-CFAR??%%%%%%%%%%%%%%%%%%%%%%%%%
M=N;P_fa=10.^（-6）;
R=16;%參考單元數(shù)
n=R/2;
L_slipper=R+1;%滑窗長度
L_move=1;%滑窗間隔
L_num=floor（（M-L_slipper）/L_move）+1;%滑窗次數(shù)
Z=zeros（1L_num）;
for?i=1:L_num
????for?j=1:L_slipper
????????Z（i）=Z（i）+Fx（（i-1）*L_move+j）;
????end
????Z（i）=Z（i）/L_slipper;
end
T=P_fa.^（-1/R）-1
S=Z.*T;
figure;plot（freq20.*log10（Fx/max（Fx））‘k‘）;
hold?on;plot（freq（n+1:M-n）20.*log10（S/max（S）））;
xlabel（‘頻率‘）;ylabel（‘?歸一化幅值/dB‘）;title（‘CA-CFAR‘）;set（gcf‘color‘‘white‘）;grid?on;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????991??2014-05-04?09:44??CACFAR_sn.m