实验随机信号重点分析.doc

实验报告 通信信号分析与处理 专 业 学 号 姓 名 王溪岩 日 期 通信信号分析与处理实验 指导书 姓 名： 王溪岩 院 （系）：电子与信息工程学院 专 业： 通信工程 学 号：j130510401 实验日期： 2016年1月 5日 成 绩： 实验一：平稳随机过程的时域特征分析 实验目的： 1、加深理解平稳随机过程数字特征的概念 2、分析平稳随机过程时域数字特征的特点 3、掌握平稳随机序列时域特征分析方法 4、掌握MATLAB分析时域特征对应的函数 实验内容： 产生高斯分布和均匀分布的随机过程，其均值为各自学号后两位的1/4，方差为各自学号的后两位。产生p=（学号后两位+1）/50的二项分布。产生信号长度分别为1000，5000，10000。分析产生的各种分布对应的均值、方差和自相关函数。 实验要求： 1、给出所求时域特征的原理公式。 2、给出所用MATLAB函数涉及的参数和使用方法。 3、所用实验均需要提供对应的完整的源程序和对应的实验结果及其必要分析。 4、不得互相抄袭。 1、实验过程与仿真 该实验主要包括函数确定、参数选择、仿真和分析几个部分，具体仿真程序和结果分析如下： 二项分布随机过程 信号产生 1）高斯分布随机过程： n=input(n=); x=0.25; o=1; m=1; R=normrnd(x,o,m,n); subplot(3,1,1); plot(R) R_a=xcorr(R); subplot(3,1,2); plot(R_a) Pf=abs(fft(R_a,2048)); subplot(3,1,3); plot(Pf) (n输入1000，5000，10000 运行结果： 结果分析：由图可看出，高斯随机分布的均值几乎在一条直线上，可看作 为恒定值，与时间无关；自相关函数是仅与时间间隔T有关的 函数，高斯随机分布为平稳过程；当n=1000时，值返回到0时 的值，此时的自相关系数最大，表明自己与本身的自相关程度最高。 2）均匀分布： m=1; n=input(n=); a=0;b=0.5; R=unifrnd(a,b,m,n); R_a=xcorr(R); subplot(3,1,1); plot(R);title(均匀随机分布); Pf=abs(fft(R_a,10000)); subplot(3,1,2); plot(R_a);title(自相关); subplot(3,1,3); plot(Pf);title(功率); 结果分析：自相关系数在时间间隔为1的时候最高。 3）二项分布 n=input(n=); m=1; p=0.02; N=1; R=binornd(N,p,m,n); subplot(3,1,1) plot(R); R_a=xcorr(R); subplot(3,1,2) plot(R_a) Pf=abs(fft(R_a,10000)); subplot(3,1,3); plot(Pf)运行结果： 结果分析：二项随机分布的值在0.5左右震荡，均值为0.5，与时间无关； 自相关函数为仅与时间间隔t有关的函数，该过程为平稳过程。 T=0时，自相关系数最高。 实验涉及函数及用法 1）normrnd函数：生成服从正态分布的随机数 语法：R＝normrnd(MU,SIGMA,m,n) 说明：R＝normrnd(MU,SIGMA,m,n)： 生成m×n形式的正态分布 的随机数矩阵。MU参数代表均值，DELTA参数代表标准差 参数：MU=2/4=0.5;SIGMA=2;m=1,n=1000。 2）unifrnd函数：是生成(连续)均匀分布的随机数 语法：R = unifrnd(A,B) 生成被A和B指定上下端点[A,B]的连续均匀分布的 随机数组R。 R = unifrnd(A,B,m,n,...) 或 R =unifrnd(A,B,[m,n,...] 返回 m*n*...数组。  参数：学号=2，E=2/4=0.5,即A,B区间取[0,1], 随机数长度为5000，