CN101038232A - 一种非高斯随机振动激励信号生成方法及其装置 - Google Patents
一种非高斯随机振动激励信号生成方法及其装置 Download PDFInfo
- Publication number
- CN101038232A CN101038232A CN 200710063654 CN200710063654A CN101038232A CN 101038232 A CN101038232 A CN 101038232A CN 200710063654 CN200710063654 CN 200710063654 CN 200710063654 A CN200710063654 A CN 200710063654A CN 101038232 A CN101038232 A CN 101038232A
- Authority
- CN
- China
- Prior art keywords
- gaussian
- distribution
- function
- random
- signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Abstract
在进行产品的可靠性试验和环境试验的时候,有些动态环境的时间历程具有非高斯分布特性。本发明给出了一种生成非高斯随机振动激励信号的方法,即由给定谱型的服从高斯分布随机信号概率的分布函数与非高斯分布随机信号概率的分布函数求得一个单调递增的非线性函数,并利用这个函数可以将高斯分布的时间序列转化成非高斯分布的时间序列,并且能保证变换后的时间序列的谱型与原信号的谱型保持一致。可以在传统的电动振动台上生成具有给定功率谱密度、斜度和峭度的非高斯随机振动激励信号,可以更真实的模拟产品实际的使用环境,解决在可靠性强化试验中使用全轴振动台时出现的谱型不可控、低频能量不足等问题。
Description
技术领域
本发明属于振动试验领域,特别涉及非高斯振动激励信号的生成方法,主要应用于振动台试验。
背景技术
在进行产品的可靠性试验和环境试验的时候,许多物理系统的行为以及动态环境的时间历程是服从非高斯分布的。由于非高斯分布振动的峰值水平通常会很高,一般来说,具有非高斯振动响应的结构比振动响应为高斯分布的结构容易被破坏。但振动台所模拟的服从非高斯的随机振动激励信号要更严酷,更符合产品的实际使用环境。可以使产品在环境应力筛选中更快的析出缺陷,从而大大的节约试验时间,降低试验成本。当今的数字随机振动控制系统只提供服从高斯(正态)分布的随机振动激励信号,新型的全轴振动台可以产生非高斯分布的随机振动激励信号,但是它所生成的信号存在很多问题,如:全轴振动台低频能量不足,而且在频谱上存在很多低谷。另外,全轴振动台的频谱是由设计决定的,谱形不可控。然而,由于当今的数字随机振动控制系统只提供服从高斯分布的随机振动信号,工程上通常假定所模拟的随机振动是符合高斯分布的。为了更好的模拟产品的外场使用环境,有必要在传统的振动台上实现非高斯分布的随机振动。
发明内容
本发明的目的是提供一种方法及装置可以在传统的振动台上实现非高斯分布的随机振动。即提出了一种利用单调递增的非线性函数产生给定谱密度、斜度和峭度的零均值非高斯随机过程的方法,可用于在电动振动台上得到所需的时间序列,并应用于产品的随机振动试验。
为了实现上述目的,本发明提供了一种新型的非高斯随机振动激励信号生成方法,包括以下步骤:
首先,在信号生成模块中设定所求非高斯分布随机过程的斜度、峭度和功率谱以及高斯分布的均值和方差;
第二步:由设定的高斯分布的均值和方差,再根据给定的功率谱密度条件,按常规的随机振动信号生成方法,得到一个服从高斯分布的时间序列;
第三步:根据产品实际工作时的测试数据或产品任务剖面的具体要求,得到所需要非高斯分布随机信号的时间序列的概率密度函数;通过对所述概率密度函数进行积分得到概率分布函数;
第四步,利用非高斯概率分布函数与已知高斯分布随机信号概率分布函数之间关系导出单调递增非线性函数y=g(x);
第五步:利用y=g(x)与服从高斯分布的时间序列相乘将其转化成服从非高斯分布的时间序列,从而得到一个给定概率分布、斜度、峭度的零均值非高斯随机信号。
第六步,将所生成的非高斯随机信号通过D/A转换器变成连续的模拟量,通过低通滤波器平滑后输给功率放大器,最后经功率放大后驱动电动振动台,实现非高斯分布的随机振动试验。
实现本发明的装置主要包括:传统的电动振动台,该电动振动台包括功率放大器,激振器、试件、传感器,和非高斯随机信号生成模块、D/A转换器、低通滤波器以及功率放大器。非高斯随机信号生成模块是由下述模块实现的:
录入模块:设定并录入所求非高斯分布随机过程的斜度、峭度和功率谱以及高斯分布的均值和方差;
服从高斯分布的时间序列模块:由设定的高斯分布的均值和方差,再根据给定的功率谱密度条件,按常规的随机振动信号生成方法,得到一个服从高斯分布的时间序列;
非高斯概率分布函数模块:根据产品实际工作时的测试数据或产品任务剖面的具体要求,得到所需非高斯分布随机信号的时间序列的概率密度函数;通过对所述概率密度函数进行积分得到概率分布函数;
单调递增函数生成模块:利用非高斯概率分布函数与已知高斯分布随机信号概率分布函数之间关系导出单调递增非线性函数y=g(x);
转化模块:利用y=g(x)与服从高斯分布的时间序列相乘将其转化成服从非高斯分布的时间序列,从而得到一个给定概率分布、斜度、峭度的零均值非高斯随机信号。
非高斯随机信号生成模块生成非高斯信号之后,其输出的信号通过D/A转换器变成连续的模拟量,通过低通滤波器平滑后输给功率放大器,最后经功率放大后驱动振动台。
根据本发明可以得到非高斯随机振动控制系统,并在传统的电动振动台上实现给定谱型(频谱可控)、斜度、峭度的非高斯随机振动;解决当前应用于可靠性强化试验的全轴振动台频谱不可控,低频能量不足等问题。利用本发明可将高斯分布的时间序列转化成非高斯分布的时间序列,并且能保证变换后的时间序列的谱型与原信号的谱型保持一致。
附图说明
图1为典型随机振动试验的参考谱;
图2为非高斯随机振动信号生成框图;
图3为非高斯随机振动信号生成装置框图;
图4为根据本发明的方法生成的斜度为0,峭度为10的非高斯分布时间序列。
具体实施方式
首先,如图1、2所示,在给定参考谱如图1所示的条件下,并在给定高斯分布的均值和方差的条件下,按照传统的随机振动信号生成方法得到高斯分布时间序列,生成功率谱密度为P的服从高斯分布时间序列{xi}。
第二步,根据产品工作的实际数据可以得到所求的非高斯随机信号的概率密度,对这个概率密度进行积分,可得到其概率分布函数fY(y)
第三步,在高斯分布与非高斯分布的概率分布函数之间可以导出一个函数关系,如下式所示:
其中,X为服从高斯分布的随机变量;
Y为服从非高斯分布的随机变量;
fY(y)为随机变量Y的概率分布函数;
fX(x)为随机变量X的概率分布函数;
在得到dx/dy之后可以求得dy/dx,将其进行积分,则得到y=g(x)。
y=g(x)为单调递增非线性的函数,即
所求得的函数y=g(x)一般都是单调递增的,并且比较光滑,可以保留原始高斯波形的各种信息,从而能保证转换后谱形变化不大。
第四步:利用y=g(x)与服从高斯分布的时间序列相乘将其转化成服从非高斯分布的时间序列,从而得到一个给定概率分布、斜度、峭度的零均值非高斯随机信号。即通过设定fY(y)可以得到任意斜度S3和峭度K4。斜度通常用于评价随机过程的对称程度,峭度通常用来度量随机过程的峰值程度。斜度和峭度的定义由下式给出,
E(yn)为随机过程Y的n阶中心矩
对于高斯分布的时间历程的S3=0,K4=3,而对于非高斯分布时间历程K4>3。
由于此种方法不会改变原高斯信号的谱密度,所以不必进行功率谱修正。
再如图3所示本发明的装置主要包括:传统的电动振动台5,该电动振动台5包括功率放大器,激振器、试件、传感器,和非高斯随机信号生成模块1、D/A转换器2、低通滤波器3以及功率放大器4。
非高斯随机信号生成模块1是由下述模块实现的:
录入模块6:设定并录入所求非高斯分布随机过程的斜度、峭度和功率谱以及高斯分布的均值和方差;
服从高斯分布的时间序列模块7:由设定的高斯分布的均值和方差,再根据给定的功率谱密度条件,按常规的随机振动信号生成方法,得到一个服从高斯分布的时间序列;
非高斯概率分布函数模块8:根据产品实际工作时的测试数据或产品任务剖面的具体要求,得到所需非高斯分布随机信号的时间序列的概率密度函数;通过对所述概率密度函数进行积分得到概率分布函数;
单调递增函数生成模块9:利用非高斯概率分布函数与已知高斯分布随机信号概率分布函数之间关系导出单调递增非线性函数y=g(x);
转化模块10:利用y=g(x)与服从高斯分布的时间序列相乘将其转化成服从非高斯分布的时间序列,从而得到一个给定概率分布、斜度、峭度的零均值非高斯随机信号。
非高斯随机信号生成模块1生成非高斯信号之后,其输出的信号通过D/A转换器2变成连续的模拟量,通过低通滤波器3平滑后输给功率放大器4,最后经功率放大后驱动振动台5。
利用本发明的装置及本发明的方法生成的斜度为0,峭度为10的非高斯分布时间序列如图4所示,可见本发明的效果。
Claims (5)
1.一种非高斯随机振动激励信号的生成方法,包括以下步骤:
首先,在信号生成模块中设定所求非高斯分布随机过程的斜度、峭度和功率谱以及高斯分布的均值和方差;
第二步,由设定的高斯分布随机过程的均值和方差,再根据给定的功率谱密度条件,按常规的随机振动信号生成方法,得到一个服从高斯分布的时间序列;
其特征在于:
第三步,根据产品实际工作时的测试数据或产品任务剖面的具体要求,得到所需要非高斯分布随机信号的时间序列的概率密度函数,通过对所述概率密度函数进行积分得到概率分布函数;
第四步,利用非高斯概率分布函数与已知高斯分布随机信号概率分布函数之间关系导出单调递增非线性函数y=g(x);
第五步,利用y=g(x)与服从高斯分布的时间序列相乘将其转化成服从非高斯分布的时间序列,从而得到一个给定概率分布、斜度、峭度的零均值非高斯随机信号;
第六步,将所生成的非高斯随机信号通过D/A转换器变成连续的模拟量,再通过低通滤波器平滑后输给功率放大器,最后经功率放大后驱动电动振动台,实现非高斯分布的随机振动试验。
2.根据权利要求1所述的一种非高斯随机振动激励信号的生成方法,其特征在于单调递增函数y=g(x)可根据所求得的非高斯概率分布函数和给定的高斯分布的概率分布函数得到,
其中,X为服从高斯分布的随机变量;
Y为服从非高斯分布的随机变量;
fY(y)为随机变量Y的概率分布函数;
fX(x)为随机变量X的概率分布函数;
在得到dx/dy之后可以求得dy/dx,将其进行积分,则可得到y=g(x)。
y=g(x)为单调递增非线性的函数,即
3.一种非高斯随机振动激励信号生成方法的装置,主要包括传统的电动振动台,该电动振动台包括功率放大器,激振器、试件、传感器,其特征在于,还包括非高斯随机信号生成模块、D/A转换器、低通滤波器以及功率放大器。
4.根据权利要求3所述的一种非高斯随机振动激励信号生成装置,其特征在于,非高斯随机信号生成模块是由下述模块实现的:
录入模块:设定并录入所求非高斯分布随机过程的斜度、峭度和功率谱以及高斯分布的均值和方差;
服从高斯分布的时间序列模块:由设定高斯分布的均值和方差,根据给定的功率谱密度条件,按常规的随机振动信号生成方法,得到一个服从高斯分布的时间序列;
非高斯概率分布函数模块:根据产品实际工作时的测试数据或产品任务剖面的具体要求,得到所需非高斯分布随机信号的时间序列的概率密度函数;通过对所述概率密度函数进行积分得到概率分布函数;
单调递增函数生成模块:利用非高斯概率分布函数与已知高斯分布随机信号概率分布函数之间关系导出单调递增非线性函数y=g(x);
转化模块:利用y=g(x)与服从高斯分布的时间序列相乘将其转化成服从非高斯分布的时间序列,从而得到一个给定概率分布、斜度、峭度的零均值非高斯随机信号。
5.根据权利要求3所述的一种非高斯随机振动激励信号生成装置,其特征在于,非高斯随机信号生成模块生成非高斯信号之后,其输出的信号通过D/A转换器变成连续的模拟量,通过低通滤波器平滑后输给功率放大器,最后经功率放大后驱动振动台。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2007100636549A CN100460845C (zh) | 2007-02-07 | 2007-02-07 | 一种非高斯随机振动激励信号生成方法及其装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2007100636549A CN100460845C (zh) | 2007-02-07 | 2007-02-07 | 一种非高斯随机振动激励信号生成方法及其装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101038232A true CN101038232A (zh) | 2007-09-19 |
CN100460845C CN100460845C (zh) | 2009-02-11 |
Family
ID=38889280
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2007100636549A Expired - Fee Related CN100460845C (zh) | 2007-02-07 | 2007-02-07 | 一种非高斯随机振动激励信号生成方法及其装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100460845C (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106570280A (zh) * | 2016-11-08 | 2017-04-19 | 上海宇航系统工程研究所 | 基于函数单调性的卫星适配器可靠度区间推断方法 |
JP2018021781A (ja) * | 2016-08-02 | 2018-02-08 | 地方独立行政法人大阪産業技術研究所 | 非ガウス性振動制御装置 |
CN108490459A (zh) * | 2018-02-12 | 2018-09-04 | 千寻位置网络有限公司 | 精度与风险均衡应用于gnss位置服务的方法及系统 |
CN108562417A (zh) * | 2018-07-12 | 2018-09-21 | 北京航天希尔测试技术有限公司 | 一种激振力频段可调式气动激振器 |
CN110967184A (zh) * | 2019-12-03 | 2020-04-07 | 合肥工业大学 | 基于振动信号分布特征识别的变速箱故障检测方法和系统 |
CN111666063A (zh) * | 2020-04-22 | 2020-09-15 | 南京惟心光电系统有限公司 | 一种基于随机计算的增函数实现装置 |
CN112284518A (zh) * | 2020-09-30 | 2021-01-29 | 南京航空航天大学 | 非高斯多维随机振动波形生成方法 |
CN112784218A (zh) * | 2021-02-03 | 2021-05-11 | 中国工程物理研究院总体工程研究所 | 一种非高斯随机振动信号的峭度估计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4991107A (en) * | 1988-12-15 | 1991-02-05 | Schlumberger Technologies Limited | Vibration control system |
-
2007
- 2007-02-07 CN CNB2007100636549A patent/CN100460845C/zh not_active Expired - Fee Related
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2018021781A (ja) * | 2016-08-02 | 2018-02-08 | 地方独立行政法人大阪産業技術研究所 | 非ガウス性振動制御装置 |
CN106570280A (zh) * | 2016-11-08 | 2017-04-19 | 上海宇航系统工程研究所 | 基于函数单调性的卫星适配器可靠度区间推断方法 |
CN108490459A (zh) * | 2018-02-12 | 2018-09-04 | 千寻位置网络有限公司 | 精度与风险均衡应用于gnss位置服务的方法及系统 |
CN108562417A (zh) * | 2018-07-12 | 2018-09-21 | 北京航天希尔测试技术有限公司 | 一种激振力频段可调式气动激振器 |
CN108562417B (zh) * | 2018-07-12 | 2024-02-27 | 北京航天希尔测试技术有限公司 | 一种激振力频段可调式气动激振器 |
CN110967184A (zh) * | 2019-12-03 | 2020-04-07 | 合肥工业大学 | 基于振动信号分布特征识别的变速箱故障检测方法和系统 |
CN111666063A (zh) * | 2020-04-22 | 2020-09-15 | 南京惟心光电系统有限公司 | 一种基于随机计算的增函数实现装置 |
CN111666063B (zh) * | 2020-04-22 | 2023-09-26 | 南京惟心光电系统有限公司 | 一种基于随机计算的增函数实现装置 |
CN112284518A (zh) * | 2020-09-30 | 2021-01-29 | 南京航空航天大学 | 非高斯多维随机振动波形生成方法 |
CN112784218A (zh) * | 2021-02-03 | 2021-05-11 | 中国工程物理研究院总体工程研究所 | 一种非高斯随机振动信号的峭度估计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN100460845C (zh) | 2009-02-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101038232A (zh) | 一种非高斯随机振动激励信号生成方法及其装置 | |
CN110086403B (zh) | 一种马达驱动信号设置方法、电子设备及存储介质 | |
CN103346692B (zh) | 振动主动控制中压电作动器迟滞非线性频域补偿方法 | |
US8942930B2 (en) | System and method for generation and control of mechanical vibration | |
EP2206231B1 (en) | Signal generation apparatus and method | |
CN1714502A (zh) | 经脉冲宽度调制的噪声整形器 | |
CN111985495A (zh) | 模型部署方法、装置、系统及存储介质 | |
CN1788413A (zh) | 功率放大器预失真 | |
CN103067009A (zh) | 高精度数模转换器测试系统 | |
CN1602535A (zh) | 获得未来电压和/或电流变化的方法 | |
CN101039118A (zh) | 一种数字模拟转换系统及方法 | |
CN1298178C (zh) | 多载波自动增益校准装置和方法 | |
CN101030195A (zh) | 一种补偿示波器时基抖动的方法 | |
US8049650B2 (en) | Method for testing a high-speed digital to analog converter based on an undersampling technique | |
CN110826197A (zh) | 一种基于改进Cholesky分解闭合解的风速场模拟方法 | |
CN101127575A (zh) | 一种均匀分布随机数发生器及均匀分布随机数产生方法 | |
CN115060605A (zh) | 基于加速疲劳试验的振动疲劳寿命高精度预测方法 | |
DE112021004250T5 (de) | System und verfahren zum kompensieren von rückkopplungsverzögerungen in digitalen klasse-d-modulatoren | |
CN1536785A (zh) | 一种占空比可调高速光归零码产生方法和装置 | |
CN1650524A (zh) | 具备纠错功能的数模转换器 | |
CN111079250A (zh) | 一种电子产品疲劳寿命评估、评估模型建立方法及装置 | |
CN101043204A (zh) | 全数字d类功放中pwm结构的非线性误差校正方法 | |
CN102195650B (zh) | 模数转换系统和方法 | |
CN1272633A (zh) | 一种用于液晶显示器的驱动方案 | |
Wisland et al. | A non-feedback Δ-Σ modulator for digital-to-analog conversion using digital frequency modulation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20090211 Termination date: 20130207 |