• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 15
    • 下载费用:30 金币  

    重庆时时彩平台有赠送彩金的吗: 基于时频谱图分解的SAR时变窄带干扰抑制方法.pdf

    关 键 词:
    基于 频谱 分解 SAR 窄带 干扰 抑制 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201310239985.9

    申请日:

    2013.06.17

    公开号:

    CN103323819A

    公开日:

    2013.09.25

    当前法律状态:

    终止

    有效性:

    无权

    法律详情: 未缴年费专利权终止IPC(主分类):G01S 7/36申请日:20130617授权公告日:20141217终止日期:20150617|||授权|||实质审查的生效IPC(主分类):G01S 7/36申请日:20130617|||公开
    IPC分类号: G01S7/36; G01S13/90 主分类号: G01S7/36
    申请人: 西安电子科技大学
    发明人: 周峰; 陶明亮; 白雪茹; 刘妍
    地址: 710071 陕西省西安市太白南路2号
    优先权:
    专利代理机构: 广东秉德律师事务所 44291 代理人: 杨焕军
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201310239985.9

    授权公告号:

    |||103323819B||||||

    法律状态公告日:

    2016.08.03|||2014.12.17|||2013.10.30|||2013.09.25

    法律状态类型:

    专利权的终止|||授权|||实质审查的生效|||公开

    摘要

    基于时频谱图分解的SAR时变窄带干扰抑制方法,包括以下步骤:获取回波数据;判断回波数据是否存在窄带干扰,若存在干扰执行下一步;根据存在窄带干扰的回波数据的列向量获得混合信号时频谱图,并进行奇异值分解,获得对应的待分离信号矩阵;对待分离信号矩阵进行分离,获得干扰信号矩阵;根据最小二乘准则利用干扰信号矩阵中的干扰信号分量重构干扰信号时频谱图:从混合信号时频谱图中剔除干扰信号得到有用信号时频谱图,并对有用信号时频谱图进行逆短时傅里叶变换,得到剔除干扰后的回波数据的列向量;判断是否遍历完所有回波数据的列向量,若遍历完则得到剔除干扰后的回波数据。本发明具有运算量较小,效率较高并且没有数据缺失的优点。

    权利要求书

    权利要求书
    1.   基于时频谱图分解的SAR时变窄带干扰抑制方法,其特征在于,包括以下步骤:
    步骤1、获取回波数据;
    雷达以脉冲重复频率发射并接收脉冲,得到以距离为行向量、以方位为列向量的回波数据;
    步骤2、判断回波数据是否存在窄带干扰,若存在干扰执行步骤3,否则执行步骤7;
    步骤3、根据存在窄带干扰的回波数据的列向量获得混合信号时频谱图S,并对混合信号时频谱图S进行奇异值分解,获得对应的待分离信号矩阵Y;
    对回波数据的列向量进行短时傅里叶变换,得到混合信号时频谱图S,并按照下式对混合信号时频谱图S进行奇异值分解:
    S=UΛVH
    上式中的U为混合信号时频谱图的左特征向量矩阵,Λ为奇异值对角矩阵,V为混合信号时频谱图的右特征向量矩阵,(·)H表示矩阵共轭转置操作;
    将分解得到的混合信号时频谱图的左特征向量矩阵U的前L个特征向量作为行向量,构成待分离信号矩阵Y,其中的L为主分量个数,主分量个数L根据最小描述长度准则计算得到:
    <mrow><MI>L</MI><MO>=</MO><MI>arg</MI><MI>min</MI><MO>{</MO><MO>-</MO><MI>m</MI><MROW><MO>(</MO><MI>n</MI><MO>-</MO><MI>L</MI><MO>)</MO></MROW><MI>log</MI><MROW><MO>(</MO><MROW><MO>(</MO><MI>n</MI><MO>-</MO><MI>L</MI><MO>)</MO></MROW><MSUP><MROW><MO>(</MO><MUNDEROVER><MROW><MI>Π</MI><MI></MI></MROW><MROW><MI>p</MI><MO>=</MO><MI>L</MI><MO>+</MO><MN>1</MN></MROW><MI>n</MI></MUNDEROVER><MSUB><MI>σ</MI><MI>p</MI></MSUB><MO>)</MO></MROW><MFRAC><MN>1</MN><MROW><MI>n</MI><MO>-</MO><MI>L</MI></MROW></MFRAC></MSUP><MSUP><MROW><MO>(</MO><MUNDEROVER><MROW><MI>Σ</MI><MI></MI></MROW><MROW><MI>p</MI><MO>=</MO><MI>L</MI><MO>+</MO><MN>1</MN></MROW><MI>n</MI></MUNDEROVER><MSUB><MI>σ</MI><MI>p</MI></MSUB><MO>)</MO></MROW><MROW><MO>-</MO><MN>1</MN></MROW></MSUP><MO>)</MO><MO>+</MO><MN>0.5</MN><MI>L</MI><MROW><MO>(</MO><MN>2</MN><MI>n</MI><MO>-</MO><MI>L</MI><MO>)</MO></MROW><MI>log</MI><MROW><MO>(</MO><MI>m</MI><MO>)</MO></MROW></MROW><MO>}</MO></MROW>]]&gt;</MATH></MATHS><BR>上式中的argmin(·)表示计算使目标函数取最小值时对应变量值的运算操作,m为混合信号时频谱图的左特征向量矩阵U的行数,n为混合信号时频谱图的左特征向量矩阵U的列数,log(·)表示求对数运算操作,表示对第L+1至第n项的求乘积运算操作,σp表示奇异值对角矩阵Λ的第p个奇异值,表示对第L+1至第n项的求和运算操作;<BR>步骤4,对待分离信号矩阵Y进行分离,获得干扰信号矩阵I,包括以下子步骤;<BR>子步骤4a、随机产生大小为L×L的解混矩阵W;<BR>子步骤4b、计算迭代矩阵z;<BR>迭代矩阵z=[z1,z2,…,zp,…,zL],其中,zp为迭代矩阵z的第p个列向量,<MATHS id=cmaths0002 num="0002"><MATH><![CDATA[<mrow><MSUB><MI>z</MI><MI>p</MI></MSUB><MO>=</MO><MI>E</MI><MO>{</MO><MI>Y</MI><MSUP><MROW><MO>(</MO><MSUBSUP><MI>w</MI><MI>p</MI><MI>H</MI></MSUBSUP><MI>Y</MI><MO>)</MO></MROW><MO>*</MO></MSUP><MSUP><MROW><MO>(</MO><MN>0.1</MN><MO>+</MO><MSUP><MROW><MO>|</MO><MSUBSUP><MI>w</MI><MI>p</MI><MI>H</MI></MSUBSUP><MI>Y</MI><MO>|</MO></MROW><MN>2</MN></MSUP><MO>)</MO></MROW><MROW><MO>-</MO><MN>1</MN></MROW></MSUP><MO>}</MO><MO>-</MO><MI>E</MI><MO>{</MO><MN>0.1</MN><MSUP><MROW><MO>(</MO><MN>0.1</MN><MO>+</MO><MSUP><MROW><MO>|</MO><MSUBSUP><MI>w</MI><MI>p</MI><MI>H</MI></MSUBSUP><MI>Y</MI><MO>|</MO></MROW><MN>2</MN></MSUP><MO>)</MO></MROW><MROW><MO>-</MO><MN>2</MN></MROW></MSUP><MO>}</MO><MSUB><MI>w</MI><MI>p</MI></MSUB><MO>,</MO></MROW>]]&gt;</MATH></MATHS>E(·)表示求期望运算操作,Y为待分离信号矩阵,wp为解混矩阵W的第p个列向量,(·)H表示矩阵共轭转置操作,(·)*表示共轭运算操作,|·|表示求绝对值运算操作;<BR>子步骤4c、对解混矩阵和迭代矩阵的所有列向量依次判断是否满足若所有列向量都满足该条件,则执行子步骤4d,否则,对不满足条件的wp值进行更新,令所述执行子步骤4b;其中,||·||表示求模值运算操作,zp为迭代矩阵z的第p个列向量,wp为解混矩阵w的第p个列向量,ε表示设定的阈值;<BR>子步骤4d、将解混矩阵w各列向量去相关得到最终解混矩阵<MATHS id=cmaths0003 num="0003"><MATH><![CDATA[<mrow><MOVER><MI>W</MI><MO>~</MO></MOVER><MO>=</MO><MI>W</MI><MSUP><MROW><MO>(</MO><MSUP><MI>W</MI><MI>H</MI></MSUP><MI>W</MI><MO>)</MO></MROW><MROW><MO>-</MO><MN>1</MN><MO>/</MO><MN>2</MN></MROW></MSUP><MO>;</MO></MROW>]]&gt;</MATH></MATHS><BR>子步骤4e、根据最终解混矩阵计算干扰信号矩阵<BR>步骤5、根据最小二乘准则,将干扰信号矩阵I中的干扰信号分量按照下式重构干扰信号时频谱图:<BR><MATHS id=cmaths0004 num="0004"><MATH><![CDATA[<mrow><MSUB><MI>S</MI><MI>k</MI></MSUB><MO>=</MO><MSUB><MI>i</MI><MI>k</MI></MSUB><MROW><MO>(</MO><MSUBSUP><MI>i</MI><MI>k</MI><MO>+</MO></MSUBSUP><MI>S</MI><MO>)</MO></MROW></MROW>]]&gt;</MATH></MATHS><BR>其中,Sk表示第k个干扰信号分量的时频谱图,ik表示干扰信号矩阵I的第k个行向量,+表示求矩阵的广义逆运算操作,S为混合信号时频谱图;<BR>步骤6、从混合信号时频谱图S中剔除干扰信号得到有用信号时频谱图S′,并对有用信号时频谱图S′进行逆短时傅里叶变换,得到剔除干扰后的回波数据的列向量;<BR>有用信号时频谱图其中,|·|表示求绝对值运算操作,表示对第1项至第L项的求和运算操作,Sk为第k个干扰信号的时频谱图,exp(·)表示求指数运算操作,j为虚数符号,Θ为对混合信号时频谱图S求得的相位角;<BR>步骤7,判断是否遍历完所有回波数据的列向量,若遍历完,则执行步骤8,否则执行步骤2;<BR>步骤8,获取剔除干扰后的回波数据;<BR>前述步骤中1≤p≤L,1≤k≤L。<BR><BR>2.&nbsp;&nbsp; 如权利要求1所述的基于时频谱图分解的SAR时变窄带干扰抑制方法,其特征在于:所述步骤2中通过一维距离像的峭度来判断回波数据是否存在干扰,判断过程为:对步骤1的回波数据的列向量进行傅里叶变换,得到一维距离像,判断一维距离像的峭度κ是否大于等于设定的幅度阈值,若是,则认为回波数据中存在窄带干扰;<BR>一维距离像的峭度<BR>上式中的E(·)表示求期望运算操作,y为一维距离像的幅度值,μ为一维距离像幅值的均值,d为一维距离像幅值的标准差。<BR><BR>3.&nbsp;&nbsp; 如权利要求2所述的基于时频谱图分解的SAR时变窄带干扰抑制方法,其特征在于:所述幅度阈值为10~200。<BR></p></div> </div> </div> <div class="zlzy"> <div class="zltitle">说明书</div> <div class="gdyy"> <div class="gdyy_show"><p>说明书基于时频谱图分解的SAR时变窄带干扰抑制方法 <BR>技术领域 <BR>本发明属于信号处理技术领域,尤其涉及一种雷达信号处理技术中的基于时频谱图分解的SAR时变窄带干扰抑制方法。 <BR>背景技术 <BR>窄带干扰是针对合成孔径雷达系统(Synthetic&nbsp;Aperture&nbsp;Radar,简称SAR)的主要干扰方式之一。所谓窄带干扰是指干扰信号带宽与合成孔径雷达发射信号带宽之比很?。ㄈ?%)的信号。窄带干扰可分为非人为干扰和人为干扰两类。非人为干扰主要包括与SAR处于同一频段的通讯设备、电视网和其它低频波段辐射设备;常见的人为干扰主要包括强脉冲干扰和频率引导“瞄准式”干扰,人为窄带干扰具有干扰功率大、针对性强、干扰实现简单、形式多变等特点。这些干扰的存在降低了SAR回波数据的信噪比,尤其是在干扰功率较大的情况下,会使得图像中出现亮线并且使图像变模糊,从而严重降低了SAR图像质量。因此,如何对窄带干扰形式进行有效的识别和抑制是目前SAR干扰处理技术的重要课题。 <BR>Feng&nbsp;Zhou等人在文献“Narrow?Band&nbsp;Interference&nbsp;Suppression&nbsp;for&nbsp;SAR&nbsp;Based&nbsp;on&nbsp;Independent&nbsp;Component&nbsp;Analysis”(IEEE&nbsp;Trans.Geosci.Remote&nbsp;Sens.,2013)中提出采用独立分量分析方法进行干扰抑制。该方法利用有用回波信号与干扰信号的统计特性差异,采用独立分量分析方法直接在信号数据域提取窄带干扰,并根据最小二乘原则获得对干扰基信号的复包络的精确估计后,重构干扰信号,从原始回波中消去干扰信号,实现时变窄带干扰抑制。但该方法存在以下不足:其利用时间平滑增加观测通道,减小了样本的数目,致使估计出的干扰信号长度缩短,降低了独立性准则的计算精度,使得干扰信号分离不彻底,从而影响了分离抑制的效果。 <BR>申请号为201210001678.2、发明名称为基于复数经验模态分解的时变窄带干扰抑制方法的中国发明专利申请公开了一种利用复数经验模态分解的时变窄带干扰抑制方法。该方法将信号分解成一系列本征模态函数,在抑制窄带干扰的同时充分保留了信号本身的非线性和非平稳特征,但是该方法对局部信号分析处理未考虑干扰信号的全局统计特性,甚至会分解出虚假的信号,从而不能有效地刻画干扰信号的时频特性,导致干扰信号分离不彻底。 <BR>发明内容 <BR>针对现有技术存在的不足,本发明的目的是提供一种基于时频谱图分解的SAR时变窄带干扰抑制方法,该方法可以准确提取窄带干扰信号的时频特征,在保留有用信息的同时对时变窄带干扰进行抑制,从而实现对雷达回波信号清晰成像。 <BR>为了实现上述目的,本发明采取如下的技术解决方案: <BR>基于时频谱图分解的SAR时变窄带干扰抑制方法,包括以下步骤: <BR>步骤1、获取回波数据; <BR>雷达以脉冲重复频率发射并接收脉冲,得到以距离为行向量、以方位为列向量的回波数据; <BR>步骤2、判断回波数据是否存在窄带干扰,若存在干扰执行步骤3,否则执行步骤7; <BR>步骤3、根据存在窄带干扰的回波数据的列向量获得混合信号时频谱图S,并对混合信号时频谱图S进行奇异值分解,获得对应的待分离信号矩阵Y; <BR>对回波数据的列向量进行短时傅里叶变换,得到混合信号时频谱图S,并按照下式对混合信号时频谱图S进行奇异值分解: <BR>S=UΛVH <BR>上式中的U为混合信号时频谱图的左特征向量矩阵,Λ为奇异值对角矩阵,V为混合信号时频谱图的右特征向量矩阵,(·)H表示矩阵共轭转置操作; <BR>将分解得到的混合信号时频谱图的左特征向量矩阵U的前L个特征向量作为行向量,构成待分离信号矩阵Y,其中的L为主分量个数,主分量个数L根据最小描述长度准则计算得到: <BR><MATHS num="0001"><MATH><![CDATA[ <mrow><MI>L</MI> <MO>=</MO> <MI>arg</MI> <MI>min</MI> <MO>{</MO> <MO>-</MO> <MI>m</MI> <MROW><MO>(</MO> <MI>n</MI> <MO>-</MO> <MI>L</MI> <MO>)</MO> </MROW><MI>log</MI> <MROW><MO>(</MO> <MROW><MO>(</MO> <MI>n</MI> <MO>-</MO> <MI>L</MI> <MO>)</MO> </MROW><MSUP><MROW><MO>(</MO> <MUNDEROVER><MROW><MI>Π</MI> <MI></MI></MROW><MROW><MI>p</MI> <MO>=</MO> <MI>L</MI> <MO>+</MO> <MN>1</MN> </MROW><MI>n</MI> </MUNDEROVER><MSUB><MI>σ</MI> <MI>p</MI> </MSUB><MO>)</MO> </MROW><MFRAC><MN>1</MN> <MROW><MI>n</MI> <MO>-</MO> <MI>L</MI> </MROW></MFRAC></MSUP><MSUP><MROW><MO>(</MO> <MUNDEROVER><MROW><MI>Σ</MI> <MI></MI></MROW><MROW><MI>p</MI> <MO>=</MO> <MI>L</MI> <MO>+</MO> <MN>1</MN> </MROW><MI>n</MI> </MUNDEROVER><MSUB><MI>σ</MI> <MI>p</MI> </MSUB><MO>)</MO> </MROW><MROW><MO>-</MO> <MN>1</MN> </MROW></MSUP><MO>)</MO> <MO>+</MO> <MN>0.5</MN> <MI>L</MI> <MROW><MO>(</MO> <MN>2</MN> <MI>n</MI> <MO>-</MO> <MI>L</MI> <MO>)</MO> </MROW><MI>log</MI> <MROW><MO>(</MO> <MI>m</MI> <MO>)</MO> </MROW></MROW><MO>}</MO> </MROW>]]&gt;</MATH></MATHS> <BR>上式中的argmin(·)表示计算使目标函数取最小值时对应变量值的运算操作,m为混合信号时频谱图的左特征向量矩阵U的行数,n为混合信号时频谱图的左特征向量矩阵U的列数,log(·)表示求对数运算操作,(·)表示对第L+1至第n项的求乘积运算操作,σp表示奇异值对角矩阵Λ的第p个奇异值,(·)表示对第L+1至第n项的求和运算操作; <BR>步骤4,对待分离信号矩阵Y进行分离,获得干扰信号矩阵I,包括以下子步骤; <BR>子步骤4a、随机产生大小为L×L的解混矩阵W; <BR>子步骤4b、计算迭代矩阵z;迭代矩阵z=[z1,z2,…,zp,…,zL],其中,zp为迭代矩阵z的第p个列向量,<MATHS num="0002"><MATH><![CDATA[ <mrow> <MSUB><MI>z</MI> <MI>p</MI> </MSUB><MO>=</MO> <MI>E</MI> <MO>{</MO> <MI>Y</MI> <MSUP><MROW><MO>(</MO> <MSUBSUP><MI>w</MI> <MI>p</MI> <MI>H</MI> </MSUBSUP><MI>Y</MI> <MO>)</MO> </MROW><MO>*</MO> </MSUP><MSUP><MROW><MO>(</MO> <MN>0.1</MN> <MO>+</MO> <MSUP><MROW><MO>|</MO> <MSUBSUP><MI>w</MI> <MI>p</MI> <MI>H</MI> </MSUBSUP><MI>Y</MI> <MO>|</MO> </MROW><MN>2</MN> </MSUP><MO>)</MO> </MROW><MROW><MO>-</MO> <MN>1</MN> </MROW></MSUP><MO>}</MO> <MO>-</MO> <MI>E</MI> <MO>{</MO> <MN>0.1</MN> <MSUP><MROW><MO>(</MO> <MN>0.1</MN> <MO>+</MO> <MSUP><MROW><MO>|</MO> <MSUBSUP><MI>w</MI> <MI>p</MI> <MI>H</MI> </MSUBSUP><MI>Y</MI> <MO>|</MO> </MROW><MN>2</MN> </MSUP><MO>)</MO> </MROW><MROW><MO>-</MO> <MN>2</MN> </MROW></MSUP><MO>}</MO> <MSUB><MI>w</MI> <MI>p</MI> </MSUB><MO>,</MO> </MROW>]]&gt;</MATH></MATHS>E(·)表示求期望运算操作,Y为待分离信号矩阵,wp为解混矩阵W的第p个列向量,(·)H表示矩阵共轭转置操作,(·)*表示共轭运算操作,|·|表示求绝对值运算操作; <BR>子步骤4c、对解混矩阵和迭代矩阵的所有列向量依次判断是否满足若所有列向量都满足该条件,则执行子步骤4d,否则,对不满足条件的wp值进行更新,令所述执行子步骤4b;其中,||·||表示求模值运算操作,zp为迭代矩阵z的第p个列向量,wp为解混矩阵w的第p个列向量,ε表示设定的阈值; <BR>子步骤4d、将解混矩阵w各列向量去相关得到最终解混矩阵<MATHS num="0003"><MATH><![CDATA[ <mrow> <MOVER><MI>W</MI> <MO>~</MO> </MOVER><MO>=</MO> <MI>W</MI> <MSUP><MROW><MO>(</MO> <MSUP><MI>W</MI> <MI>H</MI> </MSUP><MI>W</MI> <MO>)</MO> </MROW><MROW><MO>-</MO> <MN>1</MN> <MO>/</MO> <MN>2</MN> </MROW></MSUP><MO>;</MO> </MROW>]]&gt;</MATH></MATHS> <BR>子步骤4e、根据最终解混矩阵计算干扰信号矩阵 <BR>步骤5、根据最小二乘准则,将干扰信号矩阵I中的干扰信号分量按照下式重构干扰信号时频谱图: <BR><MATHS num="0004"><MATH><![CDATA[ <mrow><MSUB><MI>S</MI> <MI>k</MI> </MSUB><MO>=</MO> <MSUB><MI>i</MI> <MI>k</MI> </MSUB><MROW><MO>(</MO> <MSUBSUP><MI>i</MI> <MI>k</MI> <MO>+</MO> </MSUBSUP><MI>S</MI> <MO>)</MO> </MROW></MROW>]]&gt;</MATH></MATHS> <BR>其中,Sk表示第k个干扰信号分量的时频谱图,ik表示干扰信号矩阵I的第k个行向量,+表示求矩阵的广义逆运算操作,S为混合信号时频谱图; <BR>步骤6、从混合信号时频谱图S中剔除干扰信号得到有用信号时频谱图S′,并对有用信号时频谱图S′进行逆短时傅里叶变换,得到剔除干扰后的回波数据的列向量; <BR>有用信号时频谱图其中,|·|表示求绝对值运算操作,表示对第1项至第L项的求和运算操作,Sk为第k个干扰信号的时频谱图,exp(·)表示求指数运算操作,j为虚数符号,Θ为对混合信号时频谱图S求得的相位角; <BR>步骤7,判断是否遍历完所有回波数据的列向量,若遍历完,则执行步骤8,否则执行步骤2; <BR>步骤8,获取剔除干扰后的回波数据; <BR>前述步骤中1≤p≤L,1≤k≤L。 <BR>进一步的,所述步骤2中通过一维距离像的峭度来判断回波数据是否存在干扰,判断过程为:对步骤1的回波数据的列向量进行傅里叶变换,得到一维距离像,判断一维距离像的峭度κ是否大于等于设定的幅度阈值,若是,则认为回波数据中存在窄带干扰; <BR>一维距离像的峭度 <BR>上式中的E()表示求期望运算操作,y为一维距离像的幅度值,μ为一维距离像幅值的均值,d为一维距离像幅值的标准差。 <BR>进一步的,所述幅度阈值为10~200。 <BR>本发明方法通过对回波数据进行短时傅里叶变换增加观测通道,并对时频谱图进行奇异值分解,降低待处理数据的维数,克服了现有技术中独立分量分析方法通过时间平滑增加观测通道导致样本缺失的不足,使得本发明运算量小,效率高并且没有数据缺失;其次,本发明方法充分利用回波信号和干扰信号的时频特性差异,通过盲源分离方法将干扰信号与回波信号在时频谱图域独立的分离出来,能够很好的刻画干扰信号的时频特性,克服了现有技术中基于复数经验模态分解方法仅考虑局部统计特性而造成干扰抑制不充分的缺点,使得本发明具有信号损失小,抑制效果好的优点。 <BR>附图说明 <BR>图1为本发明方法的流程图; <BR>图2为仿真实验对混合信号进行傅里叶变换的距离维幅频特性图; <BR>图3为仿真实验对混合信号进行短时傅里叶变换的二维时频谱图; <BR>图4为仿真实验对混合信号进行短时傅里叶变换的三维时频谱图; <BR>图5为仿真实验对提取的干扰信号进行短时傅里叶变换的二维时频谱图; <BR>图6为仿真实验对提取的干扰信号进行短时傅里叶变换的三维时频谱图; <BR>图7为仿真实验对剔除干扰信号后的有用信号进行短时傅里叶变换的二维时频谱图; <BR>图8为仿真实验对剔除干扰信号后的有用信号进行短时傅里叶变换的三维时频谱图。 <BR>具体实施方式 <BR>为了让本发明的上述和其它目的、特征及优点能更明显,下文特举本发明实施例,并配合所附图示,做详细说明如下。 <BR>本发明方法的基本思路是:获得的回波数据后,首先,对回波数据进行窄带干扰定性识别;其次,对存在干扰的回波数据进行短时傅里叶变换得到相应的二、三维的混合信号时频谱图,并对其进行奇异值分解,获得对应的待分离信号矩阵;接着,利用窄带干扰信号与有用信号在时频域的差异性和独立性,采用盲源分离方法对得到的存在干扰的待分离信号矩阵进行分离,得到干扰信号矩阵;之后,采用最小二乘方法重构对应的二、三维的干扰信号分量的时频谱图,并与混合信号时频谱图进行相减,得到有用回波信号时频谱图;最后,对有用回波信号时频谱图进行短时傅里叶变换得到消去干扰信号后的有用回波信号,实现SAR时变窄带干扰抑制。 <BR>以上是本发明的核心思想,下面将结合本发明实施例中的附图,对本发明实施例的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例?;诒痉⒚髦械氖凳├?,本领域普通技术人员在没有做出创造性劳动前提下获得的所有其他实施例,都属于本发明?;さ姆段?。 <BR>在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是本发明还可以采用其它不同于在此描述的其它方式来实施,本领域技术人员可以在不违背本发明内涵的情况下做类似推广,因此本发明不受下面公开的具体实施例的限制。 <BR>参照图2,图2为本发明方法的流程图,本发明方法的步骤如下: <BR>步骤1、获取回波数据; <BR>雷达以脉冲重复频率发射并接收脉冲,得到以距离为行向量、以方位为列向量的回波数据; <BR>步骤2、判断回波数据是否存在窄带干扰; <BR>本发明可优选通过一维距离像的峭度来判断回波数据是否存在干扰,判断过程为:对回波数据的列向量进行傅里叶变换,得到一维距离像,判断一维距离像的峭度κ是否大于等于设定的幅度阈值,若是,则认为回波数据的列向量中存在窄带干扰,执行步骤3,否则执行步骤7;优选的,前述幅度阈值为10~200; <BR>一维距离像的峭度 <BR>上式中的E(·)表示求期望运算操作,y为一维距离像的幅度值,μ为一维距离像幅值的均值,d为一维距离像幅值的标准差; <BR>步骤3、根据存在窄带干扰的回波数据的列向量获得混合信号时频谱图S,并对混合信号时频谱图S进行奇异值分解,获得对应的待分离信号矩阵Y; <BR>对回波数据的列向量进行短时傅里叶变换,得到混合信号时频谱图S,并按照下式对混合信号时频谱图S进行奇异值分解: <BR>S=UΛVH <BR>上式中的U为混合信号时频谱图的左特征向量矩阵,Λ为奇异值对角矩阵,V为混合信号时频谱图的右特征向量矩阵,(·)H表示矩阵共轭转置操作; <BR>将分解得到的混合信号时频谱图的左特征向量矩阵U的前L个特征向量作为行向量,构成待分离信号矩阵Y,其中的L为主分量个数,主分量个数L根据最小描述长度准则计算得到: <BR><MATHS num="0005"><MATH><![CDATA[ <mrow><MI>L</MI> <MO>=</MO> <MI>arg</MI> <MI>min</MI> <MO>{</MO> <MO>-</MO> <MI>m</MI> <MROW><MO>(</MO> <MI>n</MI> <MO>-</MO> <MI>L</MI> <MO>)</MO> </MROW><MI>log</MI> <MROW><MO>(</MO> <MROW><MO>(</MO> <MI>n</MI> <MO>-</MO> <MI>L</MI> <MO>)</MO> </MROW><MSUP><MROW><MO>(</MO> <MUNDEROVER><MROW><MI>Π</MI> <MI></MI></MROW><MROW><MI>p</MI> <MO>=</MO> <MI>L</MI> <MO>+</MO> <MN>1</MN> </MROW><MI>n</MI> </MUNDEROVER><MSUB><MI>σ</MI> <MI>p</MI> </MSUB><MO>)</MO> </MROW><MFRAC><MN>1</MN> <MROW><MI>n</MI> <MO>-</MO> <MI>L</MI> </MROW></MFRAC></MSUP><MSUP><MROW><MO>(</MO> <MUNDEROVER><MROW><MI>Σ</MI> <MI></MI></MROW><MROW><MI>p</MI> <MO>=</MO> <MI>L</MI> <MO>+</MO> <MN>1</MN> </MROW><MI>n</MI> </MUNDEROVER><MSUB><MI>σ</MI> <MI>p</MI> </MSUB><MO>)</MO> </MROW><MROW><MO>-</MO> <MN>1</MN> </MROW></MSUP><MO>)</MO> <MO>+</MO> <MN>0.5</MN> <MI>L</MI> <MROW><MO>(</MO> <MN>2</MN> <MI>n</MI> <MO>-</MO> <MI>L</MI> <MO>)</MO> </MROW><MI>log</MI> <MROW><MO>(</MO> <MI>m</MI> <MO>)</MO> </MROW></MROW><MO>}</MO> </MROW>]]&gt;</MATH></MATHS> <BR>上式中的argmin(·)表示计算使目标函数取最小值时对应变量值的运算操作,m为混合信号时频谱图的左特征向量矩阵U的行数,n为混合信号时频谱图的左特征向量矩阵U的列数,log(·)表示求对数运算操作,表示对第L+1至第n项的求乘积运算操作,σp表示奇异值对角矩阵Λ的第p个奇异值,表示对第L+1至第n项的求和运算操作; <BR>步骤4,对待分离信号矩阵Y进行分离,获得干扰信号矩阵I; <BR>窄带干扰信号与有用信号在时频域具有差异性和独立性,利用该差异性和独立性可采用盲源分离法进行分离,得到干扰信号矩阵I,其包括以下子步骤; <BR>子步骤4a、随机产生大小为L×L的解混矩阵W; <BR>子步骤4b、计算迭代矩阵z; <BR>迭代矩阵z=[z1,z2,…,zp,…,zL],其中,zp为迭代矩阵z的第p个列向量,<MATHS num="0006"><MATH><![CDATA[ <mrow> <MSUB><MI>z</MI> <MI>p</MI> </MSUB><MO>=</MO> <MI>E</MI> <MO>{</MO> <MI>Y</MI> <MSUP><MROW><MO>(</MO> <MSUBSUP><MI>w</MI> <MI>p</MI> <MI>H</MI> </MSUBSUP><MI>Y</MI> <MO>)</MO> </MROW><MO>*</MO> </MSUP><MSUP><MROW><MO>(</MO> <MN>0.1</MN> <MO>+</MO> <MSUP><MROW><MO>|</MO> <MSUBSUP><MI>w</MI> <MI>p</MI> <MI>H</MI> </MSUBSUP><MI>Y</MI> <MO>|</MO> </MROW><MN>2</MN> </MSUP><MO>)</MO> </MROW><MROW><MO>-</MO> <MN>1</MN> </MROW></MSUP><MO>}</MO> <MO>-</MO> <MI>E</MI> <MO>{</MO> <MN>0.1</MN> <MSUP><MROW><MO>(</MO> <MN>0.1</MN> <MO>+</MO> <MSUP><MROW><MO>|</MO> <MSUBSUP><MI>w</MI> <MI>p</MI> <MI>H</MI> </MSUBSUP><MI>Y</MI> <MO>|</MO> </MROW><MN>2</MN> </MSUP><MO>)</MO> </MROW><MROW><MO>-</MO> <MN>2</MN> </MROW></MSUP><MO>}</MO> <MSUB><MI>w</MI> <MI>p</MI> </MSUB><MO>,</MO> </MROW>]]&gt;</MATH></MATHS>E(·)表示求期望运算操作,Y为待分离信号矩阵,wp为解混矩阵W的第p个列向量,(·)H表示矩阵共轭转置操作,(·)*表示共轭运算操作,|·|表示求绝对值运算操作; <BR>子步骤4c、对解混矩阵和迭代矩阵的所有列向量依次判断是否满足若所有列向量都满足该条件,则执行子步骤4d,否则,对不满足条件的wp值进行更新,令所述执行子步骤4b;其中,||·||表示求模值运算操作,zp为迭代矩阵z的第p个列向量,wp为解混矩阵w的第p个列向量,ε表示设定的阈值,一般取1e?3; <BR>子步骤4d、将解混矩阵w各列向量去相关得到最终解混矩阵<MATHS num="0007"><MATH><![CDATA[ <mrow> <MOVER><MI>W</MI> <MO>~</MO> </MOVER><MO>=</MO> <MI>W</MI> <MSUP><MROW><MO>(</MO> <MSUP><MI>W</MI> <MI>H</MI> </MSUP><MI>W</MI> <MO>)</MO> </MROW><MROW><MO>-</MO> <MN>1</MN> <MO>/</MO> <MN>2</MN> </MROW></MSUP><MO>;</MO> </MROW>]]&gt;</MATH></MATHS> <BR>子步骤4e、根据最终解混矩阵计算干扰信号矩阵 <BR>步骤5、根据最小二乘准则,将干扰信号矩阵I中的干扰信号分量按照下式重构干扰信号时频谱图: <BR><MATHS num="0008"><MATH><![CDATA[ <mrow><MSUB><MI>S</MI> <MI>k</MI> </MSUB><MO>=</MO> <MSUB><MI>i</MI> <MI>k</MI> </MSUB><MROW><MO>(</MO> <MSUBSUP><MI>i</MI> <MI>k</MI> <MO>+</MO> </MSUBSUP><MI>S</MI> <MO>)</MO> </MROW></MROW>]]&gt;</MATH></MATHS> <BR>其中,Sk表示第k个干扰信号分量的时频谱图,ik表示干扰信号矩阵I的第k个行向量,+表示求矩阵的广义逆运算操作,S为混合信号时频谱图; <BR>步骤6、从混合信号时频谱图S中剔除干扰信号得到有用信号时频谱图S′,并对有用信号时频谱图S′进行逆短时傅里叶变换,得到剔除干扰后的回波数据的列向量; <BR>有用信号时频谱图其中,|·|表示求绝对值运算操作,表示对第1项至第L项的求和运算操作,Sk为第k个干扰信号分量的时频谱图,exp(·)表示求指数运算操作,j为虚数符号,Θ为对混合信号时频谱图S求得的相位角; <BR>步骤7,判断是否遍历完所有回波数据的列向量,若遍历完则执行步骤8,否则执行步骤2; <BR>步骤8,获得剔除干扰后的回波数据; <BR>前述步骤中1≤p≤L,1≤k≤L。 <BR>本发明方法可以弥补独立分量分析法导致样本缺失、复数经验模态分解法对干扰信号分离不彻底的不足,充分利用雷达信号与干扰信号时频特性的差异性,利用时频谱图分解方法直接在信号时频域提取窄带干扰并进行抑制,能够保留有用回波信号,有效地抑制时变的窄带干扰。 <BR>本发明的效果可以通过以下的仿真实验进一步说明,仿真时采用MATLAB(R2010b)软件进行仿真。 <BR>仿真条件如下: <BR>由于雷达回波信号可以近似为线性调频信号,因此,仿真观测信号由一个线性调频信号,一个窄带干扰信号和噪声混合组成,线性调频信号的调频率为7×1012Hz/s,干扰信号载频fc为8MHz,采样频率为150MHz,观测时间为[?3.4,3.4]μs,信噪比为10dB。 <BR>仿真内容 <BR>如图2所示,对混合信号进行傅里叶变换后得到图2所示的距离维幅频特性图,其中横坐标为频率单元,纵坐标为归一化幅度值。从图2可看出,由于窄带干扰的存在,线性调频信号频谱中存在明显的带内冲激干扰,在一个小频率段中幅值出现了突变。 <BR>如图3和图4所示,对混合信号进行短时傅里叶变换后,得到图3所示的混合信号二维时频谱图和图4所示的混合信号三维时频谱图。图3的横坐标为时间单元,纵坐标为频率单元,由图3可看出,斜线为线性调频信号的时频分布,水平线为干扰分量的时频分布,干扰分量的存在影响了对线性调频信号的时频特性的辨识;图4的X坐标为频率单元,Y坐标为时间单元,Z坐标为幅度值,由图4可看出,突起为干扰分量的时频分布,底部隆起的尖峰为线性调频信号的时频分布,干扰分量能量很强,压制了线性调频信号的时频特征,从而影响了对线性调频信号的时频特性的辨识。 <BR>如图5和图6所示,采用时频谱图分解方法提取窄带干扰分量得到重构的如图5所示的干扰信号二维时频谱图和如图6所示的干扰信号三维时频谱图。图5的横坐标为时间单元,纵坐标为频率单元,由图5可看出,水平线为干扰分量的时频分布,与图3相比,提取的干扰信号,只保留了干扰信号的时频信息,没有混杂线性调频信号的时频信息;图6的X坐标为频率单元,Y坐标为时间单元,Z坐标为幅度值,由图6可看出,突起为高能量干扰信号的时频分布,与图4相比,提取的干扰信号,只保留了干扰信号的时频信息,并很好的刻画了时频谱图的幅度强度和时变特性。 <BR>如图7和图8所示,将重构的干扰信号从混合信号中剔除后得到图7所示的有用信号二维时频谱图和图8所示的有用信号三维时频谱图。图7的横坐标为时间单元,纵坐标为频率单元,由图7可看出,斜线为线性调频信号的时频分布,被恢复出来的信号在时频域中已没有与窄带干扰的交叉项,只保留了原线性调频信号的时频信息,由于高能量的干扰信号被去除,噪声的时频特征也可清晰可见;图8的X坐标为频率单元,Y坐标为时间单元,Z坐标为幅度值,由图8可看出,突出的斜线表示恢复的线性调频信号的时频变化信息,底部的随机起伏表示噪声的时频特征,高能量的干扰信号被清除,恢复的线性调频信号时频信息被很好的保留。通过仿真实验说明了本方法的有效性。 <BR>以上所述,仅是本发明的较佳实施例而已,并非对本发明做任何形式上的限制,虽然本发明已以较佳实施例揭露如上,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的技术内容做出些许更动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。</p> </div> </div> </div> </div> <div class="mt10 works-comment"> <div class="works-comment-hd"> <span class="font-tahoma">关于本文</div> <div style="line-height: 25px; padding: 10px 20px;"> 本文标题:基于时频谱图分解的SAR时变窄带干扰抑制方法.pdf<br /> 链接地址:<a href="//www.4mum.com.cn/p-5778597.html">//www.4mum.com.cn/p-5778597.html</a><br /> </div> </div> </div> <div class="boxright" id="boxright" > <div class="fr detail-aside" id="Div11" style="width:270px;"> <div class="box hot-keywords mt10" style="overflow: hidden;width: 268px; border:solid 1px #dedede;" id="relatebox0"> <div class="boxHd" > <div class="fl nt-ico mr5 ml13 ico" style="margin-top: 3px;"> </div> <h2 class="fl font-normal font16 font-yahei" style="font-size: 16px; font-weight: 100; margin-left: 0px; margin-top: 3px; font-family: 微软雅黑"> 当前资源信息</h2> </div> <div id="Div2" class="author-works-list bgF" style="overflow: hidden; padding:10px 10px; "> <table><tr><td> <dt class="author-avatar-box fl"><a class="author-avatar" title="bbs3" href="u-248.html"> <img src="//www.4mum.com.cn/FlexPaper/ZoomImage/UploadPhoto/tx_20180929235856.png" style="border-radius:50px 50px; " onerror="this.src='images/noavatar_small.gif'" alt="bbs3"></a> </dt></td><td> <div class="author-name fl w100 ellipsis"> <a href="u-248.html" target="_blank"> bbs3</a><img style="height:15px; width:20px; overflow:hidden; margin-right:10px;background:url(images/bg_index_ie6_781d95ab.png) no-repeat -185px 4px;display:none" src="images/s.gif" alt="企业认证" title="企业认证"/></div> <div class="author-level-bar"> <span class="mr5 author-grade author-grade5" title="会员等级"></span> </div></td></tr></table> <div class="fl" style="width: 240px;"> <p class="kh_cpZl"> 编号: cj20190318172904525837</p> <p class="kh_cpZl"> 类型: 共享资源</p> <p class="kh_cpZl"> 格式: PDF</p> <p class="kh_cpZl"> 大?。?7.99MB</p> <p class="kh_cpZl"> 上传时间: 2019-03-18</p> </div> </div> </div> <div class="box hot-keywords mt10" style="overflow: hidden;width: 268px; border:solid 1px #dedede;" id="relatebox"> <div class="boxHd" style="padding-bottom: 0px;"> <div class="fl keywords-ico mr5 ml13 ico"> </div> <h2 class="fl font-normal font16 font-yahei" style="font-size: 16px; font-weight: 100; margin-left: 0px; margin-top: 3px; font-family: 微软雅黑"> 相关资源</h2> </div> <div id="author-works-list" class="author-works-list bgF"> <li class="vprdn"> <img alt="枪支在位检测装置.pdf" class="pdf" src="Images/s.gif" /><a target="_parent" href="//www.4mum.com.cn/p-7326895.html" title="枪支在位检测装置.pdf">枪支在位检测装置.pdf</a> </li> <li class="vprdn"> <img alt="便捷装拆式竖向送枪支架.pdf" class="pdf" src="Images/s.gif" /><a target="_parent" href="//www.4mum.com.cn/p-7326824.html" title="便捷装拆式竖向送枪支架.pdf">便捷装拆式竖向送枪支架.pdf</a> </li> <li class="vprdn"> <img alt="一种制导弹药电性能测试仿真试验平台.pdf" class="pdf" src="Images/s.gif" /><a target="_parent" href="//www.4mum.com.cn/p-7326533.html" title="一种制导弹药电性能测试仿真试验平台.pdf">一种制导弹药电性能测试仿真试验平台.pdf</a> </li> <li class="vprdn"> <img alt="一种制导弹药舵机放大器仿真测试装置.pdf" class="pdf" src="Images/s.gif" /><a target="_parent" href="//www.4mum.com.cn/p-7326506.html" title="一种制导弹药舵机放大器仿真测试装置.pdf">一种制导弹药舵机放大器仿真测试装置.pdf</a> </li> <li class="vprdn"> <img alt="北斗卫星枪支动态定位系统.pdf" class="pdf" src="Images/s.gif" /><a target="_parent" href="//www.4mum.com.cn/p-7326502.html" title="北斗卫星枪支动态定位系统.pdf">北斗卫星枪支动态定位系统.pdf</a> </li> <li class="vprdn"> <img alt="制导弹药多路时序状态信号测量电路及其测量方法.pdf" class="pdf" src="Images/s.gif" /><a target="_parent" href="//www.4mum.com.cn/p-7326480.html" title="制导弹药多路时序状态信号测量电路及其测量方法.pdf">制导弹药多路时序状态信号测量电路及其测量方法.pdf</a> </li> <li class="vprdn"> <img alt="RF集成电路的测试方法和系统.pdf" class="pdf" src="Images/s.gif" /><a target="_parent" href="//www.4mum.com.cn/p-6421345.html" title="RF集成电路的测试方法和系统.pdf">RF集成电路的测试方法和系统.pdf</a> </li> <li class="vprdn"> <img alt="电力变压器绕组故障模拟试验方法.pdf" class="pdf" src="Images/s.gif" /><a target="_parent" href="//www.4mum.com.cn/p-6421258.html" title="电力变压器绕组故障模拟试验方法.pdf">电力变压器绕组故障模拟试验方法.pdf</a> </li> <li class="vprdn"> <img alt="磁传感器.pdf" class="pdf" src="Images/s.gif" /><a target="_parent" href="//www.4mum.com.cn/p-6421149.html" title="磁传感器.pdf">磁传感器.pdf</a> </li> <li class="vprdn"> <img alt="一种利用二次反射波测量超声波渡越时间的电路.pdf" class="pdf" src="Images/s.gif" /><a target="_parent" href="//www.4mum.com.cn/p-6421130.html" title="一种利用二次反射波测量超声波渡越时间的电路.pdf">一种利用二次反射波测量超声波渡越时间的电路.pdf</a> </li> <li class="vprdn"> <img alt="接收器与通信信号的接收方法.pdf" class="pdf" src="Images/s.gif" /><a target="_parent" href="//www.4mum.com.cn/p-6420702.html" title="接收器与通信信号的接收方法.pdf">接收器与通信信号的接收方法.pdf</a> </li></div> </div> <div class="box hot-keywords mt10" id="box3" style="overflow: hidden;width: 268px; border:solid 1px #dedede;"> <div class="boxHd" style="border: none;padding-bottom: 0px;"> <div class="fl keywords-ico mr5 ml13 ico"> </div> <h2 class="fl font-normal font16 font-yahei" style="font-size: 16px; font-weight: 100; margin-left: 0px; margin-top: 3px; font-family: 微软雅黑"> 相关搜索</h2> </div> <div class="hot-keywords-list"> <a target="_blank" href="search.html?p=0&q=%e5%9f%ba%e4%ba%8e" class="tag-item ico" title="基于" hidefocus="true"><span class="ico"><em> 基于</em></span></a> <a target="_blank" href="search.html?p=0&q=%e9%a2%91%e8%b0%b1" class="tag-item ico" title="频谱" hidefocus="true"><span class="ico"><em> 频谱</em></span></a> <a target="_blank" href="search.html?p=0&q=%e5%88%86%e8%a7%a3" class="tag-item ico" title="分解" hidefocus="true"><span class="ico"><em> 分解</em></span></a> <a target="_blank" href="search.html?p=0&q=SAR" class="tag-item ico" title="SAR" hidefocus="true"><span class="ico"><em> SAR</em></span></a> <a target="_blank" href="search.html?p=0&q=%e7%aa%84%e5%b8%a6" class="tag-item ico" title="窄带" hidefocus="true"><span class="ico"><em> 窄带</em></span></a> <a target="_blank" href="search.html?p=0&q=%e5%b9%b2%e6%89%b0" class="tag-item ico" title="干扰" hidefocus="true"><span class="ico"><em> 干扰</em></span></a> <a target="_blank" href="search.html?p=0&q=%e6%8a%91%e5%88%b6" class="tag-item ico" title="抑制" hidefocus="true"><span class="ico"><em> 抑制</em></span></a> <a target="_blank" href="search.html?p=0&q=%e6%96%b9%e6%b3%95" class="tag-item ico" title="方法" hidefocus="true"><span class="ico"><em> 方法</em></span></a> <input name="ctl00$Content$hiddenCategoryID" type="hidden" id="Content_hiddenCategoryID" value="117" /> </div> </div> </div> </div> <!--end boxright--> <!--end layout01--> </div> <!--end mainpart--> </div> <div class="aspNetHidden"><input type="hidden" name="__VIEWSTATEGENERATOR" id="__VIEWSTATEGENERATOR" value="541A68AF" /></div></form> <div style="display: none" id="Div3"> </div> <div id="jubao" class="jubao" style="display: none;"> <div class="jubaoTitle"> 举报类型</div> <ul class="vprdn"> <li class="vprdn"> <input type="radio" id="jubaoSel12" name="jubaoSel12" value="12" />&nbsp;&nbsp; <label for="jubaoSel12">广告或垃圾信息</label> </li> <li class="vprdn"> <input type="radio" id="jubaoSel13" name="jubaoSel13" value="13" />&nbsp;&nbsp; <label for="jubaoSel13">色情、淫秽、低俗信息</label> </li> <li class="vprdn"> <input type="radio" id="jubaoSel14" name="jubaoSel14" value="14" />&nbsp;&nbsp; <label for="jubaoSel14">反政府、反人类、反社会等反动信息</label> </li> <li class="vprdn"> <input type="radio" id="jubaoSel15" name="jubaoSel15" value="15" />&nbsp;&nbsp; <label for="jubaoSel15">散布赌博、暴力、凶杀、恐怖或者教唆犯罪等信息</label> </li> <li class="vprdn"> <input type="radio" id="jubaoSel16" name="jubaoSel16" value="16" />&nbsp;&nbsp; <label for="jubaoSel16">侮辱、诽谤等人身攻击信息</label> </li> <li class="vprdn"> <input type="radio" id="jubaoSel17" name="jubaoSel17" value="17" />&nbsp;&nbsp; <label for="jubaoSel17">散布谣言、扰乱社会秩序,破坏社会稳定等信息</label> </li> <li style="text-align: center; height: 50px; margin-top: 10px;"> <input type="hidden" id="hidTypeID" /> <input type="button" onclick="saveJuBao()" value="提 交" />&nbsp;&nbsp;&nbsp;&nbsp;<input type="button" onclick="$('#jubao').hide();" value="取 消" /></li> </ul> </div> <div id="jubaoResponse" class="jubao" style="display: none;"> <div class="jubaoTitle"> 提示</div> <ul class="vprdn"> <li class="vprdn"> <h4 style="text-align: center; color: rgb(83, 136, 35);"> </h4> </li> <li style="text-align: center; height: 50px; margin-top: 10px;"> <input type="button" onclick="$('#jubaoResponse').hide();" value="关 闭" /> </ul> </div></div> <div class="tempdiv" style="display:none;line-height:0px;height:0px; overflow:hidden;"> </div> <!-- JiaThis Button END --> <span id="LabelScript"></span> </div> </div> <!--foot--><div class="bg_100 foot_nav_bg" style=" min-width:1200px;"> <div class="foot_nav"> <a href="//www.4mum.com.cn/h-33.html" target="_blank" rel="nofollow">关于我们</a> - <a href="//www.4mum.com.cn/h-34.html" target="_blank" rel="nofollow">网站声明</a> - <a href="//www.4mum.com.cn/h-35.html" target="_blank" rel="nofollow"> 网站地图</a> - <a href="//www.4mum.com.cn/sitemap.html" target="_blank" rel="nofollow"> 资源地图</a> - <a href="//www.4mum.com.cn/friend.aspx" target="_blank" rel="nofollow">友情链接</a> - <a rel="nofollow">网站客服</a> - <a href="//www.4mum.com.cn/h-93.html" target="_blank" rel="nofollow">联系我们</a> </div></div><div class="bg_100 siteInner_bg" style=" min-width:1200px;"> <div class="siteInner"> <p style="text-align: center; line-height: 1.75em;"><span style="font-size: 14px;"><span style="font-size: 14px; color: rgb(102, 102, 102); font-family: 微软雅黑, Arial, &#39;Times New Roman&#39;; line-height: 20px; text-align: center;">copyright@ 2017-2018 www.4mum.com.cn网站版权所有</span><br style="color: rgb(102, 102, 102); font-family: 微软雅黑, Arial, &#39;Times New Roman&#39;; font-size: 12px; line-height: 20px; text-align: center; white-space: normal; "/><span style="font-size: 14px; color: rgb(102, 102, 102); font-family: 微软雅黑, Arial, &#39;Times New Roman&#39;; line-height: 20px; text-align: center;">经营许可证编号:粤ICP备17046363号-1&nbsp;</span></span> &nbsp;</p><p><br/></p> </div></div><!--foot end--> <!-- 代码部分begin --><div class="QQ_S" style="height: 172px;position:fixed;right: 0px;bottom: 20px; top:auto;"> <div class="Q_top" onclick="HideFoot()"> <span class="signi"></span>收起</div> <div class="Q_botm"> <div class="Q_pic"> <div class="Q_pic_hide"> <a target="_blank" title="在线客服" rel="nofollow"><span class="hide_pic"></span>在线客服</a> </div> </div> <div class="Q_anser"> <div class="Q_anser_hide"><a target="_blank" title="意见反馈"> <span class="hide_pic1"></span>意见反馈 </a> </div> </div> <div class="Q_sign"> <div class="Q_sign_hide" onclick="backToTop();"><a href="javascript:void(0)" title="返回顶部"><span class="hide_pic2"></span>返回顶部 </a></div> </div> </div></div><div class="QQ_S1"> <div class="Q_top1" onclick="ShowFoot()"> <span class="signj"></span>展开</div> <div class="Q_botm1"> <div class="Q_pic1"> <div class="Q_pic1_hide"> <a href="//www.4mum.com.cn/" title="重庆时时彩单双窍门">重庆时时彩单双窍门</a></div> </div> <div class="Q_sign1"> <div class="Q_sign1_hide" onclick="backToTop();"><a href="javascript:void(0)">返回顶部</a></div> </div> </div></div> <!-- 代码部分end --> <li><a href="//www.4mum.com.cn/77d944/69034815.html ">四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 </a> 2019-05-13</li> <li><a href="//www.4mum.com.cn/euy29/3433514.html ">银保监会新规剑指大企业多头融资和过度融资 </a> 2019-05-12</li> <li><a href="//www.4mum.com.cn/a62999/3431813.html ">韩国再提4国联合申办世界杯 中国网友无视:我们自己来 </a> 2019-05-11</li> <li><a href="//www.4mum.com.cn/4kf294/98434412.html ">中国人为什么一定要买房? </a> 2019-05-11</li> <li><a href="//www.4mum.com.cn/33r392/76278011.html ">十九大精神进校园:风正扬帆当有为&nbsp;勇做时代弄潮儿 </a> 2019-05-10</li> <li><a href="//www.4mum.com.cn/1jd106/47417610.html ">粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 </a> 2019-05-09</li> <li><a href="//www.4mum.com.cn/997397/2262689.html ">太原设禁鸣路段 设备在测试中 </a> 2019-05-09</li> <li><a href="//www.4mum.com.cn/br3559/3646538.html ">拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 </a> 2019-05-08</li> <li><a href="//www.4mum.com.cn/9h7915/8581547.html ">“港独”没出路!“梁天琦们”该醒醒了 </a> 2019-05-07</li> <li><a href="//www.4mum.com.cn/3n3511/2345366.html ">陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 </a> 2019-05-06</li> <li><a href="//www.4mum.com.cn/b0059/6332815.html ">人民日报客户端辟谣:“合成军装照”产品请放心使用 </a> 2019-05-05</li> <li><a href="//www.4mum.com.cn/v5j888/7405464.html ">【十九大·理论新视野】为什么要“建设现代化经济体系”?&nbsp;&nbsp; </a> 2019-05-04</li> <li><a href="//www.4mum.com.cn/pl1926/4722623.html ">聚焦2017年乌鲁木齐市老城区改造提升工程 </a> 2019-05-04</li> <li><a href="//www.4mum.com.cn/14g861/9039202.html ">【专家谈】上合组织——构建区域命运共同体的有力实践者 </a> 2019-05-03</li> <li><a href="//www.4mum.com.cn/7jn906/5662951.html ">【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! </a> 2019-05-03</li> <script>(function(){ var src = (document.location.protocol == "") ? "//js.passport.qihucdn.com/11.0.1.js?9ed1f3a8f9c3ff069b7b95c01474c743":"https://jspassport.ssl.qhimg.com/11.0.1.js?9ed1f3a8f9c3ff069b7b95c01474c743"; document.write('<script src="' + src + '" id="sozz"><\/script>'); })(); </script> <script> (function(){ var bp = document.createElement('script'); var curProtocol = window.location.protocol.split(':')[0]; if (curProtocol === 'https') { bp.src = 'https://zz.bdstatic.com/linksubmit/push.js'; } else { bp.src = '//push.zhanzhang.baidu.com/push.js'; } var s = document.getElementsByTagName("script")[0]; s.parentNode.insertBefore(bp, s); })(); </script> <a href="//www.0908491.live">广西11选5计算器</a> <a href="//www.888393.live">微信捕鱼游戏大全免费下载</a> <a href="//www.116103.live">六合图库app下载</a> <a href="//www.3171339.live">澳洲幸运10的开奖查询</a> <a href="//www.938819.live">双色球复式中奖计算器</a> <a href="//www.646776.live">快乐时时彩开奖记录查询结果</a> <a href="//www.510707.live">陕西快乐十分今日开奖</a> <a href="//150143.live">股票涨跌百分比怎么看</a> <a href="//www.762427.live">福彩3d和值走势图带连线</a> <a href="//www.6964028.live">双色球预测最准确人</a> <a href="//www.3354746.live">pc蛋蛋怎么玩</a> <a href="//www.5460528.live">黑龙江11选5投注技巧</a> <a href="//www.216899.live">高频彩种实时开奖结果</a> <a href="//www.443198.live">鸿运线上娱乐官网</a> <a href="//www.165626.live">91v千炮捕鱼</a> <a href="//www.642981.live">麻将来了哪种赢分快</a> </body></html>