?? ex9_14.m
字號:
clear all;
close all;
clc;
%%%%%%%%%%%產(chǎn)生瑞利噪聲%%%%%%%%%%%%%%%%%%%%%%%
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(1,n);
rayleigh_noise=sqrt(2*log2(1./u))*sigma;
%%%%%%%%%%%%產(chǎn)生目標(biāo)回波%%%%%%%%%%%%%%%%%%%
N=t/ts;
s_pc_1=[zeros(1,100),1,zeros(1,N-101)];
noise=rand(1,N);
rayleigh_clutter=[zeros(1,50),rayleigh_noise,zeros(1,N-200)];
s_pc=s_pc_1+0.1*rayleigh_clutter+0.1*noise;
figure,plot((0:ts:t-ts),s_pc),
xlabel('t(單位:s)'),title('疊加了瑞利分布雜波、熱噪聲的目標(biāo)回波');
%%%%%%%%%%%%%慢門限恒虛警處理%%%%%%%%%%%%%%%%%%%
cfar_result=zeros(1,N);
cfar_result(1,1)=s_pc(1,1);
for i=2:N
cfar_result(i)=s_pc(1,i)/mean(s_pc(1,1:i));
end
figure,plot((0:ts:t-ts),cfar_result),
xlabel('t(單位:s)'),title('采用慢門限處理結(jié)果');
%%%%%%%%%%%%%%%快門限恒虛警處理%%%%%%%%%%%%%
cfar_k_result=zeros(1,N);
cfar_k_result(1,1)=s_pc(1,1)/(sqrt(2)/pi*mean(s_pc(1,2:17)));
%第l點(diǎn)恒虛警處理時(shí)噪聲均值由其后面的l6點(diǎn)的噪聲決定
for i=2:16
%第2點(diǎn)到第16點(diǎn)的恒虛警處理的噪聲均但由其前面和后面l6點(diǎn)的噪聲共同決定
noise_mean=sqrt(2)/pi*(mean(s_pc(1,1:i-1))+mean(s_pc(1,i+1:i+16)))/2;
cfar_k_result(1,i)=s_pc(1,i)/noise_mean;
end
for i=17:N-17
%正常的數(shù)據(jù)點(diǎn)的恒虛警處理的噪聲均值由其前面和后面各26點(diǎn)的噪聲中的決定
noise_mean=sqrt(2)/pi*max(mean(s_pc(1,i-16:i-1)),mean(s_pc(1,i+1:i+16)));
cfar_k_result(1,i)=s_pc(1,i)/noise_mean;
end
for i=N-16:N-1
%例數(shù)第l6點(diǎn)到倒數(shù)第2點(diǎn)恒虛警處理的噪聲均值由其前面16和后面的噪聲共同決定
noise_mean=sqrt(2)/pi*(mean(s_pc(1,i-16:i-1))+mean(s_pc(1,i+1:N)))/2;
cfar_k_result(1,i)=s_pc(1,i)/noise_mean;
end
%最后-點(diǎn)的桓虛警處理的噪聲均值由其前面l6點(diǎn)的噪聲決定
cfar_k_result(1,N)=s_pc(1,N)/(sqrt(2)/pi*mean(s_pc(1,N-16:N-1)));
figure,plot(0:ts:t-ts,cfar_k_result),xlabel('t(單位:s)'),title('采用快門限恒虛警處理處理結(jié)果');
?? 快捷鍵說明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -