CN114510778A - 基于混合智能优化lstm的轨道不平顺预测方法 - Google Patents
基于混合智能优化lstm的轨道不平顺预测方法 Download PDFInfo
- Publication number
- CN114510778A CN114510778A CN202210021594.9A CN202210021594A CN114510778A CN 114510778 A CN114510778 A CN 114510778A CN 202210021594 A CN202210021594 A CN 202210021594A CN 114510778 A CN114510778 A CN 114510778A
- Authority
- CN
- China
- Prior art keywords
- lstm
- model
- track irregularity
- prediction
- pso
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- 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
- G06F30/25—Design optimisation, verification or simulation using particle-based methods
-
- 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
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/044—Recurrent networks, e.g. Hopfield networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- Artificial Intelligence (AREA)
- Software Systems (AREA)
- Computational Linguistics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computing Systems (AREA)
- General Health & Medical Sciences (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Automation & Control Theory (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开一种基于混合智能优化LSTM的轨道不平顺预测方法,首先将时间序列数据进行预处理,然后利用PSO算法优化LSTM模型超参数,从而确定LSTM模型的网络结构;利用GA算法优化LSTM模型初始的权值阈值,确定LSTM模型的权值阈值。最后利用确定的超参数和权值阈值,对轨道不平顺数据进行训练并预测。本发明基于LSTM‑PSO‑GA模型预测轨道不平顺数据的方法,克服了传统的预测方法预测过程中精度不高的问题,且利用PSO和GA算法对LSTM参数进行优化,避免模型陷入局部最优解的问题,提高了预测收敛速度。最终实现了对轨道不平顺数据的预测,更精确预测轨道不平顺现象。
Description
技术领域
本发明属于时间序列预测技术领域,具体涉及一种基于混合智能优化LSTM的轨道不平顺预测方法。
背景技术
我国铁路交通高速化和重载化的不断发展,轨道反复承受机车车辆的荷载,逐渐改变轨道的轨距、水平及高低等几何形位和空间位置,造成轨道不平顺。铁路轨道不平顺直接威胁列车运行安全。基于轨道不平顺的历史数据,预测轨道不平顺状态发展趋势,对轨道安全性进行检查,可为轨道维护策略提供重要的理论依据。
近年来,国内外诸多学者研究轨道不平顺预测问题。雍安泰对轨道质量状态进行预测,建立轨道质量指数(Track Quality Index,TQI)的线性回归方程,并用最小二乘法求解。曲建军采用灰色***不确定性理论,以TQI时间序列数据作为指标,建立灰色非线性轨道不平顺预测模型。徐波提出基于非等间距灰色线性加权组合预测模型,采用加权组合函数对预测残差序列进行修正,提高了预测精度。贾朝龙提出分段神经网络-ARMA递推模型,对小波分解后的低频近似序列采用神经网络递推模型,高频细节序列采用ARMA模型,对轨道不平顺状态进行预测。马子骥先利用非等间距灰色模型对预测TQI序列大致趋势,然后利用Elman神经网络对初步预测结果进行残差修正,从而更准确预测TQI序列。郭浩龙提出基于灰色理论与支持向量机的组合模型,预测TQI数据。姚亚峰提出基于灰色模型和循环神经网络的组合预测模型,预测轨道不平顺趋势。史红梅提出基于遗传算法(GeneticAlgorithm,GA)和车辆轨道耦合模型相结合的轨道不平顺估计算法,将轨道不平顺问题转换为模型参数的估计问题进行求解。然而,上述方法没有考虑数据波动变化大时干扰预测结果的情况,且没有充分挖掘轨道不平顺数据的特征信息,难以形成稳定精确的轨道不平顺预测模型。
轨道不平顺数据具有非线性、随机性及突发性。长短期记忆(Long Short-TermMemory,LSTM)循环神经网络可根据序列数据时间依赖关系,检测数据的非线性特征,并建立不同特征下非线性数据的预测模型。粒子群算法(Particle swarm optimization,PSO)可用来优化神经网络模型的超参数,GA算法可以优化神经网络模型的初始权值阈值,克服模型预测过程易陷入局部最优解及收敛速度慢的问题。因此一些混合智能算法的预测方法已经成为轨道不平顺预测的有效途径,通过多种智能算法的结合,往往可以得到更好的性能。
发明内容
本发明的目的是提供一种基于混合智能优化LSTM的轨道不平顺预测方法。解决了传统的预测方法对波动较大的轨道不平顺序列数据,预测精度不高的问题。
本发明所采用的技术方案是,
基于混合智能优化LSTM的轨道不平顺预测方法,包括以下步骤:
步骤1,收集轨道不平顺序列数据;
步骤2,获取轨道不平顺序列数据,将其转换为轨道质量指数时间序列;
步骤3,对步骤2获取的轨道质量指数时间序列数据进行预处理操作;
步骤4,根据步骤3所得数据,建立LSTM模型,同时利用PSO算法优化LSTM模型的超参数;
步骤5,根据步骤3所得数据和步骤4得到的最优超参数,建立LSTM模型,利用GA算法优化LSTM模型的权值阈值,得到LSTM-PSO-GA模型;
步骤6,利用LSTM-PSO-GA模型对未来数据进行预测。
本发明的特点还在于:
步骤3中,对序列数据采用归一化处理方法进行预处理,将原始序列数据映射到[0,1]具体方法是:先计算得序列数据的最大值和最小值,分别记为Xmax和Xmin;然后使用序列数据中的每个数据减去Xmin,再除以Xmax-Xmin。
步骤4中,LSTM模型包括输入层、隐藏层、输出层、网络训练以及网络预测、模块。输入层负责对原始轨道不平顺序列进行初步处理以满足网络输入要求,隐藏层采用LSTM细胞搭建单层循环神经网络,输出层提供预测结果网络,网络预测采用迭代的方法逐点预测。
步骤4中,利用PSO算法优化LSTM模型超参数的具体方法是:首先确定LSTM模型中待优化的超参数为时间窗口和隐藏层节点数;然后初始化参数,包括种群规模(粒子数量)M、最大迭代次数Nmax、、惯性权重w、粒子位置的边界Xmax和Xmin、粒子速度的范围Vmax和Vmin、加速因子c1和c2;随机产生一个粒子种群,每个个体包括两个超参数:时间窗口和隐藏层节点数。PSO算法迭代过程中利用全局最优个体和历史最佳个体调整当前个体的位置。PSO算法迭代结束后,获取LSTM模型的最佳超参数。
步骤5中,利用GA算法优化LSTM模型权值阈值的具体方法是:初始化种群,对待优化的权值阈值进行编码;迭代过程中若当前最大适应度值无显著变化或达到种群最大迭代次数,则停止优化;否则利用选择、交叉、变异操作,更新超参数染色体,并将新的超参数赋值给LSTM神经网络。
LSTM-PSO-GA模型包括输入层、隐藏层和输出层,其中将步骤3所得到的序列数据和步骤4中LSTM模型的预测结果作为输入层的输入,输出层为LSTM-PSO-GA模型的预测结果;隐藏层使用tanh作为激活函数。
本发明的有益效果是:
本发明不仅解决了传统的预测方法对波动较大的轨道不平顺序列数据,预测精度不高的问题。并且提出了利用PSO算法和GA算法分别优化神经网络超参数和权值阈值的时间序列数据计算方法,克服了模型预测过程中易陷入局部最优解、收敛速度慢且不稳定的问题。该预测方法能够提取轨道不平顺数据的特征变化,最终实现对轨道不平顺序列数据的高准确率预测并分析,更精确地预测轨道不平顺现象。
附图说明
图1为本发明基于混合智能优化LSTM的轨道不平顺预测方法的流程图;
图2为本发明基于混合智能优化LSTM的轨道不平顺预测方法中的GA优化LSTM权值阈值流程图;
图3为本发明基于混合智能优化LSTM的轨道不平顺预测方法中的单个LSTM细胞结构图;
图4为本发明基于混合智能优化LSTM的轨道不平顺预测方法中的轨道质量指数数据图;
图5为本发明基于混合智能优化LSTM的轨道不平顺预测方法中的轨道质量指数预测结果图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明基于混合智能优化LSTM的轨道不平顺预测方法,方法流程图如图1所示。包括以下步骤:
步骤1,收集轨道不平顺序列数据;
步骤2,获取轨道不平顺时间序列数据将其转换为轨道质量指数时间序列;不平顺时间序列数据包括:轨道质量指数、轨距标准差、高低标准差和水平标准差共四项数据;
步骤3,对步骤2获取的轨道不平顺序列数据进行预处理操作;
步骤4,根据步骤3所得轨道不平顺数据,建立LSTM网络模型,同时利用PSO算法优化该LSTM模型的超参数;
步骤5,根据步骤3所得轨道不平顺数据和步骤4所得超参数,建立LSTM网络模型,利用GA算法优化该LSTM模型的权值阈值,得到LSTM-PSO-GA模型;
步骤6,利用步骤5所得模型,利用LSTM-PSO-GA模型对未来轨道不平顺序列数据进行预测。
步骤3中,对轨道不平顺序列数据采用归一化处理方法进行预处理,将原始序列数据映射到[0,1]区间,具体方法是:先计算得序列数据的最大值和最小值,分别记为Xmax和Xmin;然后使用序列数据中的每个数据减去Xmin,再除以Xmax-Xmin。
步骤4中,用PSO算法对LSTM模型超参数进行优化的具体方法是:
步骤4.1,确定LSTM模型中待优化的超参数为时间窗口和隐藏层节点数。初始化参数:种群规模(粒子数量)M、最大迭代次数Nmax、惯性权重w、粒子位置的边界Xmax和Xmin、粒子速度的范围Vmax和Vmin、加速因子c1和c2;
步骤4.2,随机产生一个粒子种群,每个个体包括两个超参数:时间窗口和隐藏层节点数。
步骤4.3,采用PSO算法进行迭代,每次迭代中构建LSTM模型具体方法包括输入层、隐藏层、输出层、网络训练以及网络预测5个功能模块。输入层负责对原始轨道不平顺序列进行初步处理以满足网络输入要求,隐藏层搭建单层循环神经网络,输出层提供预测结果,网络预测模块采用迭代的方法逐点预测。首先在输入层中,定义已经归一化后的原始轨道不平顺序列为Fo={f1,f2,…,fn},则划分的训练集和测试集可以表示为Ftr={f1,f2,…,fm}和Fte={fm+1,fm+2,…,fn},满足约束条件m<n和m,n∈N。为了适应隐藏层输入的特点,应用数据分割的方法对Ftr进行处理,设分割长度为L,则分割后的模型为X={X1,X2,…,XL},Xp={fp,fp+1,…,fm-L+p-1},并有1≤p≤L;p,L∈N。对应的期望输出为Y={Y1,Y2,…,YL},YP={fp+1,fp+2,…,fm-L+p}。然后,将X输入隐藏层,隐藏层包含L个按前后时刻连接的同构LSTM细胞,X经过隐藏层后的输出可以表示为P={P1,P2,…,PL},Pp=LSTMforward(Xp,Cp-1,Hp-1),式中Cp-1和Hp-1分别为前一个LSTM细胞的状态和输出;LSTMforward表示LSTM向前细胞计算方法。设定细胞状态向量大小为Sstate,则Cp-1和Hp-1两个向量的大小均为Sstate。隐藏层输出P、模型输入X和理论输出Y均为维度为(m-L,L)的二维数组。选用均方误差作为误差计算公式,训练过程的损失函数可以定义为:设定损失函数最小为优化目标,不断更新网络权重,进而得到最终的隐藏层网络。
步骤4.4,PSO算法迭代过程中利用全局最优个体和历史最佳个体调整当前个体的位置。
步骤4.5,PSO算法迭代结束后,获取LSTM模型的最佳超参数。
步骤5中,利用GA算法优化该LSTM模型的权值阈值的流程如图2所示,具体方法是:
步骤5.1,根据步骤4所得最优超参数确定LSTM模型的时间窗口大小和隐藏层节点数;
步骤5.2,GA算法编码:对LSTM模型的初始权值阈值编码,模型误差作为各染色体的适应度值;
步骤5.3,产生指定数量的种群;
步骤5.4,迭代开始,利用选择、交叉、变异操作,更新超参数染色体,并将新的超参数赋值给LSTM模型;
步骤5.5,若当前最大适应度值达到种群最大迭代次数,则停止优化;
步骤5.6,计算当前的预测误差,保存当前最优超参数和对应LSTM网络模型,完成轨道不平顺时间序列数据的预测;
上述构建LSTM模型并对于现有的数据进行训练并预测;构建LSTM模型具体方法包括输入层、隐藏层、输出层、网络训练以及网络预测5个功能模块。输入层负责对原始响应时间序列进行初步处理以满足网络输入要求,隐藏层采用图3表示的LSTM细胞搭建单层循环神经网络,输出层提供预测结果网络,网络预测采用迭代的方法逐点预测。
LSTM-PSO-GA模型包括输入层、隐藏层和输出层,其中将步骤3所得到的序列数据和步骤4中LSTM模型的预测结果作为输入层的输入,输出层为LSTM-PSO-GA模型的预测结果;隐藏层使用tanh作为激活函数。
本实施例采用来源于京九线上行K449+300-K449+600区段1000天内检测数据。轨检车沿轨道里程方向每隔0.25米取一次轨道不平顺检测数据,在一次检测过程中可采集该区段内1200个检测点的各项指标数据,包含轨距、水平和高低等多项指标。将每个指标每次采集的1200个数据统计得到当日轨道质量指数数据(TQI),TQI序列如图4所示,LSTM-PSO-GA模型在TQI序列中的预测结果如图5所示,不同模型的误差对比如表1所示,我们分别采用均方根误差RMSE、平均绝对误差MAE以及平均绝对百分比误差MAPE作为评价指标,分别如公式(3)(4)(5)所示,其中RMSE为标准差,N为数据样本的个数,ypredictive为预测值,ytrue为实际值。
表1不同模型的预测误差对比
Claims (6)
1.基于混合智能优化LSTM的轨道不平顺预测方法,其特征在于,包括以下步骤:
步骤1,收集轨道不平顺序列数据;
步骤2,获取轨道不平顺序列数据,将其转换为轨道质量指数时间序列;
步骤3,对步骤2获取的轨道质量指数时间序列数据进行预处理操作;
步骤4,根据步骤3所得数据,建立LSTM模型,同时利用PSO算法优化LSTM模型的超参数;
步骤5,根据步骤3所得数据和步骤4得到的最优超参数,建立LSTM模型,利用GA算法优化LSTM模型的权值阈值,得到LSTM-PSO-GA模型;
步骤6,利用LSTM-PSO-GA模型对未来数据进行预测。
2.根据权利要求1所述的基于混合智能优化LSTM的轨道不平顺预测方法,其特征在于,所述步骤3中,对序列数据采用归一化处理方法进行预处理,将原始序列数据映射到[0,1]具体方法是:先计算得序列数据的最大值和最小值,分别记为Xmax和Xmin;然后使用序列数据中的每个数据减去Xmin,再除以Xmax-Xmin。
3.根据权利要求1所述的基于混合智能优化LSTM的轨道不平顺预测方法,其特征在于,所述步骤4中,所述LSTM模型包括输入层、隐藏层、输出层、网络训练以及网络预测、模块。输入层负责对原始轨道不平顺序列进行初步处理以满足网络输入要求,隐藏层采用LSTM细胞搭建单层循环神经网络,输出层提供预测结果网络,网络预测采用迭代的方法逐点预测。
4.根据权利要求1所述的基于混合智能优化LSTM的轨道不平顺预测方法,其特征在于,所述步骤4中,利用PSO算法优化LSTM模型超参数的具体方法是:首先确定LSTM模型中待优化的超参数为时间窗口和隐藏层节点数;然后初始化参数,包括种群规模(粒子数量)M、最大迭代次数Nmax、、惯性权重w、粒子位置的边界Xmax和Xmin、粒子速度的范围Vmax和Vmin、加速因子c1和c2;随机产生一个粒子种群,每个个体包括两个超参数:时间窗口和隐藏层节点数。PSO算法迭代过程中利用全局最优个体和历史最佳个体调整当前个体的位置。PSO算法迭代结束后,获取LSTM模型的最佳超参数。
5.根据权利要求1所述的基于混合智能优化LSTM的轨道不平顺预测方法,其特征在于,所述步骤5中,利用GA算法优化LSTM模型权值阈值的具体方法是:初始化种群,对待优化的权值阈值进行编码;迭代过程中若当前最大适应度值无显著变化或达到种群最大迭代次数,则停止优化;否则利用选择、交叉、变异操作,更新超参数染色体,并将新的超参数赋值给LSTM神经网络。
6.根据权利要求1所述的基基于混合智能优化LSTM的轨道不平顺预测方法,其特征在于,所述LSTM-PSO-GA模型包括输入层、隐藏层和输出层,其中将步骤3所得到的序列数据和步骤4中LSTM模型的预测结果作为输入层的输入,所述输出层为LSTM-PSO-GA模型的预测结果;所述隐藏层使用tanh作为激活函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210021594.9A CN114510778A (zh) | 2022-01-10 | 2022-01-10 | 基于混合智能优化lstm的轨道不平顺预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210021594.9A CN114510778A (zh) | 2022-01-10 | 2022-01-10 | 基于混合智能优化lstm的轨道不平顺预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114510778A true CN114510778A (zh) | 2022-05-17 |
Family
ID=81548907
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210021594.9A Pending CN114510778A (zh) | 2022-01-10 | 2022-01-10 | 基于混合智能优化lstm的轨道不平顺预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114510778A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116050814A (zh) * | 2023-04-03 | 2023-05-02 | 北京交通大学 | 一种有砟轨道捣固方案自适应修正方法 |
CN116307302A (zh) * | 2023-05-23 | 2023-06-23 | 西南交通大学 | 轨道不平顺动静态检测数据反演方法、***及存储介质 |
WO2024060441A1 (zh) * | 2022-09-21 | 2024-03-28 | 中建材创新科技研究院有限公司 | 一种纸面石膏板高精度切断的控制优化方法 |
-
2022
- 2022-01-10 CN CN202210021594.9A patent/CN114510778A/zh active Pending
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2024060441A1 (zh) * | 2022-09-21 | 2024-03-28 | 中建材创新科技研究院有限公司 | 一种纸面石膏板高精度切断的控制优化方法 |
CN116050814A (zh) * | 2023-04-03 | 2023-05-02 | 北京交通大学 | 一种有砟轨道捣固方案自适应修正方法 |
CN116050814B (zh) * | 2023-04-03 | 2023-08-04 | 北京交通大学 | 一种有砟轨道捣固方案自适应修正方法 |
CN116307302A (zh) * | 2023-05-23 | 2023-06-23 | 西南交通大学 | 轨道不平顺动静态检测数据反演方法、***及存储介质 |
CN116307302B (zh) * | 2023-05-23 | 2023-07-25 | 西南交通大学 | 轨道不平顺动静态检测数据反演方法、***及存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111915059B (zh) | 基于注意力机制的Seq2seq泊位占有率预测方法 | |
CN109492822B (zh) | 空气污染物浓度时空域关联预测方法 | |
CN114510778A (zh) | 基于混合智能优化lstm的轨道不平顺预测方法 | |
CN109142171B (zh) | 基于特征扩张的融合神经网络的城市pm10浓度预测方法 | |
CN108388962B (zh) | 一种风电功率预测***及方法 | |
CN112734097B (zh) | 无人驾驶列车能耗预测方法、***及存储介质 | |
CN112149230B (zh) | 一种强风铁路风致列车舒适性劣化预测方法 | |
CN112734094B (zh) | 一种智慧城市智轨车辆故障基因预测方法及*** | |
CN112966853A (zh) | 基于时空残差混合模型的城市路网短时交通流预测方法 | |
CN112149915B (zh) | 一种充电站充电负荷预测方法和*** | |
CN116757057A (zh) | 基于pso-ga-lstm模型的空气质量预测方法 | |
CN113674524A (zh) | 基于lstm-gasvr的多尺度短时交通流预测建模、预测方法及*** | |
CN110766190A (zh) | 一种配电网负荷预测方法 | |
CN116187835A (zh) | 一种基于数据驱动的台区理论线损区间估算方法及*** | |
CN116311921A (zh) | 一种基于多空间尺度时空Transformer的交通速度预测方法 | |
CN115907122A (zh) | 区域电动汽车充电负荷预测方法 | |
CN113536665B (zh) | 基于特征工程和lstm的路表温度短临预测方法及*** | |
CN112098869B (zh) | 一种基于大数据的自适应电动汽车soc估计方法 | |
CN116167465A (zh) | 基于多变量时间序列集成学习的太阳辐照度预测方法 | |
CN114463978B (zh) | 一种基于轨道交通信息处理终端的数据监测方法 | |
CN115907228A (zh) | 一种基于pso-lssvm的短期电力负荷预测分析方法 | |
CN112907971B (zh) | 基于遗传算法优化时空残差模型的城市路网短时交通流预测方法 | |
CN115056829A (zh) | 多车型连续学习的列车运动状态估计方法 | |
CN113947904A (zh) | 基于s-g滤波和深信度网络的多尺度短时交通流预测方法 | |
Wang et al. | Research on parking space prediction based on long short-term memory |
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 |