随机信号处理实验报告讲诉

随机信号处理实验报告目录一、实验要求： (3)二、实验原理： (3)2.1 随机信号的分析方法 (3)2.2 随机过程的频谱 (3)2.3 随机过程的相关函数和功率谱 (4)（1）随机信号的相关函数： (4)（2）随机信号的功率谱 (4)三、实验步骤与分析 (5)3.1实验方案 (5)3.2实验步骤与分析 (5)任务一：（s1 变量）求噪声下正弦信号的振幅和频率 (5)任务二：（s1 变量）求噪声下正弦信号的相位 (8)任务三：（s1 变量）求信号自相关函数和功率谱 (11)任务四：（s变量）求噪声下信号的振幅和频率 (14)任务五：（s变量）求信号的自相关函数和功率谱 (17)3.3实验结果与误差分析 (19)（1）实验结果 (19)（2）结果验证 (19)（3）误差分析 (21)四、实验总结和感悟 (22)1、实验总结 (22)2、实验感悟 (23)五、附低通滤波器的Matlab程序 (23)一、实验要求：（学号末尾3，7）两个数据文件，第一个文件数据中只包含一个正弦波，通过MA TLAB 仿真计算信号频谱和功率谱来估计该信号的幅度，功率，频率和相位？对第二个文件数据估计其中正弦波的幅度，功率和频率？写出报告，包含理论分析，仿真程序及说明，误差精度分析等。

第一文件调用格式load FileDat01_1 s1，数据在变量s1中；第二文件调用格式load FileDat01_2 s ，数据在变量s 中。

二、实验原理：2.1 随机信号的分析方法在信号与系统中，我们把信号分为确知信号和随机信号。

其中随机信号无确定的变化规律，需要用统计特新进行分析。

这里我们引入随机过程的概念，所谓随机过程就是随机变量的集合，每个随机变量都是随机过程的一个取样序列。

随机过程的统计特性一般采用随机过程的分布函数和概率密度来描述，他们能够对随机过程作完整的描述。

但由于在实践中难以求得，在工程技术中，一般采用描述随机过程的主要平均统计特性的几个函数，包括均值、方差、相关函数、频谱及功率谱密度等来描述它们。

2.2 随机过程的频谱信号频谱分析是采用傅立叶变换将时域信号x(t)变换为频域信号X(f)，从而帮助人们从另一个角度来了解信号的特征。

时域信号x(t)的傅氏变换为：()()2j ft X f x t e dt π+∞--∞=⎰信号的时域描述只能反映信号的幅值随时间的变化情况，除只有一个频率分量的简谐波外,一般很难明确揭示信号的频率组成和各频率分量的大小。

信号的频谱X(f)代表了信号在不同频率分量处信号成分的大小，它能够提供比时域信号波形更直观，丰富的信息。

在实际的控制系统中能够得到的是连续信号x(t)的离散采样值x(nT)，因此需要利用离散信号x(nT)来计算信号x(t)的频谱。

有限长离散信号x(n)，n=0，1，…，N-1的DFT 定义为：()()10N knN n X k x n W -==∑其中0,1,2.......1,kN =- 2jNN W eπ-=2.3 随机过程的相关函数和功率谱（1）随机信号的相关函数：信号的相关性是指客观事物变化量之间的相依关系。

对于平稳随机过程X(t)和Y(t)在两个不同时刻t 和t+τ的起伏值的关联程度，可以用相关函数表示。

