CN107341554A - 一种带有次谐频的非线性振动***输出响应预测方法 - Google Patents
一种带有次谐频的非线性振动***输出响应预测方法 Download PDFInfo
- Publication number
- CN107341554A CN107341554A CN201710429592.2A CN201710429592A CN107341554A CN 107341554 A CN107341554 A CN 107341554A CN 201710429592 A CN201710429592 A CN 201710429592A CN 107341554 A CN107341554 A CN 107341554A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msup
- harmonics
- munderover
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/02—Reservations, e.g. for tickets, services or events
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- General Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Evolutionary Computation (AREA)
- Development Economics (AREA)
- Computer Hardware Design (AREA)
- Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Geometry (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Strategic Management (AREA)
- General Business, Economics & Management (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明是一种带有次谐频的非线性振动***输出响应预测方法,其特点是,通过对带有次谐频的非线性振动***预补偿一个超谐频***建立了假想源与次谐频振动***输出的等效沃尔泰拉级数(Volterra Series,VS)模型;基于此模型,通过变幅值扫频实验辨识出等效VS模型的修正广义频响函数(Modified Generalize Frequency Response Functions,MGFRFs)最后基于等效VS模型利用已辨识出MGFRFs实现带有次谐频非线性振动***响应预测,具有科学合理,计算简单,适用性强,预测精度高等优点,利用本发明所述的方法可以快速实现带有次谐频的非线性振动***输出响应预测。
Description
技术领域
本发明涉及非线性***响应预测,尤其是一种带有次谐频的非线性振动***输出响应预测方法。
背景技术
响应预测是工程中最基本最常见的问题,其主要任务在于验算结构、产品在工作时动力响应,是否满足安全要求和其他要求。目前线性***的响应预测研究已趋于成熟,而非线性***响应预测仍未完善。随着日益增长的科技需要,实际工程中非线性问题难以避免,例如在对于舰船实施精确降噪过程中,舰船结构、隔振装置、主动控制执行器的非线性特性已是不可回避的问题。非线性振动***按多项式结构划分通常分为纯输入非线性***,纯输出非线性***以及输入输出非线性***。本发明所涉及的非线性***则是其中纯输出非线性***,典型的纯输出非线性***,如Duffing振子,Vander Pol振子等。
由于非线性微分方程尚无精确的理论解,因此非线性振动响应预测通常是在使用特定函数空间描述的非参数***模型基础之上进行的,非参数***模型有:沃尔泰拉级数(Volterra Series,VS)模型、维纳级数(Wiener Series,WS)模型等。
目前VS模型已广泛应用于非线性振动***分析及控制,对于纯输入非线性***通常可用有限阶VS模型表示,然而带次谐频的纯输出非线性***,有限的单输入单输出(Single Input Single Output,SISO)的VS模型无法表示此类***,为克服此问题,Boaghe将SISO VS模型推广到多输入多输出(Multi Input Single Output,MISO)VS模型,实现了带有次谐频的非线性***时域建模,但是MISO模型处理数据复杂而且其Volterra核函数只有在某一确定输入幅值下才有意义,无法在变幅值情况下进行输出响应预测,这给实际非线性振动***分析及控制带来困难。为解决上述非线性振动***分析及控制所遇到的问题,本发明提供一种新的带次谐频出输出非线性***输出响应的预测方法。本发明可以快速实现带有次谐频的非线性振动***输出响应预测,计算简单且预测精度高。
发明内容
本发明的目的是提供一种计算简单,适用性强,预测精度高的带有次谐频的非线性振动***输出响应预测方法,利用本发明所述的方法可以快速实现带有次谐频的非线性振动***输出响应预测。
本发明的目的是通过以下技术方案来实现的:一种带有次谐频的非线性振动***输出响应预测方法,其特征是,它包括以下内容:
1)建立等效VS模型
对于输出不含次谐频时的情形,***可直接利用VS模型便可建模,然而当次谐频产生时,直接利用原始输入数据,无法用有限的VS模型表示带次谐频的非线性振动***,于是考虑通过对带次谐频的非线性振动***预补偿一个超谐频***,建立假想源r=A1/pcos(θ/p),θ=ω(t-t0)与真实输出之间的等效的VS模型,
假想源r(t)与真实输入u(t)之间的关系表示为:
这里A为输入幅值;p值由次谐频决定,如次谐频为1/2倍基频则p=2,如次谐频为1/3倍基频则p=3,以此类推,若输出不含次谐频则p=1,即假想源等于真实输入,
2)辨识等效VS模型
假想源与真实输入等效VS数学表达为:
其中hn(τ1,…,τn)为Volterra时域核函数,其多维Fourier变换:
Hn(ω1,…,ωn)为广义频响函数(Generalized Frequency Response Functions,GFRFs);
定义修正广义频响函数(Modified GFRFs,MGFRFs)
含有次谐频的非线性振动***输出可以表示为:
其中:
这里int(·)表示取整数,
通过变幅值扫频实验,利用实验获得的输入输出数据根据LS算法便可辨识出各MGFRFs,显然计算简单易行,
3)响应预测
利用已辨识出的MGFRFs,根据式(5)~式(6)即能够实现在任意输入幅值下带有次谐频的非线性振动***输出响应预测。
本发明的一种带有次谐频的非线性振动***输出响应预测方法,具有科学合理,计算简单,适用性强,预测精度高等优点,利用本发明所述的方法可以快速实现带有次谐频的非线性振动***输出响应预测。
附图说明
图1为非线性等效方法实现框图;
图2为基于等效VS模型对带有次谐频的非线性振动***输出响应预测的准备阶段流程图;
图3为基于等效VS模型对带有次谐频的非线性振动***输出响应预测的预测阶段流程图;
图4为输入频率2rad/s的响应谱图(Response Spectrum Map,RSM);
图5为带1/2次谐频响应预测值与真实值曲线对比图。
具体实施方式
本发明的一种带有次谐频的非线性振动***输出响应预测方法,包括以下内容:
1)建立等效VS模型
对于输出不含次谐频时的情形,***可直接利用VS模型便可建模,然而当次谐频产生时,直接利用原始输入数据,无法用有限的VS模型表示带次谐频的非线性振动***,于是考虑通过对带次谐频的非线性振动***预补偿一个超谐频***,建立假想源r=A1/pcos(θ/p),θ=ω(t-t0)与真实输出之间的等效的VS模型,如图1所示,图1中:G为预补偿***,P为目标带次谐的非线性振动***,H为等效非线性***;
假想源r(t)与真实输入u(t)之间的关系表示为:
其中,A为输入幅值;p值由次谐频决定,如次谐频为1/2倍基频则p=2,如次谐频为1/3倍基频则p=3,以此类推,若输出不含次谐频则p=1,即假想源等于真实输入;
2)辨识等效VS模型
假想源与真实输入等效VS数学表达为:
其中hn(τ1,…,τn)为Volterra时域核函数,其多维Fourier变换:
Hn(ω1,…,ωn)为广义频响函数(Generalized Frequency Response Functions,GFRFs);
定义修正广义频响函数(Modified GFRFs,MGFRFs)
于是含有次谐频的非线性振动***输出可以表示为:
其中:
这里int(·)表示取整数,
通过变幅值扫频实验,利用实验获得的输入输出数据根据LS算法便可辨识出各MGFRFs,显然计算简单易行,
参照图2,基于等效VS模型对带有次谐频的非线性振动***输出响应预测的准备阶段流程为:
步骤一:对带次谐频的非线性振动***进行变幅值扫频实验,并对输入输出进行数据采集;
步骤二:将步骤一提取出的输入和输出频率数据运用FFT算法提取频率成份并进行对比,确定输出含次谐频的输入幅值范围和不含次谐频输入幅值范围以及适用的频率范围;
步骤三:根据步骤一中数据处理得到输入输出频率成份,合理确定p值;
步骤四:利用式(2)~(6)建立带有次谐频纯输出非线性***的等效VS模型;
步骤五:利用实验数据辨识利用LS算法分别辨识含次谐频输入范围的等效VS模型的MGFRFs,及不含次谐频输入幅值范围的等效VS模型的MGFRFs;
步骤六:利用已有数据对已建立的等效VS模型进行验证,判断是否满足精度要求,若满足则转步骤七,若否则返回步骤四,通过增加其中等效VS模型非线性截断阶数,提高模型精度以达到目标要求;
步骤七:输出“***建模完成,可进行输出响应预测”,
3)响应预测
利用已辨识出的MGFRFs,根据式(5)~式(6)即能够实现在任意输入幅值下带有次谐频的非线性振动***输出响应预测。
参照图3,基于等效VS模型对带有次谐频的非线性振动***输出响应预测的预测阶段流程为:
步骤一:接收需预测输入数据;
步骤二:提取输入数据频率成分及幅值;
步骤三:由准备阶段的输入频率分析范围判断所需预测的输入频率是否在分析范围以内,若是则转步骤四,若否则输出“不符合预测要求”;
步骤四:依据准备阶段数据,判断输入幅值所在区间是否存在次谐频,若是则转步骤五,若否则转步骤六;
步骤五:调用***含次谐频等效VS模型的MGFRFs,根据式(5)~(6)预测输出响应;
步骤六:调用不含次谐频等效VS模型的MGFRFs,根据式(5)~(6)预测输出响应;
步骤七:输出预测响应结果。
图4为输入频率2rad/s的响应频谱图(response spectrum map,RSM),可以从图4中看出当输入幅值在8.1至12存在1/2次谐频,利用本发明的方法对***进行响应预测,结果如图5所示,预测值与真实值曲线重合度高,可以看出本发明预测结果的精确性,验证了本发明方法的有效性。
以上具体实施方式仅用以说明本发明的技术方案,而非对其限制,所属领域的普通技术人员应该理解,参照上述实施例所作的任何形式的修改、等同变化均在本发明权利要求保护范围之内。
Claims (1)
1.一种带有次谐频的非线性振动***输出响应预测方法,其特征是,它包括以下内容:
1)建立等效沃尔塔拉级数(Volterra Series,VS)模型
对于输出不含次谐频时的情形,***可直接利用VS模型便可建模,然而当次谐频产生时,直接利用原始输入数据,无法用有限的VS模型表示带次谐频的非线性振动***,于是考虑通过对带次谐频的非线性振动***预补偿一个超谐频***,建立假想源r=A1/pcos(θ/p),θ=ω(t-t0)与真实输出之间的等效的VS模型,
假想源r(t)与真实输入u(t)之间的关系表示为:
<mrow>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>int</mi>
<mrow>
<mo>(</mo>
<mi>p</mi>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
</munderover>
<mo>{</mo>
<mo>&lsqb;</mo>
<msubsup>
<mi>C</mi>
<mrow>
<mi>p</mi>
<mo>-</mo>
<mi>k</mi>
</mrow>
<mi>k</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>C</mi>
<mrow>
<mi>p</mi>
<mo>-</mo>
<mn>1</mn>
<mo>-</mo>
<mi>k</mi>
</mrow>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>&rsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mi>k</mi>
</msup>
<msup>
<mn>2</mn>
<mrow>
<mi>p</mi>
<mo>-</mo>
<mn>1</mn>
<mo>-</mo>
<mn>2</mn>
<mi>k</mi>
</mrow>
</msup>
<msup>
<mi>A</mi>
<mrow>
<mn>2</mn>
<mi>k</mi>
<mo>/</mo>
<mi>p</mi>
</mrow>
</msup>
<msup>
<mi>r</mi>
<mrow>
<mi>p</mi>
<mo>-</mo>
<mn>2</mn>
<mi>k</mi>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>}</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
这里A为输入幅值;p值由次谐频决定,如次谐频为1/2倍基频则p=2,如次谐频为1/3倍基频则p=3,以此类推,若输出不含次谐频则p=1,即假想源等于真实输入,
2)辨识等效VS模型
假想源与真实输入等效VS数学表达为:
<mrow>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<munderover>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<mi>&infin;</mi>
</mrow>
<mrow>
<mo>+</mo>
<mi>&infin;</mi>
</mrow>
</munderover>
<mo>...</mo>
<munderover>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<mi>&infin;</mi>
</mrow>
<mrow>
<mo>+</mo>
<mi>&infin;</mi>
</mrow>
</munderover>
<msub>
<mi>h</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&tau;</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>&tau;</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<munderover>
<mo>&Pi;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<msub>
<mi>&tau;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>d&tau;</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
其中hn(τ1,···,τn)为Volterra时域核函数,其多维Fourier变换:
<mrow>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&omega;</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>&omega;</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<mi>&infin;</mi>
</mrow>
<mrow>
<mo>+</mo>
<mi>&infin;</mi>
</mrow>
</msubsup>
<mo>...</mo>
<msubsup>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<mi>&infin;</mi>
</mrow>
<mrow>
<mo>+</mo>
<mi>&infin;</mi>
</mrow>
</msubsup>
<msub>
<mi>h</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&tau;</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>&tau;</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mi>j</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&omega;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&tau;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<mo>...</mo>
<mo>+</mo>
<msub>
<mi>&omega;</mi>
<mi>n</mi>
</msub>
<msub>
<mi>&tau;</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<munderover>
<mo>&Pi;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>d&tau;</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
Hn(ω1,···,ωn)为广义频响函数(Generalized Frequency Response Functions,GFRFs);
定义修正广义频响函数(Modified GFRFs,MGFRFs)
含有次谐频的非线性振动***可以表示为:
其中:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<mi>y</mi>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>s</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>int</mi>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
</munderover>
<mfrac>
<msup>
<mi>A</mi>
<mrow>
<mn>2</mn>
<mi>s</mi>
</mrow>
</msup>
<msup>
<mn>2</mn>
<mrow>
<mi>s</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mfrac>
<msub>
<mi>R</mi>
<mrow>
<mn>2</mn>
<mi>s</mi>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mi>l</mi>
</msub>
<mo>=</mo>
<mfrac>
<msup>
<mi>A</mi>
<mi>l</mi>
</msup>
<msup>
<mn>2</mn>
<mrow>
<mi>l</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>l</mi>
<mo>,</mo>
<mi>l</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<msup>
<mi>A</mi>
<mrow>
<mi>l</mi>
<mo>+</mo>
<mn>2</mn>
</mrow>
</msup>
<msup>
<mn>2</mn>
<mrow>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msup>
</mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>l</mi>
<mo>+</mo>
<mn>2</mn>
<mo>,</mo>
<mi>l</mi>
</mrow>
</msub>
<mn>...</mn>
<mo>+</mo>
<mfrac>
<msup>
<mi>A</mi>
<mi>K</mi>
</msup>
<msup>
<mn>2</mn>
<mrow>
<mi>K</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>K</mi>
<mo>,</mo>
<mi>l</mi>
</mrow>
</msub>
<mo>,</mo>
<mi>K</mi>
<mo>=</mo>
<mn>2</mn>
<mi>int</mi>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>-</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mo>&rsqb;</mo>
<mo>+</mo>
<mi>l</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>b</mi>
<mi>l</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<msup>
<mi>A</mi>
<mi>l</mi>
</msup>
<msup>
<mn>2</mn>
<mrow>
<mi>l</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mfrac>
<msub>
<mi>I</mi>
<mrow>
<mi>l</mi>
<mo>,</mo>
<mi>l</mi>
</mrow>
</msub>
<mo>-</mo>
<mfrac>
<msup>
<mi>A</mi>
<mrow>
<mi>l</mi>
<mo>+</mo>
<mn>2</mn>
</mrow>
</msup>
<msup>
<mn>2</mn>
<mrow>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msup>
</mfrac>
<msub>
<mi>I</mi>
<mrow>
<mi>l</mi>
<mo>+</mo>
<mn>2</mn>
<mo>,</mo>
<mi>l</mi>
</mrow>
</msub>
<mn>...</mn>
<mo>-</mo>
<mfrac>
<msup>
<mi>A</mi>
<mi>K</mi>
</msup>
<msup>
<mn>2</mn>
<mrow>
<mi>K</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mfrac>
<msub>
<mi>I</mi>
<mrow>
<mi>K</mi>
<mo>,</mo>
<mi>l</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
1
这里int(·)表示取整数,
通过变幅值扫频实验,利用实验获得的输入输出数据根据LS算法便可辨识出各MGFRFs,显然计算简单易行,
3)响应预测
利用已辨识出的MGFRFs,根据式(5)~式(6)即能够实现在任意输入幅值下带有次谐频的非线性振动***输出响应预测。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710429592.2A CN107341554A (zh) | 2017-06-09 | 2017-06-09 | 一种带有次谐频的非线性振动***输出响应预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710429592.2A CN107341554A (zh) | 2017-06-09 | 2017-06-09 | 一种带有次谐频的非线性振动***输出响应预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107341554A true CN107341554A (zh) | 2017-11-10 |
Family
ID=60220550
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710429592.2A Pending CN107341554A (zh) | 2017-06-09 | 2017-06-09 | 一种带有次谐频的非线性振动***输出响应预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107341554A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110350515A (zh) * | 2019-06-20 | 2019-10-18 | 西安交通大学 | 一种适用于故障电流分析的柔性直流电网建模方法 |
-
2017
- 2017-06-09 CN CN201710429592.2A patent/CN107341554A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110350515A (zh) * | 2019-06-20 | 2019-10-18 | 西安交通大学 | 一种适用于故障电流分析的柔性直流电网建模方法 |
CN110350515B (zh) * | 2019-06-20 | 2020-11-10 | 西安交通大学 | 一种适用于故障电流分析的柔性直流电网建模方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Shao et al. | New delay-dependent stability criteria for neural networks with two additive time-varying delay components | |
Xu et al. | Stochastic stability for nonlinear systems driven by Lévy noise | |
El-Zahar | Applications of adaptive multi step differential transform method to singular perturbation problems arising in science and engineering | |
Singh et al. | On Painlevé analysis, symmetry group and conservation laws of Date–Jimbo–Kashiwara–Miwa equation | |
CN105760347A (zh) | 一种基于数据/极值联合对称延拓的hht端点效应抑制方法 | |
CN102636693A (zh) | 一种结合fft与非线性最小二乘的谐波分析算法 | |
Cazzulani et al. | Linear and non-linear systems identification for adaptive control in mechanical applications vibration suppression | |
CN107341554A (zh) | 一种带有次谐频的非线性振动***输出响应预测方法 | |
Voliansky et al. | Chaotic time-variant dynamical system | |
Dai et al. | Norm-based finite-time convergent recurrent neural network for dynamic linear inequality | |
Pintoux et al. | A direct solution to the Fokker–Planck equation for exponential Brownian functionals | |
de Morais Goulart et al. | Efficient kernel computation for Volterra filter structure evaluation | |
Giehl et al. | Edge-computing enhanced privacy protection for industrial ecosystems in the context of SMEs | |
Chen et al. | Eigenfilter design of linear-phase FIR digital filters using neural minor component analysis | |
Tavazoei | Toward searching possible oscillatory region in order space for nonlinear fractional-order systems | |
Hoffmann et al. | Linear parameter-varying control of complex mechanical systems | |
Adam et al. | Golden section, Fibonacci sequence and the time invariant Kalman and Lainiotis filters | |
Liu et al. | Time-invariant Lyapunov matrix equation solving by neural networks with simulink models | |
Li et al. | A general alternating-direction implicit Newton method for solving complex continuous-time algebraic Riccati matrix equation | |
Zhou et al. | An Improved Digital Filtering Method Based on Weighted Average of Input and Output Data | |
Hussain et al. | Harmonie detection at initialization with Kaiman filter | |
Zhan et al. | A new robust-control-oriented system identification method | |
Kusumasari et al. | Simulation of hyperbolic mean curvature flow with an obstacle in the closed curve | |
Lu et al. | On stable inversion for linear systems | |
Yao et al. | A multi-resolution pid controller based on the wavelet transform |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20171110 |
|
WD01 | Invention patent application deemed withdrawn after publication |