clear?all;
close?all;
clc;
%%%%%%%%%%%產生瑞利噪聲%%%%%%%%%%%%%%%%%%%%%%%
sigma=2;
t=1e-3;
fs=1e6;
ts=1/fs;
t1=0.05e-3:1/fs:0.2e-3-1/fs;
n=length（t1）;
rand（‘state‘0）;
u=rand（1n）;
rayleigh_noise=sqrt（2*log2（1./u））*sigma;

%%%%%%%%%%%%產生目標回波%%%%%%%%%%%%%%%%%%%
N=t/ts;
s_pc_1=[zeros（1100）1zeros（1N-101）];
noise=rand（1N）;
rayleigh_clutter=[zeros（150）rayleigh_noisezeros（1N-200）];
s_pc=s_pc_1+0.1*rayleigh_clutter+0.1*noise;
figureplot（（0:ts:t-ts）s_pc）
xlabel（‘t（單位:s）‘）title（‘疊加了瑞利分布雜波、熱噪聲的目標回波‘）;

%%%%%%%%%%%%%慢門限恒虛警處理%%%%%%%%%%%%%%%%%%%
?cfar_result=zeros（1N）;
?cfar_result（11）=s_pc（11）;
?for?i=2:N
?????cfar_result（i）=s_pc（1i）/mean（s_pc（11:i））;
?end
?figureplot（（0:ts:t-ts）cfar_result）
xlabel（‘t（單位:s）‘）title（‘采用慢門限處理結果‘）;

%%%%%%%%%%%%%%%快門限恒虛警處理%%%%%%%%%%%%%
cfar_k_result=zer