在离散情况下，信号x(n)和y(n)的相关函数定义为：∑∑-=-+=101N txy N /)t (y )t (x ),t (N R τττ τ,t=0,1，2,……N-1随机信号的自相关函数表示波形自身不同时刻的相似程度。

与波形分析、频谱分析相比，它具有能够在强噪声干扰情况下准确地识别信号周期的特点。

一般来说，信号与噪声在时域内有明显不同，信号前后是有关联的，存在相关性；而噪声在不同时刻基本上不存在关联，即不存在相关性．利用这种相关性原理，已成为从强噪声中提取弱信号的重要手段。

这种技术的理论基础是信息论和随机过程理论，这种检测方法被称为相关检测。

（2）随机信号的功率谱随机信号的功率谱密度是随机信号的各个样本在单位频带内的频谱分量消耗在一欧姆电阻上的平均功率之统计均值，是从频域描述随机信号的平均统计参量，表示X(t)的平均功率在频域上的分布。

它只反映随机信号的振幅信息，而没有反映相位信息。

随机过程的功率谱密度为：]2|)(|lim [)(2TX E x G Ti T ω∞→= －∞＜ω＜＋∞随机信号的平均功率就是随机信号的均方值，功率谱密度曲线下的总面积（即随机信号的全部功率）等于随机信号的均方值。

随机信号的功率谱与它的自相关函数构成一对傅里叶变换对。

三、 实验步骤与分析本实验利用Matlab 软件编程来实现数据文件中波形的仿真与分析，最后通过滤波器还原正弦信号，与结果进行比较。

3.1 实验方案3.2 实验步骤与分析任务一：（s1 变量）求噪声下正弦信号的振幅和频率 （1）原理：采用傅立叶变换将时域信号x(t)变换为频域信号X(f)，并作出幅频曲线进行分析，离散信号x(n)，n=0，1，…，N-1的DFT 公式如下：()()10N knNn X k x n W -==∑其中0,1,2.......1,kN =- 2j NN W eπ-=在Matlab 的编程实现时，运用的是快速算法傅里叶算法FFT ，它是DFT 的快速算法。

因为给定的数据文件中采样点数N=4096，所以取采样频率fs=4096Hz 。

（2）Matlab仿真结果及分析图1 随机信号的时域图形图2 随机信号的频域图形由时域图形可知，正弦信号被噪声“淹没”了，所以时域上看不出任何信号的特征，进行傅里叶变换，频域特征如图2所示。

已知采样频率fs=4096Hz，所以Nyquist频率为fs/2=2048Hz，傅里叶变换的数据具有对称性，整个频谱以Nyquist频率为对称轴，所以频谱分析的时候只要截取0~2048Hz范围内的频谱进行分析。

由频谱曲线可知，信号在82Hz处有一个峰值，大小为4021，所以可以得出：信号频率：=82f Hz信号绝对幅度：2=⨯=A4021 1.96N（3）附Matlab程序及说明clc; %清空clear all; %清除所有变量close all; %关闭所有窗口load('C:\Users\caolili\Desktop\FileDat01_1.mat')fs=4096; %设定采样频率N=4096; %采样点数n=0:N-1;t=n/fs; %采样时间间隔subplot(211); %两行一列第一幅图plot (n,s1); %画出时域波形xlabel('时间t(1/4096s)');ylabel('信号s1');title('原信号时域波形');grid;%进行FFT变换并做频谱图y=fft(s1,N); %进行fft变换mag=abs(y); %求幅值f=(0:length(y)-1)'*fs/length(y); %进行对应的频率转换subplot(212); %两行一列第二幅图plot(f,mag); %作频谱图xlabel('频率(Hz)');ylabel('幅值');title('信号的幅频谱图N=4096');grid;任务二：（s1 变量）求噪声下正弦信号的相位 （1）原理设观测数据为：1，,...2,1,0),()2cos()(0-=++=N n n w n f A n x φπ式中，w （n ）为已知方差2σ的高斯白噪声，正弦信号的幅度A 和频率f 0为已知。

一种估计φ的估计量为⎭⎬⎫⎩⎨⎧⎥⎦⎤⎢⎣⎡+-=∑∑-=-=∧101000)2sin()()2cos()(ln Im N n N n n f n x j n f n x ππφ定义信噪比为222σA SNR =。

由上面已经求的信号幅度A=1.96，f=82Hz ，N=4096。

○1产生服从特定概率分布的观测数据x （n ）； ○2利用估计算法计算估计量∧φ； ○3上述过程重复M 次，产生M 个∧φ的实现 ○4利用∑=∧∧=∧M ii M11θμθ确定估计量的均值。

○5利用∧=∧∧∑∧∧-=MiM122)(1θθμθσ确定估计量的方差。

○6利用直方图来确定PDF ：首先计算落入某指定区间的次数，然后再除以总的实现次数得到概率，再除以区间长度得到PDF 估计。

（2）Matlab仿真结果及分析根据利用matlab仿真得到的PDF估计如图3所示。

图3 随机信号的相位的统计特性由上面的仿真结果可知初相位约在0时具有最大概率，所以：θ=（3）附Matlab程序及说明clc; %清空clear all; %清除所有变量close all; %关闭所有窗口load('C:\Users\caolili\Desktop\FileDat01_1.mat')N=4096;A=1.96; %A为正弦信号幅值，f为其频率f=82;n=0:N-1;s0=300; %设置循环次数for m=1:s0y1=0;任务三：（s1 变量）求信号自相关函数和功率谱 （1）原理：对于噪声中信号的功率谱分析，有传统方法和现代建模方式。

本次实验中主要采用传统谱估计的自相关法，又称为间接法或BT 法。

具体步骤是先由)(n x N 估计出自相关函数)(ˆm r ，然后对)(ˆm r 求傅里叶变换得到)(n x N 的功率谱，记之为)(ˆw P BT，并以此作为对)(w P 的估计，即1,)(ˆ)(ˆ-≤=--=∑N M e m r w P jwmMMm BT。

（2）Matlab 仿真结果及分析在Matlab 中主要用C (Xn) = xcorr(xn,'unbiased')函数来计算Xn 的自相关函数，然后对其进行傅里叶变换，便得到它的功率谱。

图4 随机信号的自相关函数图4 随机信号的功率谱原信号在时域上时被噪声淹没，经过自相关后可以看出信号为正弦信号，并且由自相关函数图象可知：混合信号的平均功率：2[()](0) 2.934E X n R ==正弦信号的平均功率：()2220lim cos 1.922TT T AP A w t dt θ-→∞=+==⎰在功率谱图像上，极值点坐标为（82,27.72），正好对应正弦信号的频率为82Hz ，与上面的频谱分析一致。

对功率谱密度曲线积分也可求出信号的平均功率。

（3）附Matlab程序及说明%自相关函数clc; %清空clear all; %清除所有变量close all; %关闭所有窗口load('C:\Users\caolili\Desktop\FileDat01_1.mat')fs=4096; %设定采样频率N=4096; %采样点数n=0:N-1;t=n/fs;Lag=300; %延迟样点数[c,lags]=xcorr(s1,Lag,'unbiased'); %对原始信号进行无偏自相关估计subplot(1,2,1);plot(n,s1); %绘制原信号的时域波形xlabel('时间t(1/4096s)');ylabel('信号s1');title('带噪声的信号波形');grid on;subplot(1,2,2);plot(lags/fs,c); %绘制自相关函数图象xlabel('时间t');ylabel('Rx(t)');title('带噪声的信号自相关函数');grid on;%功率谱clc; %清空clear all; %清除所有变量close all; %关闭所有窗口load('C:\Users\caolili\Desktop\FileDat01_1.mat')fs=4096; %设定采样频率N=4096; %采样点数n=0:N-1;t=n/fs;Lag=300; %延迟样点数[c,lags]=xcorr(s1,Lag,'unbiased'); %求信号的自相关函数fy=fft(c,N); %对自相关函数做FFT变换t1=0:round(N/2-1);a=t1*fs/N;P=10*log10(fy(t1+1)); %纵坐标为相对功率谱密度，单位dB/Hz figure(gcf);plot(a,P);ylabel('功率谱密度dBw/Hz');title('信号的功率谱');grid;任务四：（s变量）求噪声下信号的振幅和频率（1）原理：同“任务一”的原理相同（2）Matlab仿真结果及分析图5 随机信号的时域波形（s变量）图6 随机信号的频谱（s变量）对于频谱图的局部放大如下图所示：图7 随机信号的局部放大的频谱（s变量）由频域分析可以发现，信号的频谱图上有两个峰值，由Matlab 计算得两个极点分别为：（82,4053） （86,8150）所以信号由两个频率相近的正弦信号组成，根据图形分析，有用信号应该是两个正弦信号相加，形如：()1122sin 2cos 2s s n n S n A f A f f f ππ⎛⎫⎛⎫=⨯+⨯ ⎪ ⎪⎝⎭⎝⎭根据“任务一”的计算方式，可以得出： 有用信号的频率：182f Hz =286f Hz =有用信号的绝对振幅：124053 1.98A N=⨯=228153 3.98A N=⨯= 对于此处采样点数N 和采样频率fs 的确定要满足频率分辨率的要求，即：max min12f f f f f N-∆=-=所以，要能有效的区分频率轴上的两个频率点f1和f2，有效数据长度必须满足以下关系式：122sf f f N<-所以此处取采样频率4096s f Hz =，采样点数4096N =，满足要求。