CN104319909B - 基于最小二乘法的感应电能传输系统负载识别方法 - Google Patents

基于最小二乘法的感应电能传输系统负载识别方法 Download PDF

Info

Publication number
CN104319909B
CN104319909B CN201410626387.1A CN201410626387A CN104319909B CN 104319909 B CN104319909 B CN 104319909B CN 201410626387 A CN201410626387 A CN 201410626387A CN 104319909 B CN104319909 B CN 104319909B
Authority
CN
China
Prior art keywords
load
primary
value
process parameter
parameter vector
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
CN201410626387.1A
Other languages
English (en)
Other versions
CN104319909A (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.)
Chongqing Huachuang Intelligent Technology Research Institute Co ltd
Wang Zhihui
Original Assignee
Chongqing University
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 Chongqing University filed Critical Chongqing University
Priority to CN201410626387.1A priority Critical patent/CN104319909B/zh
Publication of CN104319909A publication Critical patent/CN104319909A/zh
Application granted granted Critical
Publication of CN104319909B publication Critical patent/CN104319909B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/30Circuit design
    • G06F30/36Circuit design at the analogue level

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Inverter Devices (AREA)

Abstract

为了建立高效、可靠、稳定的感应电能传输系统,并且实现系统最高效率跟踪控制,本发明公开一种基于最小二乘法的感应电能传输系统负载识别方法,利用最小二乘法对系统过程参数进行辨识,通过过程参数向量及负载间的线性关系完成对负载识别,将直接对负载识别问题转换为了对系统参数辨识问题,算法简单,计算方便,对特征参数的采样要求较低,同时能够满足实际应用要求。

Description

