信号处理的技术chapter5.pptVIP

* * * * * * * * * * * * * * * * * * * * * * §5.4 对周期图作平滑后的功率谱估计 定义：韦尔奇结合巴特利特法与窗口法进行了修正，并证明了对长度N的实平稳随机序列x(n)以2:1覆盖分为L段时，不仅满足一致估计条件，且谱估计的方差最小。 2:1覆盖分段： 5.4.3 韦尔奇（Welch）修正法 ...... 0 M/2 M-1 3M/2-1 2M-1 i=2 M i=L-1 M i=1 M i=3 M i=L M M n * §5.4 对周期图作平滑后的功率谱估计 分段后的每个段序列乘上窗序列w(n)，求付氏变换： 其中，w(n)的长度为M，可为矩形、三角、升余弦窗。 如三角窗时： 或： * §5.4 对周期图作平滑后的功率谱估计 升余弦窗（哈明）时： 每段的修正周期图： (43) 与(32)式相比，多乘一个w(n)，且分母多了U，U代表窗函数能量： (44) * §5.4 对周期图作平滑后的功率谱估计 整个长度为N的实平稳随机序列的功率谱估值是上述L个修正周期图的平均： (45) 这称welch分段修正周期图平均法（可以证明此种以2：1覆盖分段法时的方差最小）。Welch已证明： 其中不分段时简化为： (46) * §5.4 对周期图作平滑后的功率谱估计 在计算机上实现时，对(45)式进行频域取样： (47) 可采用FFT降低计算量。 * §5.4 对周期图作平滑后的功率谱估计 功率谱估计小结： 1. 直接周期图估计IN(ω)，∵N→∞时var[IN(ω)]≠0, 不满足一致估计条件，一般不用（对谱估计精度要求高时）。 2.分段周期图平滑法Bxx(ω)，满足一致估计条件，但必须每段长度MN时才适应。 3.窗口平滑法Sxx(ω)，也满足一致估计条件（但仍需窗口长度MN），优点是var[Sxx(ω)]≤var[IN(ω)]且可按要求的方差大小，调节窗口长度M。 4.Welch修改法BWxx(ω)，满足一致估计条件，2:1覆盖分段时可得到比以上方法更小的方差，缺点是增加计算复杂度，但便于用FFT快速运算。 5. 参数化谱估计（略），具有不降低谱分辨率的特点，用于雷达、目标分类、医学图像处理，而经典谱估计用于一般的工程上。 * §5.5 功率谱的测量方法 1. 偏置（奇型）离散傅立叶变换（ODFT） (48) x(n)为实序列，0≤k≤N-1。 (49) * §5.5 功率谱的测量方法 2. 由ODFT取代DFT运算 将x(n)分成前后两半： ODFT的偶数谱线： * §5.5 功率谱的测量方法 其中， 称规范运算 (50) (51) * §5.5 功率谱的测量方法 (50)式表明：长度为N的实序列x(n)，其ODFT的偶数谱线可由复序列u(n)的N/2点FFT运算求得（(50)），而复序列u(n)由x(n)的前半部分作实部、后半部分作虚部的规范运算得到。ODFT[x(n)]的奇数谱线则由(49)式直接得到： (52) 由于FFT运算时N是2为基数的数，∴D[N-(2k+1)]即奇数谱线，它直接由偶数谱线D(2k)取共轭得到。（因ODFT与一般的DFT仅差半个取样间隔，可代替实信号序列的谱) * §5.5 功率谱的测量方法 3. 功率谱估值的测量方法举例 [例] 不分段的welch周期图与一般周期图的测量流程由(46)式表示的不分段的welch周期图，频域取样时为： (53) 当w(n)=1 (0≤n≤N-1)时，即一般周期图： (54) 式(53)中的 * §5.5 功率谱的测量方法 (53)式可写为 于是可由下述的流程图测量： 窗口加权 w(n) 规范运算 N/2点 ODFT R2+I2 除NU 或除N 补奇数 谱线 内存U x(n) K1 K2 x’(n)=x(n)w(n) D(2k) 或 K1断、K2合：不分段 K1合、K2断： * * * * * * * * * * * * * * * * * * * * * 第五章 功率谱估计(Power Spectrum Estimation) * 参考章节 （1）胡广书：第11章 （2）奥本海姆：第11章 （3）离散时间信号处理：第11章 * §5.1 估值与一致估计条件 1. 算术平均及与概率平均的关系 算术平均： (1a) 自相关算术平均： (2) * §5.1 估值与一致估计条件 由3-(26b)，mx＝0时： 以xnxn+m取代 ，则功率谱算术平均 ： (3) * §5.1 估值与一致估计条件 2. 从算术平均到有限长段的估值 平均的估值： (1b) 自相关估值： (4) 由式(3)，功率谱估值： (5) （这也称为周期