基于MATLAB的感应电机矢量控制系统的仿真
基于Matlab的交流感应电机矢量控制系统研究

矢量控制又称为 磁场定 向控制 , 即把交 流 电机空
坐标 轴 。若 忽 略 由反 电动 势 引起 的交 叉耦 合 , 只需 检测 出定 子 电流 的 d轴分 量 , 就可 以观 测 转 子磁 通 幅值 , 当转 子磁 通恒定 时 , 电磁 转 矩 与定 子 电流 的 q 轴 分 量成正 比 , 过 控 制 定 子 电 流 的 叮轴分 量 就 可 通 以实 现对 电磁转矩 的控 制 。 由电压方 程 的 d轴分 量 控 制转 子磁 通 , 分 量 控 制 转 矩 从 而 实 现 磁 通 和 q轴 转 矩 的解 耦 控 制 。我 们 把 这 样 选 取 的 d轴 称 为 轴 , q轴 为逆 时钟 转 9 。 即垂 直 与矢 量 的 方 而 0,
Smu ik. Be a e o h ia v ntg so r d t n lP D o to lr we u e e sn l e r n c n r 1 i ln c us ft e d s d a a e fta ii a I c n r l , s d a n w ig e n u o o to .Th e - o e e r a s n b lt n ai i a e b e e tfe y t e smu ai n rs ls o a ii a d v ld t h v e n t si d b h i lto e u t . y y i
Ab t a t sr c :Ba e n t t e tc lmo e fAC n u to tr,a mprv d AC i d to trv co o to s d o hemah mai a d lo i d c in moo ni o e n ucin moo e trc n r l
司的 Back 又将这 种一般 化 的概 念形成 系统理 论 , l he s 以磁场定 向控 制 的名称 在 S m n e e i e sR n w上 发表 , e 并 申请了专 利。随着各 国专家学 者研究 的不 断深入 , 电 子计算 机技术 、 大规模集成 电路技术 以及 电力 电子器 件技术 的发展 , 特别是数 字信号处理 来自 ( S ) D P 的应用 ,
基于MATLABSIMULINK感应电机矢量控制系统仿真

在上图间接磁场定向控制中, 主要包括以下 几个模块: ①速度调节器( ASR )、电流调节器( ACR ) 和磁通调节器( A R )。 ASR 输入为转速命令给 定值和实际转速, 输出为 q 轴电流命令值。ACR 输入为 d(q)轴电流指令值和 d(q)轴电流实际值, 输出为 d(q)轴电压。 A R 输入包含 d 轴电流和 经过弱磁环节输出的磁通指令值, 输出为 d 轴电 流指令值。 ASR 、 ACR 和 A R 内部结构类似, 本文只给出了如图 5 所示的 A R 内部仿真模型。 ②同步角计算。该模块完成同步角的计算, 输入变量包括电机实际转速、q 轴电流指令值和 磁通指令值, 输出为坐标变换的同步角速度和同 步角。内部结构如图 6 所示。 该环节完成 d 轴电 ③前馈电压解耦环节[1]。
3 系统仿真模型
3.1 感应电机仿真模型 式(1)中,当参考坐标系定向于静止坐标系 时, 结合公式(2)、 (3)利用 MATLAB 中 S-function 可以建立如图 2 所示的感应电机的仿真模型[2]。
r e dr
Lm e i r s 1 ds
e Lm i qs
(7)
图 4 矢量控制核心算法仿真模块
[2] [1]
2.1 感应电机原理 在任意两相旋转坐标系下感应电机的数学 模型如可用下列方程来描述[3]:
dids dt ( Rs 1 Lm
Ls
r
)ids iqs (
Ls Lr r
)dr ...
diqs dt ids ( Rs
Lmr
Ls Lr
本文采用 SIMULINK 软件实现的仿真模型, 采用模块化设计, 可以很方便的将各个模块连接 起来组成一个系统, 改变控制策略, 只需将其中 的功能子模块替换为相应的功能模块, 而系统的 其他模块无须做任何变动,大大缩短了建模时 间,提高了系统仿真模型的通用性。
基于MATLAB的永磁同步电机矢量控制系统仿真

3 矢量控制与坐标变换模块
矢量控制 模 块 实 现 的 就 是 PMSM 的 矢 量 控 制 算 法 , 其 中 ,
空间矢量变换 ABC- DQ 和 DQ- ABC 的变换矩阵分别为:
! CABC- DQ=
2 3
"
##sinθ
#
#
#
#
##cosθ
#
#
#
1 #
#
#
# $
!
2
sin(
θ-
2 3
π)
cos( θ- 2 π) 3
运行方式由电机电磁转矩符号决定 (为正则是电动机状态, 为
负则是发电机状态)。为了简化模型, 可以假定转子永磁磁极在
定子上产生的感应磁通是正弦分布的, 并且由于通常永磁同
步电机的气隙较大, 可以近似地忽略定电机铁心的磁饱和。
因此永磁同步电机在 d- q 轴的电压方程为:
d dt
id
=
1 Ld
ud
-
R1 Ld
π)
cos(θ+ 2 3π) Nhomakorabea$
1 !2
& & & & ’
根 据 上 面 两 式 在 MATLAB 环 境 下 可 分 别 得 到 dq/abc 和
abc/dq 坐标变换的子模块, 用以实现 PMSM 的矢量控制算法,
将电流转换为电压。
4 电流滞环型 PWM模块
电 流 滞 环 PWM 模 块 实 现 的 是 PMSM 的 滞 环 电 流 控 制 方
The S imula tion Ba s e d on Ma tla b for Ve ctor Control of P e rma ne nt Ma gne t S ynchronous Motor
基于matlab 生成C代码的电机矢量控制仿真模型研究

