gaussianpu.asv

来自「雷达系统建模与仿真的matlab程序」· ASV 代码 · 共 28 行

function y=gaussianpu(x)
%由数据量为n的高斯白噪声产生向量为n，功率谱为高斯型的高斯随机向量
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
N=0:20;
f=20;
T=1/512;
c=2*f*T*sqrt(pi)*exp(-4*f^2*pi^2*T^2*N.^2);
n=length(x);
y=zeros(1,n);
for k=1:n
    for i=50:-1:0
        if ((k-i)<=0)
            continue;
        else
            y(k)=y(k)+c(51-i)*x(k-i);
        end
    end
    for i=50:100
        if ((k-i)<=0)
            continue;
        else
            y(k)=y(k)+c(i-49)*x(k-i);
        end
    end
end
y=0.5*y;