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

    重庆时时彩大小公式: 一种高精度的传递函数采样方法.pdf

    关 键 词:
    一种 高精度 传递函数 采样 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201710081416.4

    申请日:

    2017.02.15

    公开号:

    CN106777834A

    公开日:

    2017.05.31

    当前法律状态:

    实审

    有效性:

    审中

    法律详情: 实质审查的生效IPC(主分类):G06F 17/50申请日:20170215|||公开
    IPC分类号: G06F17/50 主分类号: G06F17/50
    申请人: 华中科技大学
    发明人: 陈学东; 姜伟; 吴九林
    地址: 430074 湖北省武汉市洪山区珞喻路1037号
    优先权:
    专利代理机构: 华中科技大学专利中心 42201 代理人: 梁鹏
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201710081416.4

    授权公告号:

    |||

    法律状态公告日:

    2017.06.23|||2017.05.31

    法律状态类型:

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

    摘要

    本发明属于机械系统结构动力学振动计算与分析相关技术领域,其公开一种高精度的传递函数采样方法,所述传递函数采样方法包括以下步骤:(1)构建机械系统的动力学方程;(2)通过系统振动模态计算方法得到系统振动频率;(3)设置采样方式、采样带宽及采样点数以生成基础传递函数采样频率点序列;(4)提取处于传递函数分析带宽内的模态频率及设置每阶模态频率对应细化采样点数的总数目,以形成传递函数分析带宽内的模态频率对应的细化采样点序列;(5)将基础传递函数采样频率点序列与细化采样点序列合并以生成最终的传递函数采样频率点序列;(6)对最终的传递函数采样频率点序列进行传递函数数值计算以得到最终离散数据的传递函数图。

    权利要求书

    1.一种高精度的传递函数采样方法,其包括以下步骤:
    (1)构建机械系统的动力学方程;
    (2)通过系统振动模态计算方法得到系统振动频率;
    (3)设置采样方式、采样带宽及采样点数以生成基础传递函数采样频率点序列;
    (4)提取处于传递函数分析带宽内的模态频率及设置每阶模态频率对应细化采样点数
    的总数目,以形成传递函数分析带宽内的模态频率对应的细化采样点序列;
    (5)将基础传递函数采样频率点序列与细化采样点序列合并以生成最终的传递函数采
    样频率点序列;
    (6)对最终的传递函数采样频率点序列进行传递函数数值计算以得到最终离散数据的
    传递函数图。
    2.如权利要求1所述的高精度的传递函数采样方法,其特征在于:所述基础传递函数采
    样频率点序列与所述细化采样点序列合并后,通过按从小到大的顺序排列来生成所述的最
    终的传递函数采样频率点序列。
    3.如权利要求1所述的高精度的传递函数采样方法,其特征在于:所述基础传递函数采
    样频率点序列为ω1、ω2、ω3、...、ωk、...、ωM,
    其中
    ωstart为带宽起始频率点,ωend为带宽截止频率点,M为总采样点数。
    4.如权利要求1所述的高精度的传递函数采样方法,其特征在于:所述细化采样点序列
    为:
    ωd11、ωd12、...、ωd1j、...、ωd1p、...、ωdi1、ωdi2、...、ωdij、...、ωdip、...(i=1,
    2,...,m)
    其中,p为对应每阶模态频率点前后细化采样点的总数目,其作用是在每阶模态频率点
    前后按指定规则生成p个细化采样点,以避免因采样频率过低而产生较大的传递函数误差;
    系统N阶模态频率中的m阶模态频率处于传递函数分析带宽内。

    说明书

    一种高精度的传递函数采样方法

    技术领域

    本发明属于机械系统结构动力学振动计算与分析相关技术领域,更具体地,涉及
    一种高精度的传递函数采样方法。

    背景技术

    传递函数分析是机械系统动力学分析中的重要环节,是机械系统固有特性的表
    征,机械系统受激振动的响应特性只与激励的性质及系统激励-响应点间的传递特性(传递
    函数)有关。通过获取不同传递途径下系统传递函数的幅频特性图和相频特性图,可以揭示
    系统在不同频域带宽内的响应特性,以对系统进行实时评估,因此,传递函数分析在机械系
    统的设计、评估、结构改进以及模型修正、修改等过程中有重要的作用。

    传递函数的获取方法有两种:第一种是通过实验进行传递函数的测量获??;第二
    种是先建立系统的等效动力学模型,再通过传递函数的理论计算与分析方法获得。

    在传递函数特性曲线中(以幅频图为例),其核心的内容为传递函数中共振峰的频
    率以及共振峰对应的幅值,传递函数共振峰频率与系统模态频率间有着必然联系。目前,国
    内外通用的各种机械系统多体动力学软件(以ADAMS为例)在传递函数计算与分析过程中,
    所采取的采样方法主要有简单的线性频率分布及对数频率分布两种。用户设置传递函数显
    示带宽的开始频率、截止频率以及频率点数之后,通过选取线性等分(linear spacing of
    steps)或者对数等分(logarithmic spacing of steps)的频率采样方式以得到固定的一
    组频率序列,软件经计算及采样过程后,能够得到用户指定输入通道-输出通道间离散数据
    的传递函数曲线。但是,由于并未考虑到传递函数共振峰频率与系统模态频率的特殊关系,
    在相同采样方式下,不同的采样频率(采样点数)会对最终传递函数计算结果曲线有很大影
    响,当采样频率过低即采样点数较少时,无论是线性采样方式还是对数采样方式,均会产生
    共振频率移位、共振幅值过低的现象,甚至在某些极端情况下,部分传递函数共振峰会因为
    采样频率过低而直接被遗漏:如重叠峰值(或相近峰)出现漏采的情况,精度较低,不利于对
    系统的计算及评估。相应地,本领域存在着发展一种精度较高的传递函数采样方法的技术
    需求。

    发明内容

    针对现有技术的以上缺陷或改进需求,本发明提供了一种高精度的传递函数采样
    方法,其结合系统模态分析的数据(模态频率),在低频率条件下,通过细化采样以增加可能
    共振点(即系统模态频率点)附近采样数据点的密集程度,在总采样点数少量增加的情况下
    高精度的获取了传递函数最重要的特征、共振频率点以及共振峰值,有利于机械系统结构
    动力学振动计算与分析。

    为实现上述目的,本发明提供了一种高精度的传递函数采样方法,其包括以下步
    骤:

    (1)构建机械系统的动力学方程;

    (2)通过系统振动模态计算方法得到系统振动频率;

    (3)设置采样方式、采样带宽及采样点数以生成基础传递函数采样频率点序列;

    (4)提取处于传递函数分析带宽内的模态频率及设置每阶模态频率对应细化采样
    点数的总数目,以形成传递函数分析带宽内的模态频率对应的细化采样点序列;

    (5)将基础传递函数采样频率点序列与细化采样点序列合并以生成最终的传递函
    数采样频率点序列;

    (6)对最终的传递函数采样频率点序列进行传递函数数值计算以得到最终离散数
    据的传递函数图。

    进一步的,所述基础传递函数采样频率点序列与所述细化采样点序列合并后,通
    过按从小到大的顺序排列来生成所述的最终的传递函数采样频率点序列。

    进一步的,所述基础传递函数采样频率点序列为ω1、ω2、ω3、...、ωk、...、ωM。

    其中,

    ωstart为带宽起始频率点,ωend为带宽截止频率点,M为总采样点数。

    进一步的,所述细化采样点序列为:

    ωd11、ωd12、...、ωd1j、...、ωd1p、...、ωdi1、ωdi2、...、ωdij、...、ωdip、...(i=1,
    2,...,m)

    其中,p为对应每阶模态频率点前后细化采样点的总数目,其作用是在每阶模态频
    率点前后按指定规则生成p个细化采样点,以避免因采样频率过低而产生较大的传递函数
    误差;系统N阶模态频率中的m阶模态频率处于传递函数分析带宽内。

    总体而言,通过本发明所构思的以上技术方案与现有技术相比,本发明较佳实施
    方式提供的高精度的传递函数采样方法,其结合系统模态分析的数据(模态频率),在低频
    率条件下,通过细化采样以增加可能共振点(即系统模态频率点)附近采样数据点的密集程
    度,在总采样点数少量增加的情况下高精度的获取了传递函数最重要的特征、共振频率点
    以及共振峰值,有利于机械系统结构动力学振动计算与分析。

    附图说明

    图1是本发明较佳实施方式提供的高精度的传递函数采样方法的流程图。

    图2是图1中的高精度的传递函数采样方法涉及的三自由度弹簧阻尼-质量系统的
    示意图。

    图3(a)为对数采样方式、采样点数为101时,分别采用本发明的传递函数采样方法
    与现有方法的获得的输入Fi-输出ao传递函数幅频曲线的对比图。

    图3(b)为对数采样方式、采样点数为101时,分别采用本发明的高精度的传递函数
    采样方法与现有方法获得的输入Fi-输出ao传递函数幅频曲线的局部放大对比图(第三峰值
    区间)。

    图3(c)为线性采样方式、采样点数为101时,分别采用本发明的高精度的传递函数
    采样方法与现有方法获得的输入Fi-输出ao传递函数幅频曲线的对比图。

    图3(d)为线性采样方式、采样点数为101时,分别采用本发明的高精度的传递函数
    采样方法与现有方法获得的输入Fi-输出ao传递函数幅频曲线的局部放大对比图(第三峰值
    区间)。

    图4(a)为对数采样方式、采样点数为1001时,分别采用本发明的高精度的传递函
    数采样方法与现有方法获得的输入Fi-输出ao传递函数幅频曲线的对比图。

    图4(b)为对数采样方式、采样点数为1001时,分别采用本发明的高精度的传递函
    数采样方法与现有方法获得的输入Fi-输出ao传递函数幅频曲线的局部放大对比图(第二、
    三峰值区间)。

    图4(c)为线性采样方式、采样点数为1001时,分别采用本发明的高精度的传递函
    数采样方法与现有方法获得的输入Fi-输出ao传递函数幅频曲线的对比图。

    图4(d)为线性采样方式、采样点数为1001时,分别采用本发明的高精度的传递函
    数采样方法与现有方法获得的输入Fi-输出ao传递函数幅频曲线的局部放大对比图(第二、
    三峰值区间)。

    图5(a)为对数采样方式、采样点数分别为101、1001时,采用现有方法获得的输入
    Fi-输出ao传递函数幅频曲线的对比图。

    图5(b)为对数采样方式、采样点数分别为101、1001时,采用现有方法获得的输入
    Fi-输出ao传递函数幅频曲线的局部放大对比图(第二、三峰值区间)。

    图5(c)为线性采样方式、采样点数分别为101、1001时,采用现有方法获得的输入
    Fi-输出ao传递函数幅频曲线的对比图。

    图5(d)为线性采样方式、采样点数分别为101、1001时,采用现有方法获得的输入
    Fi-输出ao传递函数幅频曲线的局部放大对比图(第二、三峰值区间)。

    图6(a)为对数采样方式、采样点数分别为101、1001时,采用发明的高精度的传递
    函数采样方法获得的输入Fi-输出ao传递函数幅频曲线的对比图。

    图6(b)为对数采样方式、采样点数分别为101、1001时,采用本发明的高精度的传
    递函数采样方法获得的输入Fi-输出ao传递函数幅频曲线的局部放大对比图(第二、三峰值
    区间)。

    图6(c)为线性采样方式、采样点数分别为101、1001时,采用本发明的高精度的传
    递函数采样方法获得的输入Fi-输出ao传递函数幅频曲线的对比图。

    图6(d)为线性采样方式、样点数分别为101、1001时,采用本发明的高精度的传递
    函数采样方法获得的输入Fi-输出ao传递函数幅频曲线的局部放大对比图(第二、三峰值区
    间)。

    具体实施方式

    为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对
    本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并
    不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要
    彼此之间未构成冲突就可以相互组合。

    请参阅图1及图2,本发明较佳实施方式提供的高精度的传递函数采样方法,具体
    为一种可自动捕捉谐振峰值的非规则采样方法,所述传递函数采样方法可以在相对较低的
    采样频率设置下显著提高传递函数计算的输出精度,且可广泛应用于各种结构振动特性仿
    真。

    具体地,所述的高精度的传递函数采样方法主要包括以下步骤:

    步骤一,构建机械系统的动力学方程。具体地,本实施方式以三自由度弹簧阻尼-
    质量系统为例对本发明提供的高精度的传递函数采样方法进行说明,三自由度弹簧阻尼-
    质量系统的参数为:质量块质量分别为m1=1kg,m2=2kg,m3=1kg,弹簧阻尼连接参数分别
    为k1=1×106N/m,c1=1N·s/m,k2=2×104N/m,c2=1N·s/m,k3=1×104N/m,c3=1N·s/m,
    k4=0.9×106N/m,c4=1N·s/m,传递函数计算时输入通道为作用在质量块1上的扫频力信
    号Fi,输出通道为测量质量块3的振动加速度信号ao。以质量块的平动位移为广义坐标X=
    [x1 x2 x3]T,则得机械系统的振动方程如下:


    其中:


    步骤二,通过系统振动模态计算方法得到系统振动频率。具体地,系统振动频率分
    别为ωd1、ωd2、ωd3、...、ωdi、...、ωdN。

    其中,ωdi为动力学系统第i阶有阻尼振动模态频率,N为系统模态频率总阶数。本
    实施方式中,对于三自由度弹簧阻尼-质量系统,ωd1=19.33Hz、ωd2=151.83Hz、ωd3=
    160.75Hz。

    步骤三,设置常规传递函数采样方式、采样带宽及采样点数,以生成基础传递函数
    采样频率点序列。具体地,所述基础传递函数采样频率点序列(以对数等分采样方式为例)
    为ω1、ω2、ω3、...、ωk、...、ωM。

    其中:


    ωstart为带宽起始频率点,ωend为带宽截止频率点,M为总采样点数。

    本实施方式中,设置采样带宽为1Hz~1000Hz,采样点数为101,常规传递函数采样
    方式分别为对数频率采样和线性频率采样,则有:

    A:对数频率采样

    beta=(lg1000-lg1)/(101-1)=0.03,ω1=1,ωM=1000,ωk=1×10(k-1)×0.03Hz

    生成基础传递函数采样频率点序列:

    ω1=1、ω2=1.07、...、ωk=1×10(k-1)×0.03、...、ω100=933.25、ω101=1000Hz

    B:线性频率采样

    beta=(1000-1)/(101-1)=9.99,ω1=1,ωM=1000,ωk=(9.99k-8.99)Hz

    生成基础传递函数采样频率点序列:

    ω1=1、ω2=10.99、...、ωk=(9.99k-8.99)、...、ω100=990.01、ω101=1000Hz

    步骤四,设置每阶模态频率点前后细化采样点的总数目,以及提取系统N阶模态频
    率中处于传递函数分析带宽内的m阶模态频率,得到的细化采样点形成该前m阶模态频率点
    附近的细化采样点序列。本实施方式中,传递函数分析带宽为待分析的带宽。具体地,所述
    细化采样点序列为:

    ωd11、ωd12、...、ωd1j、...、ωd1p、...、ωdi1、ωdi2、...、ωdij、...、ωdip、...(i=1,
    2,...,m)。

    其中,p为对应每阶模态频率点前后细化采样点的总数目,其作用是在每阶模态频
    率点前后按指定规则生成p个细化采样点,以避免因采样频率过低而产生较大的传递函数
    误差。

    本实施方式中采样带宽1~1000Hz内包含系统前三阶固有振动频率(即m=3),分
    别为ωd1=19.33Hz、ωd2=151.83Hz、ωd3=160.75Hz。由于整体采样方式为对数采样,对于
    每个共振点,在共振频率以及其前后按规律共增加3个细化频率点(即P=3),总共生成以下
    所示的9个点数的细化采样点序列:

    ωd11=19.05、ωd12=19.33、ωd13=19.62、ωd21=148.06、ωd22=151.83、

    ωd23=155.69、ωd31=156.72、ωd32=160.75、ωd33=164.89

    即:[19.05,19.33,19.62,148.06,151.83,155.69,156.72,160.75,164.89]

    本实施方式的具体细化采样点关系为ωdi2=ωdi,

    步骤五,将基础传递函数采样频率点序列与细化采样点序列合并,并按从小到大
    的顺序排列以作为最终的传递函数采样频率点序列。具体地,

    A:对数频率采样

    [1,1.07,...,18.2,19.05,19.33,19.50,19.62,...,144.54,148.06,151.83,
    154.88,155.69,156.72,160.75,164.89,...,1000]

    B:线性频率采样

    [1,10.99,19.05,19.33,19.62,20.98,,...,148.06,150.85,151.83,155.69,
    156.72,160.75,160.84,164.89,...,1000]

    步骤六,通过底层程序??槎宰钪盏拇莺裳蛄薪写莺导扑?,
    以得到最终离散数据形式的传递函数图。请参阅图3,图3为线性采样方式,采样点数为101
    时,分别采用本发明较佳实施方式提供的高精度的传递函数采样方法及现有方法获得的多
    个不同的输入Fi-输出ao传递函数幅频曲线的对比图。

    采用本发明较佳实施方式提供的高精度的传递函数采样方法,分别设置常规传递
    函数采样方式为对数频率采样和线性频率采样,采样带宽为1Hz~1000Hz,采样点数为
    1001,得到最终传递函数图。图4为采样点数为1001时,本发明提供的方法与现有方法的输
    入Fi-输出ao传递函数幅频曲线的对比图,图4(a)、4(b)为对数采样方式,其中图4(b)为第
    二、三峰值区间局部放大对比图,图4(c)、4(d)为采用线性采样方式时,分别采用本发明的
    方法及现有方法获得输入的Fi-输出ao传递函数幅频曲线的对比图,其中图4(d)为第二、三
    峰值区间局部放大对比图。

    请参阅图5,图5(a)为对数采样方式,采样点数分别为101、1001时,采用现有方法
    获得的输入Fi-输出ao传递函数幅频曲线的对比图,图5(b)为采样点数分别为101、1001时,
    采用现有方法获得的输入Fi-输出ao传递函数幅频曲线的第二、三峰值区间局部放大对比
    图。图5(c)、5(d)为采用线性采样方式,采样点数分别为101、1001时,采用现有方法获得的
    输入Fi-输出ao传递函数幅频曲线的对比图,其中图5(d)为第二、三峰值区间局部放大对比
    图。

    请参阅图6,图6(a)、6(b)为采用对数采样方式,采样点数分别为101、1001时,采用
    本发明提供的方法获得的输入Fi-输出ao传递函数幅频曲线的对比图,其中图6(b)为第二、
    三峰值区间的局部放大对比图。图6(c)、6(d)为采用线性采样方式,采样点数分别为101、
    1001时,采用本发明提供的方法获得的输入Fi-输出ao传递函数幅频曲线对比图,其中图6
    (d)为第二、三峰值区间局部放大对比图。

    根据不同采样方式(对数/线性),不同采样点数(101/1001)下,分别采用现有采样
    方法和本发明提供的方法获得的计算结果,对传递函数频带范围内的第一、二、三峰值的频
    率和幅值误差进行统计,如表1所示。

    表1不同采样条件下传递函数误差对比表





    自图5(a)、5c)和表1可以看出,现有方法的采样精度与采样频率(或采样点数)的
    高低有直接关系,在采样频率较低的时候,采样精度差,严重时甚至出现重叠峰遗漏的情
    况,如当采样点数为101时,相近双峰误采为单峰,且峰值频率和幅值均存在很大误差。同
    时,对比图5(a)与5(c),结合表1内的峰值误差数据可以看出,现有采样方法的采样精度也
    与其采样方式直接相关,当为对数采样方式时,高频区间传递函数的采样精度较低,而当为
    线性采样方式时,则低频区间传递函数的采样精度相对较低。

    通过图6可以看出,相比于现有采样方法,采样方式以及采样频率的高低对于本发
    明较佳实施方式提供的高精度的传递函数采样方法的采样精度影响较小,无论线性采样方
    式还是对数采样方式,在低采样频率下,所述的高精度的传递函数采样方法依旧能够比较
    准确的获取共振峰的频率及幅值,从而达到频响曲线高精度获取的目的??杉?,本发明较佳
    实施方式提供的高精度的传递函数采样方法基本不受采样方式的选择以及采样频率的高
    低的影响,始终能够提供正确有效的分析结果,提高了采样精度,有利于机械系统结构动力
    学振动计算与分析。

    本发明较佳实施方式提供的高精度的传递函数采样方法,其结合系统模态分析的
    数据(模态频率),在低频率条件下,通过细化采样以增加可能共振点(即系统模态频率点)
    附近采样数据点的密集程度,在总采样点数少量增加的情况下高精度的获取了传递函数最
    重要的特征、共振频率点以及共振峰值,有利于机械系统结构动力学振动计算与分析。

    本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以
    限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含
    在本发明的?;し段е?。

    关于本文
    本文标题:一种高精度的传递函数采样方法.pdf
    链接地址://www.4mum.com.cn/p-6021051.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
  • 幸运飞艇两面盘怎么算中 pk10赛车计划手机版 广东11选五推荐 黑龙江时时号码 黑马计划软件免费版 pk10全天免费手机计划 广东11选五计划软件下载 腾讯体育女篮直播 内蒙古时时开奖走势图百度百度贴吧 如何购买大小单双最合理 pk10极速赛车计划软件 田广双色球预测专家 新疆时时走势图二星 今期36码 彩票阶梯投注方案 广东快乐十分稳赚技巧