重庆大学本科学生毕业设计(论文)基于matlab 生成C代码的电机矢量控制仿真模型研究学生:曾宇航学号:20114346指导教师:余传祥副教授专业:电气工程与自动化重庆大学电气工程学院二O一五年六月Graduation Design(Thesis) of Chongqing UniversityGenerated C code of motor vector control simulation model research based onMATLABUndergraduate: Zeng YuhangStudent Id:20114346Supervisor:Prof Yu ChuanxiangMajor: Electrical engineering and automationSchool of Electrical EngineeringChongqing UniversityJune 2015摘要电机在国民生产中占据重要地位,而传统的电机控制开发流程相较于工业技术的发展已经相对落后,本次毕业设计便是对一种前沿、高效的电机控制系统开发流程进行论述和验证。
课题选择研究对象为三相交流异步电机并采用矢量控制系统进行控制。
首先通过MATLAB/SIMULINK对所选电机进行建模,之后对矢量控制系统进行建模。
根据实验所选用的电机进行参数配置,配合矢量控制系统的数学模型完成整个仿真模型的构建,根据仿真结果不断矫正仿真模型、优化模型结构,并确定仿真模型的最优参数配置;然后根据仿真模型的控制模块并结合控制系统的硬件电路构建控制系统的C代码生成模型,生成所需的DSP可执行C 代码,将生成的可执行代码下载到以TI系列DSP为控制核心的硬件控制系统中进行硬件实验,矫正模型参数并验证生成代码的正确性。
在整个课题进展过程中根据生成代码的运行结果多次对仿真模型进行参数修正,并同步对代码生成模型进行修正,在如此反复过程中,优化了控制系统模型,使得生成代码能够更高效的运作,最后成功的完成了整个毕业设计,验证了此电机控制系统开发流程及生产可行性代码的正确性和可行性。
基于MATLAB永磁同步电动机矢量控制系统的仿真研究

