myawgn.m

来自「matlab仿真的GPS信号捕获程序」· M 代码 · 共 18 行

function [ g1,g2 ] = MyAwgn( n, div)

%****** complex Gauss white noise
%****** 均值为0
%n为向量维数
%div  %复高斯白噪声的方差
%实高斯白噪声的方差为div/2
mean = 0;%meanan
u1 = rand(1,n);%在【0，1】上均匀分布的n个数
u2 = rand(1,n);

a = sqrt(-2.*log(u1));%由随机信号分析得到
b = 2*pi.*u2;

g1 = sqrt(div/2)*a.*cos(b)+mean;%实高斯白噪声，作为实部
g2 = sqrt(div/2)*a.*sin(b)+mean;%实高斯白噪声，作为虚部