基于最小二乘法的感应电能传输系统负载识别方法
技术领域
本发明涉及感应电能传输(Inductive power transfer)技术,简称IPT技术,尤其涉及一种基于最小二乘法的感应电能传输系统负载识别方法。
背景技术
感应电能传输技术是一种借助高频电磁场将电能从电源端耦合到负载端的电能传输新技术。IPT技术理论及其关键技术的深入研究,推动了IPT技术在电动车充/供电、电子产品充/供电、生物医电以及照明系统等诸多领域的广泛应用。
由于系统负载的变化,源自能量接收端(副边电路)的反射阻抗会在能量发射端(原边电路)有着相应的变化,使得原边电路的固有频率发生漂移,从而不再与工作频率匹配。这将导致系统偏离软开关工作点,从而影响系统的功率传输能力,增加开关损耗以及EMI。此外,当负载变化,系统需根据当前负载情况调整输出能量模式,实现系统的最高效率跟踪控制。
现有技术中,文献[1-3]基于系统的能量模型,从能量守恒角度建立方程,完成对负载的辨识,此类方法仅需采集较少状态变量便可完成负载辨识。文献[4]在初始时刻给系统注入能量,随后使其自由振荡,通过推导出原边电流衰减率与负载之间的关系完成负载辨识,但不是一种稳态条件下的负载辨识策略。
参考文献:
[1]Dai Xin,Sun Yue,Tang Chunsen,et al.Dynamic parametersidentification method for inductively coupled power transfer system[C].Sustainable Energy Technologies(ICSET),2010IEEE International Conference on,Kandy,2010:1-5.
[2]戴欣,王智慧,唐春森,等.感应电能传输系统参数辨识与恒流控制[J].重庆大学学报,2011,(4).
[3]Zhi-hui Wang,Xiao Lv,Yue Sun,et al.A simple approach for loadidentification in current-fed inductive power transfer system[C].Power SystemTechnology(POWERCON),2012IEEE International Conference on,Auckland,2012:1-5.
[4]Zhi-hui Wang,Yu-Peng Li,Yue Sun,et al.Load Detection Model ofVoltage-Fed Inductive Power Transfer System[J].Power Electronics,IEEETransactions on,2013,28(11):5233-5243.
发明内容
为了建立高效、可靠、稳定的感应电能传输系统,并且实现系统最高效率跟踪控制,本发明提供一种基于最小二乘法的感应电能传输系统负载识别方法,该方法利用最小二乘法完成对系统过程参数的辨识,通过过程参数向量及负载间的线性关系完成对负载的辨识,将直接对负载辨识的问题转换为了对系统参数辨识的问题,降低负载辨识时的采样要求。
为达到上述目的,本发明所采用的具体技术方案如下:
一种基于最小二乘法的感应电能传输系统负载识别方法,其关键在于按照以下步骤进行:
步骤1:根据感应电能传输系统的拓扑结构构建系统的等效电路模型,该等效电路模型包括原边逆变器输出的方波电流源iac,原边谐振电容Cp以及原边谐振电容电压up,原边线圈Lp以及原边线圈电流ip,副边谐振电容Cs以及副边谐振电容电压us,副边线圈Ls以及副边线圈电流is,副边等效电阻R1
步骤2:选取原边逆变器输出的方波电流iac以及原边线圈电流ip作为状态变量,构建系统状态差分方程;
步骤3:根据系统状态差分方程得到系统的数据矩阵和过程参数向量;
步骤4:利用最小二乘法对过程参数向量进行估计;
步骤5:利用过程参数向量与系统参数之间的线性关系进行分解变换,得到系统的负载信息;
所述系统参数包括原边谐振电容值Cp、原边线圈电感值Lp、副边谐振电容电容值Cs、副边线圈电感值Ls、副边等效电阻值R1、原边线圈Lp与副边线圈Ls之间的互感值M、采样周期T。
针对PS结构的电流型IPT系统而言,感应电能传输系统的原边采用并联谐振回路,副边采用串联谐振回路,则步骤2所构建的系统状态差分方程为:
其中的系数c=2/T,Ip(k)表示原边线圈电流ip在第k个采样点的值,Iac(k)表示原边逆变器输出的方波电流iac在第k个采样点的值。
根据步骤2所构建的系统状态差分方程得到的系统的数据矩阵:
过程参数向量:且:
ΖL=HLθ=[Ip(k) Ip(k+1) … Ip(k+L)]T,L表示需要构建的状态方程的个数。
结合最小二乘算法,步骤4中按照对过程参数向量θ进行最小二乘估计,其中ΗT L表示数据矩阵ΗL的转置。
为了便于运算,数据矩阵HL为9×9的方阵,L的取值为9。
结合统计分布,利用过程参数向量与系统参数之间的线性关系进行分解变换,可以得到9个副边等效电阻辨识值,求其期望值作为所述副边等效电阻R1
本发明的显著效果是:利用最小二乘法对系统过程参数进行辨识,通过过程参数向量及负载间的线性关系完成对负载识别,将直接对负载识别问题转换为了对系统参数辨识问题,算法简单,计算方便,对特征参数的采样要求较低,能够满足实际应用要求。
附图说明
图1是具体实施例中PS结构的电流型IPT系统电路原理图;
图2是图1的等效电路模型;
图3是原边谐振电流及逆变电流的仿真波形图;
图4是原边谐振电流及逆变电流的实测波形图;
图5是具体实施例中负载识别结果的分布直方图;
图6是本发明的实用数据分析图。
具体实施方式
下面结合附图对本发明的具体实施方式以及工作原理作进一步详细说明。
如图1所示,PS结构的电流型IPT系统,该系统由原边电路以及副边电路两个独立部分所组成。在原边部分,直流电压源Edc提供整个系统的电能输入,Ldc为滤波电感,Edc与Ldc串联可近似视为电流源并在稳态时产生近似恒定电流idc。开关管S1~S4构成原边高频逆变电路,两组开关管S1,S4及S2,S3轮流导通将DC激励电流idc逆变成为高频方波电流iac,原边电容Cp、原边线圈Lp构成原边并联谐振电路;在副边部分,副边电容Cs、副边线圈Ls构成副边串联谐振电路,开关管D1~D4以及滤波电容C构成副边整流电路,将高频交流电压转为直流电压uo作用于负载RL上。M为耦合电感Lp、Ls之间互感,Rp、Rs分别为电感Lp、Ls串联等效电阻。
针对这样的IPT系统而言,其负载识别按照以下方式进行:
步骤1:根据感应电能传输系统的拓扑结构构建系统的等效电路模型,该等效电路模型包括原边逆变器输出的方波电流源iac,原边谐振电容Cp以及原边谐振电容电压up,原边线圈Lp以及原边线圈电流ip,副边谐振电容Cs以及副边谐振电容电压us,副边线圈Ls以及副边线圈电流is,副边等效电阻R1
如图1所示的电路结构,对于原边电路,电压源Edc产生直流电压流经滤波电感Ldc,通过高频逆变电路后输出近似方波的电流iac;对于副边电路,带有滤波电容C的整流电路与负载RL并联时可等效为负载Req,等效负载Req的阻抗由式(1)给出:
Req=8/π2RL≈0.81RL (1)
因此可得到如图2所示的等效电路。
其中,iac为原边逆变器输出方波电流,up为原边谐振电容Cp端电压,ip为流经原边线圈电流,us为原边谐振电容Cs端电压,is为流经原边线圈电流,电阻R1为电阻Req、Rs等效串联电阻,即:
R1=Req+Rs (2)
对于一个系统,总能获得有如下形式的差分方程:
z(k)+a1z(k-1)+…+amz(k-m)=b1u(k)+…+bn+1u(k-n) (3)
由于其中k存在一系列取值,因此可以进一步将式(3)改写为矩阵形式:
ΖL=ΗL·θ (4)
其中
ΖL=[z(k) z(k+1) … z(k+L)]T (5)
θ=[a1 … am b1 … bn+1]T (7)
式中ΗL为系统数据矩阵,L为数据长度。
显然,可以通过最小二乘法,得到过程参数向量θ的最小二乘估
计:
此时,辨识得到过程参数向量而过程参数向量与系统参数(如电阻、电感、电容)之间存在线性关系,通过对线性方程进行分解变换,便可获得负载信息,完成负载辨识。
针对本例而言,步骤2:选取原边逆变器输出的方波电流iac以及原边线圈电流ip作为状态变量,构建系统状态差分方程;
为了获得IPT系统的过程参数向量及数据矩阵,首先建立系统非线性微分方程。
由图2所示的等效电路图,构建如下方程:
选取原边逆变器输出方波电流iac以及原边线圈电流ip作为考虑的状态变量,由式(9)得出iac,ip关系式:
将上式进行Laplace变换,得到Iac(s)、Ip(s)间关系——Iac(s)、Ip(s)分别为iac、ip的拉氏变换:
其中,
为了获得较为精准的系统离散模型,经过分析比较,选用Tustin法对式(11)进行离散化处理。处理后得到如下关系式:
其中,Iac(z-1)、Ip(z-1)分别为iac、ip的Z变换;
式中c=2/T,T为采样周期。
由式(14)可以容易地构建出IPT系统的一组差分方程:
这里Ip(k)表示原边线圈电流ip在第k个采样点的值,Iac(k)表示原边逆变器输出的方波电流iac在第k个采样点的值。
步骤3:根据系统状态差分方程得到系统的数据矩阵和过程参数向量;本例中,系统的数据矩阵为:
过程参数向量:
且:
ΖL=HLθ=[Ip(k) Ip(k+1) … Ip(k+L)]T (20)
L表示需要构建的状态方程的个数,为了简化计算,L的取值为9,使得数据矩阵HL为9×9的方阵。
步骤4:按照公式(8),即可利用最小二乘法对过程参数向量进行估计;
步骤5:利用过程参数向量与系统参数之间的线性关系进行分解变换,得到系统的负载信息;
联合式(8)、(18)、(19)、(20)可辨识得到IPT系统的过程参数估算值由于系统中耦合电感值、互感值、补偿电容以及各电感寄生电阻值均为已知,并且检测得到原边逆变器输出方波电流iac以及原边线圈电流ip的瞬时值,因此可以容易地获得参数估算值与负载R1之间关系,进而联合式(1)、(2)获得最终结果,完成负载辨识。
在具体实施过程中,利用过程参数向量与系统参数之间的线性关系进行分解变换,可以得到9个副边等效电阻辨识值,求其期望值作为所述副边等效电阻R1
为了进一步验证本方案的可行性,下面基于Matlab/Simulink仿真平台建立系统仿真模型,并基于主电路拓扑搭建实验系统,其中原、副边谐振回路的主要参数由表1给出。
表1 IPT系统仿真及实验参数
图3、4分别给出了iac、ip的仿真波形以及实验波形。可以看出,iac为近似交变方波,ip为正弦波,与前文分析一致。图4所示波形中,由于开关管反并联二极管反向恢复时间不能忽略,因此在开关管切换时刻,二极管存在瞬时放电过程,引起逆变输出电流iac出现尖峰。
对电流iac,ip进行采样,在负载大小为30Ω、50Ω的情况下,分别进行50次辨识。
图5中(a)、(b)分别为负载值为30Ω、50Ω时50次负载辨识结果的统计图。图中横轴给出了负载辨识值的分布区间,纵轴则给出了在此分布区间中负载辨识值的统计个数。
由负载辨识结果的分布图可以得出,当负载值为30Ω时,仿真结果中误差率在5%以内的占97%,相对的,实验结果误差较大,其误差率在5%以内的辨识结果占91%;当负载值为50Ω时,仿真结果和实验结果误差率在5%以内的比例分别为94%以及88%。
为了更直观地表明辨识精度的平均水平,表2进一步给出了上文获得的负载辨识结果的数据分析。
表2 负载辨识结果分析
由表2,当负载RL为30Ω以及50Ω时,仿真得到负载辨识结果期望值的误差率分别为1.05%和1.58%,实验得到负载辨识结果期望值的误差率分别为3.45%和3.53%。因此,在多次辨识后,对得到的数据进行分析,并取其期望值作为最终负载辨识结果能显著提高结果的精准度。另外,通过图5以及表2可以看出实验结果的分布区间较宽、总体均方差较大、辨识结果误差较大,这是由于噪声干扰、实验仪器的误差以及算法存在一定的误差等原因。
显然,不同辨识次数与其得到辨识结果的误差率之间存在一定规律,在实际应用中应根据不同情况选择不同的辨识次数。图6通过仿真给出了当负载RL为30Ω以及50Ω时,辨识次数与负载辨识结果误差率间规律。
如图6,随着辨识次数的增加,误差率逐渐下降。值得注意的是,当辨识次数达到50次时,误差率下降曲线渐渐平滑。因此,在实际应用中如需较高的辨识精度,辨识次数选择50次即可满足,继续增加辨识次数不但不再使误差率有显著的下降,反而会增加计算成本,影响辨识的实时性。相反,如需较快的辨识速度,可根据实际情况降低辨识次数,当然这样也会相应的降低辨识精度。
最后需要说明的是,本实施例是以PS结构电流型的IPT系统为例,首先建立系统的非线性高阶微分方程,并完成对数据矩阵的构建,通过最小二乘法对系统参数进行求解,将直接对负载辨识的问题转换为了对系统参数的求解,最后通过仿真及实验验证了此方法的可行性。此IPT系统的负载辨识方法也可适用于其他结构类型的IPT系统,能完成副边可等效为电阻负载的IPT系统的负载辨识。