基于MATLAB永磁同步电动机矢量控制系统的仿真研究永磁同步电动机(Permanent Magnet Synchronous Motor,PMSM)是一种应用广泛的高性能电机。
在工业领域,PMSM通常采用矢量控制方法来实现精确的速度和位置控制。
本文基于MATLAB对PMSM矢量控制系统进行仿真研究,探讨其工作原理及性能。
首先,PMSM的矢量控制系统由控制器、电机和传感器三部分组成。
其中,控制器根据电机的反馈信号和期望输出来计算电机的控制信号。
传感器用于测量电机的转速、位置和电流等参数,反馈给控制器。
通过调节控制信号,控制器可以实现电机的速度和位置控制。
在PMSM的矢量控制系统中,通常采用dq轴矢量控制方法,将三相电流转换为直流参考轴和旋转参考轴的dq坐标系,进而对电机进行控制。
其次,本文利用MATLAB软件对PMSM矢量控制系统进行了仿真实验。
首先,建立了PMSM电机的数学模型,包括电机的动态方程、反电动势方程和电流方程。
然后,在MATLAB环境中编写程序,实现电机模型的数值求解和控制算法的计算。
通过调节控制参数,可以对电机的速度和位置进行精确控制,并实时监测电机的工作状态。
在仿真实验中,通过改变电机的负载情况、工作电压和控制参数等条件,分析了PMSM矢量控制系统的性能。
实验结果表明,当负载增加时,电机的转动惯量增大,控制系统的响应时间变长,但依然可以实现精确的速度和位置控制。
当电机的工作电压增加时,电机的输出功率和转速增大,但也会产生更大的电流和损耗。
当控制参数的比例增益和积分时间常数变化时,系统的稳定性和动态性能均会受到影响,需要进行合理的调节。
总结起来,本文基于MATLAB对PMSM矢量控制系统进行了仿真研究,探讨了其控制原理和性能。
通过仿真实验,可以深入理解PMSM矢量控制系统的工作原理,优化系统的参数和性能,并为实际应用提供参考。
基于MATLAB永磁同步电动机矢量控制系统的仿真研究

基于MATLAB永磁同步电动机矢量控制系统的仿真研究近年来,永磁同步电动机(Permanent Magnet Synchronous Motor,PMSM)作为一种高效率、高功率密度和高控制精度的电机,被广泛应用于工业和汽车领域。
针对永磁同步电动机的控制问题,矢量控制(Vector Control)成为了一种重要的控制策略。
本文将使用MATLAB对永磁同步电动机矢量控制系统进行仿真研究。
首先,我们需要建立永磁同步电动机的动态模型。
永磁同步电动机是一种非线性多变量系统,其数学模型可以描述为:\begin{cases}\frac{{d\theta}}{{dt}} = \Omega_m \\\frac{{d\Omega_m}}{{dt}} = \frac{1}{{J}}(T_{em} - T_{L}) \\ \frac{{di_q}}{{dt}} = \frac{1}{{L_q}}(v_q - R_s i_q -\Omega_m L_d i_d + e_f) \\\frac{{di_d}}{{dt}} = \frac{1}{{L_d}}(v_d - R_s i_d +\Omega_m L_q i_q)\end{cases}\]其中,$\theta$为转子位置,$\Omega_m$为电机机械角速度,$T_{em}$为电磁转矩,$T_{L}$为负载转矩,$i_q$和$i_d$为电流的直轴和正交轴分量,$v_q$和$v_d$为电压的直轴和正交轴分量,$R_s$为电机电阻,$L_q$和$L_d$为电机的定子轴和直轴电感,$e_f$为反电势。
接下来,我们可以使用MATLAB建立永磁同步电动机的矢量控制系统。
首先,我们需要设计控制器,其中包括速度环控制器和电流环控制器。
速度环控制器用于调节电机的机械角速度,电流环控制器用于控制电机的电流。
在速度环控制器中,我们可以选择PID控制器,其输入为速度误差,输出为电机的电压指令。
基于MATLABSimulinkSimPowerSystems的永磁同步电机矢量控制系统建模与仿真

基于MATLABSimulinkSimPowerSystems的永磁同步电机矢量控制系统建模与仿真一、本文概述随着电力电子技术和控制理论的快速发展,永磁同步电机(Permanent Magnet Synchronous Motor, PMSM)因其高效率、高功率密度和优良的调速性能,在电动汽车、风力发电、机器人和工业自动化等领域得到了广泛应用。
然而,PMSM的高性能运行依赖于先进的控制系统,其中矢量控制(Vector Control, VC)是最常用的控制策略之一。
矢量控制,也称为场向量控制,其基本思想是通过坐标变换将电机的定子电流分解为与磁场方向正交的两个分量——转矩分量和励磁分量,并分别进行控制,从而实现电机的高性能运行。
这种控制策略需要对电机的动态行为和电磁关系有深入的理解,并且要求控制系统能够快速、准确地响应各种工况变化。
MATLAB/Simulink/SimPowerSystems是MathWorks公司开发的一套强大的电力系统和电机控制系统仿真工具。
通过Simulink的图形化建模环境和SimPowerSystems的电机及电力电子元件库,用户可以方便地进行电机控制系统的建模、仿真和分析。
本文旨在介绍基于MATLAB/Simulink/SimPowerSystems的永磁同步电机矢量控制系统的建模与仿真方法。
将简要概述永磁同步电机的基本结构和运行原理,然后详细介绍矢量控制的基本原理和坐标变换方法。
接着,将通过一个具体的案例,展示如何使用Simulink和SimPowerSystems进行永磁同步电机矢量控制系统的建模和仿真,并分析仿真结果,验证控制策略的有效性。
将讨论在实际应用中可能遇到的挑战和问题,并提出相应的解决方案。
通过本文的阅读,读者可以对永磁同步电机矢量控制系统有更深入的理解,并掌握使用MATLAB/Simulink/SimPowerSystems进行电机控制系统仿真的基本方法。
感应电机矢量控制系统的仿真

