CN108536946A - 基于Landweber迭代法的运行工况传递路径分析方法 - Google Patents
基于Landweber迭代法的运行工况传递路径分析方法 Download PDFInfo
- Publication number
- CN108536946A CN108536946A CN201810283791.1A CN201810283791A CN108536946A CN 108536946 A CN108536946 A CN 108536946A CN 201810283791 A CN201810283791 A CN 201810283791A CN 108536946 A CN108536946 A CN 108536946A
- Authority
- CN
- China
- Prior art keywords
- operating condition
- test
- reference point
- transfer path
- analysis method
- 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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/10—Noise analysis or noise optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了一种基于Landweber迭代法的运行工况传递路径分析方法,基于试验运行工况设计和试验运行工况数据测量,利用观测到的参考点响应信号和目标点响应信号反求传递率函数矩阵。其包括如下步骤,步骤1,通过布置传感器采集试验工况下参考点、目标点响应信号,将观测信号分别进行快速傅里叶变换;步骤2,选择适用于不适定Landweber迭代法的迭代格式进行对传递率函数矩阵的估计;步骤3,将得到的传递率函数矩阵反乘实际工况测得的参考点响应,计算各传递路径贡献量和总贡献量。本方法在求解传递率函数矩阵过程中对于原始数据的较小扰动,不会导致所求得的解产生过大波动,拥有很好的抵抗噪声能力,求得解的拟合程度较高,反映实际情况能力优良。
Description
技术领域
本发明涉及机械振动信号处理和减振降噪领域,特别涉及基于Landweber迭代法的运行工况传递路径分析方法。
背景技术
水下兵器的声隐身性能声隐身性能是衡量其安全性和作战能力的重要指标,因此,壳体结构装备振动噪声的有效监测与控制对于提高装备性能具有重要工程意义。振动或噪声传递路径的识别是问题的关键,目前常采用传递路径分析(TPA)或者运行工况传递路径分析(OTPA)的方法来识别振动或噪声传递路径,而传统的传递路径分析方法由于其复杂的频响函数测试以及载荷识别过程导致在工程实际中难以快速有效的识别传递路径,运行工况传递路径分析方法是近年来出现的一种传递路径快速分析方法,该方法利用试验工况数据识别传递利率函数矩阵,并将其用于实际工况数据,从而得到传递路径贡献量结果,该方法简单快捷,被广泛运用于工程实际中。
现有OTPA方法均存在传递率函数矩阵拟合程度低、稳定性不足等缺陷,即解反映实际情况能力、抵抗噪声能力无法满足实际工程需求。因此,对于OTPA,需要采用更加有效的数学手段以确保矩阵良态,在高效的前提下提高精度。
发明内容
本发明的目的在于提供基于Landweber迭代法的运行工况传递路径分析方法,以解决上述问题。
为实现上述问题,本发明采用以下技术方案:
基于Landweber迭代法的运行工况传递路径分析方法,包括如下步骤:
步骤1,选取所分析机械系统激励点的临近位置作为参考点,选取待分析位置为目标点;
步骤2,根据所分析机械系统设计若干种工况试验,使机械系统在每种试验工况下运行,采集每种试验工况下参考点和目标点响应信号,得到参考点响应信号矩阵和目标点响应信号矩阵,并对其进行快速傅里叶变换得到相应的频响函数矩阵X和Y;
步骤3,设置迭代初值X0=αX*,计算迭代系数矩阵Xk=Xk-1+X0(I-XXk-1),在传递率函数矩阵Tk=XkY满足Morozov误差原则前持续迭代系数矩阵Xk,直至迭代得到的传递率函数矩阵Tk满足Morozov误差原则,停止迭代;
步骤4,根据所分析机械系统进行实际运行工况实验,得到一组实际运行工况数据,根据该组实际运行工况数据和迭代得到的传递率函数矩阵计算得到估计的目标点总贡献量和各路径贡献量
步骤5,针对各条传递路径对各路径贡献量进行排序后计算出各传递路径的贡献量占比,完成运行工况传递路径分析。
进一步的,步骤2、步骤4中均采用振动加速度传感器测量响应信号。
进一步的,步骤2中试验工况数量大于等于参考点数量,步骤4中实际运行工况数量大于等于参考点数量。
进一步的,步骤2中X(ω)=[x1(ω),x2(ω),...,xr(ω)]为全部试验工况下所有参考点响应信号,r表示试验工况数量,表示第i种试验工况下所有参考点响应信号,表示试验工况i下参考点j的振动响应信号,m表示参考点数量。
进一步的,步骤2中Y(ω)=[y1(ω),y2(ω),...,yr(ω)]为全部试验工况下所有目标点响应信号,r表示试验工况数量,表示第i种试验工况下所有目标点响应信号,表示试验工况i下目标点g的振动响应信号,n表示目标点数量。
进一步的,步骤3中迭代初值X0=αX*,α为松弛因子,且0<α<1/||X||2,X*表示参考点频率响应函数X的共轭转置。
进一步的,步骤3中计算迭代系数矩阵Xk=Xk-1+X0(I-XXk-1),I为单位矩阵,k表示迭代次数。
进一步的,步骤3中Morozov误差原则为:当迭代得到的传递率函数矩阵满足||X(ω)Tk(ω)-Y(ω)||≤τδ时停止迭代,常数τ=2,误差水平δ=10-3。
进一步的,步骤4中表示OTPA识别的参考点j对目标点i的传递路径贡献量。
与现有技术相比,本发明有以下技术效果:
本发明根据所分析机械系统设计试验工况,同时采集试验工况下参考点响应信号和目标点响应信号,针对传递率函数矩阵求解的不适定反问题特点,利用Landweber迭代法对传递率函数矩阵的求解进行迭代控制,设置Morozov误差原则为迭代终止准则,极大地提高了传递率函数矩阵求解的稳定性的同时有效地提高了其拟合程度,从而提高其抵抗噪声干扰能力和反映实际情况能力。然后采集一组实际运行工况下待分析机械系统的参考点响应信号和目标点响应信号,将迭代得到的传递率函数矩阵与实际工况下测得的参考点响应信号相乘,得到传递路径贡献量识别结果,对不同路径贡献量进行排序,得到各个传递路径的贡献量占比,完成运行工况传递路径分析。本发明克服了以往运行工况传递路径分析方法中传递率函数矩阵求解的拟合程度低和稳定性不足带来的识别误差问题,提高了传递路径贡献量识别精度。
附图说明
图1为本发明优选实例中所述的球状辐射声源仿真系统;
图2为本发明优选实例中所述的声源A源贡献量;
图3为本发明优选实例中所述的声源B源贡献量;
图4为本发明优选实例中所述的目标点总贡献量;
图5为本发明优选实例中所述的各贡献量识别精度柱状图;
图6为本发明优选实例中所述的声源A对目标点的源贡献量相对误差;
图7为本发明优选实例中所述的声源B对目标点的源贡献量相对误差;
图8为本发明优选实例中所述的特征频率处目标点总贡献量相对误差;
图9为本发明优选实例中所述的特征频率处目标点总贡献量相对误差;
图10为本发明优选实例中所述的求解稳定性比较柱状图;
具体实施方式
下面结合具体的实施例对本发明做进一步的详细说明,所述是对本发明的解释而不是限定。
基于Landweber迭代法的运行工况传递路径分析方法,包括如下步骤:
步骤1,选取所分析机械系统激励点的临近位置作为参考点,选取待分析位置为目标点;
步骤2,根据所分析机械系统设计若干种工况试验,使机械系统在每种试验工况下运行,采集每种试验工况下参考点和目标点响应信号,得到参考点响应信号矩阵和目标点响应信号矩阵,并对其进行快速傅里叶变换得到相应的频响函数矩阵X和Y;
步骤3,设置迭代初值X0=αX*,计算迭代系数矩阵Xk=Xk-1+X0(I-XXk-1),在传递率函数矩阵Tk=XkY满足Morozov误差原则前持续迭代系数矩阵Xk,直至迭代得到的传递率函数矩阵Tk满足Morozov误差原则,停止迭代;
步骤4,根据所分析机械系统进行实际运行工况实验,得到一组实际运行工况数据,根据该组实际运行工况数据和迭代得到的传递率函数矩阵计算得到估计的目标点总贡献量和各路径贡献量
步骤5,针对各条传递路径对各路径贡献量进行排序后计算出各传递路径的贡献量占比,完成运行工况传递路径分析。
步骤2、步骤4中均采用振动加速度传感器测量响应信号。
步骤2中试验工况数量大于等于参考点数量,步骤4中实际运行工况数量大于等于参考点数量。
步骤2中X(ω)=[x1(ω),x2(ω),...,xr(ω)]为全部试验工况下所有参考点响应信号,r表示试验工况数量,表示第i种试验工况下所有参考点响应信号,表示试验工况i下参考点j的振动响应信号,m表示参考点数量。
步骤2中Y(ω)=[y1(ω),y2(ω),...,yr(ω)]为全部试验工况下所有目标点响应信号,r表示试验工况数量,表示第i种试验工况下所有目标点响应信号,表示试验工况i下目标点g的振动响应信号,n表示目标点数量。
步骤3中迭代初值X0=αX*,α为松弛因子,且0<α<1/||X||2,X*表示参考点频率响应函数X的共轭转置。
步骤3中计算迭代系数矩阵Xk=Xk-1+X0(I-XXk-1),I为单位矩阵,k表示迭代次数。
步骤3中Morozov误差原则为:当迭代得到的传递率函数矩阵满足||X(ω)Tk(ω)-Y(ω)||≤τδ时停止迭代,常数τ=2,误差水平δ=10-3。
步骤4中表示OTPA识别的参考点j对目标点i的传递路径贡献量。
本发明根据所分析机械系统设计试验工况,同时采集试验工况下参考点响应信号和目标点响应信号,针对传递率函数矩阵求解的不适定反问题特点,利用Landweber迭代法对传递率函数矩阵的求解进行迭代控制,设置Morozov误差原则为迭代终止准则,极大地提高了传递率函数矩阵求解的稳定性的同时有效地提高了其拟合程度,从而提高其抵抗噪声干扰能力和反映实际情况能力。本发明在设计试验工况过程中,选取所分析机械系统临近激励点位置作为参考点,待分析位置作为目标点,参考点个数与激励源个数相同,并且试验工况数量大于等于参考点数量。将迭代得到的传递率函数矩阵和实际工况下测得的参考点响应信号相乘,得到传递路径贡献量识别结果,对不同路径贡献量进行排序,得到各个传递路径的贡献量占比,完成运行工况传递路径分析。本发明克服了以往运行工况传递路径分析方法中传递率函数矩阵求解的拟合程度低和稳定性不足带来的识别误差问题,提高了传递路径贡献量识别精度。
优选实例,为说明基于Landweber迭代法的OTPA相对于以往OTPA方法的优越性,对比基于截断奇异值分解的传统OTPA和基于Tikhonov正则化的OTPA方法,设计球状辐射声源仿真系统,参见图1所示。
设声源1与声源2的半径均为a=b=0.34m;
声源表面振动速度均为v1=v2=60m/s,且同相位;
声源1和声源2到参考点距离均为r1ref=r2ref=0.85m;
声源1和声源2到目标点距离均为r1=r2=10.25;
声源1到声源2的距离r12=13.6m。
参考点和目标点处的声压是声源1和声源2贡献量的向量和,可以表达为:
式中:noise(ω)——随机噪声。
设置分析频率为0~100Hz,频率分辨率为Δf=1Hz;
声源1和声源2分别在25Hz和50Hz输出信号;
随机噪声在0~100Hz输出信号。
截断奇异值分解(TSVD)参数设置:
截断奇异值分解方法的关键是选择合适的截断系数k。截断系数k取值过大,会导致传递率函数矩阵严重病态(不稳定);截断系数过小,会导致截断误差过大,传递率函数矩阵的拟合程度不足。
传统OTPA根据奇异值贡献率:
来确定截断系数k,在本次仿真过程中取奇异值贡献率为3%。
Tikhonov正则化参数设置:
Tikhonov正则化参数λ的取值非常关键,λ>0且λ∈(σn,σ1)。
最常用的正则化参数选择方法有广义交叉检验方法和L曲线方法,本次仿真选用L曲线法选取正则化参数,即求解[log(||Axλ-b||),log(||xλ||)]曲率极大值确定最优正则化参数λ。
Landweber迭代法终止准则设置:
如前所述,利用Landweber迭代法求解反问题的关键是要寻找一个合适的终止准则使迭代近似解在半收敛附近停止迭代。常用的迭代准则是Morozov误差原则:
为使迭代近似解在半收敛附近停止迭代,设置常数τ为τ=2,误差水平δ取δ=10-3。
仿真步骤设定:
1)基于10组试验工况数据,利用上述三种方法分别求解传递率函数矩阵;
2)再仿真1组实际工况下的数据作为校验标准,包括单独开启声源A时目标点的响应、单独开启声源B时目标点的响应以及同时开启两个声源时目标点的响应;
3)利用步骤1)中求解的传递率函数矩阵和实际工况下的参考点响应分别求出声源A对目标点的源贡献量、声源B对目标点的源贡献量以及目标点的总贡献量;
4)对比步骤2)和步骤3)的结果差异,并分别计算出相对误差;
5)按整数倍增大噪声倍数至10倍,并分别对不同噪声倍数重复上述过程。
图2-图4所示为所选取的随机噪声增大至6倍时三种方法的贡献量对比图,其中practical代表实际工况下测得的目标点响应,图2所示为声源A源贡献量对比图,图3所示为声源B源贡献量对比图,图4所示为目标点总贡献量。
图5所示为随机噪声增大至6倍时三种方法的各贡献量识别精度的精确比较柱状图,表明了本方法的贡献量识别精度明显好于另外两种对比方法。
图6-图9所示为逐步增大噪声倍数至10倍时源贡献量相对误差和目标点总贡献量相对误差变化趋势折线图,图6所示为声源A对目标点的源贡献量相对误差对比折线图,图7所示为声源B对目标点的源贡献量相对误差对比折线图,图8所示为目标点总贡献量在特征频率25Hz处相对误差对比折线图,图9所示为目标点总贡献量在特征频率50Hz处相对误差对比折线图,表明了本方法的对源贡献量的识别精度和总贡献量的识别精度在整体上要明显优于另外两种对比方法,显示识别的传递率函数矩阵具有更高的拟合程度。
为了评价本方法识别传递率函数矩阵的稳定性,分别计算三种方法的各组相对误差相邻数值间的方差并求和,即:其中δi为随机噪声增大至i倍时计算得到的贡献量相对误差。
图10所示为三种方法的贡献量求解稳定性的对比图,表明本方法计算得到的各贡献量识别相对误差相邻数值间的方差和远小于另外两种对比方法,即本方法的求解稳定性具有明显优势,显示了本方法传递率函数矩阵的抵抗噪声干扰能力明显优于另外两种对比方法,在实际工程应用中具有更好的应用价值。
Claims (9)
1.基于Landweber迭代法的运行工况传递路径分析方法,其特征在于,包括如下步骤:
步骤1,选取所分析机械系统激励点的临近位置作为参考点,选取待分析位置为目标点;
步骤2,根据所分析机械系统设计若干种工况试验,使机械系统在每种试验工况下运行,采集每种试验工况下参考点和目标点响应信号,得到参考点响应信号矩阵和目标点响应信号矩阵,并对其进行快速傅里叶变换得到相应的频响函数矩阵X和Y;
步骤3,设置迭代初值X0=αX*,计算迭代系数矩阵Xk=Xk-1+X0(I-XXk-1),在传递率函数矩阵Tk=XkY满足Morozov误差原则前持续迭代系数矩阵Xk,直至迭代得到的传递率函数矩阵Tk满足Morozov误差原则,停止迭代;
步骤4,根据所分析机械系统进行实际运行工况实验,得到一组实际运行工况数据,根据该组实际运行工况数据和迭代得到的传递率函数矩阵计算得到估计的目标点总贡献量和各路径贡献量
步骤5,针对各条传递路径对各路径贡献量进行排序后计算出各传递路径的贡献量占比,完成运行工况传递路径分析。
2.根据权利要求1所述的基于Landweber迭代法的运行工况传递路径分析方法,其特征在于,步骤2、步骤4中均采用振动加速度传感器测量响应信号。
3.根据权利要求1所述的基于Landweber迭代法的运行工况传递路径分析方法,其特征在于,步骤2中试验工况数量大于等于参考点数量,步骤4中实际运行工况数量大于等于参考点数量。
4.根据权利要求1所述的基于Landweber迭代法的运行工况传递路径分析方法,其特征在于,步骤2中X(ω)=[x1(ω),x2(ω),...,xr(ω)]为全部试验工况下所有参考点响应信号,r表示试验工况数量,表示第i种试验工况下所有参考点响应信号,表示试验工况i下参考点j的振动响应信号,m表示参考点数量。
5.根据权利要求1所述的基于Landweber迭代法的运行工况传递路径分析方法,其特征在于,步骤2中Y(ω)=[y1(ω),y2(ω),...,yr(ω)]为全部试验工况下所有目标点响应信号,r表示试验工况数量,表示第i种试验工况下所有目标点响应信号,表示试验工况i下目标点g的振动响应信号,n表示目标点数量。
6.根据权利要求1所述的基于Landweber迭代法的运行工况传递路径分析方法,其特征在于,步骤3中迭代初值X0=αX*,α为松弛因子,且0<α<1/||X||2,X*表示参考点频率响应函数X的共轭转置。
7.根据权利要求1所述的基于Landweber迭代法的运行工况传递路径分析方法,其特征在于,步骤3中计算迭代系数矩阵Xk=Xk-1+X0(I-XXk-1),I为单位矩阵,k表示迭代次数。
8.根据权利要求1所述的基于Landweber迭代法的运行工况传递路径分析方法,其特征在于,步骤3中Morozov误差原则为:当迭代得到的传递率函数矩阵满足||X(ω)Tk(ω)-Y(ω)||≤τδ时停止迭代,常数τ=2,误差水平δ=10-3。
9.根据权利要求1所述的基于Landweber迭代法的运行工况传递路径分析方法,其特征在于,步骤4中表示OTPA识别的参考点j对目标点i的传递路径贡献量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810283791.1A CN108536946A (zh) | 2018-04-02 | 2018-04-02 | 基于Landweber迭代法的运行工况传递路径分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810283791.1A CN108536946A (zh) | 2018-04-02 | 2018-04-02 | 基于Landweber迭代法的运行工况传递路径分析方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108536946A true CN108536946A (zh) | 2018-09-14 |
Family
ID=63482756
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810283791.1A Pending CN108536946A (zh) | 2018-04-02 | 2018-04-02 | 基于Landweber迭代法的运行工况传递路径分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108536946A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109492316A (zh) * | 2018-11-19 | 2019-03-19 | 西安交通大学 | 一种基于ar参数模型的otpa串扰削减方法 |
CN112202695A (zh) * | 2020-08-05 | 2021-01-08 | 重庆大学 | 一种基于Landweber迭代算法的欠采样数字预失真方法及系统 |
CN112560180A (zh) * | 2020-12-16 | 2021-03-26 | 广东工业大学 | 多点安装结构的传递路径分析方法 |
CN113761470A (zh) * | 2021-09-01 | 2021-12-07 | 国家电投集团河南电力有限公司开封发电分公司 | 基于有限参考点的结构部件全息振型测试方法 |
CN113776648A (zh) * | 2021-07-26 | 2021-12-10 | 重庆大学 | 一种基于正则化总体最小二乘法的工况传递路径分析方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106354929A (zh) * | 2016-08-29 | 2017-01-25 | 西安交通大学 | 基于刚度变化原理的承载结构载荷传递路径可视化方法 |
CN106528955A (zh) * | 2016-10-18 | 2017-03-22 | 西安交通大学 | 基于阻尼奇异值分解的运行工况传递路径分析方法 |
-
2018
- 2018-04-02 CN CN201810283791.1A patent/CN108536946A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106354929A (zh) * | 2016-08-29 | 2017-01-25 | 西安交通大学 | 基于刚度变化原理的承载结构载荷传递路径可视化方法 |
CN106528955A (zh) * | 2016-10-18 | 2017-03-22 | 西安交通大学 | 基于阻尼奇异值分解的运行工况传递路径分析方法 |
Non-Patent Citations (3)
Title |
---|
WEI CHENG等: "Tikhonov regularization-based operational transfer path analysis", 《MECHANICAL SYSTEMS AND SIGNAL PROCESSING》 * |
张军: "求解不适定问题的快速 Landweber 迭代法", 《数学杂志》 * |
成玮等: "Tikhonov正则化在运行工况传递路径分析的应用", 《振动、测试与诊断》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109492316A (zh) * | 2018-11-19 | 2019-03-19 | 西安交通大学 | 一种基于ar参数模型的otpa串扰削减方法 |
CN112202695A (zh) * | 2020-08-05 | 2021-01-08 | 重庆大学 | 一种基于Landweber迭代算法的欠采样数字预失真方法及系统 |
CN112560180A (zh) * | 2020-12-16 | 2021-03-26 | 广东工业大学 | 多点安装结构的传递路径分析方法 |
CN112560180B (zh) * | 2020-12-16 | 2022-04-22 | 广东工业大学 | 多点安装结构的传递路径分析方法 |
CN113776648A (zh) * | 2021-07-26 | 2021-12-10 | 重庆大学 | 一种基于正则化总体最小二乘法的工况传递路径分析方法 |
CN113776648B (zh) * | 2021-07-26 | 2023-09-08 | 重庆大学 | 一种基于正则化总体最小二乘法的工况传递路径分析方法 |
CN113761470A (zh) * | 2021-09-01 | 2021-12-07 | 国家电投集团河南电力有限公司开封发电分公司 | 基于有限参考点的结构部件全息振型测试方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108536946A (zh) | 基于Landweber迭代法的运行工况传递路径分析方法 | |
CN107247259B (zh) | 基于神经网络的k分布海杂波形状参数估计方法 | |
Tan et al. | Reducing data dimensionality through optimizing neural network inputs | |
CN104166804B (zh) | 一种基于时频域单源点稀疏成分分析的工作模态辨识方法 | |
CN110275148A (zh) | 一种海杂波幅度分布参数估计方法及系统 | |
CN106468770A (zh) | K分布杂波加噪声下的近最优雷达目标检测方法 | |
CN112163373A (zh) | 基于贝叶斯机器学习的雷达系统性能指标动态评估方法 | |
CN107885928A (zh) | 考虑测量误差的步进应力加速性能退化可靠性分析方法 | |
CN110245443A (zh) | 一种基于映射函数的多类型输出仿真模型可信度综合评估方法 | |
CN112147571B (zh) | 基于正则正交匹配追踪和蝙蝠算法的声源方位角估计方法 | |
CN110428185B (zh) | 基于伪量测模型的电-热互联综合能源系统抗差状态估计方法 | |
US11768180B1 (en) | Method for ultrasonic guided wave quantitative imaging in form of variable array | |
CN108846200B (zh) | 一种基于迭代法的准静态桥梁影响线识别方法 | |
CN106199537A (zh) | 基于逆高斯纹理海杂波幅度分布参数的分位点估计方法 | |
CN112214880B (zh) | 一种基于虚实融合的多阶段维修性评估方法 | |
CN107704704A (zh) | 一种基于贝叶斯理论的继电器类单机贮存可靠性评估方法 | |
CN107832789A (zh) | 基于平均影响值数据变换的特征加权k近邻故障诊断方法 | |
CN110889207B (zh) | 一种基于深度学习的体系组合模型可信度智能评估方法 | |
CN108761384B (zh) | 一种抗差的传感器网络目标定位方法 | |
CN104731762B (zh) | 基于循环移位的立方相位信号参数估计方法 | |
CN110297223A (zh) | 基于bp神经网络的脉冲激光回波数目检测的方法 | |
CN110263762A (zh) | 一种基于输出的半潜式海洋平台能量传递路径分析方法 | |
CN106353743A (zh) | 匹配于等效形状参数的近最优雷达目标检测方法 | |
CN110568074B (zh) | 基于非接触多点测振与Hilbert变换的风力机叶片裂纹定位方法 | |
CN109031328B (zh) | 含时序增益功能激光测距机的探测概率曲线簇测量方法 |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180914 |
|
RJ01 | Rejection of invention patent application after publication |