radar1.m

对MTD雷达中进行杂波滤除的滤波器组

close all
clc
clear all
N = 2e4;
NN = 100;

%方差
pn0 = 1;
pn1 = 10;
%pn1 = 1;   

A1 = rand(1,N);        
A2 = rand(1,N);

%高斯分布
ni = pn1*(-2*log(A1)).^(1/2).*cos(2*pi*A2);
nq = pn1*(-2*log(A1)).^(1/2).*sin(2*pi*A2);
n0 = pn0*(-2*log(A1)).^(1/2).*cos(2*pi*A2);
nz1 = ni-j*nq;
%包络情况
nz2 = (ni.^2+nq.^2).^(1/2);   

S = randn (1,N);   
S = 10*randn(1,N);
%瑞利分布
nRay = sqrt(-2*pn1.^2*log(1-A1));
%指数分布
nexp = -log(1-A1);

%统计直方图
 x0 = linspace(-5,5,NN);
 x1 = linspace(0,5,NN);


figure(1)
hist(ni,100)
% hold on
figure(2)
hist(S,100)
figure(3)
hist(nz2,100)
%axis([0 5 0 700])
figure(4)
hist(nRay,NN)
figure(5)
hist(nexp,NN)

%高斯分布的自相关函数和功率谱
ni_gauss = xcorr(ni,'biased');
ni_Pw = abs(fft(ni_gauss));

figure(6)
plot(ni_gauss),
figure(7)
plot(10*log10(ni_Pw/max(ni_Pw)));

xn = zeros(1,NN);
xn(50) = 1;
figure(8)
plot(abs(fft(xn,128)));

figure(9)
plot(ni)
figure(10)
plot(n0);