CN103092565B - 一种真随机信号产生方法 - Google Patents
一种真随机信号产生方法 Download PDFInfo
- Publication number
- CN103092565B CN103092565B CN201110336579.5A CN201110336579A CN103092565B CN 103092565 B CN103092565 B CN 103092565B CN 201110336579 A CN201110336579 A CN 201110336579A CN 103092565 B CN103092565 B CN 103092565B
- Authority
- CN
- China
- Prior art keywords
- sequence
- group
- random signal
- true random
- random number
- 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.)
- Active
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.2)通过位置随机数取伪随机序列对应位置的数值及此位置后个位置的数值形成一组序列;将反转随机数与对应位置随机数的一组序列相乘,得到新的一组序列;
步骤(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.2)通过位置随机数取伪随机序列对应位置的数值及此位置后个位置的数值形成一组序列。
步骤(3.3)将对应同一位置随机数的一组序列进行加窗,加窗类型可选,如半正弦窗、汉明窗、海宁窗等。
步骤(3.4)将对应不同位置随机数的不同组序列按一定比例进行搭接,得到真随机信号。优选搭接比例为50%。
实施例2
本实施例与实施例1的区别在于:步骤(3)中采用以下操作通过时域随机化产生真随机信号:
步骤(3.1)随机生成一组位置随机数和一组反转随机数。
步骤(3.2)通过位置随机数取伪随机序列对应位置的数值及此位置后个位置的数值形成一组序列;将反转随机数与对应位置随机数的一组序列相乘,得到新的一组序列。
步骤(3.3)将对应同一位置随机数的新的一组序列进行加窗,加窗类型可选,如半正弦窗、汉明窗、海宁窗等。
步骤(3.4)将对应不同位置随机数的不同组新的序列按一定比例进行搭接,得到真随机信号。优选搭接比例为50%。
Claims (6)
1.一种真随机信号产生方法,其特征在于:包括以下步骤:
步骤(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为频率谱线数;
步骤(2.2)对互谱矩阵进行Cholesky分解,得到下三角矩阵L;
步骤(2.3)生成独立白噪声谱矩阵,其与下三角矩阵L相乘得到驱动谱矩阵;
步骤(2.4)对驱动谱矩阵进行IFFT变换,得到伪随机序列。
2.根据权利要求1所述的真随机信号产生方法,其特征在于:
频率谱线数nSL为100,或200,或400,或800,或1600,或3200。
3.根据权利要求2所述的真随机信号产生方法,其特征在于:
步骤(3)具体包括以下步骤:
步骤(3.1)随机生成一组位置随机数;
步骤(3.2)通过位置随机数取伪随机序列对应位置的数值及此位置后个位置的数值形成一组序列;
步骤(3.3)将对应同一位置随机数的一组序列进行加窗;
步骤(3.4)将对应不同位置随机数的不同组序列按一定比例进行搭接,得到真随机信号。
4.根据权利要求2所述的真随机信号产生方法,其特征在于:
步骤(3)具体包括以下步骤:
步骤(3.1)随机生成一组位置随机数和一组反转随机数;
步骤(3.2)通过位置随机数取伪随机序列对应位置的数值及此位置后个位置的数值形成一组序列;将反转随机数与对应位置随机数的一组序列相乘,得到新的一组序列;
步骤(3.3)将对应同一位置随机数的新的一组序列进行加窗;
步骤(3.4)将对应不同位置随机数的不同组新的序列按一定比例进行搭接,得到真随机信号。
5.根据权利要求3或4所述的真随机信号产生方法,其特征在于:
步骤(3.3)中加窗类型为半正弦窗,或汉明窗,或海宁窗。
6.根据权利要求3或4所述的真随机信号产生方法,其特征在于:
步骤(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 CN103092565A (zh) | 2013-05-08 |
CN103092565B true 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 (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 | 南京航空航天大学 | 用于多维随机振动试验的控制装置和控制方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7116781B2 (en) * | 1993-11-18 | 2006-10-03 | Digimarc Corporation | Counteracting geometric distortions in watermarking |
-
2011
- 2011-10-31 CN CN201110336579.5A patent/CN103092565B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 |
---|
多维随机振动试验中的互谱控制技术;祝济之;《航天器环境工程》;20101031;第27卷(第5期);全文 * |
多维随机振动试验条件制定方法研究;陈立伟等;《强度与环境》;20101231;第37卷(第6期);全文 * |
多输入多输出振动试验控制系统的理论、算法及实现;贺旭东;《中国博士学位论文全文库(工程科技II辑)》;20070630(第06期);第13-54页,图2.2 * |
Also Published As
Publication number | Publication date |
---|---|
CN103092565A (zh) | 2013-05-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2013192550A3 (en) | Anti-jagged 1/jagged 2 cross-reactive antibodies, activatable anti-jagged antibodies and methods of use thereof | |
EA201692376A1 (ru) | Терапевтические конъюгаты с сульфатированными дендримерами для внутриклеточного нацеливания | |
WO2012005982A3 (en) | Reporter for rna polymerase ii termination | |
WO2012139039A3 (en) | Bisphenol compounds and methods for their use | |
WO2010000127A9 (en) | Method to induce and expand therapeutic alloantigen- specific human regulatory t cells in large-scale | |
CA2819111A1 (en) | Adjusting a dose distribution setting for a technical device for tumour therapy | |
CN103092565B (zh) | 一种真随机信号产生方法 | |
MY155287A (en) | Methods of generating hybrid/chimeric cells, and uses thereof | |
JP2008528687A5 (zh) | ||
Filippi et al. | Structure and reactivity of supramolecular systems acting as highly selective" chemo-logic" devices in the gas phase | |
Margan et al. | CONFORMATIONAL ANALYSIS FOR PROSTAGLANDIN E2. I. | |
Calza et al. | Drugs as new contaminants: photocatalytic degradation process as a tool to study environmental fate of β-blockers | |
Wörman et al. | Using a spectral approach to compare dynamic and static head driven hyporheic exchange | |
Yu | Remaking Xiamen: overseas Chinese and regional transformation in architecture and urbanism in the early 20thcentury | |
Falcinelli et al. | Molecular beam studies of the combustion reaction of oxygen atoms with propene: Primary products, branching ratios, and role of intersystem crossing | |
Emanov et al. | Seismic microzoning by the method of standing waves | |
Douglass et al. | Connections Between a High-Resolution Model and Aquarius Salinity Measurements | |
Asai et al. | Spatial distribution of toroidal flow velocity in a field-reversed configuration palsma | |
Gleyzer | Search for the dark matter signature in the lepton jet final state at the center of mass energy= 7 TeV | |
UA35329U (ru) | Способ интегрированной защиты посевов риса от сорняков | |
Zadonina et al. | Influence of model parameters on synthesized high-frequency strong-motion waveforms | |
Fitzgerald et al. | The Statistical State Dynamics of Stacked Jet Formation from Stratified Turbulence | |
Crepp | Spectra and Dynamical Masses of Brown Dwarfs and Extrasolar Planets | |
Nemeth et al. | The shape of the magnetodisk of Saturn-as revealed by in situ measurements | |
Suvorova et al. | LPM mode under long-lasting radial IMF in the subsolar region |
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 |