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

    重庆时时彩个位预测: 基于脉内干涉的太赫兹宽带雷达微多普勒解模糊方法.pdf

    关 键 词:
    基于 干涉 赫兹 宽带 雷达 多普勒 模糊 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201710167846.8

    申请日:

    2017.03.21

    公开号:

    CN106772308A

    公开日:

    2017.05.31

    当前法律状态:

    实审

    有效性:

    审中

    法律详情: 实质审查的生效IPC(主分类):G01S 7/41申请日:20170321|||公开
    IPC分类号: G01S7/41 主分类号: G01S7/41
    申请人: 中国人民解放军国防科学技术大学
    发明人: 秦玉亮; 邓彬; 游鹏; 杨琪; 王宏强; 罗成高
    地址: 410073 湖南省长沙市开福区德雅路109号
    优先权:
    专利代理机构: 北京中济纬天专利代理有限公司 11429 代理人: 陆薇薇
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201710167846.8

    授权公告号:

    |||

    法律状态公告日:

    2017.06.23|||2017.05.31

    法律状态类型:

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

    摘要

    本发明提供一种基于脉内干涉的太赫兹宽带雷达微多普勒解模糊方法,通过利用宽带雷达解线频调方式(Dechirp)接收信号的特征,基于脉内干涉实现了微多普勒解模糊,具有高效稳健等特点,且同时不增加系统复杂度。

    权利要求书

    1.一种基于脉内干涉的太赫兹宽带雷达微多普勒解模糊方法,其特征在于,包括以下
    步骤:
    步骤S100:采用宽带雷达解线频调方式接收待监测微动目标的宽带回波信号,对所得
    宽带回波信号进行残余视频相位补偿后,得到补偿宽带回波信号;
    步骤S200:逐个抽取所述补偿宽带回波信号中的信号,对所抽取的每个信号的频点进
    行频率干涉,得到多个干涉信号;
    步骤S300:对所述干涉信号进行时频分析,得到解模糊微多普勒曲线,逆Radon变换方
    法处理所述解模糊微多普勒曲线,得到待监测微动目标的微动估计参数。
    2.根据权利要求1所述的基于脉内干涉的太赫兹宽带雷达微多普勒解模糊方法,其特
    征在于,所述步骤S200中还包括对所述干涉信号进行累加的步骤。
    3.根据权利要求2所述的基于脉内干涉的太赫兹宽带雷达微多普勒解模糊方法,其特
    征在于,所述累加步骤中所用干涉差频为B/2。
    4.根据权利要求1~3中任一项所述的基于脉内干涉的太赫兹宽带雷达微多普勒解模
    糊方法,其特征在于,所述时频分析为采用线性时频分布中的短时傅里叶变换方法进行。
    5.根据权利要求4所述的基于脉内干涉的太赫兹宽带雷达微多普勒解模糊方法,其特
    征在于,所述待监测微动目标为旋转目标。

    说明书

    基于脉内干涉的太赫兹宽带雷达微多普勒解模糊方法

    技术领域

    本发明涉及雷达系统设计和雷达信号处理技术领域,具体的涉及一种基于脉内干
    涉的太赫兹宽带雷达微多普勒解模糊方法。

    背景技术

    利用雷达手段实现目标运动特征提取是一种广泛使用的方法,具有全天时、全天
    候、高精度等特点,在遥感测量、军事侦察等方面具有十分重要的应用。太赫兹(Terahertz,
    THz)频段通常指频率在0.1THz到10THz(对应波长30μm-3mm)之间的电磁波,其频率介于毫
    米波与红外光之间,处于宏观电子学向微观光子学的过渡频段,在电磁波谱中占有很特殊
    的位置,具有与其它波段不同的特殊性质。近年来,随着太赫兹源、检测和相关器件的突破,
    太赫兹雷达技术发展迅速,在高分辨成像和运动精细测量方面逐渐显现出优势。

    微动(Micro-motion)指的是目标或目标部件除质心平动以外的振动、转动等微小
    运动,最早由美国海军实验室的V.C.Chen教授提出,由微动引起的多普勒频移称为微多普
    勒(Micro-Doppler)。自然界和生活中的微动现象十分普遍,比如人体的呼吸心跳,直升机
    旋翼的旋转,车辆桥梁等的振动等等,

    然而在实际测量中,根据采样定理,雷达信号的脉冲重复频率(pulse repetition
    frequency,PRF)限定了只有-PRF/2到PRF/2之间的频率能被直接观测到。当多普勒值位于-
    PRF/2到PRF/2之外时,将发生模糊。太赫兹雷达由于其较高的载频,相比常规的微波雷达,
    对目标微动有其敏感性,这是太赫兹雷达的优势;但是,较高的载频也使得在太赫兹频段目
    标微多普勒值极易超出观测范围而模糊。

    同时利用雷达设备对这些旋转目标进行监测有十分重要的意义。但是这些旋转目
    标在运动中由于机械原因或者外部作用力的原因往往会伴随着振动干扰,对这些目标的振
    动情况同样需要进行监测,通过对振动的监测,可以为设备运行状态和安全性提供重要的
    辅助信息。

    利用雷达手段实现目标微动参数估计已经有相当广泛的研究,然而太赫兹频段的
    研究还远远不足。在微波段已有学者对微多普勒模糊问题进行了研究,主要是从雷达系统
    设计和信号处理两个方面进行考虑。

    现有技术中:

    在系统设计方面,主要是进行多通道和多重频测量。多通道是通过将大天线分裂
    为几个小天线或者将几个天线分置于不同的平台来形成多通道,利用多通道数据的联合处
    理来实现多普勒解模糊。多重频是通过需要发射多重频(通常是2~3个)脉冲串以实现模糊
    性。

    在信号处理方面,主要是通过信号域和时频域,利用目标先验信息,采用数字图像
    处理领域的逆Radon变换、Hough变换等方法对模糊的时频曲线进行校正或者补偿以实现微
    多普勒解模糊。

    上述现有微多普勒解模糊方案中,多通道和多重频测量都有其固有问题,且严重
    增加了雷达系统的复杂度,使得系统稳定性和实用性恶化。现有信号处理的手段综合来说
    具有过程复杂、计算量大、普适性差等问题,难以在实际系统中广泛应用。因此,如何快速高
    效的实现微多普勒解模糊,充分发挥太赫兹雷达微多普勒敏感性的优势并同时实现解模糊
    问题仍是一个有待研究的问题,针对这一问题,当前并没有有效的手段。

    发明内容

    本发明的目的在于提供一种基于脉内干涉的太赫兹宽带雷达微多普勒解模糊方
    法,该发明解决了宽带雷达微动目标微多普勒解模糊的技术问题。

    本发明提供一种基于脉内干涉的太赫兹宽带雷达微多普勒解模糊方法,包括以下
    步骤:

    步骤S100:采用宽带雷达解线频调方式接收待监测微动目标的宽带回波信号,对
    所得宽带回波信号进行残余视频相位补偿,得到补偿宽带回波信号;

    步骤S200:逐个抽取补偿宽带回波信号中的信号,对所抽取的每个信号的频点进
    行频率干涉,得到多个干涉信号;

    步骤S300:对干涉信号进行时频分析,得到解模糊微多普勒曲线,并利用逆Radon
    变换方法处理解模糊微多普勒曲线,得到待监测微动目标的微动估计参数。

    进一步地,步骤S200中还包括对干涉信号进行累加的步骤。

    进一步地,累加步骤中所用干涉差频为B/2。

    进一步地,时频分析为采用线性时频分布中的短时傅里叶变换方法进行。

    进一步地,待监测微动目标为旋转目标。

    本发明的技术效果:

    1、本发明提供的基于脉内干涉的太赫兹宽带雷达微多普勒解模糊方法,通过利用
    宽带雷达解线频调方式(Dechirp)接收信号的特征,利用扫频周期内部若干个频点的采样
    信号进行干涉并累加以实现微多普勒解模糊,基于脉内干涉实现了微多普勒解模糊,具有
    高效稳健等特点,且同时不增加系统复杂度。从信号处理领域解决了微动目标微多普勒模
    糊问题,实现了微动参数估计的同时减少对雷达系统复杂度的增加。

    2、本发明提供的基于脉内干涉的太赫兹宽带雷达微多普勒解模糊方法,利用宽带
    雷达解线频调方式接收方式得到微动目标宽带回波信号,对回波信号进行残余视频相位
    (residual video phase,PVP)补偿;然后抽取出宽带信号中两个不同频点的信号进行频率
    干涉,并将所有干涉信号进行累加以消除交叉项的干扰;最后对干涉后的信号进行时频分
    析,得到不模糊的微多普勒曲线,并利用图像处理领域中的逆Radon变换进行处理来估计待
    监测微动目标的微动参数。利用脉内干涉原理,将微多普勒表达式中的载频转换成了干涉
    的差频,急剧减小了微多普勒值,同时又保存了回波中的其它信息,有利于后续分析处理使
    用。

    3.本发明提供的基于脉内干涉的太赫兹宽带雷达微多普勒解模糊方法,能够快速
    稳定的实现微多普勒解模糊,且可应用于多目标场景,在充分发挥太赫兹频段优势的同时,
    解决了多普勒模糊的难题。

    具体请参考根据本发明的基于脉内干涉的太赫兹宽带雷达微多普勒解模糊方法
    提出的各种实施例的如下描述,将使得本发明的上述和其他方面显而易见。

    附图说明

    图1是本发明提供的基于脉内干涉的太赫兹宽带雷达微多普勒解模糊方法的流程
    示意图;

    图2是本发明提供的优选实施例流程示意图;

    图3是本发明提供的优选实施例中所用旋转目标结构示意图;

    图4是本发明提供的优选实施例中未经处理的角反射器回波时频分布示意图,其
    中a)为旋转目标在转速为40r/min时;b)为旋转目标在转速为60r/min时;

    图5是本发明提供的优选实施例中采用本发明提供方法处理后的角反射器回波时
    频分布示意图,其中a)为旋转目标在转速为40r/min时;b)为旋转目标在转速为60r/min时;

    图6是本发明提供的优选实施例中逆Radon变换之后的参数空间结果示意图,其中
    a)为旋转目标在转速为40r/min时;b)为旋转目标在转速为60r/min时。

    具体实施方式

    构成本申请的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实
    施例及其说明用于解释本发明,并不构成对本发明的不当限定。

    参见图1~2,本发明提供的基于脉内干涉的太赫兹宽带雷达微多普勒解模糊方
    法,包括以下步骤:

    步骤S100:采用宽带雷达解线频调方式接收待监测微动目标的宽带回波信号,对
    所得宽带回波信号进行残余视频相位补偿,得到补偿宽带回波信号;

    步骤S200:逐个抽取补偿宽带回波信号中的信号,对所抽取的每个信号的频点进
    行频率干涉,得到多个干涉信号;

    步骤S300:对干涉信号进行时频分析,得到解模糊微多普勒曲线,并利用逆Radon
    变换方法处理解模糊微多普勒曲线,得到旋转目标微动估计参数。

    此处所得解模糊微多普勒曲线是指不模糊的多普勒曲线。该方法可用于处理各类
    目标运动中产生的微动,例如旋转目标。本文中宽带回波信号是指宽带雷达系统获得的信
    号。

    优选的,待监测微动目标为旋转目标。本发明提供的方法用于处理旋转目标时效
    果最好,由于旋转目标往往具有较大的旋转半径或者角速度,因此需要一定的波束宽度才
    能覆盖目标;此外,旋转目标在微波段的微多普勒值一般在几百至几万赫兹的范围,处于一
    个既能清晰观测微动现象又不至于严重混叠的区间。因此,利用微波雷达进行旋转目标参
    数估计和成像具有一定优势。

    本发明提供方法的具体说明如下:

    该方法中所接受的信号可以通过宽带雷达系统获得,具体所用雷达可以为宽带脉
    冲雷达,也可以是宽带调频连续波雷达。所用宽带雷达系统发射的信号具有如下式所示模
    式:


    其中,为距离快时间,tm为方位慢时间,t为总时间,Tp为脉冲雷达脉冲宽度或者调
    频连续波雷达的扫频周期,fc为雷达载频,γ为调频率,j为虚数单位。

    宽带雷达解线频调方式接收的思想是将一个参考距离处的回波信号与目标信号
    进行共轭相乘,假设参考距离为Rref,参考信号可以表示为:


    其中,c为光速。假设在距离雷达R0处有K和微动散射中心,这些散射中心在雷达视
    线上的投影一般可以视作简谐运动,则简谐运动的距离Rk的变化表达式为:


    其中ak,ωk和分别为第K个散射中心的微动幅度、微动角速度和初始相位。则根
    据发射信号表达式(1),所监测到的微动回波写为:


    每个微动散射中心的微多普勒值可以表示为:


    宽带雷达解线频调方式处理后的宽带回波信号形式为:


    其中RΔk=Rk-Rref。

    通过宽带雷达解线频调方式接收得到宽带回波信号之后,根据其表达式(6)可以
    看出,其回波相位中的第一项是与待监测微动目标运动有关的分
    量,第二项和第三项为包络斜置项和残余视频相位,需要对第二
    项和第三项进行补偿。所用补偿方法在宽带雷达信号处理中讨论的比较多,详细可参见[保
    铮,邢孟道,王彤.雷达成像原理[M].电子工业出版社.2008.]基本思路是将式(6)对快时间
    进行傅里叶变换,然后乘以下式(7),最后进行傅里叶逆变换得到公式(8)所示的补偿后宽
    带回波信号表达式,所得结果即为经过补偿的第二项和第三项。


    补偿后的宽带回波信号表达式为:


    步骤S200:脉内干涉处理

    得到如式(8)的宽带回波信号之后,其脉冲内部的每个采样点即相当于一个单频
    下的回波,假设雷达系统每个脉冲的采样数目为N,就相当于有N个单频回波,这N个单频均
    位于fc-B/2与fc+B/2之间,其中B为信号带宽。其中第i个单频信号的表达式为:


    选取第i个与第l个单频信号进行干涉处理,处理之后的信号形式为:


    其中,前面K项为与K个散射中心相关的微多普勒项,因而公式(10)中的
    是需要获取的信息。
    第二项为交叉项,需要进行抑制。

    优选的,为了降低公式(10)中交叉项的影响,所述步骤S200中还包括对所述干涉
    信号进行累加的步骤。具体为:分别取i=1,2…N/2,l=N/2+1,N/2+2…N,将第i个和第l个
    信号共轭相乘得到共N/2个信号向量,并进行累加,经过处理后的信号形式为:



    经过累加处理,能有效提高所得估值的准确度。

    这样就得到了包含K个散射中心的回波信号,且这K个散射中心的微多普勒值表达
    式变为


    其中B为发射信号带宽,一般远小于信号载频fc。将式(12)与式(5)相比可以看出,
    原表达式中比较大的值fc已经被一个较小的值B/2代替,这就相当于将微多普勒值缩小了
    2fc/B倍,这样就不会超出微多普勒的不模糊范围[-PRF/2,PRF/2],因此原本模糊的多普勒
    曲线经过上述处理后不会再发生模糊。从而提高了监测精度。该步骤能消除交叉项的干扰。

    步骤S300:基于时频分析和逆Radon变换的微动参数估计。

    得到消除模糊的微动目标回波信号之后,就可以对其进行时频分析以观察微多普
    勒曲线。常用的时频分析方法有线性时频分布和二次Cohen类时频分布,优选的,为了减小
    交叉项的干扰,本发明采用线性时频分布中的短时傅里叶变换方法(Short Time Fourier
    Transform,STFT),信号s(t)的STFT表达式为:


    其中h(t)为窗函数。在得到时频分布之后,由于微动目标在雷达视线上的投影为
    简谐运动,其微多普勒为正弦形式。根据这一特点,本发明采用逆Radon变换对时频分布进
    行处理,将其转换到参数空间以实现参数估计。逆Radon变换通常被用于是字图像处理领
    域,它可以将位于图像中心的正弦曲线转换成参数空间的特显点,特显点的位置可以推出
    正弦曲线的参数。因而本发明提供方法利用该方法的这一特性,从而通过对多普勒曲线进
    行处理后得到相应的估计参数。逆Radon变换方法可以按现有步骤进行。具体的包括以下步
    骤:

    设平面上一条正弦曲线的二维表达式为:


    其中A为正弦曲线的最大值,为其初相。

    这个包含正弦曲线的图像经过逆Radon变换的表达式为:


    经过如公式(15)所示的推倒过程,能将含有正弦曲线的图像经过逆Radon的表达
    式,得到参数空间图,参数空间图上所包含的特显点分别对应正弦曲线的参数。经过上述变
    换可将正弦曲线转换为参数空间的点然后进行参数估计。即特显点到图像中心的距离为正
    弦曲线的最大值,特显点与图像中心连线与图像竖直线的夹角为正弦曲线的初相。所得特
    显点位置就相当于得到了微动参数中的微多普勒值和初相。因此,经过逆Radon变换可以将
    微动目标参数估计问题转化为逆Radon变换域参数空间特显点的位置提取问题,可以直接
    估计得到微多普勒值和初相。

    本发明采用时频分析和逆Radon变换结合的方法,对解模糊后的信号进行时频分
    析以得到完整的时频曲线,然后根据微动目标运动的简谐性,利用逆Radon变换将正弦时频
    曲线变换成参数空间的点,以实现最后的参数估计。

    本发明提供的基于脉内干涉的太赫兹宽带雷达微多普勒解模糊方法,基于脉内干
    涉的微多普勒解模糊,即将宽带雷达解线频调方式接收后的宽带回波信号视为若干单频的
    回波信号,选择合适的差频进行干涉,并将所有相同差频的干涉结果进行累加,最后进行时
    频分析和参数提取。有别于传统的基于多通道或者多重频的系统级解决方案,本发明方法
    具有实现简单,不增加系统复杂度等优势;有别于传统信号级解决方案,本方法具有计算量
    小、性能稳健等特点。

    优选的,将若干个干涉之后的信号进行累加以消除交叉项的影响,在实际中为了
    达到较好的交叉项抑制目的,选择合适的差频信号进行共轭相乘,以便将载频减小为二者
    载频之差,在实现中一般取差频信号为带宽的一半,因而选择干涉差频为B/2,这样的话,如
    果每个脉冲有N个采样,就会有N/2组干涉信号进行累加,交叉项抑制效果较好。

    以下结合具体实例对本发明提供的方法进行详细说明。

    以一个载频为220GHz的宽带雷达系统为例,以旋转电机带动角反射器作为目标,
    说明了本方法的有效性。实验中系统载频221.59GHz,带宽12.816GHz,脉冲重复周期1ms,即
    PRF为1000,每个脉冲内采样点数N为4096,观测脉冲数目8192个。两个旋转角反射器的旋转
    半径分别为16cm和24cm,两次实验中旋转角速度设置值为40r/min和60r/min。本实例以雷
    达实验系统中常用的旋转目标为例进行说明。所处理旋转目标如图3所示,是一个四面的二
    面角结构,具有类似于理想散射点的散射特性,在处理中可以视为理想点目标。

    在上述条件下,采用本发明提供的方法进行处理。转速为40r/min时两个角反射器
    的微多普勒理论值分别为990.08Hz和1485.1Hz,宽带雷达解线频调方式接收后信号的时频
    分布如图4中的a)所示。转速为60r/min时两个角反射器的微多普勒理论值分别为1485.1Hz
    和2227.7Hz,宽带雷达解线频调方式接收后信号的时频分布如图4中的b)所示。转速为40r/
    min和60r/min时,两个角反射器的微多普勒值均远超过PRF/2。

    从图4的时频分布也可以看出来,由于旋转目标的微多普勒值远远超过PRF/2而造
    成模糊,时频分布图上无法得到完整的正弦曲线。通过本发明提供的方法处理之后得到的
    干涉信号的时频分布图如图5所示??梢钥闯?,经过脉内干涉处理之后,旋转目标的微多普
    勒值被缩小了2fc/B=34.58倍,所得结果中不会发生微多普勒模糊,从其时频分布上也可
    以看出完整的正弦曲线。

    对解模糊之后的时频分布进行逆Radon变换之后的参数空间如图6所示。通过逆
    Radon变换,提取特显点与图像中心的距离,可以估计得到在40r/min下两个目标的微多普
    勒值为28.2Hz和42.5Hz,乘以缩比倍数2fc/B=34.58之后即可得到微多普勒估计值
    975.157Hz和1469.65Hz。

    可以估计得到在60r/min下两个目标的微多普勒值为42.4Hz和64.7Hz,乘以缩比
    倍数2fc/B=34.58之后即可得到微多普勒估计值1466.192Hz和2330.692Hz。

    将估计值与理论值比较可知,本发明的参数估计相对误差小于2%,验证了本发明
    方法的有效性。

    本领域技术人员将清楚本发明的范围不限制于以上讨论的示例,有可能对其进行
    若干改变和修改,而不脱离所附权利要求书限定的本发明的范围。尽管己经在附图和说明
    书中详细图示和描述了本发明,但这样的说明和描述仅是说明或示意性的,而非限制性的。
    本发明并不限于所公开的实施例。

    通过对附图,说明书和权利要求书的研究,在实施本发明时本领域技术人员可以
    理解和实现所公开的实施例的变形。在权利要求书中,术语“包括”不排除其他步骤或元素,
    而不定冠词“一个”或“一种”不排除多个。在彼此不同的从属权利要求中引用的某些措施的
    事实不意味着这些措施的组合不能被有利地使用。权利要求书中的任何参考标记不构成对
    本发明的范围的限制。

    关于本文
    本文标题:基于脉内干涉的太赫兹宽带雷达微多普勒解模糊方法.pdf
    链接地址://www.4mum.com.cn/p-6001249.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
  • 北京快乐8 pc蛋蛋预测尽享官网 赌场游戏 广西快3下载开奖结果 云南11选5选号技巧 扑克牌游戏 单机版成都麻将app 5元以下股票推荐 腾讯分分彩开奖号查询 单机捕鱼达人手机版 分分彩通用稳赚方案 福彩组六奖金多少 极速时时彩计划1期中 浙江快乐彩图表手机板 上海时时彩开奖记录 蓝洞棋牌app棋牌中心