《运动控制系统》课程设计学院:班级:姓名:学号:日期:成绩:感应电机矢量控制系统的仿真摘要:本文先分析了异步电机的数学模型和坐标变换以及矢量控制基本原理,然后利用Matlab /Simulink软件进行感应电机的矢量控制系统的仿真。
采用模块化的思想分别建立了交流异步电机模块、逆变器模块、矢量控制器模块、坐标变换模块、磁链观测器模块、速度调节模块、电流滞环PWM调节器,再进行功能模块的有机整合,构成了按转子磁场定向的异步电机矢量控制系统仿真模型。
仿真结果表明了该系统转速动态响应快、稳态静差小、抗负载扰动能力强,验证了交流电机矢量控制的可行性和有效性。
关键词:异步电机;坐标变换;矢量控制;Simulink仿真一、异步电机的动态数学模型和坐标变换异步电机的动态数学模型是一个高阶、非线性、强耦合的多变量系统,异步电机的数学模型由下述电压方程、磁链方程、转矩方程和运动方程组成。
电压方程:礠链方程:转矩方程:运动方程:异步电机的数学模型比较复杂,坐标变换的目的就是要简化数学模型。
异步电机数学模型是建立在三相静止的ABC坐标系上的,如果把它变换到两相坐标系上,由于两相坐标轴互相垂直,两相绕组之间没有磁的耦合,仅此一点,就会使数学模型简单了许多。
(1)三相--两相变换(3/2变换)在三相静止绕组A、B、C和两相静止绕组α、β之间的变换,或称三相静止坐标系和两相静止坐标系间的变换,简称 3/2 变换。
(2)两相—两相旋转变换(2s/2r变换)从两相静止坐标系到两相旋转坐标系 M、T 变换称作两相—两相旋转变换,简称 2s/2r 变换,其中 s 表示静止,r 表示旋转。
图1、异步电动机的坐标变换结构图二、感应电机矢量控制原理感应电机是指定转子之间靠电磁感应作用,在转子内感应电流以实现机电能量转换的电机。
感应电机是异步电机的一种,异步电机主要是指感应电机。
以上所讲,异步电机的动态数学模型是一个高阶、非线性、强耦合的多变量系统,通过坐标变换,可以使之降阶并化简,但并没有改变其非线性、多变量的本质。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
图 1 矢量控制过程框图
2. 2 感应电机解耦模型
感应电动机的解耦控制是从感应电动机的数学模型出
= 26. 526N ·m; 1. 5s后减负载 , 使得总的负载转矩 TL = 16.
526N·m。所得到的机械角速度 、电磁转矩以及解耦后的电
流 d轴和 q轴分量波形如图 4所示 。
进行调节并达到稳定 。1s后负载的加入使转速稍微下降但 很快就能跟踪上给定转速 ,电磁转矩相应的增加 。1. 5s后加 入减负载 ,使得转速稍有增加并很快调节到给定转速 , 电磁 转矩相应下降 ,同时很快稳定 。说明此调速性能的抗扰性能 好 。从图 4中还可以看出 ,在负载转矩变化过程中定子电流 的 q轴分量随着负载转矩的变化而相应的变化 ,而 d轴分量 基本保持不变 ,这是因为转子磁场定向的感应电机定子电流 的 d轴分量用于建立主磁通 ,而定子电流的 q轴分量产生的 磁通用于抵消转子电流产生的磁通 ,感应电机基本上实现了 完全解耦控制 。
2 感应电机矢量控制系统的仿真模型
感应电机矢量控制变频调速系统中 ,逆变器 、异步电动 机及矢量控制器是关键部分 。因此 ,对整个系统在数学建模 时 ,需要分别考虑异步电动机模型 、逆变器以及矢量控制技 术在 MATLAB中的实现 。 2. 1 异步电机的矢量控制理论
异步电机定子通过三相平衡电流 iA 、iB 、iC 和两相垂直的
电压解耦型的异步电动机模型方程如下所示 。方程式 ( 3) 、
( 5) 、( 6)表达了这种解耦关系 。转子磁场定向电压解耦型
的异步电动机模型可写为 [1 ] :
u sd
=ri s sd
+δL
s
·
i ds
-
ω m
rL
s
i
sq
+ 1
+1δrψ· r
(1)
usq
= rs isq
+δL
s
·
i sq
+ωm rδL s isd
第 23卷 第 4期 文章编号 : 1006 - 9348 ( 2006) 04 - 0326 - 04
计 算 机 机矢量控制系统的仿真
尉冰娟 ,王明渝 ,张淑国
(重庆大学高电压电工新技术教育部重点实验室 ,重庆 400044)
摘要 :该文根据感应电机矢量控制的基本原理 ,基于 MATLAB / SIMUL INK构造了电压型 SPWM 逆变器供电的感应电机按转 子磁场定向的矢量控制系统的仿真模型 。通过仿真试验验证了模型的正确性 ,并分析验证了该文所建立的调速系统具有快 速的启动性能及良好的抗扰动性能 ,基本上实现了完全解耦控制 ;另外 ,分析了感应参数变化对调速系统的影响 ,说明转子 电阻及励磁电感等参数对获得良好的调速系统的性能的重要性 。该模型可通用于笼型异步电机 ,在使用时只需输入不同的 电机参数即可 。 关键词 :感应电机 ;矢量控制 ;仿真 ;参数变化 中图分类号 : TP391. 9 文献标识码 : B
达式为 :
u dx
=
-
ω m
rδL s
i
sq
=
-
ω m
rL
′s isy
(9)
udy =ωm rL ′s isx + (Ls - L ′s )ωm r | im r |
( 10)
式中 , L ′s为定子暂态电感 L ′s = (Ls - Lm2 /L r ) , | im r |为转子励
磁电流的模 。我们以此解耦分量建立感应电机的解耦电路
图如图 3所示 。图 3中包含了坐标转换模块 ,电压控制型变
— 327 —
图 3 感应电机矢量控制系统仿真模型
频器模块 , SPWM模块 ,磁通模型 ,解耦模型以及速度和转矩 等的 PI调节器等模型 ,这些模块可由 SIMUL INK模型库分别 建立 ,然后封装成 Subsystem。
3 仿真实验
d和 q轴分量 ; Tr - 转子时间常数 ; Te - 电磁转矩 ;δ- 定子电
感因数 ,δ= 1 - Lm2 /L rLs;δr - 转子漏磁因数 ,δr = 1 - L r /Lm 。 在感应电动机的模型中 , 由于转子磁场定向控制的结
果
, 工程上可以认为
ψ r
=常数
, 这样感应电动机的电压方程
可写为 :
usd
= rs isd
+δL
s
·
i sd
-
ω m
rδL s
isq
(7)
usq
=
rs isq
+δL
s
·
i sq
+ωm rδLs
ω + 1 +mδr rψr
(8)
方程式 ( 7) 、( 8) 分别为异步电动机的交 、直轴电压方
程 。从异步电动机的交 、直轴电压方程可以看到 , 当转子磁
通恒定时 ,转子磁场定向的感应电机解耦电压分量的一般表
仿真实验 1:对图 3所示模型进行仿真
仿真所选的电机参数为 : P = 4KW , p = 2, Rs = 1. 37Ω, R r
= 1. 1Ω, Lm = 0. 14122H, Lds = 0. 004872611H, Ldr =
010079577471H,
ω e
= 303.
441665 rad / s,
平衡电流 αi 、βi 可产生等效的旋转磁场 。因而三相平衡电流 iA 、iB 、iC 和两相电流 αi 、βi 之间存在着确定的矢量变换关系 。 以产生同样的旋转磁场为准则 ,在三相定子坐标系下的定子 电流 iA 、iB 、iC 通过三相 /二相坐标变换可以等效成两相静止 坐标系下的交流电流 αi 、βi 。再通过按转子磁场定向的旋转 变换 ,可以等效成同步旋转坐标系下的直流电流 id、iq。这样 就实现了电流励磁分量与转矩分量之间的解耦 ,从而达到对 交流电机的磁链和电流分别控制的目的 。 d - q绕组相当于 直流电机的励磁绕组和电枢绕组 。由于这些被控矢量在物
发 ,应用现代控制理论中的解耦控制方法 , 通过状态反馈 , 使
得原来复杂的多变量 、非线性 、强耦合的系统解耦线性化 , 以
实现电动机的转速与转子磁链之间的动态解耦 。对于感应
型的交流电动机 ,有定子磁场 、气隙磁场和转子磁场三种磁
场定向方式 ,但以转子磁场定向构成的系统不仅可使电动机
的转矩控制和励磁控制完全解耦且较易实现 ,转子磁场定向
的机理 ,首先要获得正弦波和三角波两种信号 ,但是在 SIM 2
UL INK中没有提供三角波的模块 ,这里通过查表的方法来产
生 。然后根据自然采样法让三角波和正弦波相比较 ,在交点
处产生方波 。这样就得到图 2所示的 SPWM 功能模块 。
图 2 SPWM 功能模块
2. 4 磁通模型 转子磁场定向参照系下的磁通模型的建立方程如方程
1 引言
70年代初提出的矢量控制方法 ,通过坐标变换和磁场定 向控制 ,把交流电动机的定子电流分解成磁场定向坐标的磁 场电流分量和与之相垂直的坐标转矩电流分量 ,从而实现两 者之间的解耦 ,得到类似于直流电机的转矩模型 ,并可仿照 直流电机进行快速的转矩控制和磁通控制 ,使系统动态性能 得到显著改善 ,从而使交流电机的调速技术取得了突破性的 进展 。目前 ,运用矢量控制已成为当今交流变频调速系统的 主流 。MATLAB提供的 SIMUL INK是一个用来对动态系统 进行建模 、仿真和分析的软件包 ,它具有模块化 、可重栽 、可 封装 、面向结构图编程及可视化等特点 ,可大大提高系统仿
+ωm r1
+1δrψr
(2)
·
isd = im r + Tr im r
(3)
ω mr
=ωr
i
+ sq
Tr
i
m
r
(4)
Te
= 1
+1δrLm
im r isq
(5)
ψ r
= Lm
im r
(6)
式中 u, i - 电压和电流 ; L, r,ω - 电感 、电阻 、角速度和磁链 ;
下标 s, r, m - 定 、转子和磁场量 ;下标 d, q - 旋转坐标系的轴
来实现电机的矢量控制 。首先将电机的三相定子电流经过
3 /2变换转换为转子磁通空间相量 ,将直轴电流分量输入磁 通模型计算出磁通量 。用转速参考量和电机的测量转速差
经 P I调节器获得电机转矩参考量 。由于本设计中采用的是 电压控电压型逆变器 ,于是需要将电流参考量转化为电压参 考量 ,所以经过一个解耦电路和 P I调节器就获得了直轴和 交轴的电压参考 ,将电压参考量变换后就获得了电机需要的 定子电压 。转子磁场定向的感应电动机矢量控制系统原理
收稿日期 : 2005 - 03 - 18
— 326 —
真的效率和可靠性 。在多种矢量控制方法中 , 按转子磁场 定向的矢量控制运用较为普遍 ,本文将结合这种矢量控制和 SIMUL INK的特点 ,介绍一种感应电机按转子磁场定向的矢 量控制系统的建模仿真方法模型 ,利用仿真模型 ,进行控制 系统的仿真实验 。为同类调速系统提供有效 、可靠的研究分 析依据 。
S im ula tion of Vector Con trol System of Induction M ach ine Ba sed on M ATLAB
W E I B ing - juan,WAN G M ing - yu, ZHAN G Shu - guo
(H igh voltage and electric engineering new technology, Key laboratory of m inistry of education, ChongqingUniversity, Chongqing 400044, China)