nb_gaussnoise.m

来自「包括载波跟踪和伪码相位跟踪」· M 代码 · 共 15 行

%产生窄带高斯噪声
function [Gauss_Noise_NB,Noise_NB_Power] = NB_GaussNoise(n0,Len,Bandwidth,fs,fi)
% n0:产生的无限带宽的高斯白噪声方差即功率     
% len:产生的噪声序列长度
% Gauss_Noise_NB:产生的窄带高斯噪声
% Noise_NB_Power:窄带高斯噪声功率
Gauss_Noise = n0*randn(1,Len);      % randn产生均值为0、方差为1的随机序列
w1 = fi-Bandwidth/2;               % 以数字中频为中心的频带上下限               
w2 = fi+Bandwidth/2;
wp = [w1 w2]*2/fs;                  % 归一化数字截止频率
Nm = 15;
b = fir1(2*Nm,wp);                  % 31个系数，即30阶海明窗FIR带通滤波           l         
Gauss_Noise_NB = filtfilt(b,1,Gauss_Noise);
Noise_NB_Power = Gauss_Noise_NB*Gauss_Noise_NB'/Len;