CN102736636B - 跟踪系统中基于角度信息的前馈控制方法 - Google Patents

跟踪系统中基于角度信息的前馈控制方法 Download PDF

Info

Publication number
CN102736636B
CN102736636B CN201210192786.2A CN201210192786A CN102736636B CN 102736636 B CN102736636 B CN 102736636B CN 201210192786 A CN201210192786 A CN 201210192786A CN 102736636 B CN102736636 B CN 102736636B
Authority
CN
China
Prior art keywords
target
miss distance
signal
angle information
feed
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
Application number
CN201210192786.2A
Other languages
English (en)
Other versions
CN102736636A (zh
Inventor
唐涛
包启亮
刘儒贞
吴琼燕
李志俊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Optics and Electronics of CAS
Original Assignee
Institute of Optics and Electronics of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Institute of Optics and Electronics of CAS filed Critical Institute of Optics and Electronics of CAS
Priority to CN201210192786.2A priority Critical patent/CN102736636B/zh
Publication of CN102736636A publication Critical patent/CN102736636A/zh
Application granted granted Critical
Publication of CN102736636B publication Critical patent/CN102736636B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Feedback Control In General (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

本发明公开一种跟踪系统中基于角度信息的前馈控制方法,在控制对象上面安装有编码器和CCD探测器,CCD探测器接收编码器的角度信息和目标的轨迹信息并对角度信息和目标的轨迹信息做相减处理,提取目标脱靶量;跟踪控制器接收并将目标脱靶量生成并输出跟踪控制信号;第一相加单元对目标脱靶量和角度信息做加处理,提取出目标的轨迹信息;Kalman滤波器接收并从提取出的目标轨迹信号中获取目标的速度、加速度信号;前馈控制器将目标的速度、加速度信号生成并输出前馈控制信号;第二相加单元接收并对前馈控制信号和跟踪控制信号做相加处理,生成并输出驱动控制信号,用于驱动控制对象,实现对目标脱靶量的闭环校正。

Description

跟踪系统中基于角度信息的前馈控制方法
技术领域
本发明属于光电控制领域,具体涉及到一种基于角位置信息的速度、加速度前馈控制方法,可用于目标的高精度跟踪。
背景技术
在光电跟踪控制系统中,目标的高精度跟踪是非常重要的工作。美国林肯实验室设计的“火池”控制系统,利用角位置信息以及高精度的测距信息实现共轴跟踪,取得好的控制效果。在国内跟踪控制系统中,一般采用了基于角位置信息合成的速度前馈(预测滤波技术在光电经纬仪中的应用仿真,光电工程,Vol(8),2002)也取得了一定结果。但是该方法只采用了速度前馈,没有对角加速度补偿,其提高精度有限。此外,对控制系统的稳定性也没有详细的阐述。
发明内容
为了解决现有技术没有对角加速度补偿的技术方案,及提高系统精度不显著的技术缺陷,本发明的目的是提供一种跟踪系统中基于角位置信息的前馈控制方法。
为实现本发明的目的,本发明提供跟踪系统中基于角位置信息的前馈控制方法的技术解决方案如下步骤所述:
步骤S1:在控制对象上面安装有编码器和CCD探测器,CCD探测器接收编码器的角度信息和目标的轨迹信息,CCD探测器对编码器的角度信息和目标的轨迹信息做相减处理,提取目标脱靶量e(k);
步骤S2:跟踪控制器接收目标脱靶量e(k),并将目标脱靶量e(k)生成并输出跟踪控制信号;
步骤S3:第一相加单元对目标脱靶量和编码器角度信息做加处理,提取出目标轨迹信息;
步骤S4:Kalman滤波器接收提取出的目标轨迹信号,并利用Kalman滤波器从提取出的目标轨迹信号中获取目标的速度、加速度信号;
步骤S5:前馈控制器将目标的速度、加速度信号生成并输出前馈控制信号;
步骤S6:第二相加单元接收前馈控制信号和跟踪控制信号,并对前馈控制信号和跟踪控制信号做相加处理,生成并输出驱动控制信号,用于驱动控制对象,实现对目标脱靶量e(k)的闭环校正。
本发明相比当前控制方法的优点有:仅仅利用CCD探测器、编码器提供的角位置信息,实现速度、加速度前馈控制技术,提高光电跟踪系统的跟踪能力。可以有效地减小速度、加速度滞后误差;)由于只利用了角度信息,不需要距离信息,系统简单;采用了经典的Kalman滤波算法,不需要增加系统负担,稳定可靠。
附图说明
图1为本发明实现目标跟踪的控制装置示意图。
图2为本发明图1的频率域传递函数原理图。
图3为采用本发明方法的控制误差曲线。
具体实施方式
以下结合附图和具体实施方式说明本发明,本领域的技术人员可根据本说明书揭示的内容了解本发明的功效及优点。
图1示出一种实现目标跟踪的控制装置包括:跟踪控制器、前馈控制器、滤波器、控制对象;该控制装置组包含了由跟踪控制器形成的反馈环路以及前馈控制器形成的前馈环路。采用所述装置实现前馈(速度、加速度前馈)控制方法的步骤如下:
步骤S1:在控制对象上面安装有编码器和CCD探测器,CCD探测器接收编码器的角度信息和目标的轨迹信息,CCD探测器对编码器的角度信息和目标的轨迹信息做相减处理,提取目标脱靶量e(k);所述目标脱靶量e(k)是目标轨迹与控制对象同步运动编码器的角度信号之差。
步骤S2:跟踪控制器接收目标脱靶量e(k),并将目标脱靶量e(k)生成并输出跟踪控制信号;
步骤S3:第一相加单元对目标脱靶量和编码器角度信息做加处理,提取出目标轨迹信息;值得注意是,由于当前目标脱靶量一般滞后角位置信息。因此,需要用当前目标脱靶量信息与前面时刻的角度信息相加,获取目标当前的轨迹信息;
步骤S4:Kalman滤波器接收提取出的目标轨迹信号,并利用Kalman滤波器从提取出的目标轨迹信号中获取目标的速度、加速度信号,所述Kalman滤波器是基于一个四阶牛顿运动方程建立的Kalman滤波器;
步骤S5:前馈控制器将目标的速度、加速度信号生成并输出前馈控制信号(包括速度前馈控制信号、加速度前馈控制信号);
步骤S6:第二相加单元接收前馈控制信号和跟踪控制信号,并对前馈控制信号和跟踪控制信号做相加处理,生成并输出驱动控制信号,用于驱动控制对象,实现对目标脱靶量e(k)的闭环校正。
本发明提供基于角位置信息的前馈控制方法的技术解决方案中:
1、加速度、速度的前馈控制器与跟踪控制器的关系:
加速度前馈、速度前馈同时加在一点(跟踪控制器的输出),这就表明这前馈控制器的更新频率和跟踪控制器的频率同步。
2、速度、加速度获取:
利用CCD探测器、编码器提供的角度信息,建立一个四阶运动方程,通过Kalman滤波技术得到运动目标的角速度、角加速度。值得注意的是CCD探测器与编码器的角度信息在时刻上对齐。
3、前馈控制器的参数:
速度、加速度前馈控制器是由Kalman滤波器得到的目标速度、加速度分别乘以系数实现的。
前馈控制器是由Kalman滤波器的滤波得到目标的速度、加速度的线性函数,速度前馈的系数取1,加速度前馈的系数稍大于目标脱靶量滞后时间T。
为了说明本发明的功效,图2示出本发明的图1的原理图,下面从传递函数来分析跟踪系统的闭环性能,其中各个变量的定义如下:
e-Ts:CCD探测器特性的传递函数;
延迟因子;
Gff(s):前馈控制器的传递函数;
D(s):Kalman滤波器的传递函数;
Gc(s):跟踪控制器的传递函数;
Gp(s):控制对象的传递函数;
其中,s为拉普拉斯变换变换算子;T1为当前时刻编码器信号比理想时刻编码器角度信号延迟时间;Tff为等效带宽常数;f、c、p为上述变量的脚标记且无含义。
图2的信号流程为:由安装在控制对象上的CCD探测器实现目标的轨迹信号与控制对象上的编码器输出角度信号相减得到目标脱靶量的信号。目标脱靶量信号作为跟踪控制器(其传递函数为Gc(s))的输入。同时,利用目标脱靶量信号与编码器的角度信号相加得到目标的轨迹信号,该目标的轨迹信号流向Kalman滤波器(其传递函数为D(s))产生目标的速度、加速度信号,将它们作为前馈控制器(其传递函数为Gff(s))的输入。将跟踪控制器的输出以及前馈控制器(其传递函数为Gff(s))的输出相加得到驱动控制信号,并送至给控制对象,实现闭环过程。
原理分析:从原理图2可以得到该跟踪控制系统的传递函数Gclose(s)为:
G close ( s ) = e - Ts [ G c ( s ) + D ( s ) G ff ( s ) ] G p ( s ) 1 + G c ( s ) G p ( s ) e - Ts + ( e - Ts - e - T 1 s ) G ff ( s ) D ( s ) G p ( s ) - - - ( 1 )
跟踪控制系统的误差传递函数Error(s)为:
Error ( s ) = [ 1 - e - T 1 s D ( s ) G ff ( s ) G p ( s ) ] 1 + G c ( s ) G p ( s ) e - Ts + ( e - Ts - e - T 1 s ) D ( s ) G ff ( s ) G p ( s ) - - - ( 2 )
前馈控制器Gff(s)由下等式给出,系统性能最优
G ff ( s ) = G p - 1 ( s ) e T 1 s D - 1 ( s ) - - - ( 3 )
很显然,由公式3可以知道:要实现对误差的全补偿,尤其控制对象的补偿是非常困难的。由于目标脱靶量滞后时间T<<1,于是可以认为 e Ts = 1 + Ts + 1 2 T 2 s 2 + 1 6 T 3 s 3 + · · · ≈ 1 + Ts , 控制对象的传递函数于是前馈控制传递函数Gff(s)设计为:
Gff(s)=s+(T1+Tff)s2            (4)
等式(4)表明了该前馈控制传递函数同时包含了速度、加速度信息。速度前馈系数为1,而加速度前馈系数为(T1+Tff),实际上就是稍大于目标脱靶量滞后时间T。
采用低通滤波器描述Kalman滤波器的特性,Tff就是等效带宽常数。由公式(2)可以知道,前馈控制对跟踪控制系统性能的改善如公式(5)描述:
1 - e - T 1 s G ff ( s ) G p ( s ) = 1 - s + ( T 1 + T ff ) s 2 1 + T ff s 1 s ( 1 + T 1 s ) = T 1 T ff s 2 ( 1 + T 1 s ) ( 1 + T ff s ) - - - ( 5 )
从误差传递函数Error(s)上看,速度、加速度前馈给闭环单元增加了两次微分环节s2,减小加速度滞后误差。如果只有速度前馈,即Gff(s)=s,对闭环单元改善性能如公式(6)给出,
1 - e - T 1 s G ff ( s ) G p ( s ) = 1 - s 1 + T ff s 1 s ( 1 + T 1 s ) = T 1 T ff s + ( T 1 + T ff ) ( 1 + T 1 s ) ( 1 + T ff s ) s - - - ( 6 )
很明显,从误差传递函数上看,速度给闭环单元只增加了一次微分环节s,无法消除加速度滞后误差。因此速度、加速度前馈可以有效地提高跟踪精度。
由等式(2)可以知道,如果T≠T1,闭环单元的稳定性可能受到破坏。这也表明目标脱靶量和编码器的角度信息在时刻上要对齐,否则前馈控制器可能无法实施。
运用实例
假设某一飞行目标的轨迹如θ=atan(0.9t-10.0)描述。θ是轨迹的角度信号,t是时间信号。跟踪控制系统的工作频率为100Hz,延迟时间T=0.03。采用上面的方法,由图3可以知道最大误差为小于2′。设计方法如下,某一实际跟踪控制对象的传递函数Gp(s)为:
G p ( s ) = 1 s ( s + 160 ) × G CCD ( s ) - - - ( 7 )
其中,Gccd(s)=e-Ts表示CCD探测器的传递函数特性,对于工作频率为100Hz的CCD探测器来说,常数T约等于0.03.
在保证相位裕量不小于45°、幅值裕量不小于6dB的情况下,可以得到比例-积分的跟踪控制器Gc(s)如下,
G c ( s ) = 18.34 + 77.62 s s - - - ( 8 )
采用Tusin变换可以得到跟踪控制器Gc(s)的离散方程为 G c ( Z ) = 18.321 × z - 0.9583 z - 1 , Z为时域标识;
图2表明要得到目标的位置、速度、加速度必须通过目标脱靶量和编码器角度信号进行合成,提取出的目标轨迹信息r(k)的位置可以表示为:
r(k)=e(k)+y(k-T)                        (9)
式中e(k)为第k时刻目标脱靶量信号,y(k-T)为第k-T时刻的编码器角度信号。一般来说,目标脱靶量信号滞后编码器角度信号为T时刻,所以目标的轨迹应该由当前的目标脱靶量与前T时刻的编码器角度信号合成得到。目标的速度、加速度信息可以由提取出的目标轨迹信息r(k)微分得到,但是这样会带来噪声。为了提高精度,采用Kalman滤波求解速度、加速度信息。标准的Kalman滤波器如下描述:
X ( k ) = Φ ( k , k - 1 ) X ( k - 1 ) + Γ ( k , k - 1 ) W ( k - 1 ) Y ( k ) = C ( k ) X ( k ) + V ( k ) - - - ( 10 )
X(k)、X(k-1)分别为Kalman滤波器k、k-1时刻的n维状态向量,Φ(k,k-1)是Kalman滤波器的n×n维状态转移矩阵,V(k)是m×1维观测噪声序列,Γ(k,k-1)是n×p维噪声输入矩阵,Ck是m×n维系数输入矩阵,Y(k)是Kalman滤波器k时刻的m×1维观测序列,W(k-1)是p×1维k-1时刻过程噪声序列。其中,n、p、m表示矩阵维数,没有实际的含义。
采用一个四阶牛顿运动方程来建立Kalman滤波器模型,于是有,
Φ ( k , k - 1 ) = 1 T T 2 / 2 T 3 / 6 0 1 T T 2 / 2 0 0 1 T 0 0 0 1 , Γ ( k , k - 1 ) = T 3 / 6 T 2 / 2 T 1 , C ( k ) = 1 0 0 0 ′ - - - ( 11 )
将提取出的目标轨迹信息r(k)作为观测量。式中′是状态向量的转置标识;通过标准的Kalman滤波求解公式可以得到目标的速度Xk(2),加速度信号Xk(3)。其中,Xk(2)、Xk(3)表示估计状态X(k)第2,3个单元的数据。前馈控制器可以设计为:
Gff(k)=Xk(2)+0.035×Xk(3)            (12)
将跟踪控制器的输出、以及前馈控制器的输出同时加到图1所示的第二相加单元输出驱动控制信号,就实现了本发明的方法。
以上所述,仅为本发明中的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉该技术的人在本发明所揭露的技术范围内,可理解想到的变换或替换,都应涵盖在本发明的包含范围之内。

Claims (3)

1.一种跟踪系统中基于角度信息的前馈控制方法,其特征在于:
步骤S1:在控制对象上面安装有编码器和CCD探测器,CCD探测器接收编码器的角度信息和目标的轨迹信息,CCD探测器对编码器的角度信息和目标的轨迹信息做相减处理,提取目标脱靶量,目标脱靶量和编码器的角度信息在时刻上要对齐;由跟踪控制系统的误差传递函数为
Error ( s ) = [ 1 - e - T 1 s D ( s ) G ff ( s ) G p ( s ) ] 1 + G c ( s ) G p ( s ) e - Ts + ( e - Ts - e - T 1 s ) D ( s ) G ff ( s ) G p ( s )
可以知道,T≠T1闭环单元的稳定性受到破坏;表明目标脱靶量和编码器的角度信息在时刻上要对齐,否则前馈控制器无法实施;式中:为延迟因子;D(s)为Kalman滤波器的传递函数;Gff(s)为前馈控制器的传递函数;Gp(s)为控制对象的传递函数;Gc(s)为跟踪控制器的传递函数;e-Ts为CCD探测器特性的传递函数;T为目标脱靶量滞后时间;T1为当前时刻编码器信号比理想时刻编码器角度信号延迟时间;
步骤S2:跟踪控制器接收目标脱靶量,并将目标脱靶量生成并输出跟踪控制信号;
步骤S3:第一相加单元对目标脱靶量和编码器角度信息做加处理,提取出目标轨迹信息;
步骤S4:Kalman滤波器接收提取出的目标轨迹信号,并利用Kalman滤波器从提取出的目标轨迹信号中获取目标的速度、加速度信号;
步骤S5:前馈控制器将目标的速度、加速度信号生成并输出前馈控制信号;
步骤S6:第二相加单元接收前馈控制信号和跟踪控制信号,并对前馈控制信号和跟踪控制信号做相加处理,生成并输出驱动控制信号,用于驱动控制对象,实现对目标脱靶量的闭环校正。
2.如权利要求1所述跟踪系统中基于角度信息的前馈控制方法,其特征在于:所述Kalman滤波器是基于一个四阶牛顿运动方程建立的Kalman滤波器。
3.如权利要求1所述跟踪系统中基于角度信息的前馈控制方法,其特征在于:前馈控制器是由Kalman滤波器的滤波得到目标的速度、加速度的线性函数,速度前馈的系数取1,加速度前馈的系数稍大于目标脱靶量滞后时间。
CN201210192786.2A 2012-06-12 2012-06-12 跟踪系统中基于角度信息的前馈控制方法 Active CN102736636B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210192786.2A CN102736636B (zh) 2012-06-12 2012-06-12 跟踪系统中基于角度信息的前馈控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210192786.2A CN102736636B (zh) 2012-06-12 2012-06-12 跟踪系统中基于角度信息的前馈控制方法

Publications (2)

Publication Number Publication Date
CN102736636A CN102736636A (zh) 2012-10-17
CN102736636B true CN102736636B (zh) 2014-12-10

Family

ID=46992293

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210192786.2A Active CN102736636B (zh) 2012-06-12 2012-06-12 跟踪系统中基于角度信息的前馈控制方法

Country Status (1)

Country Link
CN (1) CN102736636B (zh)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103576677B (zh) * 2013-11-19 2015-12-02 中国科学院光电技术研究所 光电跟踪控制系统中跟踪误差的估计方法
CN104331902B (zh) * 2014-10-11 2018-10-16 深圳超多维科技有限公司 目标跟踪方法、跟踪装置和3d显示方法及显示装置
CN104535078B (zh) * 2014-12-31 2017-03-29 中国人民解放军63655部队 一种基于标志点的光电设备对飞行目标的测量方法
CN106403991A (zh) * 2015-07-31 2017-02-15 北京航天计量测试技术研究所 一种基于脱靶量偏置的多光路安装误差补偿控制方法
CN105955027B (zh) * 2016-05-30 2019-05-14 中国科学院光电技术研究所 一种基于多阶运动信息估计的前馈控制方法
CN106647257B (zh) * 2016-10-14 2020-01-03 中国科学院光电技术研究所 一种基于正交最小二乘的前馈控制方法
CN106444386B (zh) * 2016-11-29 2019-03-12 凯迈(洛阳)测控有限公司 一种光电设备跟踪脱靶量时滞补偿控制方法
CN106647725B (zh) * 2017-02-22 2019-01-11 北京爱科迪通信技术股份有限公司 一种位置滞后补偿装置和方法
CN109708667B (zh) * 2018-12-26 2021-04-13 中国久远高新技术装备公司 一种基于激光陀螺的双动态目标跟踪引导方法
CN111791241B (zh) * 2020-06-16 2022-02-08 浙江华消科技有限公司 机器人轨迹误差校正方法及装置、消防机器人
CN112684817B (zh) * 2020-12-17 2023-01-03 中国工程物理研究院应用电子学研究所 一种提高光电跟踪系统跟踪精度的方法
CN114241009B (zh) * 2021-12-24 2023-06-02 普宙科技有限公司 一种目标跟踪方法、系统、存储介质及电子设备

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101770034A (zh) * 2010-01-29 2010-07-07 中国科学院空间科学与应用研究中心 机群链路的前馈/反馈组合式载波跟踪方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02130608A (ja) * 1988-11-11 1990-05-18 Nkk Corp サーボ機構の追従制御方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101770034A (zh) * 2010-01-29 2010-07-07 中国科学院空间科学与应用研究中心 机群链路的前馈/反馈组合式载波跟踪方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
JP平2-130608A 1990.05.18 *
用卡尔曼滤波实现单轴跟踪系统的前馈复合控制;罗续成;《遥测遥控》;20050331;第26卷(第2期);1-8 *
罗续成.用卡尔曼滤波实现单轴跟踪系统的前馈复合控制.《遥测遥控》.2005,第26卷(第2期),1-8. *
陈娟等.脱靶量前馈补偿光电跟踪伺服系统的跟踪精度.《长春工业大学学报(自然科学版)》.2007,第28卷(第4期),391-393. *

Also Published As

Publication number Publication date
CN102736636A (zh) 2012-10-17

Similar Documents

Publication Publication Date Title
CN102736636B (zh) 跟踪系统中基于角度信息的前馈控制方法
Shtessel et al. Smooth second-order sliding modes: Missile guidance application
CN102280572B (zh) 压电陶瓷执行器迟滞特性的复合线性化控制方法及其实现电路
CN101708780B (zh) 用于目标姿态跟踪的刚性航天器的控制方法
CN103760900B (zh) 一种考虑控制输入约束的船舶运动控制系统
CN107168369B (zh) 一种欠驱动无人车的编队控制方法及系统
CN103831831A (zh) 具有时变时延的非线性遥操作系统位置和力跟踪控制系统
CN101571705B (zh) 位置伺服系统与方法
CN103425131B (zh) 基于非光滑控制和扰动观测的农用拖拉机导航控制方法
CN107168072B (zh) 一种基于干扰观测器的非匹配干扰系统自抗扰控制方法
CN103051274A (zh) 基于变阻尼的二自由度永磁同步电机的无源性控制方法
CN104199301A (zh) 基于改进型自抗扰控制器的直线电机轨迹跟踪装置及方法
CN103389650A (zh) 四通道遥操作系统的双边未知死区自适应控制系统
CN104777747A (zh) 基于平方根容积卡尔曼滤波的光电跟踪转台滑模控制方法
CN103529858A (zh) 位置闭环系统最小幅相差跟踪法
CN105045103A (zh) 一种基于LuGre摩擦模型伺服机械手摩擦补偿控制系统及方法
CN103576677B (zh) 光电跟踪控制系统中跟踪误差的估计方法
CN102882447A (zh) 一种光电跟踪系统水平轴双电机装置及同步驱动跟踪控制方法
Hu et al. High-order sliding-mode observer based output feedback adaptive robust control of a launching platform with backstepping
CN108107738A (zh) 变采样率非线性驱动惯性稳定跟踪控制系统及方法
CN108333944A (zh) 基于框架角自适应调整的cmg操纵方法及系统
CN102427341B (zh) 基于卡尔曼滤波的远程迭代学习控制系统传输噪声抑制方法
CN208384428U (zh) 基于反向位移补偿的稳定转台系统
CN104369875B (zh) 基于非线性轨道计算的航天器制导控制方法及系统
CN102790584B (zh) 一种伺服驱动系统及其速度环参数整定系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant