matlab编写雷达目标检测程序，并用MTI技术抑制杂波

clear; clc; load D.mat; load D_db.mat; % load D_se.mat; % load D_db_se.mat; fs=1000; [D_max,l_loc]=max(max(abs(D))); D_l_loc=D(:,l_loc); % [D_max1,h_loc]=max(abs(D_l_loc)); D_db_l_loc=D_db(:,l_loc); index=find(D_db_l_loc>=27); % D_major=D(index,l_loc); D_major=D_l_loc; N=length(D_major); for n=3:N MTI_major(n)=D_major(n)-2*D_major(n-1)+D_major(n-2); end MTI_major(1)=MTI_major(3); MTI_major(2)=MTI_major(3); D_MAJOR=fftshift(fft(D_major,N)); MTI_MAJOR=fftshift(fft(MTI_major,N)); f=linspace(-fs/2,fs/2,N); % figure() % plot(real(D_major)); figure(1) set(gcf,'color','white'); plot(f,abs(D_MAJOR)); xlabel('频率（Hz）'); ylabel('幅度'); title('直接fft'); % R_D=D_major*D_major'; % a=length(D_major); r_D=xcorr(D_major); r_MTI=xcorr(MTI_major); figure(2) set(gcf,'color','white'); plot((1-N:N-1)/1000,abs(r_D)); xlabel('derta tau(s)'); ylabel('幅度'); title('自相关分布'); R_D=fftshift(fft(r_D,length(r_D))); R_MTI=fftshift(fft(r_MTI,length(r_MTI))); figure(3) set(gcf,'color','white'); plot(linspace(-fs/2,fs/2,length(R_D)),abs(R_D)); title('功率谱密度'); xlabel('频率（Hz）'); ylabel('功率'); figure(4) set(gcf,'color','white'); plot(f,abs(MTI_MAJOR)); xlabel('频率（Hz）'); ylabel('幅度'); title('二次对消后直接fft'); figure(5) set(gcf,'color','white'); plot((1-N:N-1)/1000,abs(r_MTI)); xlabel('derta tau(s)'); ylabel('幅度'); title('二次对消后的自相关分布'); figure(6) set(gcf,'color','white'); plot(linspace(-fs/2,fs/2,length(R_MTI)),abs(R_MTI)); title('二次对消后的功率谱密度'); xlabel('频率（Hz）'); ylabel('功率'); % v=randn(1,4000); % V=xcorr(v); % VV=fftshift(fft(V)); % figure(3) % plot(abs(V)); % figure(4) % plot(abs(VV)); % tau=linspace(-4,4,2*N-1); % r_se=1/2*exp(-abs(tau)); % figure(4) % plot(abs(r_se)); % % % R_se=fftshift(fft(r_se)); % H=sqrt(R_se); % h_t=ifft(fftshift(H)); % % % % % h_t=ifft(H); % D_se=conv(D_major,h_t); % D_se=D_se(round(linspace(1,length(D_se),N))); % % D_SE=fftshift(fft(D_se,N)); % figure(5) % plot(f,abs(D_SE)); % xlabel('频率（Hz）'); % ylabel('幅度'); % title('色噪声直接fft'); % % r_se=xcorr(D_se); % figure(6) % plot(1-N:N-1,abs(r_se)); % xlabel('derta tau'); % ylabel('幅度'); % title('色噪声自相关分布'); % % R_SE=fftshift(fft(r_se,length(r_se))); % figure(7) % plot(linspace(-fs/2,fs/2,length(R_SE)),abs(R_SE)); % title('功率谱密度'); % xlabel('频率（Hz）'); % ylabel('功率');