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

Durbin-Watson测试Durbin-Watson独立性测试可以通过计算如下Durbin参数来完成：其中X[n]是一个零均值（zero-mean）随机变量[4]。下面会讲到当D值在2附近时，X[n]和X[n-1]相关性很小。*第30页，共45页，星期日，2025年，2月5日Durbin-Watson测试为了说明Durbin-Watson测试的性质，假设和相关并且是一个遍历性过程。为了简化记号，假设N足够大，使得N=N-1，式（7-19）可写成：（7-20）由假定X[n]和X[n-1]相关，可令其中X和Z不相关，ρ是X和Y的相关系数。注意，X、Y和Z有相同的方差，用б表示。将式（7-21）代入式（7-20）得：*第31页，共45页，星期日，2025年，2月5日X和Z不相关且均值为零，所以中间项等于零。由于和方差相同，所以有因为，Durbin参数在0与4之间，且有ρ=0时D=2。D2表示正相关，而D0表示ρ值为负。以下的MATLAB函数计算Durbin参数值：*第32页，共45页，星期日，2025年，2月5日*第1页，共45页，星期日，2025年，2月5日随机信号的产生与处理信道噪声、干扰和衰落等随机波形级精确地仿真这些系统，随机影响建立准确的模型其基本构建模块是随机数发生器。如何产生采样后的随机波形（信号、干扰和噪声等），以用于仿真。在仿真环境下，所有的随机过程必须用随机变量序列来表示。随机数发生器产生的是“伪随机序列”，因为尽管它们是确定的，但在具体应用中会呈现随机性。精度随应用不同：产生一个波形来表示锁相环鉴相器输入端的噪声，对一个输入SNR为50dB的系统建立噪声波形模型所要求的精度要比输入端SNR为8dB的系统高得多。*第2页，共45页，星期日，2025年，2月5日问题在（0，1）上产生均匀分布且不相关的随机数将不相关且均匀分布的随机数映射成不相关的，具有任意（需要的）概率密度函数（pdf）的随机产生不相关的、具有高斯型pdf的随机数产生相关的、具有高斯型pdf的随机数产生相关的、具有任意（需要的）pdf的随机数伪噪声（PN）序列的产生以及几种应用于随机数序列的计算方法。*第3页，共45页，星期日，2025年，2月5日平稳与遍历性过程信号、噪声和干扰而产生的样本函数通常假设为各态历经的。原因：依次处理系统中波形的时域样本，系统中的每个点上只有一个波形（样本函数）。假设仿真所处理的波形是其内在统计模型定义的总体（Emsemble）中的一个典型成员，这时各种统量如各阶矩、信噪比和误比特率，就可以当作时间平均来计算。仿真计算得到的时间平均等于总体均值。隐含的假设条件是对应的随机过程为遍历性的。*第4页，共45页，星期日，2025年，2月5日例7-1假设随机过程的样本函数表达式为是对应随机试验的样本空间中的一个输出，每一个映射为一个相位。再假定对应的随机试验就是从均匀随机数发生器输出端抽取一个数，抽取得到的结果，这里ui在（0，1）区间上均匀分布。然后映射成相位。当幅度A和频率f固定，的值便决定了波形。K-?*第5页，共45页，星期日，2025年，2月5日f=1; %frequencyofsinusoidfs=100; %samplingfrequencyt=(0:200)/fs;%timevectorfori=1:20x(:,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);endsubplot(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))*第6页，共45页，星期日，2025年，2月5日*第7页，共45页，星期日，2025年，2月5日由x(t)、y(t)和z(t)组成的所有样本函数的时间平均都等于零。大量满足0≤ti≤2的，计算出来的的总体均值近似0；Y(t）0.8751