第八讲随机过程功率谱及性质与计算
- 格式:ppt
- 大小:710.02 KB
- 文档页数:19
三、功率谱分析字体[大] [中] [小]周期信号的功率谱为其双边幅值频谱的平方|c n|2;非周期信号的功率谱为其幅值谱密度的平方|X(ω)|2=X(ω)X*(ω)。
随机信号属于时域无限信号,其频率、幅值和相位为随机变量。
因而,采用具有统计特性的功率谱估计进行谱分析(一)自功率谱密度及其估计各态历经随机信号的功率谱密度S x(ω)与自相关函数R x(τ)为傅里叶变换偶对,即为了方便,也可用在非负频率范围内(ω>0)定义的单边功率谱密度G x(ω)代替双边功率谱密度S x(ω),两者之间的关系为自功率谱估计可分为线性估计法与非线性估计法。
前者以快速变换为基础,应用较早,也称为经典谱分析法; 后者是与时序模型结合的一种新方法,又称为现代谱分析方法。
1. 周期图各态历经随机信号的均方值ψx2为信号能量的时域描述。
巴什瓦定理表明,信号能量的时域计算与频域计算相等,即由此定义自功率谱密度及其估计为:式中表12-45 典型信号的自相关、频谱、概率密度(续)X(ω)为测试数据x(t)的傅里叶变换,X(k)为N个数据x(n)的离散傅里叶变换,由FFT直接求出。
由于X(k)具有周期函数的性质,所以称由此获得的自功率谱估计为周期图。
自相关估计x′(r)的快速傅里叶变换可作为自功率谱估计的另一计算公式以上两种估计都是自功率谱S x(ω)的有偏估计,只是偏差大小不同。
两种估计在时域对数据或对自相关估计进行截断,相当于加窗处理,致使谱估计成为真实功率谱(或称为真功率谱)与窗谱W(ω)的卷积,即Ŝx(ω)=S x(ω)*W(ω)窗谱旁瓣的泄漏效应和卷积的作用使真功率谱的尖峰数值变化,邻近点的数值变大,造成谱估计的模糊与失真以上两种估计的方差较大; 相距2π/N的各点估计值互不相关,故数据点数N越大,这些点的估计值的随机起伏越严重。
为改善谱估计的估计质量,在增大数据点数的同时,采用平均化处理和窗处理方法减小谱估计的方差。
若随机过程均值非零,则功率谱在原点有一δ函数; 若随机过程均值非零,则功率谱在原点有一δ函数; 若随机过程均值非零若含有周期分量,则在相应的频率处有δ函数; 若含有周期分量,则在相应的频率处有δ函数; 相关性与功率谱的关系为:相关性越弱,功率相关性与功率谱的关系为:相关性越弱, 相关性与功率谱的关系为谱越宽平;相关性越强,功率谱越陡窄. 谱越宽平;相关性越强,功率谱越陡窄. 1 +∞ RX (0 = ∫ GX (ωdω 2π ∞ 16
例3,已知零均值平稳过程的谱密度为ω2 + 4 GX (ω = 4 ω +10ω2 + 9 求相关函数与方差. 解: 由因式分解ω4 + 10ω2 + 9 2α α τ 由公式: e 2 2 α +ω GX (ω = ω2 + 4 = 2 × 9 / 48 6 × 5 / 48 + ω2 + 1 ω2 + 9 1 τ 3 τ RX (τ = (9e + 5e 48 RX (0 = 7 24 17
功率谱密度算例例2 设随机过程Y(t = aX(t cos(ω0t 为常量, 其中a,ω0 为常量, (t的功率谱为为GX (ω , X 的功率谱密度. 求 Y (t 的功率谱密度. a 解: Y (ω = G {GX (ω ω0 + GX (ω + ω0 } 4 RY (t, t + τ = R[Y (tY (t + τ ] a2 = RX (τ [cosω0τ + cos(2ω0t + ω0τ ] 2 +∞ GX (ω = RX (τ e jωτ dτ ∞ 2 ∫ 1 +T RX (τ = lim RX (t, t +τ dt T→∞ 2T T ∫ 18
RX (τ GX (ω 2α /( a 2 + ω 2 返回 19。
计算机与信息技术学院验证性实验报告一、实验目的:1、学会用matlab自带函数编写程序;2、掌握功率谱估计的方法。
二、实验原理:功率谱:随机信号的功率谱反映的是随机信号的频率成分及各成分的相对强弱。
功率谱估计:基于有限的数据寻找信号、随机过程或系统的频率成分。
采用两种方法:圆周图法、自相关法(根据维纳-辛钦定理)、welch法。
(1)、平滑:是用加窗的办法对单一功率谱估计加以平滑,用于自相关法求功率谱,对自相关加窗,然后再求其傅里叶变换。
(2)、welch法是对长度为N的数据段x(n)分段时,允许每一段有部分的重叠,每一段数据用一个合适的窗函数进行平滑处理,求每段数据的DFT。
周期图法求各段功率谱估计,对各段功率谱求平均并归一化处理。
三、使用仪器计算机、matlab软件四、实验步骤1、周期图法进行功率谱估计;2、自相关法进行功率谱估计;3、welch法进行功率谱估计;4、对心电图进行功率谱估计。
五、实验过程及结果1、周期图法clcFs=1000;nfft=1024;n=0:1/Fs:1;x=sin(2*pi*100*n)+2*sin(2*pi*200*n)+3*sin(2*pi*400*n)+randn(size( n));X=fft(x,nfft);Pxx=periodogram(x);t=0:round(nfft/2-1);f=t*Fs/nfft;P=10*log10(Pxx(t+1));plot(f,P)2、自相关法clcFs=1000;nfft=1024;n=0:1/Fs:1;x=sin(2*pi*100*n)+2*sin(2*pi*200*n)+3*sin(2*pi*400*n)+randn(size( n));X=xcorr(x);Pxx=fft(X,nfft);t=0:round(nfft/2-1);f=t*Fs/nfft;P=10*log10(Pxx(t+1));plot(f,P)3、welch法clcFs=1000;nfft=1024;n=0:1/Fs:1;x=sin(2*pi*100*n)+2*sin(2*pi*200*n)+3*sin(2*pi*400*n)+randn(size( n));X=fft(x,nfft);[Pxx,F]=pwelch(x,33,32,nfft,Fs);t=0:round(nfft/2-1);f=t*Fs/nfft;P=10*log10(Pxx(t+1));plot(f,P)4、对心电图进行功率谱估计clcFs=1000;nfft=1024;n=0:1/Fs:1;x=fopen('ECG');[x,cn]=fread(x,8000,'int32'); x=x';x=x(1:2:8000);X=fft(x,nfft);Pxx=abs(X).^2/length(n);t=0:round(nfft/2-1);f=t*Fs/nfft;P=10*log10(Pxx(t+1));plot(f,P)六、实验分析通过本次对随机信号功率谱估计的实验,进一步加深了对功率谱的理解,掌握了利用MATLAB编程来绘制图形的方法,通过软件的编程与运行结果,加深了对书上理论知识的理解和掌握。