• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 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
    • / 13
    • 下载费用:30 金币  

    重庆时时彩赔率: 一种基于匹配自相关算法的雷达信号TOA估计方法.pdf

    关 键 词:
    一种 基于 匹配 相关 算法 雷达 信号 TOA 估计 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201710034262.3

    申请日:

    2017.01.18

    公开号:

    CN106772271A

    公开日:

    2017.05.31

    当前法律状态:

    实审

    有效性:

    审中

    法律详情: 实质审查的生效IPC(主分类):G01S 7/35申请日:20170118|||公开
    IPC分类号: G01S7/35 主分类号: G01S7/35
    申请人: 哈尔滨工业大学
    发明人: 王勇; 马淑歌; 谢俊好; 李绍滨; 李高鹏; 张庆祥
    地址: 150001 黑龙江省哈尔滨市南岗区西大直街92号
    优先权:
    专利代理机构: 哈尔滨市松花江专利商标事务所 23109 代理人: 杨立超
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201710034262.3

    授权公告号:

    |||

    法律状态公告日:

    2017.06.23|||2017.05.31

    法律状态类型:

    实质审查的生效|||公开

    摘要

    一种基于匹配自相关算法的雷达信号TOA估计方法,本发明涉及基于匹配自相关算法的雷达信号TOA估计方法。本发明的目的是为了解决现有雷达信号估计必须建立在信号参数已知的条件下的问题。具体过程为:一、确定雷达信号的自相关函数和双检测门限;通过双检测门限检测得出雷达信号起、止时间的粗估计;二、根据步骤一得出的雷达信号起、止时间的粗估计值对雷达信号未知的参数进行估计,得出发射的雷达信号模型;三、通过发射的雷达信号模型得出匹配滤波的响应;四、根据三得到的匹配滤波的响应对接收的雷达信号进行匹配滤波,得到匹配后的雷达信号,找到匹配后雷达信号的峰值点即为雷达信号TOA。本发明用于雷达信号领域。

    权利要求书

    1.一种基于匹配自相关算法的雷达信号TOA估计方法,其特征在于:一种基于匹配自相
    关算法的雷达信号TOA估计方法具体过程为:
    步骤一、确定雷达信号的自相关函数和双检测门限;通过双检测门限检测得出雷达信
    号起、止时间的粗估计;
    步骤二、根据步骤一得出的雷达信号起、止时间的粗估计值对雷达信号未知的参数进
    行估计,得出发射的雷达信号模型;
    步骤三、通过发射的雷达信号模型得出匹配滤波的响应;
    步骤四、根据步骤三得到的匹配滤波的响应对接收的雷达信号进行匹配滤波,得到匹
    配后的雷达信号,找到匹配后雷达信号的峰值点即为雷达信号TOA;
    所述TOA为雷达信号到达时间。
    2.根据权利要求1所述一种基于匹配自相关算法的雷达信号TOA估计方法,其特征在
    于:所述步骤一中确定雷达信号的自相关函数和双检测门限;通过双检测门限检测得出雷
    达信号起、止时间的粗估计;
    具体过程为:
    步骤一一、确定雷达信号的自相关函数,假设雷达接收到的信号为m(t)=s(t)+n(t),
    其中,m(t)为雷达接收到的信号,s(t)为雷达信号,n(t)为高斯白噪声;
    采用正交采样的方法把雷达接收到的信号m(t)变成复数信号

    其中,A为复数信号幅度,为初始相位,f为载频,L(i)为复噪声,m(i)为复数信号,j为
    虚数单位,j2=-1,i为采样点个数,取值为正整数,△t为采样时间;
    自相关函数为
    <mrow> <mi>k</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <mi>m</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mi>j</mi> <mo>)</mo> </mrow> <msup> <mi>m</mi> <mo>*</mo> </msup> <mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
    将(1)带入(2),得出自相关函数k(i)=NA2e-j2πf△t+L' (3)
    其中,N为样本数,取值为正整数,m*(i+j+1)为m(i+j+1)的共轭函数,NA2e-j2πf△t为雷达
    复数信号部分,L'为自相关后的噪声部分,

    其中,L*(i+j+1)为L(i+j+1)的共轭函数,L(i+j)为复噪声;
    步骤一二、将自相关函数公式(3)通过滑动递推得出
    k(i+1)=k(i)+m(i+1+N)·m*(i+2+N)-m(i+1)·m*(i+2) (5)
    步骤一三、确定检测门限,高斯白噪声的标准差的无偏估计
    <mrow> <mover> <mi>&sigma;</mi> <mo>^</mo> </mover> <mo>=</mo> <msqrt> <mrow> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msubsup> <mi>n</mi> <mi>i</mi> <mn>2</mn> </msubsup> <mo>/</mo> <mrow> <mo>(</mo> <mi>N</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </msqrt> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
    其中,ni为复变换前第i时刻的噪声值,自相关运算后雷达复数信号服从N(0,4Nσ4)分
    布,自相关运算后雷达复数信号模值服从均值为μ,方差为的瑞利分布,并且
    由高斯白噪声的瑞利分布确定第一检测门限为:
    <mrow> <msub> <mi>V</mi> <mi>T</mi> </msub> <mo>=</mo> <mi>&mu;</mi> <mo>+</mo> <msubsup> <mi>a&sigma;</mi> <mn>1</mn> <mn>2</mn> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
    其中,a为检测门限系数,根据第一检测门限,将大于第一检测门限的自相关函数k(i)
    判断为有雷达复数信号,反之,则判断为无雷达复数信号;
    在第一检测门限的基础上进行第二检测门限,即在连续l个数据中有a个超过第一检测
    门限时,则认为有雷达复数信号;其起点为这l个数据中第一次连续有b个点超过第一检测
    门限的时刻;反之,在连续l个数据中有a个未超过第一检测门限时,则认为没有雷达复数信
    号;其终止点为这l个数据中第一次连续有b个点未超过第一检测门限的时刻;a、b为正整
    数,且a大于b。
    3.根据权利要求2所述一种基于匹配自相关算法的雷达信号TOA估计方法,其特征在
    于:所述检测门限系数a取值为2-10。
    4.根据权利要求3所述一种基于匹配自相关算法的雷达信号TOA估计方法,其特征在
    于:所述步骤三中通过发射的雷达信号模型得出匹配滤波的响应;具体过程为:
    步骤三一、确定雷达接收的信号能量,设输入雷达信号为si(t),其频谱为
    <mrow> <mi>F</mi> <mrow> <mo>(</mo> <mi>w</mi> <mo>)</mo> </mrow> <mo>=</mo> <msubsup> <mo>&Integral;</mo> <mrow> <mo>-</mo> <mi>&infin;</mi> </mrow> <mrow> <mo>+</mo> <mi>&infin;</mi> </mrow> </msubsup> <msub> <mi>s</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mi>w</mi> <mi>t</mi> </mrow> </msup> <msub> <mi>d</mi> <mi>t</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>
    则雷达接收信号频谱的傅里叶反变换表述为

    其中,H(w)是雷达接收信号的匹配滤波器的频率响应,雷达接收信号的能量为|s0(t)
    |2,w为频率,t为时间;
    步骤三二、确定雷达接收信号的噪声功率,设雷达接收信号的匹配滤波器输入信号的
    实高斯白噪声功率谱密度为N0/2,则雷达接收信号的匹配滤波器输出的平均功率为

    其中,N0为复高斯白噪声功率谱密度;
    步骤三三、计算雷达接收信号的匹配滤波器的频域响应,雷达接收信号的匹配滤波器
    的输出信噪比td时刻为
    <mrow> <msub> <mi>SNR</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>t</mi> <mi>d</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mfrac> <mrow> <mo>|</mo> <msubsup> <mo>&Integral;</mo> <mrow> <mo>-</mo> <mi>&infin;</mi> </mrow> <mrow> <mo>+</mo> <mi>&infin;</mi> </mrow> </msubsup> <mi>H</mi> <mrow> <mo>(</mo> <mi>w</mi> <mo>)</mo> </mrow> <mi>F</mi> <mrow> <mo>(</mo> <mi>w</mi> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <msub> <mi>jwt</mi> <mi>d</mi> </msub> </mrow> </msup> <msub> <mi>d</mi> <mi>w</mi> </msub> <mo>|</mo> </mrow> <mrow> <mfrac> <msub> <mi>N</mi> <mn>0</mn> </msub> <mn>2</mn> </mfrac> <msubsup> <mo>&Integral;</mo> <mrow> <mo>-</mo> <mi>&infin;</mi> </mrow> <mrow> <mo>+</mo> <mi>&infin;</mi> </mrow> </msubsup> <msup> <mrow> <mo>|</mo> <mi>H</mi> <mrow> <mo>(</mo> <mi>w</mi> <mo>)</mo> </mrow> <mo>|</mo> </mrow> <mn>2</mn> </msup> <msub> <mi>d</mi> <mi>w</mi> </msub> </mrow> </mfrac> <mn>2</mn> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow>
    在最大输出信噪比准则下,得出雷达接收信号的匹配滤波器的频域响应为

    其中F*(w)为F(w)的共轭函数,K为系数,取值为正整数。
    5.根据权利要求4所述一种基于匹配自相关算法的雷达信号TOA估计方法,其特征在
    于:所述步骤四中根据步骤三得到的匹配滤波的响应对接收的雷达信号进行匹配滤波,得
    到匹配后的雷达信号,找到匹配后雷达信号的峰值点即为雷达信号TOA;具体过程为:
    将公式(12)和公式(8)相乘G=H(w)×F(w),对G进行傅里叶反变换,得出峰值点即为雷
    达信号TOA。

    说明书

    一种基于匹配自相关算法的雷达信号TOA估计方法

    技术领域

    本发明涉及基于匹配自相关算法的雷达信号TOA估计方法。

    背景技术

    经典的信号TOA估计算法是基于信号参数已知条件下的匹配滤波器,而传统的匹
    配滤波器是在信噪比最大准则下推导出的一种维纳滤波器,对于输入信号,与其匹配的最
    佳滤波器是与其共轭反向并有一定时间延迟的一个滤波器。对于不同类型的输入信号,只
    要知道信号参数,便可以设计出最佳的匹配滤波器。但是这样的信号TOA估计算法必须建立
    在信号参数已知的条件下,而对于信号参数部分已知的情况则不适用。因此在信号参数部
    分已知条件下,找出一种不完全依赖信号参数的方法,在低信噪比条件下估计信号的TOA显
    得尤为重要。

    发明内容

    本发明的目的是为了解决现有信号估计必须建立在信号参数已知的条件下的问
    题,而提出一种基于匹配自相关算法的信号TOA估计方法。

    一种基于匹配自相关算法的雷达信号TOA估计方法具体过程为:

    步骤一、确定雷达信号的自相关函数和双检测门限;通过双检测门限检测得出雷
    达信号起、止时间的粗估计;

    步骤二、根据步骤一得出的雷达信号起、止时间的粗估计值对雷达信号未知的参
    数进行估计,得出发射的雷达信号模型;

    步骤三、通过发射的雷达信号模型得出匹配滤波的响应;

    步骤四、根据步骤三得到的匹配滤波的响应对接收的雷达信号进行匹配滤波,得
    到匹配后的雷达信号,找到匹配后雷达信号的峰值点即为雷达信号TOA;

    所述TOA为雷达信号到达时间。

    本发明的有益效果为:

    本方法提出一种基于匹配相关算法的雷达信号TOA估计方法,相比于传统算法,本
    算法不需要信号所有的先验知识,算法简单,可实时检测信号。首先需要做一次自相关运
    算,并通过采用双门限检测,提高检测概率,同时粗略估计出信号的截止时间,从而估计信
    号的未知参数,之后再做一次匹配滤波运算,最后便可以精确估计出信号的TOA。如图6可精
    确得出TOA估计为7200;如表1可得出在参数完全未知时,现有技术利用自相关算法估计信
    号TOA测量精度为17.5ns,在参数部分未知时,本发明利用匹配自相关算法估计信号TOA测
    量精度为0ns。

    附图说明

    图1为本发明接收雷达复数信号示意图;

    图2为本发明接收雷达复数信号自相关处理示意图;

    图3为本发明发射雷达复数信号频谱示意图;

    图4为本发明接收雷达复数信号频谱示意图;

    图5为本发明输出雷达复数信号频谱示意图;

    图6为本发明输出雷达复数信号时域示意图。

    具体实施方式

    具体实施方式一:本实施方式的一种基于匹配自相关算法的雷达信号TOA估计方
    法具体过程为:

    步骤一、确定雷达信号的自相关函数和双检测门限;通过双检测门限检测得出雷
    达信号起、止时间的粗估计;

    步骤二、根据步骤一得出的雷达信号起、止时间的粗估计值对雷达信号未知的参
    数进行估计,得出发射的雷达信号模型;

    步骤三、通过发射的雷达信号模型得出匹配滤波的响应;

    步骤四、根据步骤三得到的匹配滤波的响应对接收的雷达信号进行匹配滤波,得
    到匹配后的雷达信号,找到匹配后雷达信号的峰值点即为雷达信号TOA;

    所述TOA为雷达信号到达时间(B雷达发送的信号到达C雷达的时间点)。

    具体实施方式二:本实施方式与具体实施方式一不同的是:所述步骤一中确定雷
    达信号的自相关函数和双检测门限;通过双检测门限检测得出雷达信号起、止时间的粗估
    计;具体过程为:

    步骤一一、确定雷达信号的自相关函数,假设雷达接收到的信号为m(t)=s(t)+n
    (t),

    其中,m(t)为雷达接收到的信号,s(t)为雷达信号,n(t)为高斯白噪声;

    采用正交采样的方法把雷达接收到的信号m(t)变成复数信号


    其中,A为复数信号幅度,为初始相位,f为载频,L(i)为复噪声,m(i)为复数信号,
    j为虚数单位,j2=-1,i为采样点个数(连续信号离散后的点),取值为正整数,△t为采样时
    间;

    自相关函数为


    将(1)带入(2),得出自相关函数k(i)=NA2e-j2πf△t+L' (3)

    其中,N为样本数,取值为正整数,(与雷达信号长度有关),m*(i+j+1)为m(i+j+1)
    的共轭函数,NA2e-j2πf△t为雷达复数信号部分,L'为自相关后的噪声部分,


    其中,L*(i+j+1)为L(i+j+1)的共轭函数,L(i+j)为复噪声;

    步骤一二、将自相关函数公式(3)通过滑动递推得出

    k(i+1)=k(i)+m(i+1+N)·m*(i+2+N)-m(i+1)·m*(i+2) (5)

    从上式可以看出,每计算一个新的k(i)仅需要进行两次复数乘法和两次复数加
    法,因此自相关算法可实时检测信号。

    步骤一三、确定检测门限,为了提高检测性能,本文采用双门限法,第一检测门限
    随噪声变化而变化。高斯白噪声的标准差的无偏估计


    其中,ni为复变换前第i时刻的噪声值,自相关运算后雷达复数信号服从N(0,4Nσ4)
    分布,自相关运算后雷达复数信号模值服从均值为μ,方差为的瑞利分布,并且
    实际中σ由确定。在一定的虚警概率的前提下,可由
    高斯白噪声的瑞利分布确定第一检测门限为:


    其中,a为检测门限系数。

    根据前面估计的第一检测门限,将大于第一检测门限的自相关函数k(i)判断为有
    雷达复数信号,反之,则判断为无雷达复数信号;

    由于噪声是随机的,它可能大于门限,为了提高检测性能,在第一检测门限的基础
    上进行第二检测门限,即在连续l个数据中有a个超过第一检测门限时,则认为有雷达复数
    信号;其起点为这l个数据中第一次连续有b个点超过第一检测门限的时刻。反之,在连续l
    个数据中有a个未超过第一检测门限时,则认为没有雷达复数信号;其终止点为这l个数据
    中第一次连续有b个点未超过第一检测门限的时刻;a、b为正整数,且a大于b;

    信号TOA粗估计的修改,由于粗估计是在相关累加之后进行的,因此需要对粗估计
    的到达时间进行修改。并且这一修改值与信噪比有关,当k(i)首次大于第一检测门限时,并
    不意味雷达复数信号就在i时刻到达。从自相关公式(2)看出,k(i)利用了从i时刻到i+N时
    刻的所有样本,并且k(i)大于第一检测门限时,并不要求所有被利用的样本均含有雷达复
    数信号。当信噪比较低时,需要较多的样本含有雷达复数信号才能使k(i)大于第一检测门
    限,因此雷达复数信号的起始时刻接近i,反之,雷达复数信号的起始时刻则接近i+N。因此
    需要根据信号的信噪比对信号估计的到达时间进行修正。

    其它步骤及参数与具体实施方式一相同。

    具体实施方式三:本实施方式与具体实施方式一或二不同的是:所述检测门限系
    数a由要求的虚警概率确定,一般情况下,a的取值为2-10。

    其它步骤及参数与具体实施方式一或二相同。

    具体实施方式四:本实施方式与具体实施方式一至三之一不同的是:所述步骤三
    中通过发射的雷达信号模型得出匹配滤波的响应;具体过程为:

    匹配滤波是信号检测算法中的传统方法,通常根据信号的先验知识,在信噪比最
    大准则下得出。对于输入信号,与其匹配的最佳滤波器是与其共轭反向并有一定时间延迟
    的滤波器。具体步骤如下:

    步骤三一、确定雷达接收的信号能量,设输入雷达信号为si(t),其频谱为


    则雷达接收信号频谱的傅里叶反变换表述为


    其中,H(w)是雷达接收信号的匹配滤波器的频率响应,雷达接收信号的能量为|s0
    (t)|2,w为频率,t为时间;

    步骤三二、确定雷达接收信号的噪声功率,设雷达接收信号的匹配滤波器输入信
    号的实高斯白噪声功率谱密度为N0/2,则雷达接收信号的匹配滤波器输出的平均功率为


    其中,N0为复高斯白噪声功率谱密度;

    步骤三三、计算雷达接收信号的匹配滤波器的频域响应,雷达接收信号的匹配滤
    波器的输出信噪比td时刻为


    在最大输出信噪比准则下,得出雷达接收信号的匹配滤波器的频域响应为


    其中F*(w)为F(w)的共轭函数,K为系数,取值为正整数。由此便可以根据信号的先
    验知识设计出匹配滤波器。

    其它步骤及参数与具体实施方式一至三之一相同。

    具体实施方式五:本实施方式与具体实施方式一至四之一不同的是:所述步骤四
    中根据步骤三得到的匹配滤波的响应对接收的雷达信号进行匹配滤波,得到匹配后的雷达
    信号,找到匹配后雷达信号的峰值点即为雷达信号TOA;具体过程为:

    将公式(12)和公式(8)相乘G=H(w)×F(w),对G进行傅里叶反变换,得出峰值点即
    为雷达信号TOA。

    其它步骤及参数与具体实施方式一至四之一相同。

    采用以下实施例验证本发明的有益效果:

    实施例一:

    本实施例一种基于匹配相关算法的雷达信号TOA估计方法具体是按照以下步骤制
    备的:

    实验采用载频已知、调频斜率未知的线性调频信号,仿真环境:信号的采样频率
    2.4GHz,调制带宽100GHz,输入信噪比为-3dB,载频100MHz;脉宽3μs,上升沿下降沿脉宽
    300ns。

    可以得出图1-图6的结果。

    从仿真结果可以看出,在低信噪比下,接收信号淹没在在噪声中,经过自相关运
    算,信号从输出波形凸现出来,由此便可以粗略估计信号的起、止时间。通过构造匹配滤波
    器,接收信号经匹配滤波器,信号在输出波形中呈现一个尖峰,以此便可以更加精确的估计
    信号的TOA。

    在上述仿真环境下,分别用匹配自相关算法和自相关算法对其蒙特卡罗仿真,比
    较估计性能,仿真结果如表1所示:

    表1两种TOA估计方法结果的对比表



    仿真分析:经过1000次蒙特卡洛仿真结果可知,在低信噪比条件下,匹配自相关算
    法的测量精度明显优于自相关算法,证明了该发明方法的优越性。因此对信号参数部分已
    知的线性调频信号采用匹配自相关算法更优。

    本发明还可有其它多种实施例,在不背离本发明精神及其实质的情况下,本领域
    技术人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于
    本发明所附的权利要求的?;し段?。

    关于本文
    本文标题:一种基于匹配自相关算法的雷达信号TOA估计方法.pdf
    链接地址://www.4mum.com.cn/p-6001338.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 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
  • 奇趣腾讯分分彩计划软件 时时彩龙虎免费计划 龙虎和马 快乐时时和澳洲一样的吗 吉林时时开奖视频 篮彩竞彩分析预测 必赢客北京pk拾精准版 骰子单双最合理的下法 彩名堂pk10计划软件,在线软件 安徽时时预测软件下载 百盈快三 口诀 双色球彩计划app 时时彩真正稳赚方法 双色球开奖结果下软件 微信大小单双平台 分分彩通用稳赚方案方法技巧总结