6.窄带随机过程产生 - 随机信号分析实验报告.docVIP

计算机与信息工程学院综合性实验报告 专业：通信工程 年级/班级：2011级 第 3学年第 1学期 课程名称 随机信号分析 指导教师 杨育婕 本组成员 学号姓名 实验地点 计算机学院111 实验时间 周二 7-8节 项目名称 窄带随机过程的产生 实验类型 综合性 一、实验目的 1、基于随机过程的莱斯表达式产生窄带随机过程。 2、掌握窄带随机过程的特性，包括均值（数学期望）、方差、概率密度函数、相关函数及功率谱密度等。 3、掌握窄带随机过程的分析方法。 二、实验仪器或设备 一台计算机 2、ATLAB r2013a 三、实验内容及实验原理 基于随机过程的莱斯表达式 (3.1) 实验过程框图如下： 理想低通滤波器如图所示： 图1 理想低通滤波器 (3.2) 设白噪声的物理谱 ,则系统输出的物理谱为 (3.3) 输出的自相关函数为: (3.4) 可知输出的自相关函数是一个振荡函数。、及的均值，并绘出随机过程各个随机过程的自相关函数，功率谱密度图形。 四、MATLAB实验程序 function random(p,R,C) %产生一个p个点的随机过程 %--------------------------高斯窄带随机过程代码--------------------------% n=1:p; w=linspace(-pi,pi,p); wn=1/2*pi*R*C; [b,a]=butter(1,wn,low); %产生低通滤波器 Xt=randn(1,p); %产生p个点均值为0方差为1的随机数，即高斯白噪声 at=filter(b,a,Xt); %让高斯白噪声通过低通滤波器 y_at=at.*cos(w.*n); %产生随机过程a（t） y_bt=at.*sin(w.*n); %产生随机过程b（t） yt=y_at-y_bt; %产生一个p个点的高斯窄带随机过程 subplot(211) plot(yt) title(高斯窄带随机过程y(t)) subplot(212) pdf_ft=ksdensity(yt) ; plot(pdf_ft) title(y(t)的概率密度图) disp(均值如下) E_Xt=mean(y_at) E_at=mean(y_at) E_bt=mean(y_bt) E_ft=mean(yt) %-----------------------自相关函数代码如下--------------------------% figure(2) R_Xt=xcorr(Xt); %高斯白噪声X(t)的自相关函数 R_at=xcorr(at); %限带白噪声的自相关函数 R_y_at=xcorr(y_at); %随机过程a(t).coswt的自相关函数 R_y_bt=xcorr(y_bt); %随机过程b(t).coswt的自相关函数 R_ft=xcorr(yt); subplot(2,2,1); plot(R_Xt);title(高斯白噪声的自相关函数R_Xt); %并绘制图形 subplot(2,2,2) plot(R_at);title(限带白噪声的自相关函数R_a_bx); %并绘制图形 subplot(2,2,3) plot(R_y_bt);title(随机过程b(t)的自相关函数R_y_bt); subplot(2,2,4) plot(R_ft);title(高斯窄带随机过程y(t)的自相关函数R_yt); %------------------------功率谱密度代码如下---------------------------% figure(3) subplot(1,2,1) periodogram(Xt); title(高斯白噪声功率谱密度S_Xt); subplot(1,2,2) periodogram(at); title(限带白噪声功率谱密度S_a_bt); figure(4) subplot(3,1,1) periodogram(y_at); title(随机过程a(t).coswt概率密度概率密度S_y_at); subplot(3,1,2) periodogram(y_bt); t