通信系统的计算机模拟第9讲.pptVIP

通信系统的计算机模拟 第九讲 随机信号的产生与处理 信道噪声、干扰和衰落等随机 波形级精确地仿真这些系统，随机影响建立准确的模型 其基本构建模块是随机数发生器。如何产生采样后的随机波形（信号、干扰和噪声等），以用于仿真。在仿真环境下，所有的随机过程必须用随机变量序列来表示。 随机数发生器产生的是“伪随机序列”，因为尽管它们是确定的，但在具体应用中会呈现随机性。 精度随应用不同：产生一个波形来表示锁相环鉴相器输入端的噪声，对一个输入SNR为50dB的系统建立噪声波形模型所要求的精度要比输入端SNR为8dB的系统高得多。 问题 在（0，1）上产生均匀分布且不相关的随机数 将不相关且均匀分布的随机数映射成不相关的，具有任意（需要的）概率密度函数（pdf）的随机 产生不相关的、具有高斯型pdf的随机数 产生相关的、具有高斯型pdf的随机数 产生相关的、具有任意（需要的）pdf的随机数 伪噪声（PN）序列的产生以及几种应用于随机数序列的计算方法。 平稳与遍历性过程 信号、噪声和干扰而产生的样本函数通常假设为各态历经的。 原因：依次处理系统中波形的时域样本，系统中的每个点上只有一个波形（样本函数）。 假设仿真所处理的波形是其内在统计模型定义的总体（Emsemble）中的一个典型成员，这时各种统量如各阶矩、信噪比和误比特率，就可以当作时间平均来计算。仿真计算得到的时间平均等于总体均值。隐含的假设条件是对应的随机过程为遍历性的。 例7-1 假设随机过程的样本函数表达式为 f = 1; % frequency of sinusoid fs = 100; % sampling frequency t = (0:200)/fs; % time vector for i=1:20 x(:,i) = cos(2*pi*f*t+rand(1)*2*pi); y(:,i) = cos(2*pi*f*t+rand(1)*pi/2); z(:,i) = (1+rand(1))*cos(2*pi*f*t); end subplot(3,1,1); plot(t,x,k); ylabel(x(t)) subplot(3,1,2); plot(t,y,k); ylabel(y(t)) subplot(3,1,3); plot(t,z,k); ylabel(z(t)) 由x(t)、y(t)和z(t)组成的所有样本函数的时间平均都等于零。 大量满足0≤ti≤2的，计算出来的的总体均值近似0； Y(t） 0.875 1.875 0.375 1.375 0.125 0.625 1.125 1.625 -周期平稳 Z（t） 例7-2 在后面的工作中，我们经常要用到数字调制器的模型。这些调制器的基本构建模块是函数random-binary，它产生电平值为+1或-1的二进制波形，产生的比特数以及每比特的采样数是该函数的参数。 function [x, bits] = random_binary(nbits,nsamples) % This function genrates a random binary waveform of length nbits % sampled at a rate of nsamples/bit. x = zeros(1,nbits*nsamples); bits = round(rand(1,nbits)); for m=1:nbits for n=1:nsamples index = (m-1)*nsamples + n; x(1,index) = (-1)^bits(m); end end 函数random-binary的应用 函数random-binary可以仿真多个数字调制器，例如，可用如下MATLAB语句仿真一个QPSK调制器： x=random_binary(nbits, nsamples)+i*random_binary (nbits, nsaples); 产生一个长10比特的QPSK信号，采样频率为每比特8个采样点 nbits = 10; nsamples = 8; x = random_binary(nbits,nsamples)+i*random_binary(nbits,nsamples); xd = real(x); xq = imag(x); subplot(2,1,1) stem(xd,.); grid; axis([0 80 -1.5 1.5]); xlabel(Sample Index); ylabel(xd) subplot(2,1,2) stem(xq,