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

    重庆时时彩群发器: 一种针对多径信号的波束形成器设计方法.pdf

    关 键 词:
    一种 针对 信号 波束 形成 设计 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201410168758.6

    申请日:

    2014.04.24

    公开号:

    CN103941234A

    公开日:

    2014.07.23

    当前法律状态:

    授权

    有效性:

    有权

    法律详情: 授权|||实质审查的生效IPC(主分类):G01S 7/02申请日:20140424|||公开
    IPC分类号: G01S7/02; G01S7/36 主分类号: G01S7/02
    申请人: 清华大学
    发明人: 汤俊; 王成; 刘雅琪
    地址: 100084 北京市海淀区北京市100084-82信箱
    优先权:
    专利代理机构: 北京众合诚成知识产权代理有限公司 11246 代理人: 朱琨
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201410168758.6

    授权公告号:

    ||||||

    法律状态公告日:

    2016.10.05|||2014.08.20|||2014.07.23

    法律状态类型:

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

    摘要

    本发明属于雷达和通信领域,尤其是涉及一种针对多径信号的波束形成器设计方法。首先采用空间平滑技术,降低多径信号间的相关性,从而得到满秩的信号空域协方差矩阵;接着配合已知的多径信号部分来向信息,构建最小方差无畸变响应波束形成器;然后依据不同情况选择合适方案,能够获得信号与各个传感器接收数据的互相关矢量,在此基础上设计出接收与多径环境相匹配的最小均方误差波束形成器。相比于传统的波束形成器,本发明解决了多径信号下,当前波束形成器对有用信号利用不充分,阵列孔径损失的问题,提高了阵列有效孔径,且对多径信号利用更充分,为移动通信,相控阵雷达等场合中的多径信号波束形成器设计提供了可靠的方法。

    权利要求书

    权利要求书
    1.  一种针对多径信号的波束形成器设计方法,其特征在于,包括以下步骤:
    步骤1、信号处理系统依据不同来向入射信号个数将传感器阵列划分为多个子阵列,其中信号处理系统为相控阵雷达和/或通信终端;
    步骤2、通过空间平滑处理,由各子阵列协方差矩阵得到满秩的信号空域矩阵;
    步骤3、依据已知多径信号的部分来向信息,构建基于子阵列l的最小方差无畸变响应波束形成器,其中子阵列l是指某传感器阵列;
    步骤4、存储该波束形成器的输出数据和对应的各传感器接收数据;
    步骤5、判断系统噪声功率是否已知且在采样时间内恒定,如果是则计算最小方差无畸变响应波束形成器输出数据与各个传感器接收数据的互相关矢量,该矢量再减去子阵列l的最小方差无畸变响应波束形成器权值与噪声功率的乘积,结果用于构建基于全阵列的最小均方误差波束形成器;如果否,则执行步骤6;
    步骤6、判断系统传感器数量是否超过信号个数4倍,如果是则计算最小方差无畸变响应波束形成器输出数据与子阵列l之外的传感器接收数据互相关矢量,结果用于构建基于子阵列l之外传感器阵列的最小均方误差波束形成器;如果否,则执行步骤7;
    步骤7、将各个传感器数据经过一定延时后与最小方差无畸变响应波束形成器输出数据进行互相关计算,结果用于构建基于全阵列的最小均方误差波束形成器;而如果多径信号本身具有周期性,那么选择延时量接近多径信号周期。

    2.  根据权利要求1所述方法,其特征在于,所述步骤1中,第l个子阵列的k时刻传感器接收数据矢量表示为
    Xl(k)=[xl(k),…,xl+N-1(k)]T,l=1,2,…M-N+1
    其中,M为传感器总数,N为子阵列含传感器个数,xl(k)表示为传感器l在k时刻接收的数据,xl+N-1(k)表示为第l+N-1个传感器在k时刻接收的数据,(·)T表示转置操作,k=1,2,…K,K为处理器一次缓存的数据量。

    3.  根据权利要求1所述方法,其特征在于,所述步骤2中,子阵列协方差矩阵通过下式计算
    Rl=1KΣk=1KXl(k)XlH(k)]]>
    则前向空间平滑后恢复的协方差矩阵表示为
    Rsmooth=1LΣl=1LRl]]>
    其中,K为处理器一次缓存的数据量,Xl(k)表示第l个子阵列的k时刻传感器接收数据矢量,(·)H表示共轭转置操作,L为子阵列个数。

    4.  根据权利要求1所述方法,其特征在于,所述步骤3中构建的基于子阵列l的波束形成器权向量为
    wMVDR=Rsmooth-1al(θ1)alH(θ1)Rsmooth-1al(θ1)]]>
    其中,θ1为已知多径信号的部分来向信息,al(θ1)为已知多径信号的部分来向信息对应的子阵列l的阵列流形,Rsmooth为前向空间平滑后恢复的协方差矩阵,得到权向量后,MVDR波束形成器输出为
    s^1=wMVDRHXl(k)]]>
    最后信号处理系统缓存接收数据向量X(k)=[x1(k),…,xM(k)]T和MVDR-SS输出数据其中Xl(k)表示第l个子阵列的k时刻传感器接收数据矢量,k=1,2,…K。

    5.  根据权利要求1所述方法,其特征在于,所述步骤5采用噪声消去方案设计MMSE波束形成器,首先按如下公式计算接收数据和MVDR-SS输出数据的互相关矢量

    噪声消去方案的MMSE波束形成器的权向量由如下公式计算
    wMMSE=(1KΣk=1KX(k)XH(k))-1rxd]]>
    计算得权向量后,即由如下公式计算多径信号的波束形成器输出
    s^2(k)=wMMSEHX(k)]]>
    其中,X(k)为最后信号处理系统缓存接收数据向量,为MVDR波束形成器输出,(·)*表示取共轭操作,wMVDR为基于子阵列l的波束形成器权向量,表示已知的噪声能量。

    6.  根据权利要求1所述方法,其特征在于,所述步骤6采用空间分集方案设计MMSE波束形成器,划分子阵使得各个子阵的传感器个数N大于信号个数,小于子阵列个数L的三分之一,按如下公式计算接收数据和MVDR-SS输出数据的互相关矢量
    rxd=1KΣk=1KXM-N(k)s^1*(k)]]>
    之后,空间分集方案的MMSE波束形成器权向量按如下公式计算
    wMMSE=(1KΣk=1KXM-N(k)XM-NH(k))-1rxd]]>
    计算得权向量后,即由如下公式计算多径信号的波束形成器输出
    s^2(k)=wMMSEHXM-N(k)]]>
    其中,为MVDR波束形成器输出,(·)*表示取共轭操作,XM-N(k)=[xN+1(k),…,xM(k)]T为除子阵列l外的传感器接收数据向量,xN+1(k) 表示为第N+1个传感器在k时刻接收的数据,xM(k)表示为第M个传感器在k时刻接收的数据。

    7.  根据权利要求1所述方法,其特征在于,所述步骤7发明采用时间分集方案设计MMSE波束形成器,按如下公式计算接收数据和MVDR-SS输出数据的互相关矢量
    rxd=1KΣk=1K-τX(k+τ)s^1*(k)]]>
    之后,时间分集方案的MMSE波束形成器权向量按如下公式计算
    wMMSE=(1KΣk=1K-τX(k+τ)XH(k+τ))-1rxd]]>
    计算得权向量后,即由如下公式计算多径信号的波束形成器输出
    s^2(k)=wMMSEHX(k)]]>
    其中,为MVDR波束形成器输出,(·)*表示取共轭操作,X(k+τ)表示经过延时τ后的全阵列接收向量,X(k)为最后信号处理系统缓存接收数据向量。

    说明书

    说明书一种针对多径信号的波束形成器设计方法
    技术领域
    本发明属于雷达和通信领域,尤其是涉及一种针对多径信号的波束形成器设计方法。
    背景技术
    阵列信号处理是信号处理的一大分支,广泛应用在雷达、声呐、射电天文探测、医学检测和移动通信等领域,它是指将多个传感器按一定方式布置在空间不同位置上,形成传感器阵列,再对传感器组进行数据处理的技术。波束形成技术是阵列信号处理的重要方面。其实质是通过对各个传感器的数据加权,来尽可能地实现增强期望信号,抑制干扰的目的。对接收的空间信号进行处理时,与传统的单个定向传感器相比,使用波束形成技术的传感器阵列具有波束控制灵活,干扰抑制能力强,和空间分辨能力高的优点。而自适应波束形成要求根据信号环境变化自适应地改变加权因子。
    自适应波束形成技术具有非常重要的理论意义和实用价值。传统的自适应波束形成器设计方法都假设干扰与接收信号之间相互独立,然而在各类信号传输过程中,多径效应往往使得传感器阵列接收到的有用信号呈现相关性,这使得常规的自适应波束形成方法性能严重下降。以空间平滑为代表的各类去相关技术,能够有效降低信号间的相关性,然而空间平滑会降低传感器阵列的有效孔径,且在有用信号来向信息不全的情况下,采用该类技术的波束形成器会将有用信号当作干扰进行抑制。
    针对多径信号,如何在不损失阵列孔径的前提下,设计出高效、稳健的波束形成器,是亟待解决的现实问题。
    发明内容
    本发明的目的在于,针对多径信号,如何在不损失阵列孔径的前提下设计出高效、稳健的波束形成器,提出了一种针对多径信号的波束形成器设计方法,其特征在于,包括以下步骤:
    步骤1、信号处理系统依据不同来向入射信号个数将传感器阵列划分为多个子阵列,其中信号处理系统为相控阵雷达和/或通信终端;
    步骤2、通过空间平滑处理,由各子阵列协方差矩阵得到满秩的信号空域矩阵;
    步骤3、依据已知多径信号的部分来向信息,构建基于子阵列l的最小方差无畸变响应波束形成器,其中子阵列l是指某传感器阵列;
    步骤4、存储该波束形成器的输出数据和对应的各传感器接收数据;
    步骤5、判断系统噪声功率是否已知且在采样时间内恒定,如果是则计算最小方差无畸变响应波束形成器输出数据与各个传感器接收数据的互相关矢量,该矢量再减去子阵列l的最小方差无畸变响应波束形成器权值与噪声功率的乘积,结果用于构建基于全阵列的最小均方误差波束形成器;如果否,则执行步骤6;
    步骤6、判断系统传感器数量是否超过信号个数4倍,如果是则计算最小方差无畸变响应波束形成器输出数据与子阵列l之外的传感器接收数据互相关矢量,结果用于构建基于子阵列l之外传感器阵列的最小均方误差波束形成器;如果否,则执行步骤7;
    步骤7、将各个传感器数据经过一定延时后与最小方差无畸变响应波束形成器输出数据进行互相关计算,结果用于构建基于全阵列的最小均方误差波束形成器;而如果多径信号本身具有周期性,那么选择延时量接近多径信号周期。
    所述步骤1中,第l个子阵列的k时刻传感器接收数据矢量表示为
    Xl(k)=[xl(k),…,xl+N-1(k)]T,l=1,2,…M-N+1
    其中,M为传感器总数,N为子阵列含传感器个数,xl(k)表示为传感器l在k时刻接收的数据,xl+N-1(k)表示为第l+N-1个传感器在k时刻接收的数据,(·)T表示转置操作,k=1,2,…K,K为处理器一次缓存的数据量。
    所述步骤2中,子阵列协方差矩阵通过下式计算
    Rl=1KΣk=1KXl(k)XlH(k)]]>
    则前向空间平滑后恢复的协方差矩阵表示为
    Rsmooth=1LΣl=1LRl]]>
    其中,K为处理器一次缓存的数据量,Xl(k)表示第l个子阵列的k时刻传感器接收数据矢量,(·)H表示共轭转置操作,L为子阵列个数。
    所述步骤3中构建的基于子阵列l的波束形成器权向量为
    wMVDR=Rsmooth-1al(θ1)alH(θ1)Rsmooth-1al(θ1)]]>
    其中,θ1为已知多径信号的部分来向信息,al(θ1)为已知多径信号的部分来向信息对应的子阵列l的阵列流形,Rsmooth为前向空间平滑后恢复的协方差矩阵,得到权向量后,MVDR波束形成器输出为
    s^1=wMVDRHXl(k)]]>
    最后信号处理系统缓存接收数据向量X(k)=[x1(k),…,xM(k)]T和MVDR-SS输出数据其中Xl(k)表示第l个子阵列的k时刻传感器接收数据矢量,k=1,2,…K。
    所述步骤5采用噪声消去方案设计MMSE波束形成器,首先按如下公式计算接收数据和MVDR-SS输出数据的互相关矢量

    噪声消去方案的MMSE波束形成器的权向量由如下公式计算
    wMMSE=(1KΣk=1KX(k)XH(k))-1rxd]]>
    计算得权向量后,即由如下公式计算多径信号的波束形成器输出
    s^2(k)=wMMSEHX(k)]]>
    其中,X(k)为最后信号处理系统缓存接收数据向量,为MVDR波束形成器输出,(·)*表示取共轭操作,wMVDR为基于子阵列l的波束形成器权向量,表示已知的噪声能量。
    所述步骤6采用空间分集方案设计MMSE波束形成器,划分子阵使得各个子阵的传感器个数N大于信号个数,小于子阵列个数L的三分之一,按如下公式计算接收数据和MVDR-SS输出数据的互相关矢量
    rxd=1KΣk=1KXM-N(k)s^1*(k)]]>
    之后,空间分集方案的MMSE波束形成器权向量按如下公式计算
    wMMSE=(1KΣk=1KXM-N(k)XM-NH(k))-1rxd]]>
    计算得权向量后,即由如下公式计算多径信号的波束形成器输出
    s^2(k)=wMMSEHXM-N(k)]]>
    其中,为MVDR波束形成器输出,(·)*表示取共轭操作,XM-N(k)=[xN+1(k),…,xM(k)]T为除子阵列l外的传感器接收数据向量,xN+1(k)表示为第N+1个传感器在k时刻接收的数据,xM(k)表示为第M个传感器在k时刻接收的数据。
    所述步骤7发明采用时间分集方案设计MMSE波束形成器,按如下公式计算接收数据和MVDR-SS输出数据的互相关矢量
    rxd=1KΣk=1K-τX(k+τ)s^1*(k)]]>
    之后,时间分集方案的MMSE波束形成器权向量按如下公式计算
    wMMSE=(1KΣk=1K-τX(k+τ)XH(k+τ))-1rxd]]>
    计算得权向量后,即由如下公式计算多径信号的波束形成器输出
    s^2(k)=wMMSEHX(k)]]>
    其中,为MVDR波束形成器输出,(·)*表示取共轭操作,X(k+τ)表示经过延时τ后的全阵列接收向量,X(k)为最后信号处理系统缓存接收数据向量。
    本发明的有益效果在于,相比于传统的单纯利用空间平滑技术的波束形成器,本发明解决了多径信号下,当前波束形成器对有用信号利用不充分,阵列孔径损失的问题,提高了阵列有效孔径,且对多径信号利用更充分。它适用于无相关干扰的场合,为移动通信,相控阵雷达等场合中的多径信号波束形成器设计提供了可靠的方法。
    附图说明
    图1是本设计方法的总体步骤图;
    图2是本设计方法的空间平滑子阵列划分示意图;
    图3是本设计方法的噪声消去方案示意图;
    图4是本设计方法的空间分集方案示意图;
    图5是本设法方法的时间分集方案示意图。
    具体实施方式
    下面结合附图,对优选实施例作详细说明。
    本发明提出的一种针对多径信号的波束形成器设计方法,如图1所示,包括以下步骤:
    步骤1、信号处理系统(如相控阵雷达、通信终端等)依据不同来向入射信号个数,划分传感器阵列为多个子阵列;
    步骤2、通过空间平滑处理,由各子阵列协方差矩阵得到满秩的信号空域矩阵;
    步骤3、依据已知多径信号的部分来向信息,构建基于子阵列l的最小方差无畸变响应(MVDR)波束形成器,其中子阵列l是指某传感器阵列;
    步骤4、存储该波束形成器的输出数据和对应的各传感器接收数据;
    步骤5、判断系统噪声功率是否已知且在采样时间内恒定,如果是则计算最小方差无畸变响应波束形成器输出数据与各个传感器接收数据的互相关矢量, 该矢量再减去子阵列l的最小方差无畸变响应波束形成器权值与噪声功率的乘积,结果用于构建基于全阵列的最小均方误差(MMSE)波束形成器;如果否,则执行步骤6;
    步骤6、判断系统传感器数量是否超过信号个数4倍,如果是则计算最小方差无畸变响应波束形成器输出数据与子阵列l之外的传感器接收数据互相关矢量,结果用于构建基于子阵列l之外传感器阵列的最小均方误差波束形成器;如果否,则执行步骤7;
    步骤7、将各个传感器数据经过一定延时后与最小方差无畸变响应波束形成器输出数据进行互相关计算,结果用于构建基于全阵列的最小均方误差波束形成器;而如果多径信号本身具有周期性,那么选择延时量接近多径信号周期。
    下面以等距线阵和前向平滑为例说明具体操作,空间平滑子阵列划分如图2所示。传感器总数为M,子阵列含传感器个数为N。用xl(t)表示传感器l的接收数据,则第l个子阵列的k时刻传感器接收数据矢量表示为
    Xl(k)=[xl(k),…,xl+N-1(k)]T,l=1,2,…M-N+1
    其中(·)T表示转置操作??占淦交蠡值淖诱罅懈鍪齃和子阵传感器个数N都不小于相关的信号个数。所以要依据系统接收信号数量选择划分的子阵列个数,处理器一次缓存的数据量为K,则子阵列协方差矩阵通过下式计算
    Rl=1KΣk=1KXl(k)XlH(k)]]>
    其中(·)H表示共轭转置操作,前向空间平滑后恢复的协方差矩阵表示为
    Rsmooth=1LΣl=1LRl]]>
    其中,步骤3中将已知多径信号的部分来向信息记为θ1,其对应的子阵列l的阵列流形为al(θ1),则构建的基于子阵列l的波束形成器权向量为
    wMVDR=Rsmooth-1al(θ1)alH(θ1)Rsmooth-1al(θ1)]]>
    为方便表达,将上述先空间平滑再构建出的最小方差无畸变波束形成器记作MVDR-SS。得到权向量后,波束形成器输出由下式计算
    s^1=wMVDRHXl(k)]]>
    最后信号处理系统缓存接收数据向量X(k)=[x1(k),…,xM(k)]T和MVDR-SS输出数据其中k=1,2,…K。
    需要说明的是,上述步骤1以等距线阵和前向空间平滑为例介绍了具体步骤。而本发明并不局限于等距线阵,任何阵列,只要可以划分为多个具有相同阵列流行子阵列即可采用空间平滑处理,且发明对具体使用的平滑算法也无限制。
    其中步骤5中如果系统噪声功率已知且在采样时间内恒定,即噪声功率在一段时间内浮动不超过总功率的10%即为恒定,本发明采用如图3所示的噪声消去方案设计MMSE波束形成器。首先,按如下公式计算接收数据和MVDR-SS输出数据的互相关矢量

    其中(·)*表示取共轭操作,表示已知的噪声能量。噪声消去方案的MMSE波束形成器的权向量由如下公式计算
    wMMSE=(1KΣk=1KX(k)XH(k))-1rxd]]>
    计算得权向量后,即由如下公式计算多径信号的波束形成器输出
    s^2(k)=wMMSEHX(k)]]>
    其中步骤6中如果系统传感器数量超过信号个数4倍时,本发明采用如图4所示的空间分集方案设计MMSE波束形成器?;肿诱笫沟酶鞲鲎诱蟮拇衅鞲鍪齆大于信号个数,小于子阵列个数L的三分之一。
    定义除子阵列l外的传感器接收数据向量为XM-N(k)=[xN+1(k),…,xM(k)]T,首先,按如下公式计算接收数据和MVDR-SS输出数据的互相关矢量
    rxd=1KΣk=1KXM-N(k)s^1*(k)]]>
    之后,空间分集方案的MMSE波束形成器权向量按如下公式计算
    wMMSE=(1KΣk=1KXM-N(k)XM-NH(k))-1rxd]]>
    计算得权向量后,即由如下公式计算多径信号的波束形成器输出
    s^2(k)=wMMSEHXM-N(k)]]>
    步骤7中若系统不满足上述要求时,本发明采用如图5所示的时间分集方案设计MMSE波束形成器。将存储的接收数据向量延时后与MVDR-SS输出数据相关,如果已知多径信号具有周期性,那么选择周期作为延时量。
    按如下公式计算接收数据和MVDR-SS输出数据的互相关矢量
    rxd=1KΣk=1K-τX(k+τ)s^1*(k)]]>
    其中,X(k+τ)表示经过延时τ后的全阵列接收向量;
    之后,时间分集方案的MMSE波束形成器权向量按如下公式计算
    wMMSE=(1KΣk=1K-τX(k+τ)XH(k+τ))-1rxd]]>
    计算得权向量后,即由如下公式计算多径信号的波束形成器输出
    s^2(k)=wMMSEHX(k)]]>
    以上所述,仅为本发明较佳的具体实施方式,但本发明的?;し段Р⒉痪窒抻诖?,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的?;し段е?。因此,本发明的?;し段вΩ靡匀ɡ蟮谋;し段?。

    关于本文
    本文标题:一种针对多径信号的波束形成器设计方法.pdf
    链接地址://www.4mum.com.cn/p-6143304.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
  • 后三组选包胆技巧时时 彩 卖什么稳赚不赔 时时彩二星直选 无错36码特2019 好运来飞艇计划软件推荐软件 买彩票如何能稳赚不赔 网上赌龙虎的技巧 百盈彩票是黑平台吗 默认版块三肖六码公开网址 八码滚雪球倍投公式 定位胆大小单双诀窍 老虎机的规律 幸运飞艇人工计划聊天室 功夫时时彩计划手机客户端 分分pk拾 黄金计划软件官网