CN103092565A - 一种真随机信号产生方法 - Google Patents
一种真随机信号产生方法 Download PDFInfo
- Publication number
- CN103092565A CN103092565A CN2011103365795A CN201110336579A CN103092565A CN 103092565 A CN103092565 A CN 103092565A CN 2011103365795 A CN2011103365795 A CN 2011103365795A CN 201110336579 A CN201110336579 A CN 201110336579A CN 103092565 A CN103092565 A CN 103092565A
- Authority
- CN
- China
- Prior art keywords
- centerdot
- true random
- signal generating
- group
- 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
Landscapes
- Complex Calculations (AREA)
Abstract
本发明属于结构强度与环境可靠性试验技术领域,具体涉及一种真随机信号产生方法。本发明的方法包括以下步骤:参数设置;获取伪随机序列;通过时域随机化产生真随机信号。本发明的方法解决了现有真随机信号产生方法不能实现不同通道间真随机信号相干控制的技术问题。本发明的方法能够实现不同通道间真随机信号相干控制,产生的真随机信号符合正态分布,在包括结构强度与环境可靠性试验的众多技术领域,具有广阔应用前景。
Description
技术领域
本发明属于结构强度与环境可靠性试验技术领域,具体涉及一种真随机信号产生方法。
背景技术
真随机信号是指功率谱符合高斯分布的随机电压信号,广泛应用于各领域控制系统及检定系统中。现有技术中,上述控制系统和检定系统产生多通道的真随机信号,但各通道的真随机信号之间只能实现全独立或全相关,不能使通道间保持一定的相干和相位关系,即不能实现不同通道间真随机信号相干控制,而不能实现不同通道间相干控制的真随机信号,在结构强度与环境可靠性试验领域,应用范围十分有限。
发明内容
本发明需要解决的技术问题是现有的真随机信号产生方法不能实现不同通道间真随机信号的相干控制。
本发明的技术方案如下所述:
一种真随机信号产生方法,包括以下步骤:
步骤(1)参数设置;
步骤(2)获取伪随机序列;
步骤(3)通过时域随机化产生真随机信号。
步骤(1)中通过控制计算机进行参数设置,需要设置的参数包括:频率范围、参考谱形、通道数目、通道间相干系数、通道间相干相位、频率谱线数及量级。
步骤(2)具体包括以下步骤:
步骤(2.1)通过下式将参考谱形、通道间相干系数及通道间相干相位转化成互谱矩阵:
其中:
CPM表示互谱矩阵;
PSD表示参考谱形;
Δf表示频率谱线间隔,Δf=频率范围/频率谱线数;
Coe表示通道间相干系数;
Phase表示通道间相干相位;
l=1,2,……,nGroupCount,nGroupCount为通道数目;
k=1,2,……,nGroupCount,nGroupCount为通道数目;
j=1,2,3,……,nSL,nSL为频率谱线数,nSL优选为100,或200,或400,或800,或1600,或3200;
步骤(2.2)对互谱矩阵进行Cholesky分解,得到下三角矩阵L。
步骤(2.3)生成独立白噪声谱矩阵,其与下三角矩阵L相乘得到驱动谱矩阵。
步骤(2.4)对驱动谱矩阵进行IFFT变换,得到伪随机序列。
步骤(3)具体包括以下步骤:
步骤(3.1)随机生成一组位置随机数;
步骤(3.2)通过位置随机数取伪随机序列对应位置的数值及此位置后个位置的数值形成一组序列;
步骤(3.3)将对应同一位置随机数的一组序列进行加窗;
步骤(3.4)将对应不同位置随机数的不同组序列按一定比例进行搭接,得到真随机信号。
作为优选方案,步骤(3)还可以具体包括以下步骤:
步骤(3.1)随机生成一组位置随机数和一组反转随机数;
步骤(3.3)将对应同一位置随机数的新的一组序列进行加窗;
步骤(3.4)将对应不同位置随机数的不同组新的序列按一定比例进行搭接,得到真随机信号。
上述两种步骤(3)的操作中,步骤(3.3)中加窗类型可以为半正弦窗,或汉明窗,或海宁窗;搭接比例可以为50%
本发明的有益效果为:
本发明的真随机信号产生方法能够实现不同通道间真随机信号相干控制,产生的真随机信号符合正态分布,在包括结构强度与环境可靠性试验的众多技术领域,具有广阔应用前景。
具体实施方式
下面结合实施例对本发明的一种真随机信号产生方法进行详细说明。
实施例1
本实施例的真随机信号产生方法具体包括以下步骤:
步骤(1)参数设置
通过控制计算机进行参数设置,需要设置的参数包括:频率范围、参考谱形、通道数目、通道间相干系数、通道间相干相位、频率谱线数及量级。
步骤(2)获取伪随机序列
步骤(2.1)通过下式将参考谱形、通道间相干系数及通道间相干相位转化成互谱矩阵:
其中:
CPM表示互谱矩阵;
PSD表示参考谱形;
Δf表示频率谱线间隔,Δf=频率范围/频率谱线数;
Coe表示通道间相干系数;
Phase表示通道间相干相位;
l=1,2,……,nGroupCount,nGroupCount为通道数目;
k=1,2,……,nGroupCount,nGroupCount为通道数目;
j=1,2,3,……,nSL,nSL为频率谱线数,可以为100、200、400、800、1600、3200。
步骤(2.2)对互谱矩阵进行Cholesky分解,得到下三角矩阵L。
步骤(2.3)生成独立白噪声谱矩阵,其与下三角矩阵L相乘得到驱动谱矩阵。
步骤(2.4)对驱动谱矩阵进行IFFT变换,得到伪随机序列。
步骤(3)通过时域随机化产生真随机信号
步骤(3.1)随机生成一组位置随机数,随机生成位置随机数的方法为本领域技术人员公知常识,本实施例中通过TMS320C6713型DSP芯片实现。
步骤(3.3)将对应同一位置随机数的一组序列进行加窗,加窗类型可选,如半正弦窗、汉明窗、海宁窗等。
步骤(3.4)将对应不同位置随机数的不同组序列按一定比例进行搭接,得到真随机信号。优选搭接比例为50%。
实施例2
本实施例与实施例1的区别在于:步骤(3)中采用以下操作通过时域随机化产生真随机信号:
步骤(3.1)随机生成一组位置随机数和一组反转随机数。
步骤(3.3)将对应同一位置随机数的新的一组序列进行加窗,加窗类型可选,如半正弦窗、汉明窗、海宁窗等。
步骤(3.4)将对应不同位置随机数的不同组新的序列按一定比例进行搭接,得到真随机信号。优选搭接比例为50%。
Claims (8)
1.一种真随机信号产生方法,其特征在于:包括以下步骤:
步骤(1)参数设置;
步骤(2)获取伪随机序列;
步骤(3)通过时域随机化产生真随机信号。
2.根据权利要求1所述的真随机信号产生方法,其特征在于:
步骤(1)中通过控制计算机进行参数设置,需要设置的参数包括:频率范围、参考谱形、通道数目、通道间相干系数、通道间相干相位、频率谱线数及量级。
3.根据权利要求1所述的真随机信号产生方法,其特征在于:
步骤(2)具体包括以下步骤:
步骤(2.1)通过下式将参考谱形、通道间相干系数及通道间相干相位转化成互谱矩阵:
其中:
CPM表示互谱矩阵;
PSD表示参考谱形;
Δf表示频率谱线间隔,Δf=频率范围/频率谱线数;
Coe表示通道间相干系数;
Phase表示通道间相干相位;
l=1,2,……,nGroupCount,nGroupCount为通道数目;
k=1,2,……,nGroupCount,nGroupCount为通道数目;
j=1,2,3,……,nSL,nSL为频率谱线数;
步骤(2.2)对互谱矩阵进行Cholesky分解,得到下三角矩阵L。
步骤(2.3)生成独立白噪声谱矩阵,其与下三角矩阵L相乘得到驱动谱矩阵。
步骤(2.4)对驱动谱矩阵进行IFFT变换,得到伪随机序列。
4.根据权利要求3所述的真随机信号产生方法,其特征在于:
频率谱线数nSL为100,或200,或400,或800,或1600,或3200。
7.根据权利要求5或6所述的真随机信号产生方法,其特征在于:
步骤(3.3)中加窗类型为半正弦窗,或汉明窗,或海宁窗。
8.根据权利要求5或6所述的真随机信号产生方法,其特征在于:
步骤(3.4)中搭接比例为50%。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110336579.5A CN103092565B (zh) | 2011-10-31 | 2011-10-31 | 一种真随机信号产生方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110336579.5A CN103092565B (zh) | 2011-10-31 | 2011-10-31 | 一种真随机信号产生方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103092565A true CN103092565A (zh) | 2013-05-08 |
CN103092565B CN103092565B (zh) | 2016-08-03 |
Family
ID=48205185
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110336579.5A Active CN103092565B (zh) | 2011-10-31 | 2011-10-31 | 一种真随机信号产生方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103092565B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020118831A1 (en) * | 1993-11-18 | 2002-08-29 | Rhoads Geoffrey B. | Counteracting geometric distortions in watermarking |
EP1341004A1 (fr) * | 2002-03-01 | 2003-09-03 | Institut Francais Du Petrole | Méthode et dispositif de prospection sismique par émission simultanée de signaux sismiques à base de séquences pseudo aléatoires |
CN101697084A (zh) * | 2009-10-19 | 2010-04-21 | 大连海事大学 | 一种基于rls滤波器的电液伺服系统随机振动控制方法 |
CN102122151A (zh) * | 2010-12-10 | 2011-07-13 | 南京航空航天大学 | 用于多维随机振动试验的控制装置和控制方法 |
-
2011
- 2011-10-31 CN CN201110336579.5A patent/CN103092565B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020118831A1 (en) * | 1993-11-18 | 2002-08-29 | Rhoads Geoffrey B. | Counteracting geometric distortions in watermarking |
EP1341004A1 (fr) * | 2002-03-01 | 2003-09-03 | Institut Francais Du Petrole | Méthode et dispositif de prospection sismique par émission simultanée de signaux sismiques à base de séquences pseudo aléatoires |
CN101697084A (zh) * | 2009-10-19 | 2010-04-21 | 大连海事大学 | 一种基于rls滤波器的电液伺服系统随机振动控制方法 |
CN102122151A (zh) * | 2010-12-10 | 2011-07-13 | 南京航空航天大学 | 用于多维随机振动试验的控制装置和控制方法 |
Non-Patent Citations (3)
Title |
---|
祝济之: "多维随机振动试验中的互谱控制技术", 《航天器环境工程》 * |
贺旭东: "多输入多输出振动试验控制系统的理论、算法及实现", 《中国博士学位论文全文库(工程科技II辑)》 * |
陈立伟等: "多维随机振动试验条件制定方法研究", 《强度与环境》 * |
Also Published As
Publication number | Publication date |
---|---|
CN103092565B (zh) | 2016-08-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
ATE519134T1 (de) | Erzeugung eines absorptionsparametermodells | |
WO2013081854A3 (en) | A graphene nanomesh based charge sensor | |
CN103092224B (zh) | 多维随机振动控制方法 | |
CN101832849B (zh) | 基于三参量控制的振动台软启动控制方法 | |
CN103076194B (zh) | 实时混合模拟试验效果的频域评价方法 | |
CN102392785B (zh) | 一种风力机最大功率跟踪曲线的线性化方法 | |
CN104156504B (zh) | 一种发电机励磁系统参数可辨识性判定方法 | |
CN103092565A (zh) | 一种真随机信号产生方法 | |
CN107465465A (zh) | 一种高斯信道仿真方法及其仿真系统 | |
CN104077479A (zh) | 一种基于守恒迎风格式获取参量阵声场空间分布的方法 | |
CN107292034A (zh) | 基于Hyperworks的路面激励车内振动噪声仿真分析方法和系统 | |
CN104792447B (zh) | 一种大型船舶隔振装置动态耦合多载荷识别方法 | |
CN104215964A (zh) | 一种多列等差频率原波相互作用形成参量阵的声场获取方法 | |
CN103065044A (zh) | 在分形海面背景下畸形波的模拟方法 | |
CN102564787B (zh) | 基于空运行激励的数控机床模态比例因子获取方法 | |
鲁俊 et al. | Numerical investigation of vertical turbulent jets in different types of waves | |
CN105372493A (zh) | 基于三条dft复数谱线的信号幅值和相位测量方法 | |
CN104343628A (zh) | 一种含有死区补偿器的风电机组变桨距控制方法 | |
CN104360322B (zh) | 基于四阶非对称乘积型核函数的qfm信号参数估计方法 | |
Liu et al. | Evaluation of rotary damping sprinkler hydraulic performance in sprinkler irrigation system | |
Yu et al. | Research on Random Vibration Control Algorithm Based-on EV Model | |
Kwon | CO 2 Emission and Productivity of Fossil-fueled Power Plants: A Luenberger Indicator Approach | |
WO2012175243A3 (de) | Bestimmung der läuferposition einer synchronmaschine | |
Chi et al. | Eddy covariance measurements of the forest floor CO2 exchange in two contrasting forest stands in boreal Sweden | |
WO2013170133A3 (en) | Computer-generated investment index |
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 |