MATLAB 求功率谱

功率谱

概念

  功率谱(power spectrum)是功率谱密度函数的简称。它定义为单位频带内的信号功率,表示信号功率随着频率的变化情况,即信号功率在频域的分布状况。功率谱表示了信号功率随着频率的变化关系。
  功率谱常用于功率信号(区别于能量信号)的表述于分析,其曲线一般以频率为横坐标,功率为纵坐标。周期性连续星号 x ( t ) x(t) x(t)的频谱可以表示为离散的非周期序列 x n x_n xn,它的幅度频谱的平方 ∣ x n ∣ 2 |x_n|^2 xn2所排成的序列,就称之为该信号的功率谱

定义

  功率信号 f ( t ) f(t) f(t) 在时间段 t ∈ [ − T / 2 , T / 2 ] t\in[-T/2,T/2] t[T/2,T/2] 上的平均功率可以表示为:

(1) P = ∫ − T / 2 T / 2 f ( t ) 2 d t P=\int_{-T/2}^{T/2} f(t)^2{\rm d}t \tag{1} P=T/2T/2f(t)2dt(1)
  如果 f ( t ) f(t) f(t)在时间段 t ∈ [ − T / 2 , T / 2 ] t\in[-T/2,T/2] t[T/2,T/2]上可以用 f T ( t ) f_T(t) fT(t)表示,且 f T ( t ) f_T(t) fT(t)的傅里叶变换为 F T ( ω ) = F

3种经典功率谱估计方法的MATLABA代码-功率谱代码.doc 3种MATLAB的经典谱估计方法 希望对大家有用~ 附件所有代码: 直接法: 直接法又称周期图法,它是把随机序列x的N个观测数据视为一能量有限的序列,直接计算x的离散傅立叶变换,得X,然后再取其幅值的平方,并除以N,作为序列x真实功率谱的估计。 Matlab代码: clear; Fs=1000; %采样频率 n=0:1/Fs:1; %产生含有噪声的序列 xn=cos 3*cos randn); window=boxcar); %矩形窗 nfft=1024; [Pxx,f]=periodogram; %直接法 plot); 改进的直接法: 对于直接法的功率谱估计,当数据长度N太大时,谱曲线起伏加剧,若N太小,谱的分辨率又不好,因此需要改进。 1. Bartlett法 Bartlett平均周期图的方法是将N点的有限长序列x分段周期图再平均。 Matlab代码: clear; Fs=1000; n=0:1/Fs:1; xn=cos 3*cos randn); nfft=1024; window=boxcar); %矩形窗 noverlap=0; %数据无重叠 p=0.9; %置信概率 [Pxx,Pxxc]=psd; index=0:round; k=index*Fs/nfft; plot_Pxx=10*log10); plot_Pxxc=10*log10); figure plot; pause; figure plot; 2. Welch法 Welch法对Bartlett法进行了两方面的修正,一是选择适当的窗函数w,并再周期图计算前直接加进去,加窗的优点是无论什么样的窗函数均可使谱估计非负。二是在分段时,可使各段之间有重叠,这样会使方差减小。 Matlab代码: clear; Fs=1000; n=0:1/Fs:1; xn=cos 3*cos randn); nfft=1024; window=boxcar; %矩形窗 window1=hamming; %汉明窗 window2=blackman; %blackman窗 noverlap=20; %数据无重叠 range='half'; %频率间隔为[0 Fs/2],只计算一半的频率 [Pxx,f]=pwelch; [Pxx1,f]=pwelch; [Pxx2,f]=pwelch; plot_Pxx=10*log10; plot_Pxx1=10*log10; plot_Pxx2=10*log10; figure plot; pause; figure plot; pause; figure plot; 复制代码
评论 15
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值