基于斜投影的极化滤波方法，它涉及极化信号处理及抗干扰领域。当目标信号和干扰的极化状态不正交时，传统极化滤波对目标信号幅度和相位有影响，零相移极化滤波通过变换和补偿来消除目标信号幅度和相位的影响。此外，在目标信号和干扰的极化角相同、极化角相差不同或椭圆率角相同、方位角不同时，现有的极化滤波在抑制干扰的同时也抑制了目标信号。本发明专利技术要解决的技术问题是完善和发展极化滤波的方法，能实现目标信号的准确恢复，并具有计算量小、简单、实用的优点。它主要包含四步：第一步：估计目标和干扰的极化状态；第二步：分别建立目标信号和干扰的极化子空间；第三步：根据建立的子空间，计算出本发明专利技术的滤波算子；第四步：提取目标信号。本方法实现了现有的极化滤波方法的补充和完善。
【技术实现步骤摘要】
本专利技术涉及一种极化滤波的算法实现，更具体的说，本专利技术涉及一种利用斜投影技术和极化信号处理技术实施极化滤波的方法。
技术介绍
作为时域信号处理、频域信号处理和空域信号处理的补充形式，极化域信号处理技术越来越受到人们的关注，在雷达和通信系统中得到了较广泛的应用。在抗干扰技术中，极化抗干扰有不同于频域抗干扰和空域抗干扰的显著特点，它是利用干扰信号和目标信号在极化状态的差别来抑制干扰而提取信号的，因此不受干扰频率和空间入射方向等因素的限制。只要干扰和目标的极化状态存在差别，即使两者的其他特性都相同，也能通过极化滤波技术将其分离，或者对干扰进行抑制。 传统极化滤波方法虽然能够较好的抑制干扰，但是由于它是建立在正交投影的基础上，因此当目标和干扰的极化状态非正交时，抑制干扰的同时会对目标信号的幅度和相位产生影响，期望通过直接的补偿消除这些影响是十分困难的。 专利技术人在《电子学报》2004年9月第32卷第9期的文章“零相移瞬时极化滤波器”和《IEEE Transaction on Aerospace and Electronic Systems》2007年10月第43卷第4期的文章“Null phase-shift polarization filtering for high-frequencyradar”2篇文章中提出了通过线性极化矢量变换(LPVT)将目标信号的极化角变换成90度，然后进行极化滤波，进而采用幅度/相位补偿(A/PCD)恢复目标信号的幅度和相位信息，实现了对目标信号的准确还原。文章中利用高频雷达的实测数据证明了该方法的有效性。这2篇文章提出的方法在目前的极化滤波方法中是比较有影响力的算法，也是一种比较成熟的算法。 此外，在目标信号和干扰的极化角相同、极化角相差不同或者椭圆率角相同、方位角不同的条件下，现有的极化滤波(传统极化滤波和零相移极化滤波)在抑制干扰的同时，也抑制了目标信号，无法完成极化滤波的目的。因此，进一步完善和发展极化滤波的理论和极化滤波的设计方法是十分有必要的。
技术实现思路
本专利技术要解决的技术问题是完善和发展现有的极化滤波的设计方法，为实现上述目标，本专利技术提供了一种，并分别从笛卡尔坐标系和椭圆坐标系给出了相应的滤波算子。 本专利技术若采用笛卡尔坐标系描述时，这种极化状态表示为极化角和极化角相差，为了实现上述专利技术目的，本专利技术在迪笛卡尔坐标系中的步骤为 第一步确定目标信号和干扰的极化角和极化角相差； 第二步根据斜投影算子的理论模型，由第一步中得出的目标信号和干扰的极化参数，建立目标信号极化子空间S和干扰极化子空间I，分别为 T为矩阵的转置，上述公式中εs和εi分别为确定出的目标信号的极化角和干扰的极化角，δs和δi分别为确定出的目标信号极化角相差和干扰极化角相差。建立这两个极化子空间的目的在于计算相应的斜投影算子； 第三步首先计算由第二步中得到的目标信号和干扰的极化子空间对应的斜投影算子，计算公式如下 其中E为单位矩阵，H代表矩阵的共轭转置。 通过如下运算得到的滤波算子Q 其中 γ＝(cosεssinεi)2+(sinεscosεi)2+2sinεscosεssinεicosεicos(δi-δs) 第四步设定接收信号表示为y(t)，通过下式提取出目标信号s(t) s(t)＝y(t)TQ 本专利技术涉及的算法若采用椭圆坐标系描述时，这种极化状态表示为极化的椭圆率角和方位角。为了实现上述专利技术目的，本专利技术在椭圆坐标系中的步骤为 第一步确定目标信号和干扰的椭圆率角和方位角； 第二步根据斜投影算子的理论模型，由第一步中得出的目标信号和干扰的极化参数，建立目标信号极化子空间S和干扰极化子空间I，分别为 上述公式中τs和τi分别为确定出的目标信号的椭圆率角和干扰的椭圆率角， 和 分别为确定出的目标信号方位角和干扰方位角。建立这两个极化子空间的目的在于计算相应的斜投影算子； 第三步首先计算由第二步中得到的目标信号和干扰的极化子空间对应的斜投影算子，计算公式如下 其中E为单位矩阵，H为矩阵的共轭转置。 通过如下运算得到的滤波算子Q 其中 第四步设定接收信号表示为y(t)，通过下式提取出目标信号s(t) s(t)＝y(t)TQ 与现有的极化滤波方法相比，本专利技术具有以下优点1、不需要对接收信号进行极化矢量变换；2、不需要进行补偿即可保持目标信号的幅度和相位特性；3、在目标信号和干扰的极化角相同、极化角相差不同或者椭圆率角相同、方位角不同的条件下也能进行极化滤波，方便的实现目标信号和干扰的分离，或者干扰的抑制。实现了目标信号的准确恢复，并具有计算量小、简单、实用的优点。 附图说明 图1是本专利技术的步骤流程； 图2是不同方法对目标信号相位影响的比较； 图3是本专利技术方法处理前后信号的时域波形。 具体实施例方式 本专利技术能在笛卡尔坐标系和椭圆坐标系两种坐标系条件下实施，由图1可知本专利技术公开的方法的具体实施方式包含四个步骤。两种坐标系下实施的步骤相同，但对应的具体实施方式不同。 实施方式1本专利技术若采用笛卡尔坐标系描述时，这种极化状态表示为极化角和极化角相差，在笛卡尔坐标系条件下的具体实施方式为 第一步确定目标信号和干扰的极化角和极化角相差； 第二步根据斜投影算子的理论模型，由第一步中得出的目标信号和干扰的极化参数，建立目标信号极化子空间S和干扰极化子空间I，分别为 T为矩阵的转置，上述公式中εs和εi分别为确定出的目标信号的极化角和干扰的极化角，δs和δi分别为确定出的目标信号极化角相差和干扰极化角相差。建立这两个极化子空间的目的在于计算相应的斜投影算子； 第三步首先计算由第二步中得到的目标信号和干扰的极化子空间对应的斜投影算子，计算公式如下 其中E为单位矩阵，H为矩阵的共轭转置； 通过如下运算得到的滤波算子Q 其中 γ＝(cosεssinεi)2+(sinεscosεi)2+2sinεscosεssinεi cosεicos(δi-δs) 第四步设定接收信号表示为y(t)，则可通过下式提取出目标信号s(t) s(t)＝y(t)TQ 实施方式2本专利技术若采用椭圆坐标系描述时，这种极化状态表示为极化的椭圆率角和方位角，在椭圆坐标系条件下的具体实施方式为 第一步确定目标信号和干扰的椭圆率角和方位角； 第二步根据斜投影算子的理论模型，由第一步中得出的目标信号和干扰的极化参数，建立目标信号极化子空间S和干扰极化子空间I，分别为 上述公式中τs和τi分别为确定出的目标信号的椭圆率角和干扰的椭圆率角， 和 分别为确定出的目标信号方位角和干扰方位角。建立这两个极化子空间的目的在于计算相应的斜投影算子； 第三步首先计算由第二步中得到的目标信号和干扰的极化子空间对应的斜投影算子，本文档来自技高网...
【技术保护点】
基于斜投影的极化滤波方法，其特征在于在笛卡尔坐标系下，包含四个步骤： 第一步：确定目标信号和干扰的极化角和极化角相差； 第二步：根据斜投影算子的理论模型，由第一步中得出的目标信号和干扰的极化参数，建立目标信号极化子空间Ｓ和干扰极化子空间Ｉ，分别为： Ｓ＝［ｃｏｓε↓［ｓ］，ｓｉｎε↓［ｓ］ｅ↑［ｊδ↓［ｓ］］↑［Ｔ］， Ｉ＝［ｃｏｓε↓［ｉ］，ｓｉｎε↓［ｉ］ｅ↑［ｊδ↓［ｉ］］］↑［Ｔ］； Ｔ为矩阵的转置，上述公式中ε↓［ｓ］和ε↓［ｉ］分别为确定出的目标信号的极化角和干扰的极化角，δ↓［ｓ］和δ↓［ｉ］分别为确定出的目标信号极化角相差和干扰极化角相差； 第三步：首先计算由第二步中得到的目标信号和干扰的极化子空间对应的斜投影算子，计算公式如下： Ｅ↓［ＳＩ］＝Ｓ（Ｓ↑［Ｈ］Ｐ↓［Ｉ］↑［⊥］Ｓ）↑［－１］Ｓ↑［Ｈ］Ｐ↓［Ｉ］↑［⊥］ 其中Ｐ↓［Ｉ］↑［⊥］＝Ｅ－Ｐ↓［Ｉ］＝Ｅ－Ｉ（Ｉ↑［Ｈ］Ｉ）↑［－１］Ｉ↑［Ｈ］，Ｅ为单位矩阵，Ｈ代表矩阵的共轭转置； 通过如下运算得到基于斜投影的极化滤波方法的滤波算子Ｑ： Ｑ＝Ｅ↓［ＳＩ］↑［Ｔ］Ｒ＝Ｅ↓［ＳＩ］↑［Ｔ］＊＝１／γ［α，β］↑［Ｔ］ 其中 γ＝（ｃｏｓε↓［ｓ］ｓｉｎε↓［ｉ］）↑［２］＋（ｓｉｎε↓［ｓ］ｃｏｓε↓［ｉ］）↑［２］＋２ｓｉｎε↓［ｓ］ｃｏｓε↓［ｓ］ｓｉｎε↓［ｉ］ｃｏｓε↓［ｉ］ｃｏｓ（δ↓［ｉ］－δ↓［ｓ］） α＝ｃｏｓ↑［３］ε↓［ｓ］ｓｉｎ↑［２］ε↓［ｉ］－ｓｉｎε↓［ｓ］ｃｏｓ↑［２］ε↓［ｓ］ｅ↑［ｊ（δ↓［ｉ］－δ↓［ｓ］）］＋ｓｉｎ↑［２］ε↓［ｓ］ｃｏｓε↓［ｓ］ｃｏｓ↑［２］ε↓［ｉ］ｅ↑［ｊδ↓［ｓ］］－ｓｉｎ↑［３］ε↓［ｓ］ｓｉｎε↓［ｉ］ｃｏｓε↓［ｉ］ｅ↑［ｊδ↓［ｉ］］ β＝ｓｉｎ↑［３］ε↓［ｓ］ｃｏｓ↑［２］ε↓［ｉ］－ｓｉｎ↑［２］ε↓［ｓ］ｃｏｓε↓［ｓ］ｓｉｎε↓［ｉ］ｃｏｓε↓［ｉ］ｅ↑［ｊ（δ↓［ｓ］－δ↓［ｉ］）］＋ｓｉｎε↓［ｓ］ｃｏｓ↑［２］ε↓［ｓ］ｃｏｓ↑［２］ε↓［ｉ］ｅ↑［－ｊδ↓［ｓ］］－ｃｏｓ↑［３］ε↓［ｓ］ｓｉｎε↓［ｉ］ｃｏｓε↓［ｉ］ｅ↑［－ｊδ↓［ｉ］］ 第四步骤：设定接收信号表示为ｙ（ｔ），通过下式提取出目标信号ｓ（ｔ）： ｓ（ｔ）＝ｙ（ｔ）↑［Ｔ］Ｑ。...
【技术特征摘要】
【专利技术属性】
技术研发人员：毛兴鹏，张钦宇，
申请(专利权)人：哈尔滨工业大学威海，
类型：发明
国别省市：37[中国|山东]
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。