快速约束极值子空间跟踪算法.pdfVIP

第 34卷第 4期 华 北 水 利 水 电 学 院 学 报 Vo1．34NO．4 2013年 8月 JournalofNorthChinaInstituteofWaterConservancyandHydroelectricPower Aug．2013 DOI：10．3969／j．issn．1002—5634．2013．04．031 快速约束极值子空间跟踪算法 魏志强，牛慧贤 (华北水利水 电大学 数学与信息科学学院，河南 郑州 450045) 摘 要 ：子空间跟踪算法是许多工程计算问题的核心．Valizadeh等人基于约束优化 问题提 出了一种快速投影逼近 子空间跟踪算法 ．在指 出快速投影逼近子空间跟踪算法不足的同时，应用矩阵求逆引理给 出了一种新 的快 速子空间跟踪算法——快速约束极值子空间跟踪算法．仿真实验表明，所提算法是收敛且稳定的，其性能 优于或相当于几种常见的快速子空间跟踪算法． 关键词：子空间跟踪算法 ；约束投影逼近；快速算法 中图分类号 ：029 文献标识码 ：A 文章编号 ：1002—5634(2013)04—0121—04 在许多工程计算问题 中，比如数据压缩 、数字 个主特征 向量张成 的子空间range(W(t))，其 中信 滤波 、参数估计 、模式识别 、神经分析等，特别是在无 号子空 间权重矩 阵 W(t)是一个 N ×K 阶矩 阵． 线通信领域中，子空间的估计与跟踪都是非常重要 W(t)是如下约束极值问题 的…．传统的提取子空间信息 的方法是数据 自相关 ．，(W(t))= 矩阵的特征值分解 (EVD)或者是数据矩 阵的奇异 ∑ ()一w()W(i一1)“(￡)II (1) 值分解 (SVD)．这两类方法 的性 能都是非常优越 的，但是其计算复杂度都太大 ，高达 0(N3)(其 中Ⅳ W “(t)W(t)=IK (2) 是数据向量的维数)，不适合工程问题 的实时实现． 的解 ，式中 ：卢(0卢1)为遗忘因子；J为K×K阶 因此 ，如何快速地实现子空 间的估计与跟踪是一个 单位矩阵．求解此约束极值问题可得 热门问题 ，吸引了许多学者 的研究并提 出了各种解 W(t)=C (t)(C (t)“C (t))一。 (3) 决方案 ¨ ．根据计算的复杂度 ，可以将子空间跟踪 符号(·) 表示矩阵平方根逆 ，而 算法分为0( )、0(Ⅳ )和 0(NK)3类 (其 中 是 J，(i)=W (i一1)“()，i=1，2，…，t (4) 信号子空间的维数)．目前 ，研究的热点与兴趣集 中 C ()：∑ t-iX()J，()“ (5) 在第 3类算法上 ，称之为快速子 空 间跟踪算法． Valizadeh等… 给出了一种称之为约束投影逼近的 基于式 (3)，(4)，(5)计算 W(￡)的计算过程运 子空间跟踪算法 ，但其推导有不 明确 的地方．鉴于 算量很大 ，达到 0(NK)+0(K)数量级，不利于工 此，笔者在文献 [1]的基础上 ，给出了一种新的快速 程实现．以下是不同于文献[1]，类似于文献[5]，考 子空间跟踪算法． 虑新的减小计算复杂度 的方法． 1 约束投影逼近子空间跟踪 2 快速约束极值子空间跟踪 在此 ，简要叙述文献 [1]的相关结果．考虑一个 记 (￡)= (C ()“C (￡)) (6) Ⅳ维 的 随机 向