Claims (5)

1.一种基于最小二乘法的感应电能传输系统负载识别方法,其特征在于按照以下步骤进行:
步骤1:根据感应电能传输系统的拓扑结构构建系统的等效电路模型,该等效电路模型包括原边逆变器输出的方波电流源iac,原边谐振电容Cp以及原边谐振电容电压up,原边线圈Lp以及原边线圈电流ip,副边谐振电容Cs以及副边谐振电容电压us,副边线圈Ls以及副边线圈电流is,副边等效电阻R1
步骤2:选取原边逆变器输出的方波电流iac以及原边线圈电流ip作为状态变量,构建系统状态差分方程;
步骤3:根据系统状态差分方程得到系统的数据矩阵和过程参数向量;
步骤4:利用最小二乘法对过程参数向量进行估计;
步骤5:利用过程参数向量与系统参数之间的线性关系进行分解变换,得到系统的负载信息;
所述系统参数包括原边谐振电容Cp的电容值、原边线圈Lp的电感值、副边谐振电容电容Cs的电容值、副边线圈Ls的电感值、副边等效电阻R1的电阻值、原边线圈Lp与副边线圈Ls之间的互感值M、采样周期T;
当感应电能传输系统的原边采用并联谐振回路,副边采用串联谐振回路时,步骤2所构建的系统状态差分方程为:
其中的系数c=2/T,Ip(k)表示原边线圈电流ip在第k个采样点的值,Iac(k)表示原边逆变器输出的方波电流iac在第k个采样点的值,Rp表示原边线圈串联等效电阻的电阻值。
2.根据权利要求1所述的基于最小二乘法的感应电能传输系统负载识别方法,其特征在于:根据步骤2所构建的系统状态差分方程得到的系统的数据矩阵:
H L = - I p ( k - 1 ) ... - I p ( k - 4 ) I a c ( k ) ... I a c ( k - 4 ) - I p ( k ) ... - I p ( k - 3 ) I a c ( k + 1 ) ... I a c ( k - 3 ) . . . . . . . . . . . . - I p ( k + L - 1 ) ... - I p ( k + L - 4 ) I a c ( k + L ) ... I a c ( k + L - 4 )
过程参数向量:且:
ΖL=HLθ=[Ip(k) Ip(k+1) … Ip(k+L)]T,L表示需要构建的状态方程的个数。
3.根据权利要求2所述的基于最小二乘法的感应电能传输系统负载识别方法,其特征在于:步骤4中按照对过程参数向量θ进行最小二乘估计,其中ΗT L表示数据矩阵ΗL的转置,为过程参数向量θ的最小二乘估计值。
4.根据权利要求2所述的基于最小二乘法的感应电能传输系统负载识别方法,其特征在于:数据矩阵HL为9×9的方阵,L的取值为9。
5.根据权利要求4所述的基于最小二乘法的感应电能传输系统负载识别方法,其特征在于:利用过程参数向量与系统参数之间的线性关系进行分解变换,可以得到9个副边等效电阻辨识值,求其期望值作为所述副边等效电阻R1
CN201410626387.1A 2014-11-06 2014-11-06 基于最小二乘法的感应电能传输系统负载识别方法 Active CN104319909B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410626387.1A CN104319909B (zh) 2014-11-06 2014-11-06 基于最小二乘法的感应电能传输系统负载识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410626387.1A CN104319909B (zh) 2014-11-06 2014-11-06 基于最小二乘法的感应电能传输系统负载识别方法

