xggszb.m

来自「雷达噪声信号的仿真」· M 代码 · 共 41 行

%
%产生服从高斯分布的随机数
%
%clear;
%clc;
N=1024;
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=conv(c,x);


subplot(211);
plot(y);
subplot(212);
pyulear(y,8,'twosided');