通信系统仿真原理与无线应用.pptVIP

*5.5.2确定任意的功率谱密度或自相关函数思路：对一组不相关的样本进行适当的滤波，从而使之具有目标功率谱密度。不相关样本的功率谱密度在仿真带宽|f|fs/2上是常数，方差为随机数（噪声）的方差和采样频率的关系对于独立样本，其功率谱密度为常数K=N0/2Watts/Hz为得到目标功率谱密度，所需的H(f)为*随机过程通过线性系统独立样本的PSD相关随机序列的产生*举例1：设计具有给定传递函数的滤波器方法：确定传递函数在最小均方误差意义下的最佳拟合。可通过MATLAB函数yulewalk求解Yule-Walker方程来完成。实现：利用yulewalk求出滤波器的系数bk和ak，使得传递函数是给定传递函数的最小均方误差拟合。*假定给定的功率谱密度具有S(f)=K/f的形式（闪烁噪声），用来建模振荡器的相位噪声。要产生这种噪声，必须使用一个滤波器，使其传递函数为。因为当f?0时,H(f)??,故传递函数定义为利用MATALAB实现H(f)，其中频率用Nyquist频率做归一化处理，即有f0=fN/20。*f=0:200;%frequencypointsfn=200;%NyquistrateF=f./fn;%frequencyvectorM=abs(1./sqrt(F));%normalizedfrequencyresponseM=[zeros(1,11),M(11:200)];size(M)%boundfromzerofrequency[b1,a1]=yulewalk(3,F,M);%generateorder=3filter[b2,a2]=yulewalk(20,F,M);%generateorder=20filter[h1,w1]=freqz(b1,a1);%generate3-rdorderH(f)[h2,w2]=freqz(b2,a2);%generate20-rdorderH(f)subplot(2,1,1)plot(F,M,:,w1/pi,abs(h1))xlabel(归一化频率)ylabel(幅度响应)subplot(2,1,2)plot(F,M,:,w2/pi,abs(h2))xlabel(归一化频率)ylabel(幅度响应)*图中虚线表示期望特性，实线表示近似特性上图为3阶近似，下图为20阶近似*举例2:用于无线移动系统仿真的随机过程所求滤波器的传递函数为*%Generateandtesttheimpulseresponseofthefilter.clearall;clc;fd=100;%maximumdopplerimpw=jakes_filter(fd);%calltoJakesfilterfs=16*fd;ts=1/fs;time=[1*ts:ts:128*ts];subplot(2,2,1)stem(time,impw,.);grid;xlabel(Time);ylabel(ImpulseResponse);%----------------------------------------------------------%SquaretheFFTandcheckthepowertransferfunction[hf]=linear_fft(impw,128,ts);%generateH(f)forfiltersubplot(2,2,2)plot(f,abs(h.*h));grid;xlabel(Frequency);ylabel(PSD)*%PutGaussiannoisethroughandchecktheoutputpsdx=randn(1,1024);%generateGaussianinputy=filter(impw,1,x);%filterGaussianinput[output_psdff]=