CN111498150B - 一种星载大惯量旋转载荷角动量和力矩补偿方法及系统 - Google Patents
一种星载大惯量旋转载荷角动量和力矩补偿方法及系统 Download PDFInfo
- Publication number
- CN111498150B CN111498150B CN202010229891.3A CN202010229891A CN111498150B CN 111498150 B CN111498150 B CN 111498150B CN 202010229891 A CN202010229891 A CN 202010229891A CN 111498150 B CN111498150 B CN 111498150B
- Authority
- CN
- China
- Prior art keywords
- balance wheel
- speed
- rotating speed
- real
- rotary table
- 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
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/28—Guiding or controlling apparatus, e.g. for attitude control using inertia or gyro effect
- B64G1/285—Guiding or controlling apparatus, e.g. for attitude control using inertia or gyro effect using momentum wheels
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B11/00—Automatic controllers
- G05B11/01—Automatic controllers electric
- G05B11/36—Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential
- G05B11/42—Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential for obtaining a characteristic which is both proportional and time-dependent, e.g. P. I., P. I. D.
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Radar, Positioning & Navigation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Testing Of Balance (AREA)
Abstract
本发明属于星载载荷应用技术领域,具体涉及一种星载大惯量旋转载荷角动量和力矩补偿方法及系统。解决现有补偿方式无法满足大惯量高转速的载荷补偿要求的问题,系统包括第一速度转换单元、平衡轮系统、第二速度转换单元及转台系统;方法过程如下:先将卫星给定的期望转速转换为平衡轮的目标转速,平衡轮开始向目标转速加速,利用平衡轮输出的霍尔信号实时计算平衡轮的转速,将实时平衡轮转速转换为转台转速给定,并反向旋转,从而保证在加速、匀速、减速各过程中转台转速均实时与平衡轮转速进行匹配和跟踪,保证残余角动量和残余力矩补偿效果。
Description
技术领域
本发明属于星载载荷应用技术领域,具体涉及一种星载大惯量旋转载荷角动量和力矩补偿方法及系统。
背景技术
卫星的有效载荷承担着整个卫星的任务与使命,卫星有效载荷涉及光学、无线电、微波等多个领域,承担着包括对地成像、通信、资源勘察,气象观测等多种任务。常见的有效载荷包括红外扫描仪、合成孔径雷达、可见光电荷耦合器件等,卫星的有效载荷因不同的任务而有所差异。
随着卫星技术的快速发展以及对卫星功能的需求,卫星载荷的种类呈现出多样性的趋势。部分卫星携带了大惯量旋转载荷,大惯量旋转载荷一般搭载在转台上,转台驱动大惯量旋转载荷旋转的过程中会有反作用力作用到卫星上。由于卫星工作在空间自由空间中,会使卫星产生相对运动,影响卫星姿态,严重时会影响卫星正常工作。
一般可以利用卫星的动量轮来补偿旋转载荷产生的角动量和力矩,但是动量轮的输出角动量和响应能力有限,在某些场合,无法满足大惯量高转速的载荷补偿要求。
发明内容
为了解决现有补偿方式无法满足大惯量高转速的载荷补偿要求的问题,本发明提供星载大惯量旋转载荷角动量和力矩补偿方法及系统。
本发明的技术解决方案是提供一种星载大惯量旋转载荷角动量和力矩补偿方法,包括以下步骤:
步骤1、引入平衡轮,所述平衡轮搭载小惯量旋转部件,将卫星给定的期望转台工作转速转换为平衡轮的目标转速ωP*;
步骤2、驱动平衡轮向目标转速加速;
步骤3、实时采集平衡轮位置信号并根据位置信号计算平衡轮的实时转速;
步骤4、将平衡轮的实时转速转换为搭载大惯量载荷的转台实时目标转速ωZ **;计算公式如下:
其中,JZ为转台的负载惯量;JP为平衡轮的负载惯量;ωP为步骤3获取的平衡轮的实时转速;
步骤5、驱动转台按照步骤4转换后的转台实时目标转速旋转且旋转方向与平衡轮旋转方向相反。
进一步地,为了保证运行速度的准确性,步骤3与步骤4之间还包括平衡轮速度闭环控制的步骤,即检测并判断平衡轮的实时转速,当平衡轮的实时转速不等于当前目标实时转速时,返回步骤2,通过控制使得平衡轮在当前目标实时转速下旋转。
步骤5之后还包括转台速度闭环控制的步骤,即检测并判断转台实时转速,当转台实时转速不等于步骤4转换后的当前转台实时目标转速时,返回步骤5,通过控制使得转台在转换后当前实时目标转速下旋转。
进一步地,步骤3中,平衡轮转速为1/Δt(rpm),其中Δt为平衡轮旋转一周,采集平衡轮位置信号所需要的时间。
本发明还提供一种星载大惯量旋转载荷角动量和力矩补偿系统,其特殊之处在于:包括第一速度转换单元、平衡轮系统、第二速度转换单元及转台系统;
上述第一速度转换单元用于将卫星给定的期望转台工作转速转换为平衡轮的目标转速,并反馈至平衡轮系统;
上述平衡轮系统包括平衡轮驱动电路、平衡轮单元及平衡轮速度计算单元;所述平衡轮单元包括平衡轮、用于驱动平衡轮旋转的高速电机及用于采集平衡轮速度信息的霍尔传感器;
上述平衡轮驱动电路的输入端与第一速度转换单元的输出端连接,平衡轮驱动电路的输出端接平衡轮单元,用于控制高速电机驱动平衡轮向目标转速加速;上述平衡轮速度计算单元用于根据霍尔传感器的输出计算平衡轮的实时转速并将速度数据反馈至第二速度转换单元;
上述第二速度转换单元用于将平衡轮的实时转速转化为转台实时目标转速ωZ **,并反馈至转台系统,计算公式如下:
其中,JZ为转台的负载惯量;JP为平衡轮的负载惯量;ωP为获取的平衡轮的实时转速;
上述转台系统包括转台驱动电路、转台单元及转台速度计算单元;上述转台单元包括转台、用于驱动转台旋转的力矩电机和用于实时测量转台位置数据的位置测量部件;
上述转台驱动电路的输入端与第二速度转换单元的输出端连接,转台驱动电路的输出端与转台单元连接,用于驱动力矩电机以与平衡轮相反的旋转方向驱动转台按照转换后的转台实时目标转速旋转;上述转台速度计算单元用于根据位置测量部件的输出计算转台速度。
进一步地,上述平衡轮系统还包括连接在第一速度转换单元与平衡轮驱动电路之间的第一PID控制器,用于实现平衡轮转速的闭环控制;当平衡轮的实时速度不等于当前目标实时转速时,通过控制使得平衡轮在当前目标实时转速下旋转。
上述转台系统还包括连接在第二速度转换单元与转台驱动电路之间的第二PID控制器,用于实现转台转速的闭环控制;当转台实时转速不等于转换后的转台实时目标转速时,通过控制使得转台在转换后的实时目标转速下旋转。保证了运行时速度的准确性和快速性。
进一步地,上述平衡轮驱动电路及转台驱动电路均包括三相桥电路。
本发明的有益效果是:
1、本发明,为了达到更好的补偿效果,先将卫星给定的期望转速转换为平衡轮的目标转速,平衡轮开始向目标转速加速,利用平衡轮输出的霍尔信号实时计算平衡轮的转速,将实时平衡轮转速转换为转台转速给定,并反向旋转,从而保证在加速、匀速、减速各过程中转台转速均实时与平衡轮转速进行匹配和跟踪,保证残余角动量和残余力矩补偿效果。
2、本发明转台和平衡轮均利用闭环控制很好的保证速度控制精度及响应速度;提高软件算法的控制频率,保证补偿的快速性,针对星载大惯量旋转载荷角动量和力矩有很好的补偿效果。
附图说明
图1为本发明星载大惯量旋转载荷角动量和力矩补偿系统框图;
具体实施方式
以下结合附图及具体实施例对本发明做进一步地描述。
本发明为了更好的补偿旋转载荷产生的角动量和力矩,利用平衡轮实现补偿。为了减小平衡轮的体积和重量,其搭载小惯量旋转部件,工作在高转速。转台搭载大惯量载荷工作在低转速,机动能力强,可以更好的在运行过程中精确快速的跟踪匹配平衡轮转速。
本发明先将卫星给定的期望转台工作转速转换为平衡轮的目标转速;具体的转化公式如下:其中JZ为转台负载惯量;JP为平衡轮负载惯量;ωZ*为卫星给定的期望转台工作转速。然后驱动平衡轮开始向目标转速加速;实时采集平衡轮位置信号并根据位置信号计算平衡轮的实时转速;将实时平衡轮转速转换为转台转速给定,其旋转方向与平衡轮旋转方向相反。从而保证在加速、匀速、减速各过程中转台转速均实时与平衡轮转速进行匹配和跟踪,保证残余角动量和残余力矩补偿效果。
具体将平衡轮转速转换为转台实时目标转速ωZ **的公式如下:
其中,JZ为转台的负载惯量;JP为平衡轮的负载惯量;ωP为获取的平衡轮的实时转速;提高计算频率,保证计算控制过程的快速性,进一步提高补偿效果。
对平衡轮及转台均进行速度闭环控制,针对平衡轮,即检测并判断平衡轮的实时速度,当平衡轮的实时速度不等于当前目标转速时,通过控制使得平衡轮在当前目标转速下旋转。对于转台,即检测并判断转台的实时速度,当转台的实时速度不等于转换后的转台实时目标转速时,通过控制使得转台在转换后的实时目标转速下旋转。从而保证运行时速度的准确性和快速性。
如图1所示,本实施例星载大惯量旋转载荷角动量和力矩补偿系统,包括第一速度转换单元、平衡轮系统、第二速度转换单元及转台系统。第一速度转换单元用于将卫星给定的期望转台工作转速转换为平衡轮的目标转速,并反馈至平衡轮系统;其中,转台转速至平衡轮转速的转换系数,其中,JZ为转台的负载惯量;JP为平衡轮的负载惯量;卫星给定的期望转台工作转速为ωZ*;因此,转化后平衡轮的目标转速为
平衡轮系统包括依次连接的第一PID控制器、平衡轮驱动电路、平衡轮单元及平衡轮速度计算单元;平衡轮速度计算单元的输出分别与第一PID控制器的输入及第二速度转换单元的输入连接。补偿用的平衡轮单元,其内部包含了平衡轮、高速电机及霍尔传感器,其中高速电机用于驱动平衡轮高速旋转,霍尔传感器用于采集平衡轮速度信息。第一PID控制器接收第一速度转换单元输出的平衡轮的目标转速,并通过平衡轮驱动电路控制高速电机驱动平衡轮向目标转速加速,平衡轮速度计算单元根据霍尔传感器的输出计算平衡轮的实时速度并将速度数据反馈至第一PID控制器及第二速度转换单元。第一PID控制器检测并判断平衡轮的实时速度,当平衡轮的实时速度不等于当前目标转速时,通过控制使得平衡轮在当前目标转速下旋转,保证控制的精确性。本实施例平衡轮速度计算单元根据如下过程计算平衡轮转速:平衡轮内部安装了6个霍尔传感器,高速电机每旋转360度电角度产生6个脉冲,电机共有8对极,因此平衡轮每旋转一周共旋转8*360度电角度,因此平衡轮每旋转一周产生48个脉冲信号,采集产生48个脉冲所需要的时间Δt(min),则平衡轮转速为1/Δt(rpm),进行单位换算为从而得到平衡轮实时转速。
第二速度转换单元将平衡轮的实时速度转化为转台实时目标转速ωZ **,并反馈至转台系统,计算公式如下:
其中,JZ为转台的负载惯量;JP为平衡轮的负载惯量;ωP为获取的平衡轮的实时转速;
转台系统包括依次连接的第二PID控制器、转台驱动电路、转台单元及转台速度计算单元。转台速度计算单元的输出与第二PID控制器的输入连接。转台单元包括搭载大惯量载荷的转台,用于实时测量转台位置数据的位置测量部件及用于驱动转台旋转的力矩电机。
第二PID控制器接收第二速度转换单元输出的转换后转台的目标转速,并控制通过转台驱动电路控制力矩电机驱动转台与平衡轮相反的旋转方向驱动转台按照转换后的转台实时目标转速旋转;保证转台和平衡轮旋转方向相反,转速大小满足上述的关系式,就可以很好的补偿运行过程中产生的角动量和力矩。转台速度计算单元根据位置测量部件的输出计算转台转速,并将转台的实时转速反馈至第二PID控制器,当转台的实时速度不等于转换后的转台实时目标转速时,通过控制使得转台在转换后的实时目标转速下旋转。保证控制精度,避免转台运行时转速存在误差,提高转台转速的响应能力。
Claims (6)
1.一种星载大惯量旋转载荷角动量和力矩补偿方法,其特征在于,包括以下步骤:
步骤1、引入平衡轮,所述平衡轮搭载小惯量旋转部件,将卫星给定的期望转台工作转速转换为平衡轮的目标转速ωP*;
步骤2、驱动平衡轮向目标转速加速;
步骤3、实时采集平衡轮位置信号并根据位置信号计算平衡轮的实时转速;
平衡轮速度闭环控制,检测并判断平衡轮的实时转速,当平衡轮的实时转速不等于当前目标实时转速时,返回步骤2,通过控制使得平衡轮在当前目标实时转速下旋转;
步骤4、将平衡轮的实时转速转换为搭载大惯量载荷的转台实时目标转速ωZ **;计算公式如下:
其中,JZ为转台的负载惯量;JP为平衡轮的负载惯量;ωP为步骤3获取的平衡轮的实时转速;
步骤5、驱动转台按照步骤4转换后的转台实时目标转速旋转且旋转方向与平衡轮旋转方向相反;
步骤6、转台速度闭环控制,检测并判断转台实时转速,当转台实时转速不等于步骤4转换后的当前转台实时目标转速时,返回步骤5,通过控制使得转台在转换后当前实时目标转速下旋转。
2.根据权利要求1所述的星载大惯量旋转载荷角动量和力矩补偿方法,其特征在于:步骤3中,平衡轮转速为1/Δt(rpm),其中Δt为平衡轮旋转一周,采集平衡轮位置信号所需要的时间。
3.一种星载大惯量旋转载荷角动量和力矩补偿系统,其特征在于:包括第一速度转换单元、平衡轮系统、第二速度转换单元及转台系统;
所述第一速度转换单元用于将卫星给定的期望转台工作转速转换为平衡轮的目标转速,并反馈至平衡轮系统;
所述平衡轮系统包括平衡轮驱动电路、平衡轮单元及平衡轮速度计算单元;所述平衡轮单元包括平衡轮、用于驱动平衡轮旋转的高速电机及用于采集平衡轮速度信息的霍尔传感器;
所述平衡轮驱动电路的输入端与第一速度转换单元的输出端连接,平衡轮驱动电路的输出端接平衡轮单元,用于控制高速电机驱动平衡轮向目标转速加速;所述平衡轮速度计算单元用于根据霍尔传感器的输出计算平衡轮的实时转速并将速度数据反馈至第二速度转换单元;
所述第二速度转换单元用于将平衡轮的实时转速转化为转台实时目标转速ωZ **,并反馈至转台系统,计算公式如下:
其中,JZ为转台的负载惯量;JP为平衡轮的负载惯量;ωP为获取的平衡轮的实时转速;
所述转台系统包括转台驱动电路、转台单元及转台速度计算单元;所述转台单元包括转台、用于驱动转台旋转的力矩电机和用于实时测量转台位置数据的位置测量部件;
所述转台驱动电路的输入端与第二速度转换单元的输出端连接,转台驱动电路的输出端与转台单元连接,用于驱动力矩电机以与平衡轮相反的旋转方向驱动转台按照转换后的转台实时目标转速旋转;所述转台速度计算单元用于根据位置测量部件的输出计算转台速度。
4.根据权利要求3所述的星载大惯量旋转载荷角动量和力矩补偿系统,其特征在于:所述平衡轮系统还包括连接在第一速度转换单元与平衡轮驱动电路之间的第一PID控制器,用于实现平衡轮转速的闭环控制;当平衡轮的实时速度不等于当前目标实时转速时,通过控制使得平衡轮在当前目标实时转速下旋转。
5.根据权利要求4所述的星载大惯量旋转载荷角动量和力矩补偿系统,其特征在于:所述转台系统还包括连接在第二速度转换单元与转台驱动电路之间的第二PID控制器,用于实现转台转速的闭环控制;当转台实时转速不等于转换后的转台实时目标转速时,通过控制使得转台在转换后的实时目标转速下旋转。
6.根据权利要求5所述的星载大惯量旋转载荷角动量和力矩补偿系统,其特征在于:所述平衡轮驱动电路及转台驱动电路均包括三相桥电路。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010229891.3A CN111498150B (zh) | 2020-03-27 | 2020-03-27 | 一种星载大惯量旋转载荷角动量和力矩补偿方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010229891.3A CN111498150B (zh) | 2020-03-27 | 2020-03-27 | 一种星载大惯量旋转载荷角动量和力矩补偿方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111498150A CN111498150A (zh) | 2020-08-07 |
CN111498150B true CN111498150B (zh) | 2021-10-15 |
Family
ID=71872636
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010229891.3A Active CN111498150B (zh) | 2020-03-27 | 2020-03-27 | 一种星载大惯量旋转载荷角动量和力矩补偿方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111498150B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113532427B (zh) * | 2021-07-14 | 2022-08-19 | 贵州航天林泉电机有限公司 | 一种基于位置规划的卫星转台路径规划方法 |
CN118567282B (zh) * | 2024-07-31 | 2024-10-11 | 中国科学院西安光学精密机械研究所 | 一种转台及高转速平衡轮间的全数字协同控制方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103092208A (zh) * | 2013-01-09 | 2013-05-08 | 哈尔滨工业大学 | 基于sgcmg和rw的航天器高精度快速姿态机动方法 |
CN106828980A (zh) * | 2016-11-21 | 2017-06-13 | 上海卫星工程研究所 | 星上有效载荷扫描机构大干扰力矩实时补偿方法 |
CN109131942A (zh) * | 2018-10-29 | 2019-01-04 | 上海微小卫星工程中心 | 一种用于对卫星的旋转体进行消旋控制的方法及相应消旋装置 |
CN110104217A (zh) * | 2019-03-29 | 2019-08-09 | 上海卫星工程研究所 | 卫星姿态控制与大角动量补偿复用飞轮的构形与控制方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2937954B1 (fr) * | 2008-10-31 | 2011-07-29 | Thales Sa | Procede et systeme de desaturation des roues d'inertie d'un engin spatial |
US9963248B2 (en) * | 2016-02-04 | 2018-05-08 | The Boeing Company | Spin stabilization of a spacecraft for an orbit maneuver |
US10175700B2 (en) * | 2017-03-22 | 2019-01-08 | The Boeing Company | Methods and apparatus to minimize command dynamics of a satellite |
US10239570B2 (en) * | 2017-05-23 | 2019-03-26 | Stephen J. Lesko | Device and method for performing tilt compensation by rotating arms |
CN110697086B (zh) * | 2019-10-11 | 2021-08-03 | 上海航天控制技术研究所 | 以单控制力矩陀螺和两飞轮实现卫星三轴稳定控制的方法 |
-
2020
- 2020-03-27 CN CN202010229891.3A patent/CN111498150B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103092208A (zh) * | 2013-01-09 | 2013-05-08 | 哈尔滨工业大学 | 基于sgcmg和rw的航天器高精度快速姿态机动方法 |
CN106828980A (zh) * | 2016-11-21 | 2017-06-13 | 上海卫星工程研究所 | 星上有效载荷扫描机构大干扰力矩实时补偿方法 |
CN109131942A (zh) * | 2018-10-29 | 2019-01-04 | 上海微小卫星工程中心 | 一种用于对卫星的旋转体进行消旋控制的方法及相应消旋装置 |
CN110104217A (zh) * | 2019-03-29 | 2019-08-09 | 上海卫星工程研究所 | 卫星姿态控制与大角动量补偿复用飞轮的构形与控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111498150A (zh) | 2020-08-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110174899B (zh) | 一种基于敏捷卫星的高精度成像姿态指向控制方法 | |
Russell et al. | Wind tunnel and hover performance test results for multicopter uas vehicles | |
CN111498150B (zh) | 一种星载大惯量旋转载荷角动量和力矩补偿方法及系统 | |
CN105007012A (zh) | 一种车载转台控制系统及方法 | |
CN104932510A (zh) | 一种针对飞轮低速摩擦的挠性航天器姿态控制系统与方法 | |
CN101872198A (zh) | 车载摄像稳定平台 | |
Sanchez et al. | Hovering flight improvement of a quad-rotor mini UAV using brushless DC motors | |
CN107147341A (zh) | 用于控制电机转动的矢量控制方法及控制系统及稳定云台 | |
CN107193292A (zh) | 机载天线高速伺服运动平台 | |
CN107357318A (zh) | 稳定云台转动的控制方法及控制系统及稳定云台 | |
CN114148536A (zh) | 一种用于光电吊舱的线位移补偿控制方法 | |
CN104201458B (zh) | 星载雷达对航天器平台扰动实时补偿解耦方法及补偿解耦系统 | |
CN114459479B (zh) | 旋转载体姿态、位置测量装置与方法 | |
CN113867423A (zh) | 一种伺服转台方位角度跟踪控制方法及装置 | |
US6860451B1 (en) | Spacecraft spin axis reorientation method | |
CN111547177B (zh) | 自平衡机器人控制方法、系统、自平衡机器人及介质 | |
CN107219863A (zh) | 稳定云台转动的控制方法及控制系统及稳定云台 | |
CN114966694B (zh) | 适合于fmcw信号的距离多普勒成像算法 | |
CN102799186B (zh) | 远程控制设备的驱动控制器 | |
US20220388692A1 (en) | Attitude determination and control system based on a quaternion kalman filter and method thereof | |
CN111891402B (zh) | 一种基于自主机动的火星探测对地天线指向恢复方法 | |
CN104765375B (zh) | 一种飞行控制器及其辅助飞行控制方法 | |
JP2001325245A (ja) | ハイブリッドフィルタ | |
RU2368547C1 (ru) | Способ управления положением солнечной батареи космического аппарата | |
JP4963429B2 (ja) | 追尾アンテナ装置、姿勢角推定方法およびそのプログラム |
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 |