esprit算法研究应用.doc

课程设计汇报 试验名称： ESPRIT算法研究 试验日期： 姓 名： 学 号： 哈尔滨工业大学（威海） 设计任务 实现空间谱估量算法，并考察算法性能。 方案设计 由均匀线阵形式，确定阵列导向矢量； 由阵列导向矢量，对接收信号进行建模拟真； 由ESPRIT算法实现信号DOA估量； 考察算法性能和信噪比，采样率，观察时间等参数关系。 设计原理 3.1空间谱估量数学模型 空间谱估量就是利用空间阵列实现空间信号参数估量一项专门技术。整个空间谱估量系统应该由三部分组成：空间信号入射、空间阵列接收及参数估量。对应地可分为三个空间，即目标空间、观察空间及估量空间，也就是说空间谱估量系统由这三个空间组成，其框图见图1。 图1 空间谱估量系统结构 对于上述系统结构，作以下几点说明。 (1)目标空间是一个由信号源参数和复杂环境参数张成空间。对于空间谱估量系统，就是利用特定部分方法从这个复杂目标空间中估量出信号未知参数。 (2)观察空间是利用空间按一定方法排列阵元，来接收目标空间辐射信号。因为环境复杂性，所以接收数据中包含信号特征(方位、距离、极化等)和空间环境特征(噪声、杂波、干扰等)。另外因为空间阵元影响，接收数据中一样也含有空间阵列一些特征(互耦、通道不一致、频带不一致等)。这里观察空间是一个多维空间，即系统接收数据是由多个通道组成，而传统时域处理方法通常只有一个通道。尤其需要指出是：通道和阵元并不是一一对应，通道是由空间一个、多个或全部阵元合成(可用加权或不加权)，当然空间某个特定阵元可包含在不一样通道内。 (3)估量空间是利用空间谱估量技术(包含阵列信号处理中部分技术，如阵列校正、空域滤波等技术)从复杂观察数据中提取信号特征参数。 从系统框图中能够清楚看出，估量空间相当于是对目标空间一个重构过程，这个重构精度由众多原因决定，如环境复杂性、空间阵元间互耦、通道不一致、频带不一致等。 3.2 阵列信号处理 首先，考虑N个远场窄带信号入射到空间某阵列上，阵列天线由M个阵元组成，这里假设阵元数等于通道数，即各阵元接收到信号后经过各自传输信道送四处理器，也就是说处理器接收来自M个通道数据。 (3.2-1) 式中，是接收信号幅度，是接收信号相位，是接收信号频率。在窄带远场信号源假设下，有 (3.2-2) 依据式(3.2-1)和式(3.2-2)，显然有下式成立： (3.2-3) 则能够得到第L个阵元接收信号为 (3.2-4) 式中，为第L个阵元对第i个信号增益，表示第L个阵元在t时刻噪声，表示第i个信号抵达第L个阵元时相对参考阵元时延。 将M个阵元在特定时刻接收信号排列成一个列矢量，可得 (3.2-5) 在理想情况下，假设阵列中各阵元是各向同性且不存在通道不一致、互耦等原因影响，则式(3.2-4)中增益能够省略(即归一化1)，在此假设下式(3.2-5)能够简化为 (3.2-6) 将式(3.2-6)写成矢量形式以下： (3.2-7) 式中，为阵列维快拍数据矢量，为阵列维噪声数据矢量，为空间信号维矢量，A为空间阵列维流型矩阵(导向矢量阵)，且 (3.2-8) 其中导向矢量 (3.2-9) 式中，c为光速，为波长。 由上述知识可知，一旦知道阵元间延迟表示式τ，就很轻易得出待定空间阵列导向矢量或阵列流型。下面推导一下空间阵元间延迟表示式。假设空间任意两个阵元，其中一个为参考阵元(在原点)，另一个阵元坐标为(x ,y, z)，两阵元几何关系见图，图中“×”表示阵元。 图2 空间任意两阵元几何关系 由几何关系能够推导出两阵元波程差为 (3.2-10) 这里波程差其实就是在x轴上两阵元间延迟、在y轴上两阵元间延迟和在z轴上两阵元间延迟之和。 依据式(3.2-10)结论，下面给出实际环境中常见多个阵列及阵元间相互延迟表示式。 (1)平面阵　设阵元位置为，以原点为参考点，另假设信号入射参数为，分别表示方位角和俯仰角，其中方位角表示和x轴夹角。 (2)线阵设 阵元位置为，以原点为参考点，另假设信号入射参数为，表示方位角，其中方位角表示和y轴夹角（即和线阵法线夹角），则有 (3.2-11) (3)均匀圆阵 设以均匀圆阵圆心为参考点，则有 (3.2-12) 其中方位角表示和x轴夹角，r为圆半径。 3.3旋转不变子空间算法原理 3.3.1信号模型 算法介绍前，首先对信号进行建模。为了推导分析方便，将波达方向数学模型做以下理想状态假设： 1) 阵列形式为