matlab正弦波高斯白噪声均匀白噪声功率谱密度自相关函数.doc

现代通信原理 作业一 姓名：张英伟 学号：133320085208036 班级：13级理工部3班 利用matlab完成： 产生正弦波信号、均匀白噪声以及高斯白噪声并分别将两种噪声叠加到正弦波信号上，绘出波形。 分别求取均匀白噪声序列和高斯白噪声序列的自相关及功率谱密度，绘出波形。 白噪声区别及产生方法 1、定义： 均匀白噪声：噪声的幅度分布服从均匀分布，功率谱密度在整个频域内均匀分布的噪声在整个频域内均匀分布的噪声是区间在[0,1]的随机数=a+(b-(a))*rand(m,n)............（公式1） randn函数默认产生均值是0、方差是1的随机序列，所以可以用其来产生均值为0、方差为1的正态分布白噪声，即N(0，12)。利用公式： z1=a+b*randn(1,n).................（公式2） 可以产生均值为a，方差为b2 高斯白噪声，即N(a，b2)。 二、自相关函数与功率谱密度之间的关系 1、功率谱密度：每单位波携带的功率，这被称为信号的功率谱密度自相关函数描述随机信号X(t)在任意两个不同时刻t1，t2的取值之间的相关程度维纳-辛钦定理由于平均值不为零的信号不是平方可积的，所以在这种情况下就没有傅里叶变换。幸运的是维纳-辛钦定理提供了一个简单的替换方法，如果信号可以看作是平稳随机过程，那么功率谱密度就是信号自相关函数的傅里叶变换。 平稳随机过程是在固定时间和位置的概率分布与所有时间和位置的概率分布相同的随机过程x=(0:0.01:2); %采样频率100Hz y1=sin(10*pi*x); %产生频率5Hz的sin函数 plot(x,y1,b); 高斯白噪声+正弦波 z1=0.1*randn(1,201); %产生方差N（0，0.12）高斯白噪声 （b=0.01/0.1/1） plot(x,z1,b); y2=y1+z1; %叠加高斯白噪声的正弦波 plot(x,y2,b); 均匀白噪声+正弦波 z2=-.3+.6*rand(1,201); %产生-0.3到0.3的均匀白噪声 plot(x,z2,b); y3=y1+z2; %叠加均匀白噪声的正弦波 plot(x,y3,b); 4、高斯白噪声序列自相关函数及功率谱密度 z1=0.1*randn(1,201); %产生方差N（0，0.12）高斯白噪声 [r1,lags]=xcorr(z1); %自相关函数的估计 plot(lags,r1); f1=fft(r1); f2=fftshift(f1); %频谱校正 l1=(0:length(f2)-1)*200/length(f2)-100; %功率谱密度x轴 y4=abs(f2); plot(l1,y4); 5、均匀白噪声序列自相关函数及功率谱密度 z2=-.3+.6*rand(1,201); %产生-0.3到0.3的均匀白噪声 [r2,lags]=xcorr(z2); %自相关函数的估计 plot(lags,r2); f3=fft(r2); f4=fftshift(f3); %频谱校正 l2=(0:length(f4)-1)*200/length(f4)-100; %功率谱密度x轴 y5=abs(f4); plot(l2,y5);