Publications (2)

Publication Number Publication Date
CN104319909A CN104319909A (zh) 2015-01-28
CN104319909B true CN104319909B (zh) 2017-02-15

Family

ID=52375108

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410626387.1A Active CN104319909B (zh) 2014-11-06 2014-11-06 基于最小二乘法的感应电能传输系统负载识别方法

Country Status (1)

Country Link
CN (1) CN104319909B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108337919B (zh) 2015-06-01 2021-10-29 香港大学 用于标识无线充电系统中的线圈未对准/互耦的快速方法
CN105141046B (zh) * 2015-10-26 2017-08-15 重庆大学 Sp型icpt系统负载与互感识别方法
CN106816965B (zh) * 2017-01-18 2019-08-02 同济大学 一种谐振式无线充电系统中线圈自感的辨识方法
CN108199495B (zh) * 2018-01-16 2021-01-15 厦门大学 一种用于无线电能传输的双向能量主动注入式发射装置
CN108879992B (zh) * 2018-06-01 2021-11-05 上海大学 一种谐振式无线充电系统中一次侧电容与电感参数的辨识方法
CN110365120B (zh) * 2019-07-24 2021-10-08 南京航空航天大学 一种从原边辨识参数的非接触电能传输装置

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100045114A1 (en) * 2008-08-20 2010-02-25 Sample Alanson P Adaptive wireless power transfer apparatus and method thereof
CN101707371A (zh) * 2009-11-09 2010-05-12 重庆大学 小扰动状况下的电力系统负荷模型等值参数的辨识方法
CN101719726B (zh) * 2009-12-04 2011-10-26 重庆大学 非接触能量信号同步传输的方法及其装置
CN102969802B (zh) * 2012-11-01 2014-07-16 重庆大学 电流型无线供电系统负载自适应控制方法

