统计信号分析及处理实验报告.docVIP

一、实验目的和要求 了解随机过程特征估计的基本概念和方法，学会运用Matlab 函数对随机过程进行特征 估计,并且通过实验了解不同估计方法所估计出来结果之间的差异。 二、实验环境、内容和方法 内容：实验原理 设随机序列 X(n)、Y(n)为各态历经过程，样本分别为x(n)、y(n)（n=0,1,....N-1）。 均值的估计 2、方差的估计 方差估计有两种情况，如果均值 mx 已知，则 如果均值未知，那么 相关函数估计 4、功率谱估计 功率谱的估计有几种方法， （1）自相关法 先求相关函数的估计 然后对估计的相关函数做傅立叶变换， （2）周期图法 先对序列 x(n)做傅立叶变换， 则功率谱估计为 周期图法是一种非参数谱估计方法，另外还有一种修正的周期图方法，也叫Welch 法， MATLAB 有周期图和Welch 法的谱估计函数。 （3）现代谱估计技术 现代谱估计主要有参数谱估计和子空间谱估计。参数谱估计法是假定待估计功率谱的信 号是白噪声驱动线性系统的输出，常用的方法有基于最大墒估计的伯格算法和Yuler-Walk 自回归(AR)方法，这些方法是估计线性系统的参数，通常会得到比经典谱估计方法更好的估 计。子空间法也称为高分辨率谱估计或超分辨率谱估计，常用的方法有MUSIC 法和特征矢 量法，这些方法是根据相关矩阵的特征分析或特征分解得到对信号的频率分量的估计，特别 适合于线谱（即正弦信号）的估计，是低信噪比环境下检测正弦信号的有效方法。 三、实验过程描述 第一题 1）、产生一组均值为1，方差为4 的正态分布的随机序列（1000 个样本），估计该序列的均值 与方差。 画出图像的代码如下： n=1:1:1000; w(n)=1+randn(1,1000)*sqrt(4); plot(n,w(n)); mean(w(n)) ans = 1.0711 var(w(n)) ans = 4.0549 3）运行以上代码，得到如下图: 第二题 1）按如下模型产生一组随机序列： x(n)=0.8x(n-1)+w(n)，其中w(n)为均值为1，方差为4 的正态分布白噪声序列。估计过程的自相关函数与功率谱。 编写代码如下: for n=1:1:1000; if n2; x(n)=rand(1); else; x(n)=0.8*x(n-1)+randn(1)*sqrt(4)+1; end; end; n=1:1:1000; figure(1); plot(n,(x(n)));title(随机数列x（n）); Fs=1000; window=boxcar(length(x(n))); [Pxx,w]=periodogram(x(n),window,n,Fs); c=xcorr(x(n)); figure(2); subplot(1,2,1);plot(c);title(估计过程的自相关函数) subplot(1,2,2);plot(w,10*log10(Pxx));title(估计过程的功率谱) 运行以上代码，得到如下图: 第三题 1）设信号为x(n)=sin(2πf1n)+2cos(2πf2n)+w(n)，n=1,2,....,N，其中f1=0.05，f2=0.12，w(n)为 正态白噪声，试在N=356 和1024 点时，分别产生随机序列x(n)、画出x(n)的波形并估计x(n) 的相关函数和功率谱。 2）编写代码如下: f1=0.05; f2=0.12; n=1:365; x(n)=sin(2*pi*f1*n)+2*cos(2*pi*f2*n)+randn(1); Fs=1000; window=boxcar(length(x(n))); [Pxx,w]=periodogram(x(n),window,n,Fs); subplot(3,2,1);plot(n,x(n));title(N=365点时产生的随机序列x（n）) c=xcorr(x(n)) subplot(3,2,3);plot(c);title(N=365点时产生的随机序列x（n）的自相关函数) subplot(3,2,5);plot(w,10*log10(Pxx));title(N=365点时产生的随机序列x（n）的功率谱) n=1:1024; x(n)=sin(2*pi*f1*n)+2*cos(2*pi*f2*n)+randn(1); subplot(3,2,2);plot(n,x(n));title(N=1024点时产生的随机序列x（n）) c=xcorr(x(n)); subplot(3,2,4);plot(c);title(N=1024点时产生的随机序列x（n）的自相关函数) wi