为了正常的体验网站,请在浏览器设置里面开启Javascript功能!

复高斯噪声下基于功率谱的频率估计与MATLAB仿真

2017-11-21 14页 doc 114KB 110阅读

用户头像

is_482581

暂无简介

举报
复高斯噪声下基于功率谱的频率估计与MATLAB仿真复高斯噪声下基于功率谱的频率估计与MATLAB仿真 摘要:频率估计是信息科学在信号处理领域的一个重要组成部分,在许多实际信号处理应用中具有非常重要的地位,且 己广泛应用于雷达、地震、生物医学、语音处理和故障诊断等领域。本文对几种常用的频率估计方法进行了回顾,简单介绍 了其中两种频率估计方法的原理;并应用最大似然、MUSIC算法对复高斯白噪声下信号的频率进行估计,并给出了计算机 matlaB仿真结果,并对各种方法估计结果进行了比较。 关键词:频率估计;最大似然;复高斯白噪声 Frequency estimation bas...
复高斯噪声下基于功率谱的频率估计与MATLAB仿真
复高斯噪声下基于功率谱的频率估计与MATLAB仿真 摘要:频率估计是信息科学在信号处理领域的一个重要组成部分,在许多实际信号处理应用中具有非常重要的地位,且 己广泛应用于雷达、地震、生物医学、语音处理和故障诊断等领域。本文对几种常用的频率估计方法进行了回顾,简单介绍 了其中两种频率估计方法的原理;并应用最大似然、MUSIC算法对复高斯白噪声下信号的频率进行估计,并给出了计算机 matlaB仿真结果,并对各种方法估计结果进行了比较。 关键词:频率估计;最大似然;复高斯白噪声 Frequency estimation based on power spectrum in WGN and MATLAB simulation Du Yang 201321010408 (School of Communication and Information Engineering, University of Electronic Science and Technology of China Chengdu 610054) Abstract: frequency estimation is a significant component in the field of signal processing in communication science, which is rather vital in the case of processing actual signal. It has been widely applied in the field of radar, earthquake, biomedicine, speech processing, and fault diagnosis n. In this paper, some common methods of frequency estimation will be reviewed at first. Then, the principals of two of these methods will be introduced in detail. Subsequently, the writer will use maximum likelihood method and MUSIC method to estimate the frequency of signal that is in the circumstance of complex white Gauss noise and receive some matlaB simulation results. Lastly, we will compare the consequences of different ways. Keywords: frequency estimation; maximum likelihood; WGN 1 引言 1.1 研究背景 频率是参量估计中的一个重要物理量。信号频率估计是指通过对信号采样值的计算和变换,估计出淹 没于噪声中的信号的频率。要高精度和高速度测量信号频率必须要有好的频率估计方法,在通信对抗侦察 中,一种好的频率估计方法,在电子对抗中能占据更大优势。信号频率估计也是信号参数估计的一个重要 物理量,参数估计是信号检测和信号征的重点。因此对信号的频率估计以及各种适用于应用的算法 研究受到普遍关注和重视,是当前信号处理技术中一个十分活跃的课题。频率估计方法繁多,本文就针对 [1]复白噪声下的信号,采用2种方法估计其频率,并得出结论,对实际应用提供一定帮助。 1.2 频率估计的发展 频率估计常常用“谱分析”的概念来表示,在傅里叶分析中,其谱表示谐波成分的幅度随频率的变化 情况,因此是“频谱”;对于周期信号,其能量是无限的,常研究功率在频域上的分布,即“功率谱”。因 此,有时又将频率估计技术称为谱估计技术。 1822年,法国工程师傅里叶提出了谐波分析的理论,此理论奠定了信号分析和功率谱估计的理论基础;19世纪末,Schuster提出了周期图的概念,至今仍被沿用;1958年,Blackman和Tukey提出了自相关谱估计,简称为BT法,它利用有限长数据,估计自相关函数,再对该自相关函数做傅里叶变换,即得谱估计,以上方法称为经典法。 上世纪60年代末70年代初,学者们提出了许多旨在提高频率分辨率的谱估计,这些方法被称为现代法,分为两大类,一类是参数模型法,一类是非参模型法。 参数模型法主要有自回归(AR)模型、滑动平均( MA)模型和自回归滑动平均(ARMA)模型等,这些方法具有较好的频率分别能力,运算速度也较快,但是,该法的性能受参数的选取、初相位等因素的影响,因而其广泛性受到限制。 非参数模型法主要有最大似然法、多信号分类法(MUSIC)和子空间旋转不变法(ESPRIT)等。 最大似然法频率估计是通过求解似然函数的最大值来进行频率估计,此法对信噪比门限要求较低,也是最为准确的频率估计方法,而且,此法不需要事先对噪声进行学习,故应用比较简便。便L} - b}。但是,因为在求解最大似然函数的最大值时要求解一个高度非线性化的方程,其计算量非常巨大,有时甚至是无法求解的,故此法不适用实时计算。 多信号分类法(MUSIC法)在20世纪70年代末实现了向现代超分辨率频率估计技术的飞跃,但是它们需要在整个频域内进行谱峰搜索,实时性差. Roy, Paulraj和Kailath于1986年提出的ESPRIT算法,其含义就是利用旋转不变子空间估计信号参数,现己成为现代信号处理中一种主要方法,并得到了广泛的应用。ESPRIT算法估计信号参数时要求由采样数据构成两个或两个以上相同的子矩阵,利用了协方差矩阵的信号子空间的旋转不变性,因为它不需要进 [2]行谱峰搜索,因而计算量比较小,但是ESPRIT算法的估计误差通常大于MUSIC算法。 本文将应用最大似然、MUSIC算法对复高斯白噪声下信号的频率进行估计,并给出了计算机matlaB仿真结果,并对各种方法估计结果进行了比较。 2 基本理论Equation Chapter 2 Section 1 2.1 频率估计的数学模型 淹没于高斯白噪声中的正弦信号的频率估计问题是频谱分析的重要内容,正弦信号频率估计是指通过对信号采样值的计算和变换,估计出淹没于噪声中的信号的频率的过程。频率估计不仅在理论上,而且在通信、雷达、声纳等领域的实际应用中,都有着非常重要的研究价值。频率估计目前国内外己经提出了不少方法,主要分为时域、频域及时频分析算法等。下面讨论复信号和实信号模型,其中复信号模型将作为 [3]后续章节的信号参考。 [4]2.1.1 复信号模型 在数学上,人们把信号描述为一个或若干个自变量的函数形式,如指数信号、正余弦信号、复指数信号、高斯信号等。通常,复指数信号可表示为 stftAe, (2.1) ,, 其中,A为信号的振幅,指数因子 。在现实中人们可以siiRR,,,,,,,,,,,1,, 观察到正弦和余弦信号,但是并不能真正观察到复指数信号,引入复指数信号一方面是为了模型在数学上 易于处理,另一方面是复指数信号的数学表示较为广泛。当=0,即s为实数时,复指数信号就是一般的, 指数信号;当=0,即s=0。时,复指数信号是直流信号。 ,,= 假设在观察时间0标准
差为1的高斯白噪声,对两个randn函数进行复组合后乘以系数m,就组成了均值为0,标准差为m的复 2,=m高斯白噪声,即,改变m相当于改变了信噪比。 s 利用周期图法估计频率和功率谱密度的函数:[Pxx ,f ] = periodogram(x,window, nfft, fs, `range') , x为输人数据向量,nfft是FFT的长度,fs是取样频率,向量window为用于计算输人信号x的改进周期图的窗函数系数,默认值为矩形窗(boxcar),`range’若取值为twosided则计算双边功率谱,若为onesided则计算单边功率谱,此函数返回功率谱密度向量Pxx和归一化频率向量f , Pxx:和f有相同的长度,等于nfft。 下面列出利用周期图法估计频率和功率谱密度的函数的程序 %最大似然法利用周期图估计频率:mle. m randn('state',0);n=0:99;%取样点数为100 m=0.05;%设置噪声的标准差 Fs=1000;%设置采样频率为1000 x=exp(i*2*pi*0.48*n)+exp(i*(2*pi*0.50*n+pi/4))+m*(randn(size(n))+i*randn(size(n)));%定义输入的数据向量,加号前面的是信号,后面的是复高斯信号 [pxx,f]=periodogram(x,[],512,Fs,'twosided');%利用周期图进行最大似然估计,空向量[]为默认的矩形窗 P=10*log10(pxx);%取出pxx,并且对数化 plot(f,P);%图形显示 title('最大似然法利用周期图估计频率N=99,fl=0.48,f2=0.5,m=0.05');%加标题 其仿真图形为 3.2 频率的MUSIC 算法的仿真实现 MATLAB信号处理工具箱提供一种功率谱估计函数pmusic,可以实现多信号分类法。该函数将数据自相关矩阵看成由信号自相关矩阵和噪声自相关矩阵两部分组成,即数据自相关矩阵R包含有两个子空间信息:信号子空间和噪声子空间。这样,矩阵特征值向量也可分为两个子空间:信号子空间和噪声子空间。 为了求得功率谱估计,函数pmusic计算信号子空间和噪声子空间的特征值向量函数,使得在周期信号频率处函数值最大,功率谱估计出现峰值,而在其他频率处函数值最小。 其调用格式为:[Pxx,f]=pmusic(x,p,Nfft,Fs, `corr', `range') 式中,x为输入信号的向量或矩阵;输人参数p有两种情况:一种情况是一个整数标量,表示信号子空间维数标识信号子空间的维数为P,即单独利用信号子空间来估计频率。另一种情况是一个两个元素向量。P(l)指定信号子空间的最大维数。P的第二个元素p(2)和信号相关矩阵最小特征值的乘积表示一个阈值。如果特征值在阈值之下,则分配到噪声子空间,否则分配到信号子空间。‘corr’指定输入为一个相关矩阵,而不是信号数据矩阵。‘range’若为‘whole’则在频率区间[0,fs」上计算伪谱,这是x为复数时的默认值,若为‘half’则在指定的频率范围内计算伪谱,这是x为实数时的默认值。函数返回信号的伪谱S和频率向量f 。 下面列出以MUSIC算法估计信号频率的程序 %以MUSIC算法估计信号频率:music.m randn('state',1); fs=1000;n=0:99;M=2;m=0.01; %设置噪声的标准差,设置采样频率为1000,取样点数为100,M为信号空间维数,m为噪声方差 x=exp(i*2*pi*0.48*n)+exp(i*(2*pi*0.50*n+pi/4))+m*(randn(size(n))+i*randn(size(n))); %定义输入的数据向量,加号前面的是信号,后面的是复高斯信号 [X,R]=corrmtx(x,M,'mod'); %返回M+1阶的x的自相关矩阵 pmusic(R,M,'whole',512,fs,'corr');%图形显示 图(3- 2)以MUSIC算法估计信号频率 图(3- 3)以MUSIC算法估计信号频率 4 结果分析 通过对上面两幅图片的分析对比,当正弦波的频率能用周期图加以分辨时,则根据周期图谱峰的位置来估计正弦波的频率,这是一种最简便的方法,同时可避免最大似然估计遇到的困难。但在在加入一定噪声以后,正弦波频率很接近的情况下,以周期图为基础的最大似然估计方法就受采样点数的限制,周期图就无法分辨。因而采用频率分辨率更高的现代谱估计方法进行估计。可以看出多信号分类法是一种比较理想的方法。 参考文献 [1] 王立乾,赵国庆,郑文秀(基于现代谱估计的PSK信号频率估计方法[J](现代电子技术,2003(23): 44-47( [2] 张晓菲,刘振兴,陈栋(频率估计综述[J](信息技术,2011(11):58-62( [3] 赵媛媛(信号模型中若干参数估计问题的研究[D](合肥:中国科学技术大学,2011( [4] 邹建彬,张炜(复高斯白噪声下基于多种方法的频率估计[J](现代电子技术,2006(17):1-3+6( [5] 杨萃(噪声环境下频率估计算法研究[D](广州:华南理工大学,2010( [6] 汪琼,陈万培(基于MUSIC的含噪声信号频率估计[J](国外电子测量技术,2010(02):20-22( [7] 陈明(基于解相干信号的空间谱估计算法的研究[D](成都:电子科技大学,2011(
/
本文档为【复高斯噪声下基于功率谱的频率估计与MATLAB仿真】,请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑, 图片更改请在作品中右键图片并更换,文字修改请直接点击文字进行修改,也可以新增和删除文档中的内容。
[版权声明] 本站所有资料为用户分享产生,若发现您的权利被侵害,请联系客服邮件isharekefu@iask.cn,我们尽快处理。 本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用。 网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。

历史搜索

    清空历史搜索