Also Published As

Publication number Publication date
CN104319909A (zh) 2015-01-28

Similar Documents

Publication Publication Date Title
CN104319909B (zh) 基于最小二乘法的感应电能传输系统负载识别方法
CN111106676B (zh) Lcc-s型mc-wpt系统的磁耦合机构参数多目标优化方法
Aldhaher et al. High-input-voltage high-frequency class E rectifiers for resonant inductive links
CN110350673A (zh) 一种无线电能传输系统在最大效率跟踪下的阻抗匹配网络优化方法
CN102969801B (zh) 电流型无线供电系统负载识别方法
CN110350674B (zh) 基于动态耦合系数识别的wpt系统的最大效率跟踪方法
EP3031129A1 (en) Methods for parameter identification, load monitoring and output power control in wireless power transfer systems
Ni et al. Design and comparison of parallel and series resonant topology in wireless power transfer
Bloom et al. Design considerations for wireless electric vehicle charging
CN107276390A (zh) 双lcc谐振补偿无线充电系统的电磁干扰预测电路及方法
Zheng et al. Primary control strategy of magnetic resonant wireless power transfer based on steady-state load identification method
CN108649804B (zh) 基于磁耦合无线电能传输的直流变压器系统
Xiao et al. Multiple parameters estimation based on transmitter side information in wireless power transfer system
Hou et al. Output current characterization of parallel-series/series compensated resonant converter for contactless power transfer
Ma et al. Modeling and analysis of wireless power transfer system via unified full-load discrete-time model
Chen et al. Event-triggered model predictive control for series–series resonant ICPT systems in electric vehicles: A data-driven modeling method
Liu et al. A compact Class E rectifier for megahertz wireless power transfer
CN109756124B (zh) 一种用于无线电能传输的电流馈电式半桥谐振拓扑结构
Chow et al. Online regulation of receiver-side power and estimation of mutual inductance in wireless inductive link based on transmitter-side electrical information
Chung et al. Output voltage control for series-series compensated wireless power transfer system without direct feedback from measurement or communication
Zhang et al. Maximum efficiency point tracking control method for series–series compensated wireless power transfer system
Kiratipongvoot et al. Design a high-frequency-fed unity power-factor AC-DC power converter for wireless power transfer application
Xin et al. Generalized hamiltonian energy modeling method for wireless power transfer system
Jiang et al. Optimization of coils for magnetically coupled resonant wireless power transfer system based on maximum output power
CN112564311A (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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20200825

Address after: 402760 No. 92 Donglin Avenue, Biquan Street, Bishan District, Chongqing (No. 52 Factory Building)

Patentee after: Chongqing Huachuang Intelligent Technology Research Institute Co.,Ltd.

Address before: No.5-2, no.66-1, shabin Road, Shapingba District, Chongqing

Patentee before: Wang Zhihui

Effective date of registration: 20200825

Address after: No.5-2, no.66-1, shabin Road, Shapingba District, Chongqing

Patentee after: Wang Zhihui

Address before: 400044 Shapingba street, Shapingba District, Chongqing, No. 174

Patentee before: Chongqing University