www.pudn.com > leidaxitong.rar > xggszb.m

function y=xggszb1(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('滤波后的信号频谱图')