CN114714911B - 一种基于多域稀疏表达的电动客车高压负载系统故障诊断方法 - Google Patents
一种基于多域稀疏表达的电动客车高压负载系统故障诊断方法 Download PDFInfo
- Publication number
- CN114714911B CN114714911B CN202210374371.0A CN202210374371A CN114714911B CN 114714911 B CN114714911 B CN 114714911B CN 202210374371 A CN202210374371 A CN 202210374371A CN 114714911 B CN114714911 B CN 114714911B
- Authority
- CN
- China
- Prior art keywords
- domain
- vector
- sparse
- geometric
- time
- 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
- 238000003745 diagnosis Methods 0.000 title claims abstract description 47
- 238000000034 method Methods 0.000 title claims abstract description 35
- 239000013598 vector Substances 0.000 claims abstract description 122
- 239000011159 matrix material Substances 0.000 claims abstract description 54
- 230000004044 response Effects 0.000 claims abstract description 49
- 238000012360 testing method Methods 0.000 claims abstract description 46
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 15
- SAZUGELZHZOXHB-UHFFFAOYSA-N acecarbromal Chemical compound CCC(Br)(CC)C(=O)NC(=O)NC(C)=O SAZUGELZHZOXHB-UHFFFAOYSA-N 0.000 claims abstract description 12
- 230000008569 process Effects 0.000 claims abstract description 8
- 238000012545 processing Methods 0.000 claims abstract description 3
- 238000001228 spectrum Methods 0.000 claims description 21
- 230000005284 excitation Effects 0.000 claims description 15
- 238000005070 sampling Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 238000010438 heat treatment Methods 0.000 claims description 3
- 238000012804 iterative process Methods 0.000 claims description 3
- 239000003550 marker Substances 0.000 claims description 3
- 238000004804 winding Methods 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 2
- 230000007704 transition Effects 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 description 4
- 239000003990 capacitor Substances 0.000 description 3
- 238000012423 maintenance Methods 0.000 description 3
- 238000003909 pattern recognition Methods 0.000 description 3
- 230000015556 catabolic process Effects 0.000 description 2
- 238000006731 degradation reaction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000005485 electric heating Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000002372 labelling Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000002405 diagnostic procedure Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60L—PROPULSION OF ELECTRICALLY-PROPELLED VEHICLES; SUPPLYING ELECTRIC POWER FOR AUXILIARY EQUIPMENT OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRODYNAMIC BRAKE SYSTEMS FOR VEHICLES IN GENERAL; MAGNETIC SUSPENSION OR LEVITATION FOR VEHICLES; MONITORING OPERATING VARIABLES OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRIC SAFETY DEVICES FOR ELECTRICALLY-PROPELLED VEHICLES
- B60L3/00—Electric devices on electrically-propelled vehicles for safety purposes; Monitoring operating variables, e.g. speed, deceleration or energy consumption
- B60L3/0023—Detecting, eliminating, remedying or compensating for drive train abnormalities, e.g. failures within the drive train
- B60L3/0061—Detecting, eliminating, remedying or compensating for drive train abnormalities, e.g. failures within the drive train relating to electrical machines
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60L—PROPULSION OF ELECTRICALLY-PROPELLED VEHICLES; SUPPLYING ELECTRIC POWER FOR AUXILIARY EQUIPMENT OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRODYNAMIC BRAKE SYSTEMS FOR VEHICLES IN GENERAL; MAGNETIC SUSPENSION OR LEVITATION FOR VEHICLES; MONITORING OPERATING VARIABLES OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRIC SAFETY DEVICES FOR ELECTRICALLY-PROPELLED VEHICLES
- B60L3/00—Electric devices on electrically-propelled vehicles for safety purposes; Monitoring operating variables, e.g. speed, deceleration or energy consumption
- B60L3/0023—Detecting, eliminating, remedying or compensating for drive train abnormalities, e.g. failures within the drive train
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60L—PROPULSION OF ELECTRICALLY-PROPELLED VEHICLES; SUPPLYING ELECTRIC POWER FOR AUXILIARY EQUIPMENT OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRODYNAMIC BRAKE SYSTEMS FOR VEHICLES IN GENERAL; MAGNETIC SUSPENSION OR LEVITATION FOR VEHICLES; MONITORING OPERATING VARIABLES OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRIC SAFETY DEVICES FOR ELECTRICALLY-PROPELLED VEHICLES
- B60L3/00—Electric devices on electrically-propelled vehicles for safety purposes; Monitoring operating variables, e.g. speed, deceleration or energy consumption
- B60L3/0023—Detecting, eliminating, remedying or compensating for drive train abnormalities, e.g. failures within the drive train
- B60L3/003—Detecting, eliminating, remedying or compensating for drive train abnormalities, e.g. failures within the drive train relating to inverters
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60L—PROPULSION OF ELECTRICALLY-PROPELLED VEHICLES; SUPPLYING ELECTRIC POWER FOR AUXILIARY EQUIPMENT OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRODYNAMIC BRAKE SYSTEMS FOR VEHICLES IN GENERAL; MAGNETIC SUSPENSION OR LEVITATION FOR VEHICLES; MONITORING OPERATING VARIABLES OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRIC SAFETY DEVICES FOR ELECTRICALLY-PROPELLED VEHICLES
- B60L2200/00—Type of vehicles
- B60L2200/18—Buses
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
- Y04S10/52—Outage or fault management, e.g. fault detection or location
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Sustainable Development (AREA)
- Sustainable Energy (AREA)
- Power Engineering (AREA)
- Transportation (AREA)
- Mechanical Engineering (AREA)
- Tests Of Circuit Breakers, Generators, And Electric Motors (AREA)
- Testing Of Short-Circuits, Discontinuities, Leakage, Or Incorrect Line Connections (AREA)
Abstract
本发明涉及一种基于多域稀疏表达的电动客车高压负载系统故障诊断方法,属于电动客车高压负载系统级故障诊断领域。该方法分为诊断前准备和诊断过程,其中诊断前先对电动客车进行多次测试,获取电流响应信号并分别构造时域、频域及几何域下的字典矩阵;诊断过程中,先收集测试样本的电流响应信号,再用标记样本字典矩阵分别在时域、频域及几何域下稀疏表达电流响应信号,采用批量匹配追踪算法得到时域、频域及几何域下的稀疏向量,融合三种域下的稀疏向量,再对其处理并识别测试样本故障模式。本发明用于电动客车的高压负载系统故障诊断,及时发现和消除电动客车高压系统的潜在风险。
Description
技术领域
本发明属于电动客车高压负载系统级故障诊断领域,涉及一种基于多域稀疏表达的电动客车高压负载系统故障诊断方法。
背景技术
为了缓解能源危机和减少环境污染,电动客车近年来取得了快速发展,并在许多城市批量推广。然而,由于电动客车的动力电池电压平台最高达750V,电量多达500KWh,甚至更多,无论是电压平台还是存储电量都要远高于电动乘用车,并且作为载客公共交通工具,电动客车载客量也远超过私家电动乘用车,因此电动客车的高压系统安全非常重要,一旦出现高压系统故障,将严重威胁公共安全和人身安全。因此,要及时发现和消除电动客车高压系统的潜在风险,当然电动客车高压负载故障诊断系统就成为电动客车的重要组成部分。
电动客车的发展虽然已有了很大的进步,但在关键技术开发,尤其是高压负载系统故障诊断方面仍需要取得进一步突破。目前高压系统故障诊断的相关工作一直集中在部件级故障诊断,这些方法普遍用于电动乘用车故障诊断,如驱动电机故障诊断方法包括机理分析、信号分析和智能故障诊断,采用基于残差的方法对功率变换器进行故障诊断,采用故障树分析技术对充电装置进行故障诊断等。然而,对于电动客车高压负载系统而言,系统级故障诊断和定位比部件级故障诊断更有实用价值,这取决于其应用情况和维修服务系统。首先,作为公共交通系统的关键组成部分,电动客车的故障必须尽快得到解决,以便快速运营,以免给城市交通带来压力,一般来说,为了实现这一目标,采用故障组件快速替换,而不是采用故障组件现场修复,因此电动客车系统级故障诊断与定位尤为重要。
针对电动客车系统级故障诊断与定位的研究较少,这为城市客车的智能维修提供了障碍。为解决这一问题,提高电动客车维修效率,提出了一种基于激励响应分析和多域稀疏表达的系统级故障诊断与定位方案。
发明内容
有鉴于此,本发明的目的在于提供一种基于多域稀疏表达的电动客车高压负载系统故障诊断方法,以快速准确地诊断出电动客车高压负载系统的故障原因。
为达到上述目的,本发明提供如下技术方案:
一种基于多域稀疏表达的电动客车高压负载系统故障诊断方法,考虑到电动客车的行驶方案,本发明的故障诊断与定位主要在出发间隔或停车时间内进行诊断,具体包括以下步骤:
S1:在进行故障诊断前,对电动客车进行多次测试,获取其正常模式或不同故障模式下由标准激励信号引起的电流响应信号,将响应信号定义为标记数据样本,利用标记的数据样本构造时域、频域和几何域下的字典矩阵;
S2:故障诊断时,向高压负载系统注入一个标准的激励信号,并收集由激励信号引起的响应电流;分别在时域、频域和几何域上用标记样本稀疏表达电流响应信号,采用批量匹配追踪算法计算出对应的时域稀疏向量、频域稀疏向量和几何域稀疏向量;
S3:将时域稀疏向量、频域稀疏向量和几何域稀疏向量融合为联合稀疏向量,再对该联合稀疏向量进行处理,识别测试样本的模式。
进一步,在步骤S1中,具体包括:电动客车高压负载系统的故障模式包括:电池加热膜内部短路、驱动电机绕组匝间短路、DC/DC输出端子短路、母线电容性能下降(母线电容减少)以及PTC加热器性能下降(电阻增加);基于电动客车电压平台,将标准激励信号设计为多级跃迁信号,分为三个阶段:第1阶段,时间为0s~0.2s,输出电压0V;第2阶段,时间为0.2s~1s;第3阶段,时间为1s~2s。
进一步,在步骤S1中,构造时域、频域和几何域下的字典矩阵:
(1)构造时域字典矩阵:根据获取的电流响应信号,直接表示为:
DT=[DT1 DT2 ... DTi ... DTp]N×pq
其中,DTi是由ith模式的标记样本组成的矩阵,i=1,2,...,p,p表示模式数,q表示从一个模式中收集到的标记样本数,并且:
DTi=[xi,1 xi,2 ... xi,j ... xi,q]
其中,xi,j是模式i的jth标记数据样本,j=1,2,...,q。
(2)构造频域字典矩阵:
首先将原始信号转换为频谱:由于激励信号第0s~0.2s时输出电压为0V,响应信号在这段时间内的数据点为0,因此不考虑输出电压为0时的采样点,而采用后续的数据点来构造频域字典矩阵,设截断信号中的数据点个数为N,;
假设截断信号为xtr(n),n=1,2,...,N,频谱为y(k),k=0,1,...,N-1,则:
由于采样点xtr(n)是实数,因此:
频谱是对称的,只有一半频谱有效,因此采用y(k)中的前半段数据点来构造频域字典矩阵,字典中原子的长度为R,得到频域字典矩阵为:
DF=[DF1 DF2 ... DFi ... DFp]R×pq
其中,DFi是由ith模式的频谱组成的矩阵,i=1,2,...,p,p表示模式数,q表示从一个模式中收集到的标记样本数,:
DFi=[yi,1 yi,2 ... yi,j ... yi,q],其中yi,j是模式i的jth标记频谱,j=1,2,...,q。
(3)构造几何域字典矩阵:
为了提高诊断精度,从响应信号中提取14个几何特征,这些特征包含6个信号特征和8个动态性能指标;从响应信号的前节段(t=0.2s~1s)提取4个指标,从响应信号的后节段(t=1s~2s)提取4个指标;
1)前段脉搏强度:PSA=xi,j|t=0.2000s;
2)后段脉搏强度:PSP=xi,j|t=1.0000s;
3)前段的最小曲率半径:xi,j|t=0.2101s~t=0.2500s的最小曲率半径,记为MCR;
4)前段最小曲率半径位置:最小曲率圆的圆心在t轴上的位置,记为PMCR;
5)有效反应平均值:
6)前段值稳定:
7)前段高峰时间:定义前段的峰值为PVA=max(xi,j|t=0.2001s~t=0.9999s),响应信号首先在时间t=tPVA达到前段的峰值前段的峰值时间为
PTA=tPVA-0.2000s;
8)前段上升时间:在t=tSVA,t∈(0.2001s,tPVA)时刻,响应信号首先到达前段的稳定值:
那么前段的上升时间为:
RTA=tSVA-0.2000s;
9)前段超调量:
10)后段稳定值:
11)后段谷时间:将后段的谷值定义为VVP=min(xi,j|t=1.0001s~2.0000s),响应信号在t=1.0001s之后首先在t=tVVP时间到达后段谷值:
那么后段的谷时间为:
VTP=tVVP-1.0000s;
12)后段下降时间:在t=tSVP,t∈(1.0001s,tVVP)时刻,响应信号首先达到后段的稳定值:
那么后段下降时间为:
FTP=tSVP-1.0000s;
13)后段超调量:
14)前段和后段的稳定值差:
SVD=SVA-SVP;
将上述14个几何特征按如下顺序组成几何特征向量:
其中gi,j是模式i的jth标记几何特征向量,j=1,2,...,q,再用这些从不同模式中提取的几何特征向量构造几何域字典矩阵:
DG=[DG1 DG2 ... DGi ... DGp]14×pq,
其中,DGi由从ith模式中提取的几何特征向量组成,i=1,2,...,p,p表示模式数,q表示从一个模式中收集到的标记样本数,
DGi=[gi,1 gi,2 ... gi,j ... gi,q]。
进一步,在步骤S2中,所述用标记样本稀疏表达电流响应信号,在时域、频域以及几何域中,测试样本分别用过完备字典表示,其中过完备字典矩阵为不同模式下标记样本在时域、频域和几何域中构造的字典矩阵;
在时域中,直接用时域字典矩阵稀疏表达测试样本,得到测试样本的时域稀疏向量;
在频域中,先通过离散傅里叶变换将测试样本变换为电流响应频谱,再用频域字典矩阵稀疏表达测试样本的频谱,得到频域稀疏向量;
在几何域中,先从测试样本中提取出几何特征构成几何特征向量,再用几何域字典矩阵稀疏表达测试样本的几何特征向量,得到几何域稀疏向量。
进一步,所述稀疏表达具体为:假设测试样本为zk,k=x,y,g,zx,18000×1为时域响应信号,zy,9000×1为测试样本的频谱,zg,14×1为测试样本的几何特征向量;Dk,k=x,y,g为过完备字典矩阵,Dx,18000×pq为时域字典矩阵,Dy,9000×pq为频域字典矩阵,Dg,14×pg为几何域字典矩阵;稀疏表达表示为:
zk=Dk·θk,pq×1,k=x,y,g
其中,θk,pq×1,k=x,y,g为稀疏向量,由批量匹配追踪算法获得。
进一步,所述批量匹配追踪算法,具体如下:
算法输入:字典矩阵:Dk;测试样品:zk;图案数量:p;每批支持向量数nSV;迭代次数:nIter;
算法输出:稀疏向量:θk,pq×1;
中间变量:迭代计数器:times=1,2,...,nIter;不同向量之间的余弦相似度:Product;非零元素在θk中的位置:θk中的非零元素:/>为测试样本选择的支持向量:MSV;残差向量:Res;
运算过程:
a:变量初始化:times=1;最初的θk:θk,0=(0 0 ... 0)T pq×1;θk中非零元素的位置索引:θk,Pos,0=[ ];初始残差向量:Res0=zk;
b:计算:计算残差向量Restimes-1与字典矩阵中每个列向量kcol的余弦相似度,其中k=x,y,g,col=1,2,...,pq:
Productcol=(kcol T·Restimes-1)/(||kcol||·||Restimes-1||)
然后选择Productcol的nSV最大值,这些值的位置索引为:
位置索引向量θk,Pos通过迭代过程进行更新和扩展:
θk,Pos=θk,Pos∪θk,Pos,times
基于位置索引向量,选出对表示测试样本贡献最大的一批支持向量:
同时,将字典矩阵Dk中列索引为θk,Pos的列向量设置为空向量:
之后用最小二乘法计算θk中的元素对应于MSV,times中的支持向量:
θk,Element,times=(MSV,times T·MSV,times)-1·MSV,times T·zk
再计算新的残差向量Restimes:
Restimes=zk-MSV,times·θElement,times
c:迭代:将上述b过程执行nIter次,得到θk中的非零元素θk,Element及其位置θk,Pos,稀疏向量为:
θk,pq×1:θk,θk,Pos=θk(θk,Pos)=θk,Element。
进一步,在步骤S3中,首先融合时域稀疏向量θx,pq×1、频域稀疏向量θy,pq×1和几何域稀疏向量θg,pq×1:
其中,θpq×1是融合的稀疏向量,θij是ith模式下jth测试样本的稀疏向量,i=1,2,...,p,j=1,2,...,q,p表示模式数,q表示从一个模式中收集到的标记样本数;
从θpq×1中提取一个权重分布向量,记为W:
Wp×1=(w1 w2 ... wp)T
其中,wi是与ith模式对应的融合稀疏向量中的元素的总和:
与测试样本具有相同模式的标记样本向量对稀疏表达的贡献最大,因此,权重分布向量W中最大元素的索引即为测试样本的模式或故障模式:
Pattern=i,s.t.max(wi),i=1,2,...,p。
本发明的有益效果在于:本发明能及时发现和消除电动客车高压系统的潜在风险,避免因高压系统故障而引起安全事故,同时避免给公共安全和个人安全带来严重影响。本发明具有准确率高、实用性强和部署简单方便等特点,对于工程应用具有重要意义。
本发明的其他优点、目标和特征在某种程度上将在随后的说明书中进行阐述,并且在某种程度上,基于对下文的考察研究对本领域技术人员而言将是显而易见的,或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书来实现和获得。
附图说明
为了使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作优选的详细描述,其中:
图1为高压负载故障诊断方案;
图2为高压负载故障诊断系统;
图3为稀疏向量融合与模式识别;
图4为基于多域稀疏表达的数据分析;
图5为响应信号的几何特征;
图6为批量匹配追踪算法的稀疏表达;
图7为稀疏向量的模式识别过程。
具体实施方式
以下通过特定的具体实例说明本发明的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本发明的精神下进行各种修饰或改变。需要说明的是,以下实施例中所提供的图示仅以示意方式说明本发明的基本构想,在不冲突的情况下,以下实施例及实施例中的特征可以相互组合。
其中,附图仅用于示例性说明,表示的仅是示意图,而非实物图,不能理解为对本发明的限制;为了更好地说明本发明的实施例,附图某些部件会有省略、放大或缩小,并不代表实际产品的尺寸;对本领域技术人员来说,附图中某些公知结构及其说明可能省略是可以理解的。
本发明实施例的附图中相同或相似的标号对应相同或相似的部件;在本发明的描述中,需要理解的是,若有术语“上”、“下”、“左”、“右”、“前”、“后”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此附图中描述位置关系的用语仅用于示例性说明,不能理解为对本发明的限制,对于本领域的普通技术人员而言,可以根据具体情况理解上述术语的具体含义。
请参阅图1~图7,为一种基于多域稀疏表达的电动客车高压负载系统故障诊断方法。电动客车高压系统由动力电池系统和高压负载系统组成。动力电池系统由相对成熟的电池管理系统(BMS)管理和保护,本发明的诊断主要集中在高压负载系统。高压负载系统由驱动电机、电池电加热膜、电除霜器、转向电机、空压机、DC/DC、母线电容和PTC加热器等组成。在电动客车的高电压负载中,根据历史维护数据,发现电加热膜、驱动电机、DC/DC、母线电容、PTC加热器等负载故障率较高。同时,这些负载可在快速维修服务站快速更换。因此,本发明旨在快速准确地诊断出电动客车高压负载系统的故障原因。
图1所示为本发明的电动客车高压负载系统故障诊断基本结构。电动客车在断电状态下,在出发间隔或停车期间,向高压负载系统注入一个标准的激励信号,由激励信号引起的响应电流通过高压负载系统。为了收集响应电流信号,电流传感器串联到电路上。然后,利用高速数据采集装置采集这些响应电流信号。
图2所示为激励信号的产生和响应信号的采集:高压负载故障诊断系统与负载供电电路并联,部署在电动客车高压配电单元(PDU)中。电动客车驱动过程中K1、K2开关闭合,K3/K4/K5/K6开关打开;在停车时间内,K1、K2开关打开,K3/K4/K5/K6开关关闭,诊断系统接入高压系统。诊断系统由信号发生器和数据采集模块组成。基于电动客车电压平台,将标准激励信号设计为多级阶跃信号,如表1所示:
表1标准激励信号
阶段 | 1 | 2 | 3 |
时间 | 0~0.2 | 0.2~1 | 1~2 |
电压(V) | 0 | 1000 | 576 |
图3和图4所示为数据处理方法:数据采集模块的采样频率为10KS/s,从而在诊断过程中采集数据点。但由于响应信号的第1-2000个数据点为零,因此不考虑这2000个数据点,而采用第2001-20000个数据点构造时域字典矩阵,则截断信号中的数据点个数为N=18000。随后,利用图3和图4所示的算法对采集到的数据进行分析。
假设从一个模式中收集到的标记样本数为q,模式数为p,则时域字典矩阵表示为:
DT=[DT1 DT2 ... DTi ... DTp]N×pq
其中,DTi是由ith模式的标记样本组成的矩阵,i=1,2,...,p,并且:
DTi=[xi,1 xi,2 ... xi,j ... xi,q]
其中,xi,j是模式i的jth标记数据样本,j=1,2,...,q。
如表2所示,选择了5种故障模式,此外,正常状态也是一种模式,因此,在本发明中,p=6。
表2高压负载系统的故障模式
故障部件 | 故障模式 |
电池加热膜 | 内部短路 |
驱动电机 | 绕组匝间短路 |
DC/DC输出端子 | 输出端子短路 |
母线电容 | 性能下降(母线电容减少) |
PTC | 性能下降(电阻增加) |
数据转换和特征提取:
时频转换:在频域稀疏表达之前,将原始信号变换为频谱。由于响应信号的幅值在0~0.2s为零,所以第1-2000个采样点不进行变换。因此,截断信号中的数据点数是18000。
假设截断信号为xtr(n),n=1,2,...,N,频谱为y(k),并且k=0,1,...,N-1,这里n=1,2,...,N,那么:
考虑到采样点xtr(n)是实数:
频谱是对称的,因此只有一半的频谱有效,这里用y(k)中的第1-9000个数据点来构造频域字典矩阵,字典中原子的长度是R=9000。
因此,频域字典矩阵为:
DF=[DF1 DF2 ... DFi ... DFp]R×pq
其中,DFi是频域字典矩阵,由ith模式的频谱组成,i=1,2,...,p:
DFi=[yi,1 yi,2 ... yi,j ... yi,q]
这里,yi,j是模式i的jth标记频谱,j=1,2,...,q。
图5所示为响应信号的几何特征提取:为了提高诊断精度,从响应信号中提取了14个几何特征,这些特征包含6个信号特征和8个动态性能指标;从响应信号的前节段(t=0.2s~1s)提取4个指标,从响应信号的后节段(t=1s~2s)提取4个指标。
(1)前段脉搏强度:
PSA=xi,j|t=0.2000s
(2)后段脉搏强度:
PSP=xi,j|t=1.0000s
(3)前段的最小曲率半径:
xi,j|t=0.2101s~t=0.2500s的最小曲率半径,记为MCR。
(4)前段最小曲率半径位置:
圆心(最小曲率圆)在t轴上的位置,记为PMCR。
(5)有效反应平均值:
(6)前段值稳定:
(7)前段高峰时间:
定义前段的峰值为:
PVA=max(xi,jt=0.2001s~0.9999s)
响应信号首先在时间t=tPVA达到前段的峰值
前段的峰值时间为:
PTA=tPVA-0.2000s
(8)前段上升时间:
在t=tSVA,t∈(0.2001s,tPVA)时刻,响应信号首先达到前段的稳定值:
前段的上升时间为:
RTA=tSVA-0.2000s
(9)前段超调量:
(10)后段稳定值:
(11)后段谷时间:将后段的谷值定义为:
VVP=min(xi,j|t=1.0001s~2.0000s)
响应信号在t=1.0001s之后首先在t=tVVP时间到达后段谷值:
后段的谷时间为:
VTP=tVVP-1.0000s
(12)后段下降时间:
在t=tSVP,t∈(1.0001s,tVVP)时刻,响应信号首先达到后段的稳定值:
后段下降时间为:
FTP=tSVP-1.0000s
(13)后段超调量:
(14)前段和后段的稳定值差:
SVD=SVA-SVP
几何特征向量由这14个特征按如下顺序组成,gi,j是模式i的jth标记几何特征向量:
然后,用这些从不同模式中提取的几何特征向量构造几何域字典矩阵:
DG=[DG1 DG2 ... DGi ... DGp]14×pq
其中,DGi由从ith模式中提取的几何特征向量组成,i=1,2,...,p:
DGi=[gi,1 gi,2 ... gi,j ... gi,q]
图6所示为基于批量匹配追踪算法的稀疏表达:在时域、频域和几何域,测试样本分别用过完备字典表示,过完备字典由标记样本在时域、频域及几何域下的字典矩阵组成。
假设测试样本为zk,k=x,y,g,zx,18000×1为时域响应信号,zy,9000×1为测试样本的频谱,zg,14×1为测试样本的几何特征向量。Dk,k=x,y,g为过完备字典,Dx,18000×pq为时域字典矩阵,Dy,9000×pq为频域字典矩阵,Dg,14×pg为几何域字典矩阵。然后,稀疏表达表示为:
zk=Dk·θk,pq×1,k=x,y,g
其中,θk,pq×1,k=x,y,g为稀疏向量,未知且待确定。在本发明中,θk是通过批量匹配追踪算法获得的:
输入:
字典矩阵:Dk;
测试样品:zk;
图案数量:p;
每批支持向量数:nSV;;
迭代次数:nIter;
输出:
稀疏向量:θk,pq×1;
中间变量:
迭代计数器:times=1,2,...,nIter;
不同向量之间的余弦相似度:Product;
非零元素在θk中的位置:
θk中的非零元素:
为测试样本选择的支持向量:MSV;
残差向量:Res;
程序:
a.变量初始化:
迭代计数器:times=1;
最初的θk:θk,0=(0 0 ... 0)T pq×1;
θk中非零元素的位置索引:θk,Pos,0=[ ];
初始残差向量:Res0=zk;
b.计算
计算残差向量Restimes-1与字典矩阵kcol中每个列向量的余弦相似度,其中k=x,y,g,col=1,2,...,pq:
Productcol=(kcol T·Restimes-1)/(||kcol||·||Restimes-1||)
然后,选择Productcol的nSV最大值,这些值的位置索引为:
位置索引向量通过迭代过程进行更新和扩展:
θk,Pos=θk,Pos∪θk,Pos,times
基于位置索引向量,选出对表示测试样本贡献最大的一批支持向量:
同时,将字典矩阵Dk中列索引为的列向量设置为空向量:
之后,用最小二乘法计算中的元素对应于MSV,times中的支持向量:
然后,计算新的残差向量Restimes:
c.迭代
等式到等式中描述的过程执行nIter次。
最后得到中的非零元素/>及其位置/>稀疏向量/>为:
/>
图7所示为模式识别与故障诊断:
首先融合时域稀疏向量频域稀疏向量/>和几何域稀疏向量/>
这里,是融合的稀疏向量。然后,从融合后的向量中提取一个权重分布向量,该权重分布向量记为W:
Wp×1=(w1 w2 … wp)T
其中,wi是与ith模式对应的融合稀疏向量中元素的总和:
基于稀疏表达理论,与测试样本具有相同模式的标记样本向量对稀疏表达的贡献最大,因此,权重分布向量W中最大元素的索引是测试的模式或(故障)模式样本:
Pattern=i,s.t.max(wi),i=1,2,…,p。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (9)
1.一种基于多域稀疏表达的电动客车高压负载系统故障诊断方法,其特征在于:包括以下步骤:
S1:在进行故障诊断前,对电动客车进行多次测试,获取其正常模式或不同故障模式下由标准激励信号引起的电流响应信号,将响应信号定义为标记数据样本,利用标记的数据样本构造时域、频域和几何域下的字典矩阵;
S2:故障诊断时,向高压负载系统注入一个标准的激励信号,并收集由激励信号引起的响应电流;分别在时域、频域和几何域上用标记样本稀疏表达电流响应信号;采用批量匹配追踪算法计算出对应的时域稀疏向量、频域稀疏向量和几何域稀疏向量;
S3:将时域稀疏向量、频域稀疏向量和几何域稀疏向量融合为联合稀疏向量,再对该联合稀疏向量进行处理,识别测试样本的模式。
2.根据权利要求1所述的一种基于多域稀疏表达的电动客车高压负载系统故障诊断方法,其特征在于:电动客车高压负载系统的故障模式包括:电池加热膜内部短路、驱动电机绕组匝间短路、DC/DC输出端子短路、母线电容减少以及PTC加热器电阻增加;标准激励信号设计为多级跃迁信号,分为三个阶段:第1阶段,时间为0s~0.2s,输出电压0V;第2阶段,时间为0.2s~1s;第3阶段,时间为1s~2s。
3.根据权利要求1所述的一种基于多域稀疏表达的电动客车高压负载系统故障诊断方法,其特征在于:在步骤S1中,构造时域字典矩阵:根据获取的电流响应信号,直接表示为:
DT=[DT1 DT2...DTi...DTp]N×pq
其中,DTi是由ith模式的标记样本组成的矩阵,i=1,2,...,p,其中p表示模式数,q表示从一个模式中收集到的标记样本数,并且:
DTi=[xi,1 xi,2...xi,j...xi,q]
其中,xi,j是模式i的jth标记数据样本,j=1,2,...,q。
4.根据权利要求1所述的一种基于多域稀疏表达的电动客车高压负载系统故障诊断方法,其特征在于:在步骤S1中,构造频域字典矩阵:
首先将原始信号转换为频谱:不考虑输出电压为0时的采样点,截断信号中的数据点数为N;
假设截断信号为xtr(n),n=1,2,...,N,频谱为y(k),k=0,1,...,N-1,则:
由于采样点xtr(n)是实数,因此:
采用y(k)的前半段数据点构造频域字典矩阵,字典中原子的长度为R,得到频域字典矩阵为:
DF=[DF1 DF2...DFi...DFp]R×pq
其中,DFi是由ith模式的频谱组成的矩阵,i=1,2,...,p,其中p表示模式数,q表示从一个模式中收集到的标记样本数,
DFi=[yi,1 yi,2...yi,j...yi,q],
其中yi,j是模式i的jth标记频谱,j=1,2,...,q。
5.根据权利要求1所述的一种基于多域稀疏表达的电动客车高压负载系统故障诊断方法,其特征在于:在步骤S1中,构造几何域字典矩阵:首先从响应信号中提取14个几何特征:
前段脉搏强度:PSA=xi,j|t=0.2000s;
后段脉搏强度:PSP=xi,j|t=1.0000s;
前段的最小曲率半径:xi,j|t=0.2101s~t=0.2500s的最小曲率半径,记为MCR;
前段最小曲率半径位置:最小曲率圆的圆心在t轴上的位置,记为PMCR;
有效反应平均值:
前段值稳定:
前段高峰时间:定义前段的峰值为PVA=max(xi,j|t=0.2001s~t=0.9999s),响应信号首先在时间t=tPVA达到前段的峰值前段的峰值时间为
PTA=tPVA-0.2000s;
前段上升时间:在t=tSVA,t∈(0.2001s,tPVA)时刻,响应信号首先到达前段的稳定值:
那么前段的上升时间为:
RTA=tSVA-0.2000s;
前段超调量:
后段稳定值:
后段谷时间:将后段的谷值定义为VVP=min(xi,j|t=1.0001s~2.0000s),响应信号在t=1.0001s之后首先在t=tVVP时间到达后段谷值:
那么后段的谷时间为:
VTP=tVVP-1.0000s;
后段下降时间:在t=tSVP,t∈(1.0001s,tVVP)时刻,响应信号首先达到后段的稳定值:
那么后段下降时间为:
FTP=tSVP-1.0000s;
后段超调量:
前段和后段的稳定值差:
SVD=SVA-SVP;
将上述14个几何特征按如下顺序组成几何特征向量:
其中gi,j是模式i的jth标记几何特征向量,j=1,2,...,q,再用这些从不同模式中提取的几何特征向量构造几何域字典矩阵:
DG=[DG1 DG2...DGi...DGp]14×pq,
其中,DGi由从ith模式中提取的几何特征向量组成,i=1,2,...,p,其中p表示模式数,q表示从一个模式中收集到的标记样本数,
DGi=[gi,1 gi,2...gi,j...gi,q]。
6.根据权利要求1所述的一种基于多域稀疏表达的电动客车高压负载系统故障诊断方法,其特征在于:在步骤S2中,所述用标记样本稀疏表达电流响应信号,在时域、频域以及几何域中,测试样本分别用过完备字典矩阵表示,其中过完备字典矩阵为不同模式下标记样本在时域、频域和几何域中构造的字典矩阵;
在时域中,直接用时域字典矩阵稀疏表达测试样本,得到测试样本的时域稀疏向量;
在频域中,先通过离散傅里叶变换将测试样本变换为电流响应频谱,再用频域字典矩阵稀疏表达测试样本的频谱,得到频域稀疏向量;
在几何域中,先从测试样本中提取出几何特征构成几何特征向量,再用几何域字典矩阵稀疏表达测试样本的几何特征向量,得到几何域稀疏向量。
7.根据权利要求6所述的一种基于多域稀疏表达的电动客车高压负载系统故障诊断方法,其特征在于:所述稀疏表达具体为:
假设测试样本为zk,k=x,y,g,zx,18000×1为时域响应信号,zy,9000×1为测试样本的频谱,zg,14×1为测试样本的几何特征向量;Dk,k=x,y,g为过完备字典矩阵,Dx,18000×pq为时域字典矩阵,Dy,9000×pq为频域字典矩阵,Dg,14×pg为几何域字典矩阵;稀疏表达表示为:
zk=Dk·θk,pq×1,k=x,y,g
其中,θk,pq×1,k=x,y,g为稀疏向量,由批量匹配追踪算法获得。
8.根据权利要求7所述的一种基于多域稀疏表达的电动客车高压负载系统故障诊断方法,其特征在于:所述批量匹配追踪算法,其过程为:
算法输入:字典矩阵:Dk;测试样品:zk;图案数量:p;每批支持向量数nSV;迭代次数:nIter;
算法输出:稀疏向量:θk,pq×1;
中间变量:迭代计数器:times=1,2,...,nIter;不同向量之间的余弦相似度:Product;非零元素在θk中的位置:θk中的非零元素:/>为测试样本选择的支持向量:MSV;残差向量:Res;
运算过程:
a:变量初始化:times=1;最初的θk:θk,0=(0 0...0)T pq×1;θk中非零元素的位置索引:θk,Pos,0=[];初始残差向量:Res0=zk;
b:计算:计算残差向量Restimes-1与字典矩阵中每个列向量kcol的余弦相似度,其中k=x,y,g,col=1,2,...,pq:
Productcol=(kcol T·Restimes-1)/(||kcol||·||Restimes-1||)
然后选择Productcol的nSV最大值,这些值的位置索引为:
位置索引向量θk,Pos通过迭代过程进行更新和扩展:
θk,Pos=θk,Pos∪θk,Pos,times
基于位置索引向量,选出对表示测试样本贡献最大的一批支持向量:
同时,将字典矩阵Dk中列索引为θk,Pos的列向量设置为空向量:
之后用最小二乘法计算θk中的元素对应于MSV,times中的支持向量:
θk,Element,times=(MSV,times T·MSV,times)-1·MSV,times T·zk
再计算新的残差向量Restimes:
Restimes=zk-MSV,times·θElement,times
c:迭代:将上述b过程执行nIter次,得到θk中的非零元素θk,Element及其位置θk,Pos,稀疏向量为:
θk,pq×1:θk,θk,Pos=θk(θk,Pos)=θk,Element。
9.根据权利要求1所述的一种基于多域稀疏表达的电动客车高压负载系统故障诊断系统,其特征在于:在步骤S3中,首先融合时域稀疏向量θx,pq×1、频域稀疏向量θy,pq×1和几何域稀疏向量θg,pq×1:
其中,θpq×1是融合的稀疏向量,θij是ith模式下jth测试样本的稀疏向量,i=1,2,...,p,j=1,2,...,q,p表示模式数,q表示从一个模式中收集到的标记样本数;
从θpq×1中提取一个权重分布向量,记为W:
Wp×1=(w1 w2...wp)T
其中,wi是与ith模式对应的融合稀疏向量中的元素的总和:
与测试样本具有相同模式的标记样本向量对稀疏表达的贡献最大,因此,权重分布向量W中最大元素的索引即为测试样本的模式或故障模式:
Pattern=i,s.t.max(wi),i=1,2,...,p。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210374371.0A CN114714911B (zh) | 2022-04-11 | 2022-04-11 | 一种基于多域稀疏表达的电动客车高压负载系统故障诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210374371.0A CN114714911B (zh) | 2022-04-11 | 2022-04-11 | 一种基于多域稀疏表达的电动客车高压负载系统故障诊断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114714911A CN114714911A (zh) | 2022-07-08 |
CN114714911B true CN114714911B (zh) | 2024-05-28 |
Family
ID=82243110
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210374371.0A Active CN114714911B (zh) | 2022-04-11 | 2022-04-11 | 一种基于多域稀疏表达的电动客车高压负载系统故障诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114714911B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116361727A (zh) * | 2023-03-28 | 2023-06-30 | 重庆大学 | 一种基于音频特征和SRC-Adaboost的电池换电系统驱动齿轮故障诊断方法 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102435934A (zh) * | 2011-09-20 | 2012-05-02 | 北京工业大学 | 一种随机采样的模拟电路压缩传感测量和信号重构方法 |
CN103728130A (zh) * | 2013-10-10 | 2014-04-16 | 西安交通大学 | 一种基于稀疏分解的风力发电机组故障特征提取方法 |
CN105738109A (zh) * | 2016-02-22 | 2016-07-06 | 重庆大学 | 基于稀疏表示与集成学习的轴承故障分类诊断方法 |
CN106828106A (zh) * | 2016-12-29 | 2017-06-13 | 北京交通大学 | 一种面向动车组牵引电机的故障预警方法 |
CN109061537A (zh) * | 2018-08-23 | 2018-12-21 | 重庆大学 | 基于观测器的电动车辆锂离子电池传感器故障诊断方法 |
CN109975697A (zh) * | 2019-03-04 | 2019-07-05 | 河南理工大学 | 一种基于原子稀疏分解的高压断路器机械故障诊断方法 |
CN110954761A (zh) * | 2019-11-04 | 2020-04-03 | 南昌大学 | 基于信号稀疏表示的npc三电平逆变器故障诊断方法 |
CN112965001A (zh) * | 2021-02-09 | 2021-06-15 | 重庆大学 | 一种基于实车数据的动力电池组故障诊断方法 |
CN113140054A (zh) * | 2020-01-02 | 2021-07-20 | 株洲中车时代电气股份有限公司 | 一种列车冷却系统的异常检测方法及装置 |
WO2022011754A1 (zh) * | 2020-07-16 | 2022-01-20 | 苏州大学 | 一种基于自适应流形嵌入动态分布对齐的故障诊断方法 |
CN114217256A (zh) * | 2021-11-25 | 2022-03-22 | 华中科技大学 | 基于自适应稀疏滤波的电机驱动系统变流器故障诊断方法 |
-
2022
- 2022-04-11 CN CN202210374371.0A patent/CN114714911B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102435934A (zh) * | 2011-09-20 | 2012-05-02 | 北京工业大学 | 一种随机采样的模拟电路压缩传感测量和信号重构方法 |
CN103728130A (zh) * | 2013-10-10 | 2014-04-16 | 西安交通大学 | 一种基于稀疏分解的风力发电机组故障特征提取方法 |
CN105738109A (zh) * | 2016-02-22 | 2016-07-06 | 重庆大学 | 基于稀疏表示与集成学习的轴承故障分类诊断方法 |
CN106828106A (zh) * | 2016-12-29 | 2017-06-13 | 北京交通大学 | 一种面向动车组牵引电机的故障预警方法 |
CN109061537A (zh) * | 2018-08-23 | 2018-12-21 | 重庆大学 | 基于观测器的电动车辆锂离子电池传感器故障诊断方法 |
CN109975697A (zh) * | 2019-03-04 | 2019-07-05 | 河南理工大学 | 一种基于原子稀疏分解的高压断路器机械故障诊断方法 |
CN110954761A (zh) * | 2019-11-04 | 2020-04-03 | 南昌大学 | 基于信号稀疏表示的npc三电平逆变器故障诊断方法 |
CN113140054A (zh) * | 2020-01-02 | 2021-07-20 | 株洲中车时代电气股份有限公司 | 一种列车冷却系统的异常检测方法及装置 |
WO2022011754A1 (zh) * | 2020-07-16 | 2022-01-20 | 苏州大学 | 一种基于自适应流形嵌入动态分布对齐的故障诊断方法 |
CN112965001A (zh) * | 2021-02-09 | 2021-06-15 | 重庆大学 | 一种基于实车数据的动力电池组故障诊断方法 |
CN114217256A (zh) * | 2021-11-25 | 2022-03-22 | 华中科技大学 | 基于自适应稀疏滤波的电机驱动系统变流器故障诊断方法 |
Non-Patent Citations (2)
Title |
---|
Fault Diagnosis of Electric City Bus High-Voltage Load System Based on Multidomain Sparse Representation;Xianglong You;《IEEE TRANSACTIONS ON TRANSPORTATION ELECTRIFICATION》;20240331;第10卷(第1期);1207-1221 * |
大数据下机械智能故障诊断的机遇与挑战;雷亚国;贾峰;孔德同;林京;邢赛博;;机械工程学报;20170822(05);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114714911A (zh) | 2022-07-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Guo et al. | A data-driven remaining capacity estimation approach for lithium-ion batteries based on charging health feature extraction | |
Wang et al. | Probability based remaining capacity estimation using data-driven and neural network model | |
CN101131417B (zh) | 电池监测系统 | |
CN114714911B (zh) | 一种基于多域稀疏表达的电动客车高压负载系统故障诊断方法 | |
CN111007417A (zh) | 基于不一致性评估的电池组soh和rul预测方法及系统 | |
Lu et al. | Analysis of the key factors affecting the energy efficiency of batteries in electric vehicle | |
CN108303649A (zh) | 一种电池健康状态识别方法 | |
CN110703112A (zh) | 一种基于局部充电数据的电池组状态的在线估计方法 | |
CN113064939A (zh) | 一种新能源车辆三电系统安全特征数据库构建方法 | |
CN112965001B (zh) | 一种基于实车数据的动力电池组故障诊断方法 | |
CN113608140A (zh) | 一种电池故障诊断方法及系统 | |
Khalfi et al. | An electric circuit model for a lithium-ion battery cell based on automotive drive cycles measurements | |
Amirkhani et al. | Electric vehicles driving range and energy consumption investigation: A comparative study of machine learning techniques | |
Yin et al. | Voltage fault diagnosis of power batteries based on boxplots and gini impurity for electric vehicles | |
CN115879048A (zh) | 一种基于wrfmda模型的串联电弧故障识别方法和系统 | |
Chen et al. | Meta rule-based energy management strategy for battery/supercapacitor hybrid electric vehicles | |
CN113687251A (zh) | 一种基于双模型的锂离子电池组电压异常故障诊断方法 | |
CN116911161A (zh) | 一种结合数据增强的深度学习暂态电压稳定评估方法 | |
CN116736170A (zh) | 基于半监督学习和k均值聚类算法的动力电池SOH估计方法 | |
CN113419172A (zh) | 基于gcn和vmd-ed的新能源汽车三相异步电机故障识别 | |
Zhang et al. | Development of the driving cycle for Dalian city | |
Nair et al. | Predicting the Impact of Road Conditions on Battery Health Via Transfer Learning | |
CN113420494A (zh) | 一种超级电容贝叶斯概率融合建模方法 | |
Guo et al. | Early diagnosis of battery faults through an unsupervised health scoring method for real-world applications | |
CN114814391A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |