CN119740195A 基于多源信息的航空变流器故障特征提取及诊断方法

(12)发明专利申请(10)申请公布号CN119740195A(21)申请号202510255475.3(71)申请人南京航空航天大学(54)发明名称(57)摘要和频谱获取及特征提取，实现多源信息的融23以下方程组得到：{2=12;r(xi,x;)+μ=γ(x;,P)i=1,…,,这里条插值分别得到上包络线emax(t)和下包络线emin(t),计算上下包络线的均值到满足条件的IMF₁(t);然后，从原始序列中减去IMF₁(t),得到剩余序列4对于粗粒化后的序列Zj,计算其样本熵；给定嵌入维数n和容限r,定义向量类似地，计算n+1维向量的Cn+1(r);样本熵对于时间序列x(t),在计算自相关时，考虑不同延迟T下的子序列x₁={x(1),x(2),…,x(N-T)}和X₂={x(T+1子序列X₁={x(1),x(2),…,x(N-T)}和y₂={y(t+1),y(T+2),…,y(N)}过DTW算法找到最优路径P,计算DTW距离DTW(x₁,y₂);改进后的5⑨n,k,T,w(t)=An,k,其中An,k是归一化常数，6对于每次迭代m计算残差rm-1(t)与所有可能的时频原子φn,K,T,u(t)的内积x(t)≈ZM=1(rm-1,4nm,km,Tn,0m>·φnm,km,Tm,wm(时频原子的频率支撑范围或相关能量分布来确定；然后，带宽B计算为分别得到上包络线emax(t)和下包络线emin(t),计算上下包络线的均值到得到满足条件的IMF₁(t);然后，从原始信号中减去IMF₁(t),得到剩余序列希尔伯特变换求包络的过程为：对每个IMFi(t)进行希尔伯特变换，得到Zi(t)=IMF₂(t)+jy₁(t),其幅值即为IMF₂的瞬时幅78Coiflet小波Coif3进行处理：第j层分解的近似系数Cj,k计算公式为：9 [0017]其中，软阈值部分γ[sgn(dj,k)(|d;,k|-λ)]能平滑噪声，硬阈值部分[0020]异常判定的过程为：为每个传感器邻域设定一个正常信息熵阈值范围IMF₂(t),IMF₃(t),…,IMFn(t)和残余项rn(t),即类似地，计算n+1维向量的Cn+1(r);样本熵[0026]对于时间序列x(t),在计算自相关时，考虑不同延迟T下的子序列路径P的约束条件为i₁=1,j₁=1,iz=N-子序列X₁={x(1),x(2),…,x(N-t)}和y₂={y(T+1),y(t+2),…,y(N)},通过DTW算法找到最优路径P,计算DTW距离DTW(x₁,y₂);改进后的互相关函数9i,j,k=∑m,n,1Wm,n,ifi+m-1,j+n-1,k+l-1+b,其中Wm,n,I是卷积核元素，b是偏置，希尔伯特变换H[uk(t)],将其与Uk(t)构成解析信号uk(t)+jH[uk(t)],通过指数项 变分问题为,约束条件为,其中An,k是归一化常数，[0040]更新残差Ym(t)=rm-1(t)-(rm-1,4nm,kmuTm,m>·4nmkm,Tm,①m(t);x(t)≈ZM=1(rm-1,4nm,km,Tn,0m>·φnm,km,Tm,wm(通过时频原子的频率支撑范围或相关能量分布来确定；然后，带宽B计算为[0046]希尔伯特变换求包络的过程为：对每个Zi(t)=IMFi(t)+jyi(t),其幅值即为IMF的瞬时幅Yi,j,k=∑m,n,IWm,n,lXi+m-1,j+n-1,k+l-1+b,附图说明 [0073]第j层分解的近似系数Cj,k计算公式为：,其中，,hj.k(n)是与所选小波基对应的低通滤波器系数， [0074]第j层分解的细节系数dj,k计算公式为：,其中[0076]S1.3、考虑到航空机载变流器运行状态的动态变化，结合变流器温度变化率为混合参数(0<γ<1),处理后到满足条件的IMF₁(t)。然后，从原始序列中减去IMF₁(t),得到剩余序列类似地，计算n+1维向量的Cn+1(r)。样本熵决这一问题。对于时间序列x(t),在计算自相关时，考虑不同延迟T下的子序列 子序列X₁={x(1),x(2),…,x(N-T)}和y₂={y(T+1),y(T+2),…,y(N)},通过DTW算法找到最优路径P,计算DTW距离DTW(x₁,y₂)。改进后的互相关函数,其中An,k是归一化常数，通过时频原子的频率支撑范