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

    重庆时时彩过年停几天: 一种火箭橇试验连续振动传递谱确定方法.pdf

    关 键 词:
    一种 火箭 试验 连续 振动 传递 确定 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201610423010.5

    申请日:

    2016.06.14

    公开号:

    CN106404159A

    公开日:

    2017.02.15

    当前法律状态:

    授权

    有效性:

    有权

    法律详情: 授权|||实质审查的生效IPC(主分类):G01H 17/00申请日:20160614|||公开
    IPC分类号: G01H17/00 主分类号: G01H17/00
    申请人: 北京航天控制仪器研究所
    发明人: 魏宗康; 董小燕; 魏珍虹
    地址: 100854 北京市海淀区北京142信箱403分箱
    优先权:
    专利代理机构: 中国航天科技专利中心 11009 代理人: 臧春喜
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201610423010.5

    授权公告号:

    ||||||

    法律状态公告日:

    2019.03.12|||2017.03.15|||2017.02.15

    法律状态类型:

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

    摘要

    一种火箭橇试验连续振动传递谱确定方法,属于力学环境试验技术领域。步骤如下:(1)通过传感器测量出设备X和设备Y在运行时的加速度值序列aIx(k)、aIy(k);(2)把aIx(k)、aIy(k)作为输入量,经过高通滤波后求得振动值序列adx(k),ady(k);(3)基于AR模型求解adx(k)和ady(k)的振动功率谱的连续函数Px(f)、Py(f);(4)修正振动功率谱连续函数Px(f)、Py(f)的幅值,得到P′x(f)、P′y(f);(5)计算设备X至设备Y的连续振动传递谱该方法能够直接获取高频振动信息,有效去除了噪声,实现了随频率连续变化的各级振动传递特性的量化描述,精度高。

    权利要求书

    1.一种火箭橇试验连续振动传递谱确定方法,其特征在于步骤如下:
    (1)在火箭橇试验过程中,分别通过传感器采集设备X和设备Y运行时的加速度值,得到
    设备X的加速度值序列aIx(k)和设备Y的加速度值序列aIy(k),k为采样时刻,k=1,2,…,N,
    所述设备X和设备Y的加速度值均包含随机振动信息和过载信息,所述传感器采样频率fs满
    足4KHz<fs<10KHz,采样周期为Ts=1/fs;
    (2)采用高通滤波器对加速度值序列aIx(k)进行滤波,得到设备X的振动值序列adx(k);
    采用高通滤波器对加速度值序列aIy(k)进行滤波,得到设备Y的振动值序列ady(k);
    (3)基于振动AR模型计算adx(k)的振动功率谱连续函数Px(f)以及ady(k)的振动功率谱
    连续函数Py(f);
    (4)分别对步骤(3)得到的振动功率谱进行幅值修正,得到修正后的振动功率谱连续函
    数P′x(f)和P′y(f);
    (5)利用公式确定设备X至设备Y的连续振动传递谱Φxy(f),所述设备X
    至设备Y的连续振动传递谱为橇体至被试产品的连续振动传递谱、橇体至减振平台的连续
    振动传递谱或减振平台至被试产品的连续振动传递谱。
    2.根据权利要求1所述的一种火箭橇试验连续振动传递谱确定方法,其特征在于:所述
    步骤(2)中高通滤波器对加速度值序列进行滤波得到振动值序列的实现方法为:
    (2.1)定义高通滤波器的转折频率为fz,阻尼系数为ξ,则高通滤波器的系数为
    <mrow> <msub> <mi>b</mi> <mn>0</mn> </msub> <mo>=</mo> <mfrac> <mrow> <mn>4</mn> <mo>/</mo> <msup> <mrow> <mo>(</mo> <mn>2</mn> <msub> <mi>&pi;f</mi> <mi>z</mi> </msub> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mn>4</mn> <mi>&xi;</mi> <mo>/</mo> <mrow> <mo>(</mo> <mn>2</mn> <msub> <mi>&pi;f</mi> <mi>z</mi> </msub> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mn>4</mn> <mo>/</mo> <msup> <mrow> <mo>(</mo> <mn>2</mn> <msub> <mi>&pi;f</mi> <mi>z</mi> </msub> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mn>4</mn> <mi>&xi;</mi> <mo>/</mo> <mrow> <mo>(</mo> <mn>2</mn> <msub> <mi>&pi;f</mi> <mi>z</mi> </msub> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mn>1</mn> </mrow> </mfrac> <mo>;</mo> </mrow>
    <mrow> <msub> <mi>n</mi> <mn>1</mn> </msub> <mo>=</mo> <mo>-</mo> <mfrac> <mn>2</mn> <mrow> <mn>1</mn> <mo>+</mo> <mn>2</mn> <msub> <mi>&pi;f</mi> <mi>z</mi> </msub> <msub> <mi>&xi;T</mi> <mi>s</mi> </msub> </mrow> </mfrac> <mo>;</mo> </mrow>
    <mrow> <msub> <mi>n</mi> <mn>2</mn> </msub> <mo>=</mo> <mfrac> <mrow> <mn>1</mn> <mo>-</mo> <mn>2</mn> <msub> <mi>&pi;f</mi> <mi>z</mi> </msub> <msub> <mi>&xi;T</mi> <mi>s</mi> </msub> </mrow> <mrow> <mn>1</mn> <mo>+</mo> <mn>2</mn> <msub> <mi>&pi;f</mi> <mi>z</mi> </msub> <msub> <mi>&xi;T</mi> <mi>s</mi> </msub> </mrow> </mfrac> <mo>;</mo> </mrow>
    <mrow> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>=</mo> <mfrac> <mrow> <mn>2</mn> <mo>-</mo> <mn>8</mn> <mo>/</mo> <msup> <mrow> <mo>(</mo> <mn>2</mn> <msub> <mi>&pi;f</mi> <mi>z</mi> </msub> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> <mrow> <mn>4</mn> <mo>/</mo> <msup> <mrow> <mo>(</mo> <mn>2</mn> <msub> <mi>&pi;f</mi> <mi>z</mi> </msub> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mn>4</mn> <mi>&xi;</mi> <mo>/</mo> <mrow> <mo>(</mo> <mn>2</mn> <msub> <mi>&pi;f</mi> <mi>z</mi> </msub> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mn>1</mn> </mrow> </mfrac> <mo>;</mo> </mrow>
    <mrow> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>=</mo> <mfrac> <mrow> <mn>4</mn> <mo>/</mo> <msup> <mrow> <mo>(</mo> <mn>2</mn> <msub> <mi>&pi;f</mi> <mi>z</mi> </msub> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <mn>4</mn> <mi>&xi;</mi> <mo>/</mo> <mrow> <mo>(</mo> <mn>2</mn> <msub> <mi>&pi;f</mi> <mi>z</mi> </msub> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mn>1</mn> </mrow> <mrow> <mn>4</mn> <mo>/</mo> <msup> <mrow> <mo>(</mo> <mn>2</mn> <msub> <mi>&pi;f</mi> <mi>z</mi> </msub> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mn>4</mn> <mi>&xi;</mi> <mo>/</mo> <mrow> <mo>(</mo> <mn>2</mn> <msub> <mi>&pi;f</mi> <mi>z</mi> </msub> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mn>1</mn> </mrow> </mfrac> <mo>;</mo> </mrow>
    (2.2)k=1时,高通滤波器的状态变量x1(1)=0、x2(1)=0;
    (2.3)高通滤波器利用公式ad(k)=(n1-d1)x1(k)+(n2-d2)x2(k)+b0aI(k)对第k个时刻的
    加速度值aI(k)进行滤波得到振动值ad(k),进入步骤(2.4);
    (2.4)k的值加1,判断k是否小于等于N时,如果是,则利用公式
    更新x1(k)和x2(k),返回步骤(2.3);否则,高
    通滤波器滤波结束,得到振动值序列。
    3.根据权利要求1所述的一种火箭橇试验连续振动传递谱确定方法,其特征在于:所述
    步骤(3)的实现方式为:
    (3.1)设adx(k)的自相关序列为Rx(i),ady(k)的自相关序列为Ry(i),i=0,1,2,…,N-1,
    Rx(i)和Ry(i)的计算公式如下:
    <mrow> <msub> <mi>R</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mi>N</mi> <mo>-</mo> <mi>i</mi> </mrow> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>l</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mi>i</mi> </mrow> </munderover> <msub> <mi>a</mi> <mrow> <mi>d</mi> <mi>x</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>l</mi> <mo>)</mo> </mrow> <msub> <mi>a</mi> <mrow> <mi>d</mi> <mi>x</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>l</mi> <mo>+</mo> <mi>i</mi> <mo>)</mo> </mrow> </mrow>
    <mrow> <msub> <mi>R</mi> <mi>y</mi> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mi>N</mi> <mo>-</mo> <mi>i</mi> </mrow> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>l</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mi>i</mi> </mrow> </munderover> <msub> <mi>a</mi> <mrow> <mi>d</mi> <mi>y</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>l</mi> <mo>)</mo> </mrow> <msub> <mi>a</mi> <mrow> <mi>d</mi> <mi>y</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>l</mi> <mo>+</mo> <mi>i</mi> <mo>)</mo> </mrow> </mrow>
    (3.2)取q<N,根据自相关序列Rx(i)计算adx(k)的q阶AR模型的各项系数a1、a2、…、aq,根
    据自相关序列Ry(i)计算ady(k)的q阶AR模型的各项系数b1、b2、…、bq,计算公式如下:
    <mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>a</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>a</mi> <mi>2</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>a</mi> <mi>q</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mo>-</mo> <msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>R</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <msub> <mi>R</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <mrow> <msub> <mi>R</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mi>q</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>R</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <msub> <mi>R</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <mrow> <msub> <mi>R</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mi>q</mi> <mo>-</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mrow></mrow> </mtd> <mtd> <mrow></mrow> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mrow></mrow> </mtd> <mtd> <mrow></mrow> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mrow></mrow> </mtd> <mtd> <mrow></mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>R</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mi>q</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <msub> <mi>R</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mi>q</mi> <mo>-</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <mrow> <msub> <mi>R</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>R</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>R</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>R</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mi>q</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
    σx2=Rx(0)
    <mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>b</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>b</mi> <mi>2</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>b</mi> <mi>q</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mo>-</mo> <msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>R</mi> <mi>y</mi> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <msub> <mi>R</mi> <mi>y</mi> </msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <mrow> <msub> <mi>R</mi> <mi>y</mi> </msub> <mrow> <mo>(</mo> <mi>q</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>R</mi> <mi>y</mi> </msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <msub> <mi>R</mi> <mi>y</mi> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <mrow> <msub> <mi>R</mi> <mi>y</mi> </msub> <mrow> <mo>(</mo> <mi>q</mi> <mo>-</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mrow></mrow> </mtd> <mtd> <mrow></mrow> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mrow></mrow> </mtd> <mtd> <mrow></mrow> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mrow></mrow> </mtd> <mtd> <mrow></mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>R</mi> <mi>y</mi> </msub> <mrow> <mo>(</mo> <mi>q</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <msub> <mi>R</mi> <mi>y</mi> </msub> <mrow> <mo>(</mo> <mi>q</mi> <mo>-</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <mrow> <msub> <mi>R</mi> <mi>y</mi> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>R</mi> <mi>y</mi> </msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>R</mi> <mi>y</mi> </msub> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>R</mi> <mi>y</mi> </msub> <mrow> <mo>(</mo> <mi>q</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
    σy2=Ry(0);
    (3.3)设adx(k)为受方差为1的白噪声序列u1(n)激励后的结果,则基于振动AR模型计算
    adx(k)的振动功率谱连续函数Px(f)的公式为:
    <mrow> <msub> <mi>P</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mi>f</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <mo>&CenterDot;</mo> <mfrac> <msubsup> <mi>&sigma;</mi> <mi>x</mi> <mn>2</mn> </msubsup> <mrow> <mo>|</mo> <mn>1</mn> <mo>+</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>q</mi> </munderover> <mrow> <msub> <mi>a</mi> <mi>m</mi> </msub> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mn>2</mn> <mi>&pi;</mi> <mi>f</mi> <mi>m</mi> </mrow> </msup> </mrow> <msup> <mo>|</mo> <mn>2</mn> </msup> </mrow> </mfrac> </mrow>
    其中q为振动AR模型的阶次,f为频率,j为实数单位,j2=-1;
    (3.4)设ady(k)为受方差为1的白噪声序列u2(n)激励后的结果,则基于振动AR模型计算
    ady(k)的振动功率谱连续函数Py(f)的公式为:
    <mrow> <msub> <mi>P</mi> <mi>y</mi> </msub> <mrow> <mo>(</mo> <mi>f</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <mo>&CenterDot;</mo> <mfrac> <msubsup> <mi>&sigma;</mi> <mi>y</mi> <mn>2</mn> </msubsup> <mrow> <mo>|</mo> <mn>1</mn> <mo>+</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>q</mi> </munderover> <mrow> <msub> <mi>b</mi> <mi>m</mi> </msub> <msup> <mi>z</mi> <mrow> <mo>-</mo> <mi>j</mi> <mn>2</mn> <mi>&pi;</mi> <mi>f</mi> <mi>m</mi> </mrow> </msup> </mrow> <msup> <mo>|</mo> <mn>2</mn> </msup> </mrow> </mfrac> <mo>.</mo> </mrow>
    4.根据权利要求3所述的一种火箭橇试验连续振动传递谱确定方法,其特征在于:所述
    步骤(4)的实现方式为:
    (4.1)取fp=(p-1)/T,其中p=1,2,…,N/2,代入步骤(3.3)和(3.4)可分别求得adx(k)
    的振动功率谱的功率谱序列Px(fp)以及ady(k)的振动功率谱序列Py(fp),T为N个采样周期的
    持续时间;
    (4.2)根据adx(k)得到设备X的振动功率谱密度序列Φx(fp),根据ady(k)得到设备Y的振
    动功率谱密度序列Φy(fp);
    (4.3)利用公式计算设备X的振动功率谱连续函数的幅值修正系数cx,
    利用公式计算设备Y的振动功率谱连续函数的幅值修正系数cy;
    (4.4)设备X修正后的振动功率谱连续函数P′x(f)=cxPx(f),设备Y修正后的振动功率
    谱连续函数P′y(f)=cyPy(f)。
    5.根据权利要求4所述的一种火箭橇试验连续振动传递谱确定方法,其特征在于:所述
    步骤(4.2)中根据adx(k)得到设备X的振动功率谱密度序列Φx(fp)的实现方式为:
    (i)基于正频域,利用如下公式计算设备X振动数据的线谱序列cx(fp):
    <mrow> <msub> <mi>c</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>p</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>2</mn> <mi>N</mi> </mfrac> <mi>F</mi> <mrow> <mo>(</mo> <mi>p</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>2</mn> <mi>N</mi> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>a</mi> <mrow> <mi>d</mi> <mi>x</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mi>p</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mn>2</mn> <mi>&pi;</mi> <mo>/</mo> <mi>N</mi> </mrow> </msup> <mo>,</mo> <mi>g</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>N</mi> <mo>/</mo> <mn>2</mn> </mrow>
    其中fp为频率,fp=(p-1)/T,p=1,2,…,N/2;
    (ii)根据线谱序列计算设备X的振动功率谱密度序列Φx(fp):
    <mrow> <msub> <mi>&Phi;</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>p</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msub> <mi>c</mi> <mi>x</mi> </msub> <msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>p</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mi>T</mi> <mo>,</mo> <mi>p</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>N</mi> <mo>/</mo> <mn>2</mn> </mrow>
    其中F(p)为离散傅立叶变换;
    根据ady(k)得到设备Y的振动功率谱密度序列Φy(fp)的实现方式为:
    (i)利用如下公式计算设备Y的振动数据的线谱序列cy(fp):
    <mrow> <msub> <mi>c</mi> <mi>y</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>p</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>2</mn> <mi>N</mi> </mfrac> <mi>F</mi> <mrow> <mo>(</mo> <mi>p</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>2</mn> <mi>N</mi> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>a</mi> <mrow> <mi>d</mi> <mi>y</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mi>p</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mn>2</mn> <mi>&pi;</mi> <mo>/</mo> <mi>N</mi> </mrow> </msup> <mo>,</mo> <mi>p</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>N</mi> <mo>/</mo> <mn>2</mn> </mrow>
    (ii)根据线谱序列计算设备Y的振动功率谱密度序列Φy(fp):
    <mrow> <msub> <mi>&Phi;</mi> <mi>y</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>p</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msub> <mi>c</mi> <mi>y</mi> </msub> <msup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>p</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mi>T</mi> <mo>,</mo> <mi>p</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>N</mi> <mo>/</mo> <mn>2.</mn> </mrow>

    说明书

    一种火箭橇试验连续振动传递谱确定方法

    技术领域

    本发明涉及一种火箭橇试验连续振动传递谱确定方法,属于力学环境试验技术领
    域。

    背景技术

    火箭橇是用火箭发动机作为推力,在轨道上高速运行的滑橇。由于火箭橇在运行
    过程中可产生大过载,因此,不同类型的产品都安装在橇体上开展火箭橇试验。

    火箭橇试验过程中,被试产品根据试验要求有两种安装方式,第一种方式是被试
    产品直接安装于橇体上;第二种方式是被试产品安装于减振平板上,而减振平板通过金属
    减振器安装于橇体上。通过振动测试系统实时测量和记录橇体和被试产品(第二种方式的
    话还包括减振平板)的振动和冲击信号,待试验结束后读取记录存储的数据并进行数据处
    理。

    振动测试系统由传感器、数据采集记录系统、触发装置和信号传输电缆等组成。测
    试系统中传感器、数据采集系统和触发装置通过信号传输电缆连接,各部件在橇体上的安
    装可采用压条、螺栓紧固的方式,具体安装部位包括:

    a)在橇体、减振平板和被试产品壳体上分别安装不少于3个输入轴正交的传感器,
    分别测量安装位置处的航向加速度、竖向加速度和侧向加速度;

    b)数据采集系统、触发装置直接安装于橇体上;

    c)信号传输电缆根据布线方式采用捆扎成束就近与橇体固连。

    目前,在数据处理过程中存在以下问题:

    (1)由于传感器测量的加速度信息包括了低频过载信息和高频振动信息,为得到
    振动信息,首先采用低通滤波器对加速度信息进行滤波得到过载信息,然后,加速度信息减
    去过载信息得到振动信息。而缺少实现高频振动信息的直接获取方法;

    (2)对各级振动信息只开展了周期功率谱分析,特点是噪声大,缺少有效的降噪处
    理方法;

    (3)只是定性分析各级的振动量级大小、有无谐振峰值等,不能对橇体至减振平
    台、减振平台至被试产品、橇体至被试产品的振动传递特性进行分析;

    (4)振动传递只能实现离散频率点,而不能实现连续频率段的分析,因此,分辨率
    较低。

    发明内容

    本发明解决的技术问题是:克服现有技术的不足,提供一种火箭橇试验连续振动
    传递谱确定方法,能够直接获取高频振动信息,有效去除了噪声,实现了各级振动传递特性
    的连续的量化描述,精度高。

    本发明的技术方案是:一种火箭橇试验连续振动传递谱确定方法,步骤如下:

    (1)在火箭橇试验过程中,分别通过传感器采集设备X和设备Y运行时的加速度值,
    得到设备X的加速度值序列aIx(k)和设备Y的加速度值序列aIy(k),k为采样时刻,k=1,
    2,…,N,所述设备X和设备Y的加速度值均包含随机振动信息和过载信息,所述传感器采样
    频率fs满足4KHz<fs<10KHz,采样周期为Ts=1/fs;

    (2)采用高通滤波器对加速度值序列aIx(k)进行滤波,得到设备X的振动值序列adx
    (k);采用高通滤波器对加速度值序列aIy(k)进行滤波,得到设备Y的振动值序列ady(k);

    (3)基于振动AR模型计算adx(k)的振动功率谱连续函数Px(f)以及ady(k)的振动功
    率谱连续函数Py(f);

    (4)分别对步骤(3)得到的振动功率谱进行幅值修正,得到修正后的振动功率谱连
    续函数P′x(f)和P′y(f);

    (5)利用公式确定设备X至设备Y的连续振动传递谱Φxy(f),所述设
    备X至设备Y的连续振动传递谱为橇体至被试产品的连续振动传递谱、橇体至减振平台的连
    续振动传递谱或减振平台至被试产品的连续振动传递谱。

    所述步骤(2)中高通滤波器对加速度值序列进行滤波得到振动值序列的实现方法
    为:

    (2.1)定义高通滤波器的转折频率为fz,阻尼系数为ξ,则高通滤波器的系数为






    (2.2)k=1时,高通滤波器的状态变量x1(1)=0、x2(1)=0;

    (2.3)高通滤波器利用公式ad(k)=(n1-d1)x1(k)+(n2-d2)x2(k)+b0aI(k)对第k个时
    刻的加速度值aI(k)进行滤波得到振动值ad(k),进入步骤(2.4);

    (2.4)k的值加1,判断k是否小于等于N时,如果是,则利用公式
    更新x1(k)和x2(k),返回步骤(2.3);否则,高
    通滤波器滤波结束,得到振动值序列。

    所述步骤(3)的实现方式为:

    (3.1)设adx(k)的自相关序列为Rx(i),ady(k)的自相关序列为Ry(i),i=0,1,2,…,
    N-1,Rx(i)和Ry(i)的计算公式如下:



    (3.2)取q<N,根据自相关序列Rx(i)计算adx(k)的q阶AR模型的各项系数a1、a2、…、
    aq,根据自相关序列Ry(i)计算ady(k)的q阶AR模型的各项系数b1、b2、…、bq,计算公式如下:


    σx2=Rx(0)


    σy2=Ry(0);

    (3.3)设adx(k)为受方差为1的白噪声序列u1(n)激励后的结果,则基于振动AR模型
    计算adx(k)的振动功率谱连续函数Px(f)的公式为:


    其中q为振动AR模型的阶次,f为频率,j为实数单位,j2=-1;

    (3.4)设ady(k)为受方差为1的白噪声序列u2(n)激励后的结果,则基于振动AR模型
    计算ady(k)的振动功率谱连续函数Py(f)的公式为:


    所述步骤(4)的实现方式为:

    (4.1)取fp=(p-1)/T,其中p=1,2,…,N/2,代入步骤(3.3)和(3.4)可分别求得adx
    (k)的振动功率谱的功率谱序列Px(fp)以及ady(k)的振动功率谱序列Py(fp),T为N个采样周
    期的持续时间;

    (4.2)根据adx(k)得到设备X的振动功率谱密度序列Φx(fp),根据ady(k)得到设备Y
    的振动功率谱密度序列Φy(fp);

    (4.3)利用公式计算设备X的振动功率谱连续函数的幅值修正系
    数cx,利用公式计算设备Y的振动功率谱连续函数的幅值修正系数cy;

    (4.4)设备X修正后的振动功率谱连续函数Px′(f)=cxPx(f),设备Y修正后的振动
    功率谱连续函数Py′(f)=cyPy(f)。

    所述步骤(4.2)中根据adx(k)得到设备X的振动功率谱密度序列Φx(fp)的实现方
    式为:

    (i)基于正频域,利用如下公式计算设备X振动数据的线谱序列cx(fp):


    其中fp为频率,fp=(p-1)/T,p=1,2,…,N/2;

    (ii)根据线谱序列计算设备X的振动功率谱密度序列Φx(fp):


    其中F(p)为离散傅立叶变换;

    根据ady(k)得到设备Y的振动功率谱密度序列Φy(fp)的实现方式为:

    (i)利用如下公式计算设备Y的振动数据的线谱序列cy(fp):


    (ii)根据线谱序列计算设备Y的振动功率谱密度序列Φy(fp):


    本发明与现有技术相比的有益效果是:

    (1)本发明采用基于AR模型谱估计的方法,给出了一种精确描述火箭橇试验连续
    振动传递谱的计算方法,克服了以前各级之间的振动传递定性描述或只实现离散频率点,
    而不能实现连续频率段分析的不足,实现了对橇体至减振平台、减振平台至被试产品、橇体
    至被试产品随频率连续变化的振动传递特性的量化描述,精度高。

    (2)本发明传感器采集的加速度值通过高通滤波器滤波后,实现了对高频振动信
    息的直接提取,相对于先提取过载再求差的方法,提高了数据的转换效率。

    (3)本发明采用振动AR模型及其幅值修正方法得到功率谱连续函数,既可以精确
    给出功率谱在各频率点的峰值和谷值信息,又可精确给出幅值信息,实现了火箭橇试验中
    无分辨率限制的振动传递谱精确描述。

    (4)相对传统的周期功率谱具有噪声大的特点,本发明通过振动AR模型及其幅值
    修正方法对振动数据进行降噪处理,有效减小了振动的噪声。

    附图说明

    图1本方法的流程图;

    图2在火箭橇试验中被试品基于周期谱的功率谱;

    图3在火箭橇试验中被试品基于AR模型修正前的初始功率谱;

    图4在火箭橇试验中被试品基于AR模型修正后的精确功率谱;

    图5为火箭橇试验中橇体的振动功率谱密度;

    图6为火箭橇试验中减振平板的振动功率谱密度;

    图7为火箭橇试验中被试产品上的振动功率谱密度;

    图8为火箭橇试验中减振平板在橇体振动激励下的振动传递谱分布;

    图9为火箭橇试验中被试产品相对减振平板振动激励下的振动传递谱分布;

    图10为火箭橇试验中被试产品在橇体振动激励下的振动传递谱分布。

    具体实施方式

    下面结合附图1对本发明做详细说明。

    火箭橇试验连续振动传递谱包括橇体至被试产品的连续振动传递谱、橇体至减振
    平台的连续振动传递谱或减振平台至被试产品的连续振动传递谱。因此,如图1所示,本发
    明提出一种火箭橇试验过程中,设备X至设备Y的连续振动传递谱确定方法,步骤如下:

    1、在火箭橇试验过程中,分别通过传感器采集设备X和设备Y运行时的加速度值,
    得到设备X的加速度值序列aIx(k)和设备Y的加速度值序列aIy(k),k为采样时刻,k=1,
    2,…,N,所述设备X和设备Y的加速度值均包含随机振动信息和过载信息,所述传感器采样
    频率fs满足4KHz<fs<10KHz,采样周期为Ts=1/fs。

    以橇体至被试品的振动传递谱为例,安装于橇体的传感器在每个采样周期采样橇
    体运行时的激振源加速度值,得到橇体运行时的激振源加速度值序列aIx(k);同时,安装于
    被试品的传感器在每个采样周期采样被试品运行时的加速度值,得到被试品运行时的加速
    度值序列aIy(k)。

    2、采用高通滤波器对加速度值序列aIx(k)进行滤波,得到设备X的振动值序列adx
    (k);采用高通滤波器对加速度值序列aIy(k)进行滤波,得到设备Y的振动值序列ady(k),实
    现方法为:

    (2.1)定义高通滤波器的转折频率为fz,阻尼系数为ξ,则高通滤波器的系数为






    (2.2)k=1时,高通滤波器的状态变量x1(1)=0、x2(1)=0;

    (2.3)高通滤波器利用公式ad(k)=(n1-d1)x1(k)+(n2-d2)x2(k)+b0aI(k)对第k个时
    刻的加速度值aI(k)进行滤波得到振动值ad(k),进入步骤(2.4);

    (2.4)k的值加1,判断k是否小于等于N时,如果是,则利用公式
    更新x1(k)和x2(k),返回步骤(2.3);否则,高
    通滤波器滤波结束,得到振动值序列。

    3、计算adx(k)基于AR模型的振动功率谱连续函数Px(f)以及ady(k)基于AR模型的振
    动功率谱连续函数Py(f),方法为:

    (3.1)设adx(k)的自相关序列为Rx(i),ady(k)的自相关序列为Ry(i),i=0,1,2,…,
    N-1,Rx(i)和Ry(i)的计算公式如下:



    (3.2)取q<N,根据自相关序列Rx(i)计算adx(k)的q阶AR模型的各项系数a1、a2、…、
    aq,根据自相关序列Ry(i)计算ady(k)的q阶AR模型的各项系数b1、b2、…、bq,计算公式如下:


    σx2=Rx(0)


    σy2=Ry(0);

    (3.3)设adx(k)为受方差为1的白噪声序列u1(n)激励后的结果,则基于振动AR模型
    (自回归模型,auto-regressive)计算adx(k)的振动功率谱连续函数Px(f)的公式为:


    其中q为设备X振动AR模型的阶次,f为频率,j为实数单位,j2=-1。

    (3.4)设ady(k)为受方差为1的白噪声序列u2(n)激励后的结果,则基于振动AR模型
    计算ady(k)的振动功率谱连续函数Py(f)的公式为:


    4、对步骤(3)得到的功率谱进行幅值修正,得到修正后的功率谱连续函数P′x(f)
    和P′y(f),方法为:

    (4.1)取fp=(p-1)/T,其中p=1,2,…,N/2,可分别求得设备X振动AR模型的功率
    谱序列Px(fp)以及Y振动AR模型的功率谱序列Py(fp);

    (4.2)根据adx(k)得到设备X的振动功率谱密度序列Φx(fp),根据ady(k)得到设备Y
    的振动值功率谱密度序列Φy(fp);

    根据adx(k)得到设备X的振动功率谱密度序列Φx(fp)的实现方式为:

    (a)基于正频域,利用如下公式计算设备X振动数据的线谱序列cx(fp):


    其中fp为频率,fp=(p-1)/T,p=1,2,…,N/2;

    (b)根据线谱序列计算设备X的振动功率谱密度序列Φx(fp):


    其中T为N个采样周期的持续时间;F(p)为离散傅立叶变换;

    根据ady(k)得到设备Y的振动功率谱密度序列Φy(fp)的实现方式为:

    (a)利用如下公式计算设备Y的振动数据的线谱序列cy(fp):


    (b)根据线谱序列计算设备Y的振动功率谱密度序列Φy(fp):


    (4.3)利用公式计算设备X的振动功率谱连续函数的幅值修正系
    数cx,利用公式计算设备Y的振动功率谱连续函数的幅值修正系数cy;

    (4.4)设备X修正后的振动功率谱连续函数Px′(f)=cxPx(f),设备Y修正后的振动
    功率谱连续函数Py′(f)=cyPy(f)。

    5、根据步骤(4)得到的功率谱连续函数P′x(f)、P′y(f),利用公式
    确定设备X至设备Y的一个连续振动传递谱Φxy(f),其中设备X至设备Y的连续振动传递谱为
    橇体至被试产品的连续振动传递谱、橇体至减振平台的连续振动传递谱或减振平台至被试
    产品的连续振动传递谱。

    本发明为了得到连续的振动传递谱,在火箭橇试验采用AR模型谱估计。由于AR模
    型是一个有理分式,利用平滑特性可实现对连续频率段的分析,同时在阶次较高时,各频率
    点的波峰和波谷也较为准确。但AR模型的缺点是整体幅值受随机信号是否平稳等因素影响
    而相对真实值有偏离,因此,本发明进一步提出了对AR模型进行幅值修正的方法,以适应火
    箭橇试验的非平稳数据处理要求。从而既可以精确给出功率谱在各频率点的峰值和谷值信
    息,又可精确给出幅值信息,实现了火箭橇试验中,橇体至减振平台、减振平台至被试产品、
    橇体至被试产品无分辨率限制的振动传递谱精确描述,精度高。

    实施例1:

    以被试品为例,验证本发明提出的功率谱幅值修正方法的有效性。

    其中设定条件如下:在火箭橇试验中,传感器采样频率为9.6kHz,被试品基于周期
    谱的振动功率谱如图2所示,基于AR模型的初始功率谱如图3所示,可以看出,二者幅值有些
    许差别。采用本发明方法,功率谱幅值修正系数为c=0.723,把该值代入修正公式后,被试
    品基于AR模型的精确功率谱如图4所示。比较图3和图4可以看出,图4的功率谱曲线相对图3
    的功率谱曲线整体往下移,幅值得到修正,峰谷形状保持不变。比较图4和图2可以看出,采
    用AR模型的功率谱相对平滑,噪声得到衰减。

    实施例2:

    从火箭橇试验中根据多传感器的振动信息计算各级之间的连续振动传递谱,实现
    振动传递特性分析。

    设定条件如下:在火箭橇试验中,采用多级安装方式,即:被试产品安装于减振平
    板上,再由减振平板通过隔振器安装于橇体。传感器采样频率为9.6kHz。在开始运行至停止
    的过程,由橇体安装的航向传感器、减振平板安装的航向传感器以及惯性平台(被试产品)
    安装的航向传感器分别采集记录对应的加速度值信息,经过高通滤波器得到振动数据,并
    进一步得到功率谱密度。

    橇体的振动功率谱密度信息如图5所示、减振平板的振动功率谱密度信息如图6所
    示、惯性平台(被试产品)的振动功率谱密度信息如图7所示。根据本发明的连续振动传递谱
    计算方法,可分别计算出:

    (1)减振平板至橇体(减振平板在橇体振动激励下)的连续振动传递谱分布,如图8
    所示;

    (2)惯性平台(被试产品)至减振平板(惯性平台(被试产品)在减振平板振动激励
    下)的连续振动传递谱分布,如图9所示;

    (3)惯性平台(被试产品)至橇体(惯性平台(被试产品)在橇体振动激励下)的连续
    振动传递谱分布,如图10所示。

    从图中可以看出,惯性平台(被试产品)相对于减振平板的振动量级放大,由于减
    振平板的衰减作用有限,从而造成橇体振动激励时惯性平台(被试产品)的振动量级也放
    大。

    以上所述,仅为本发明一个具体实施方式,但本发明的?;し段Р⒉痪窒抻诖?,任
    何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都
    应涵盖在本发明的?;し段е?。

    本发明未详细说明部分属于本领域技术人员公知常识。

    关于本文
    本文标题:一种火箭橇试验连续振动传递谱确定方法.pdf
    链接地址://www.4mum.com.cn/p-6013955.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码网站多少 3d包胆多少钱 快3大小单双玩法 中国体肓彩票开奖大乐 平推也稳赚的倍投 北京pk10官网开奖视频 北京汽车pk10在哪里玩 北京pk10微信群全天计划 艇走势技巧贴吧 双色球开奖结果下软件 双色球模拟选号投注 pk10长期稳赚的方法 旺旺时时采彩分析计划软件 幸运飞艇6码倍投计划 免费快乐飞艇计划软件