xggszb1.m

有关雷达系统的仿真

function y=xggszb(N)
%
%产生相关高斯杂波
%

%
%产生服从高斯分布的随机数
%
clear;
clc;
u1=rand(1,N);  %标准正态分布的随机数
u2=rand(1,N);  %
for i=1:N
    x(i)=sqrt((-2)*log2(u1(i)))*cos(2*pi*u2(i));
end

%
%求滤波器系数
%
f=20;
f0=512;
t=1/f0;
k1=2*f*t*sqrt(pi);
k2=(-4)*(f^2)*(pi^2)*(t^2);
c(1)=k1;
for i=2:11
    c(i)=k1*exp(k2*((i-1)^2));
    c(i)=c(i)/2;
    c(i+10)=c(i);
end
c1=c;
for i=1:5
    t=c(i);
    c(i)=c(12-i);
    c(12-i)=t;
end                %最终产生的滤波器系数是偶对称的

%
%对信号进行滤波
%
y=filter(c,1,x);
z=fft(y);

subplot(131);
hist(x,50);  
title('原始信号');

subplot(132);
hist(y,50);
title('滤波后的信号');

%subplot(133);
%plot(z);
%title('滤波后的信号频谱图')