离散时间随机信号讲课资料.pptVIP

第八章 离散时间随机信号;;;;;;;例8.1计算零均值过程的相关序列和协方差序列。（Matlab上机实验）程序代码如下 :;subplot(223) [r1,lags]=xcorr(x,length(x)/2,unbiased); plot(lags,r1); xlabel(n);ylabel(R_x(n)); title(随机过程的自相关序列); grid; subplot(224) [r2,lags]=xcov(x,length(x)/2,unbiased); plot(lags,r2); xlabel(n);ylabel(Cov_x(n)); title(随机过程的协方差序列); grid;;例8.2 已知序列x(n)为一长度为N＝104的高斯随机噪声，计算它的样本均值、样本方差。（Matlab上机实验）程序代码如下：;;;;;;;;;;;;【例8.3】随机信号经过线性非移变系统后的均值（Matlab上机实验）。设滤波器的单位脉冲响应h(n)是由下式所给出的离散衰减余弦信号的前61个值所决定： h(n)=(0.95)ncos(0.1nπ),输入信号为均值为3、方差为1的Gauss白噪声序列x(n)。该序列通过滤波器后，得到输出序列y(n)。验证式(8.4.22)。 代码如下： %% clc,clear N=1e5; sigma=1; x=3+sigma*randn(1,N); n=0:60; h=0.95.^n.*cos(0.1*n*pi); y=filter(h,1,x); disp(输出信号的均值m_y); mean(y) disp(式(8.4.22)计算得到的m_y); mean(x)*sum(h);;;;;【例8.4】线性非移变系统单位脉冲响应的估计。（Matlab上机实验） 滤波器的单位脉冲响应h(n)由下式所给出的离散衰减余弦信号的前61个值所决定：h(n)=(0.95)ncos(0.1nπ), 为利用互相关序列求该滤波器的单位脉冲响应，首先产生1000个均值为零、方差为1的Gauss白噪声序列w(n)。该序列通过滤波器后，得到输出序列y(n)。然后通过下式计算互相关函数的估计;代码如下： %% clc,clear N=1e3; sigma_sq=1; w=sigma_sq*randn(1,N); n=0:60; h=0.95.^n.*cos(0.1*n*pi); y=filter(h,1,w); R_wy=xcorr(w,y,60,unbiased); h_est=R_wy/sigma_sq; figure; plot(h,.);hold on; plot(h_est(60:120));;;;;;本章小结