周期图及其改进方法功率谱估计报告.docxVIP

功率谱估计 姓 名：李立荣 学 号：201000800343 班 级：通信一班 指导老师：刘若伦 -AZ-.—. 一、前吕: Spectral estimation (谱估计)的Fl标是基于一个有限的数据集合描述一个信号 的功率(在频率上的)分布。功率谱估计在很多场合下都是有用的，包括对宽带 噪声湮没下的信号的检测。本论文主要研究了功率谱估计的最常用的方法一一周 期图法及其改进的两种方法，给出了各方法的matlab仿真结果及分析比较，并 尝试讨论给出各方法的优劣及原因。 二、估计质量的评价标准 1.无偏性 对于待估参数，不同的样本值就会得到不同的估计值，这样，要确定一个 估计量的好坏，就不能仅仅依据某次抽样的结果来衡量，而必须由大量抽样的结 果来衡量。我们希望估计量的均值(数学期望)应等于未知参数的真值，这就是 所谓无偏性(Unbiasedness)的要求。 A 定义:设来自总体X的一个样本，是总体参数°的一个估计量，若E0 = 0 , A 则称是的无偏估计量(Unbiased Estimator) o 一个估计量如果不是无偏的 A A 就称它是有偏估计量。eb-o称为估计量。的偏差。无偏估计的实际意义就是无 A 系统偏差，估计量是否无偏是评价估计量好坏的一个重要标准，若￡0工°,但 A lirn E9 -0 A 有 i ,则称°是°的渐近无偏估计。 2.有效性 比较两个无偏估计量优劣的一个重要标准就是观察它们哪一个取值更集屮 于待估参数的真值附近，即哪一个估计量的方差更小，这就是下面给出的有效性 (Effectiveness)概念。 定义： 2 2 A A 设与02=02(X|X2，.?X”)都是总体参数0的无偏估计，若 W.)W2)? 则称必比N更有效。在的所有无偏估计量中，如果存在一个估计量A ,它的 方差最小，则此估计量应当最好，并称此估计量久为。的最小方差无偏估计, 也称其为最有效的. 3.相合性 A 估计量°的无偏性和有效性都是在样本容量n固定的情况下讨论的。由于 A A 估计量°和样本容量n有关，我们自然希望当n很大时，一次抽样得出的的〃值 能以很大的概率充分接近被估参数 ,这就提出了相合性(Consistency)(一致 性)的要求。 定义： A A 设e=e(x、,X2，..xj是总体参数的估计量，如果对任意￡°都有 lim川 lim川 n-oo 0 — 0} = \ A A ,则称是的相合估计量(或一致估计量)。是的相合估 A 计就意味着依概率收敛于。?根据大数定律，无论总体X服从什么分布，只要 其k阶原点矩严存在,则对任意00 都有 1 i 1 i 〃一8 Ak= -±X\ 所以样本的k阶原点矩 n曰 始终是总体k阶原点矩从的相合估计。进 一步地，可以证明：只要相应的总体矩存在，矩估计必定是相合估计。特别地， 一 o f 2 u = x总是八氏 的相合估计，样本方差S「和样本的二阶中心矩S”都是 总体方差/的相合估计S和s又都是b的相合估计。 三、仿真及分析 1、周期图法 周期图法，它是把随机序列x(n)的N个观测数据视为一能量有限的序列，直 接计算x(n)的离散傅立叶变换，得X(k),然后再取其幅值的平方，并除以N,作 为序列x (n)真实功率谱的估计。 其程序实现及仿真结果如下： clear all;clc;close all; Fs=round(15*25.5) ; %采样频率383,可改变参数 n=0:l/Fs:l;%采样点数 密产牛含冇噪声的序列改变n的取值范围观察图形的变换 xn=cos(2*pi*40*n)+sqrt(0 ?10)*randn(size(n)); %sqrt (n) , 1/n代表信噪比，可改变参数 window=boxcar (length (xn) ) ; %矩形窗 nfft=1024; [Pxx, f ] =periodogram (xnz window, nf ft, Fs); 殆直接法 plot(f,10*logl0(Pxx)); title ( 1直接周期图法的功率谱估计 (N=383 SNR=0.1 ) 1 ) ;xlabel (频率 (f/Hz) * ) ; ylabel (* 功率/频率(power/f) 1); grid; 直接周期图法的功率谱估计 (n=383 SNR=-10)20 40 60 80 100 120 140 160 180 200频率(f/Hz)? i— —r !v? ? 直接周期图法的功率谱估计 (n=383 SNR=-10) 20 40 60 80 100 120 140 160 180 200 频率(f/Hz) ? i— —r !v? ? -- ? o 5 0 5 --1-1 (V」①Mod)褂暴、斉 5 2 ■ 5 3