psim永磁同步电动机矢量控制仿真
详细介绍了 如何用psim软件实现永磁电动机的矢量仿真Vol 16 No900系统仿真学报May 2004对永磁同步电机(PMSM)数学模型的电压方程式(4进行2.4SPWM模块abc/da坐标变换,可得dq坐标系下的电压方程式(8):正弦脉宽调制,它是以正弦波作为基准的调制波,三角0d=rid+pod-oo(8)波作为载波,当调制波与载波相交时,由他们的交点确定变4=rig+ peyd+oou频器开通的时间,从而产生等幅不等宽的脉冲波形。图8是其中SPWM变频器的控制回路,一组三相对称的正弦参考信号Og =lsig, Od=liia+OatUa、Ub、Ua与三角波参考信号U,相比较,作为三相桥臂6个功率开关元件的控制信号。式中:、口—d、q相定子磁链;a—转子磁场对定子的交链;、Ld、q相绕组申感。图5、图6中的U、U子模块的功能就是实现方程式(8),心子模块的底层结构如图5所示,U子模块的底层结构如图6所示。Kp1载波他号[Iq图8SPWM变频器控制回路图5U4子模块结构框图2.5速度控制模块速度控制模块的结构较为简单,如图9所示,输入:参考转速和实际转速的差值,输出:q相电流参考值lqrfq其中,Kp为P控制器中P(比例)的参数,KpT为P控制器中I(积分)的参数, Saturation饱和跟幅模块将输出的q相参考电流幅值限定在要求范围內。图6U子模块结构框图Ia ret23坐标变换模块坐标变换模块实现的是d旋转坐标系下的两相相电压x口Un、Ug向abc静止坐标系的三相电压Ua、Ub、U的等Kp/T效变换。与矢量控制模坎类似,d2abc模块实现的是dq两相向ab三相的变换,模块的底层结构嬗鹵7腙示,功由图9速度控制模块结构框图dq/abc电流变换方程式(9戾实现26电压逆变模块C。=Ucos6+ sine+Ucos(-120°)+Usin(-120°)+U电压逆变模块实现的是逆变器功能,输入为SPWM模U= U, cos(日+12°)+Usin(日+120°)+U块给出的逆变控制信号,输出为三相相电压。图10是电压逆变模块结构框图。该模块可采用中提供的通用逆变模块搭建,只需3对iT功率开关器件,反向并联续流二极管c根据SPWM模块给出的控制信号,控制6个开关器件顺序四导通和关断,从而产生三相相电压输出。2s0本本去J2PI/3图10电压逆变模块结构框图图7dq2abc模块结构框图C1994-2010chinaAcadcmicJOurnalElcctronicPublishingHousc.Allrightsrcscrved.http://www.cnki.nctVoL 16 No 5May 2004纪志成,等∶基于PSIM永磁同步电机矢量控制系统的仿真建模·901·3仿真结果组电感L=0.06H,转动惯量J=000179kgm2,额定转速1500rmin,极对数n=2。为了验证所设计的PMSM根据上述所建立的PMSM矢量控制系统的仿真模型,矢量控制系统仿真模型的静、动态性能,系统带负载T在PSIM6.0的仿真环境下进行了仿真,PMSM电机参数设Nm起动,得到系统转、转矩、d-q两相相电流仿真曲置为:电机功率p=500W,直流电压l=220V,定子相绕线如图114所示。组电阻R=432,定子d相绕组电感L=0027H,q相绕200013001200110015000010008v8006500003001.0000200400100t(s)t(3)图11转速响应曲线图12转矩响应曲线10016.001400-0504000.002005c01.00200300500t(s)图13d相电流波形图14q相电流波形由仿真滅形可以看出:在n=150 Or/in的参考转速下,性能威者模拟相同的实验糸件比较不同控制策喲的优劣,系统带负载启动响应快速且平稳,如两相电流波形较为理汋分析和设计永磁冋步电机控制系统提供了有效地手段和想,稳态运行时转速无静差。仿真结果证明了本文所提岀的工具,也为实际电机控制系统的设计和调试提供了新的思这种新型PMsⅥ仿真建模方法的有效性。路。结论参考文献:[I P Pillay, R Krishnan. Modeling, simulation, and analysis of本文在分析PMSM数学模型的基础上,提出了基于permanent-magnet motor drives, Part 2: The permanent-magnetPSIM的PMSM控制系统仿真模型。该控制系统采用速度环synchronous motor drive [] IEEE Trans. on Industry ApplicationsPⅠ控制和电流环矢量控制的双闭环控制方法并在1989,25(2):265-273[21 Pragasan Pillay, R Krishnan. Modeling of permanent magnet motorPSIM60 SIMCAD环境下对该控制系统进行了设计与仿真。drives[J]. IEEE Trans on Industry Electronics, 1988. 35(4): 537-541仿真结果表明:波形符合理论分析,系统能平稳运行,具有③CiP. Zhu jG.IaQ,P,ctal. Simulation of nonlinear switche较好的静、动态特性。采用该PMSM仿真模型,可以十分eluctance motor drives with PSIM[C]. Proceedings of Electrica便捷地实现、验证控制算法,只需对部分功能模块迸行替换Machines and Systems, 2001, 2: 1061-10644]纪志成薛花沈艳霞.永磁同步电机调速系统的樸糊PⅠ智能控制或修改,就可实现控制策珞的改换或改进,不仅可以节省控新方法[电工技术学报,2003,18(6:53-58制方案的设计周期,快速验证所设计的控制算法,更可以充「5Y. SJeon, H.S. Mok,G, H. Choe,cta. a new simulation model of分利用计算机仿真的优越性,通过修改系统参变量或人为加PMSM motor with rcal back EMF wavcform[C]. Procccding from入不同扰动因素来考察不同实验条件下电机系统的动、静态Computers in Power Electronics, 2000, 16-18: 217-220C1994-2010chinaAcadcmicJOurnalElcctronicPublishingHousc.Allrightsrcscrved.http://www.cnki.nct
- 2021-05-06下载
- 积分:1
北航矩阵论学习笔记
北京航空航天大学矩阵理论学习笔记,总结版,学霸总结,可以放心下载使用北京航空航天大学张京蕊工程系统工程系月录§0补充公式§1 Jordan(约当)标准形(简介)§2线性变换与矩阵.24§3欧式空间与QR分解.48§4常用矩阵分解●鲁D●●·,,,,,74§5范数与级数.81§6广义逆A..97§7直积拉直及应用105矩阵理论A笔记北京航空航天大学张京蕊工程系统工程系§0补充公式令A=(a)mxn∈C",风x)=4o+a1x0定义f(4)=a0+a1A+…+amAm,其中I=l若g(x)=bo+b1x+…+bkx,(x)g(x)=g(x)(x),则f(4)“g(A)=g(A)f(A)分块公式A10令A,A1,A2为方阵00 A(2)f(A),fx)为多项式令A=,A1,4为方阵AO(2)f(4)相似关系:A∽B,(PAP=B)则:(1)(P1AP)=P!AP,(k=0,1,2,(2)f(PAP)=PfA)P,f(x)为多项式许尔公式( schur):每个复方阼,A-(a)nxm都相似丁上三角形。共113页矩阵理论A笔记第1页北京航空航天大学张京蕊工程系统工程系即:P-1AP=其中41,,的次序可以任意指定Pf:用归纳法n=1时成立可以设为(n=1阶方阵成立对于n阶方阵A=(an)2×n设特征值为A,…,n取为对应的特征向量,记为a1≠0,A1=1ax1把a1扩展为可逆方阵Q=(a1,02,xn)22e又:g(a,a,…,.)=(Qa,Qba2,,Qan)其中Qe1,aQ0Q4=QA(a1a2,…an)2-I(Aa,,AAQ=(Qa,、+)…,(*)其中A1为(n-1所阶0人:0 A为由假设,对于A1必有(n-1)阶P,可推出PAPEg知n阶方阵A,适合A=0,则A+|=1共113页矩阵理论A笔记第2页北京航空航天大学张京蕊工程系统工程系Pf:A=0→任意特征值A=0→>=0即全体特征值为00,,00由需要P1AP=→PAP+7=1pAP+PP|=P(4+1)P=14+1→A+1=-1注(1)若AB(相似),则AB有相同特征值A,可引入记号:谱集(4)={2,2,…,λ}(全体特征值,含重复)A∽B→o()=o(B)(2)A∽B→1-A=1-B-(2-4元一2)…(-n),特征多项式PAP=B=A-A=p(1-A)P=A-B引理:若A0A2,则M-A|-|M1-4|-1-A1|2-A2→ar(4)=o(A)∪a(42k+1,Ak-2,…n1f(x2)设B,f(x)为多项式,则f(B)=o f(,)引理:若n阶方阵A的谱集(4)=1,42,…},则)的全体特社值为)2,…,),x)为多项式Pf:由许尔定理,A∽B→f(4)∽f(B)f(x)的全体特征值为(A1)(42),,()},fx)为多项式例如:4为A的特征值→x为4的特征值。(x)=x)共113页矩阵理论A笔记第3页北京航空航天大学张京蕊工程系统工程系引理:令B,f(x)=x-B|=(x-41)(x-12)….(x-n)则fB)=(B-1D(B-21)…(B-A1D=0Pf:当n=2时,B=0x2f(x)=(x-1)(x-2)000→f(B)-(B-41)(B-21)(2-元)0(00∴得证★ Cayley公式:设n阶方阵A的特征多项式为f(x)=|x-A|=a+a1x+…,+x则f4)=anl+a14+…,+4=0Pf:由许尔PAP=B=→P(4)P=fp3P)=f(B)=0(引理)定义若多项式x)使(4)=0,则称(x)为A的个零化式结论方阵A的特征多项式)=1x1-4为A的一个零化式g特征多项式fx)=x2可知:f(A)=A2+1=+I=00-1Hx)=|xI-A|=(x-)(x+i,(i=√-1,t2=-1)f(A)=(A-i)(4+i1=0也可取P=则PPAP=,对角形共113页矩阵理论A笔记第4页北京航空航天大学张京蕊工程系统工程系g:知A则A"=0Onxn由 Cayley特征多项式:f(x)=x"→f(4)=4"=0Ex 1. A=求P使得PP为对角阵,并验证 Cayley定理2.A=cd/,求fx)=x1-4验证f4)-0补充知识( schur公式、 Cayley公式)应用由A"=-(a0I+a1A+1A·AanA+a142+…+a.,A把①代入②→Am1=(-)+(+)4+…+(+)41可知:任何和(m≥n)都可写成,4,,A的线性组合任何多项式g(A),可写成lA,…,4的组合。Fg:若A|≠0,fx)=xI-A|=a0+a1x+…+x",ao=|-A|≠则A可用A的多项式表示∵a1A+a242+…+an21A-+A"--a072A(a1+a24+…+an-142+A)Aa1+…+an1A"2+A-1零化式定义:若g(x)=b+b1x+…+bnx,使得g(4)=bn+b14+…+bn4m=0,称g(x)为方阵A的零化式注:方阵A的零化式有无穷多个∴取特征多项式x)则4)=0任取式M(x),f(A(4)=0→f(x)(x)也是零化式极小式定义:在方阵A的零化式集合中,去次数最小的且首项系数为1的零化式m(x),称它为A的极小式共113页矩阵理论A笔记第5页北京航空航天大学张京蕊工程系统工程系注:极小式唯一性质:①极小式m(x)必为特征多项式fx)=|xI-A的因式。②特征多项式fx)=|x1-A的每个单因子(x-4)也是极小式的因子)f(x)=|x1-4=(x-x)(x-2)则极小式m(x)=(x-x)(x-2)y…(x-,),且1≤l1≤m1,1≤l2≤m2,…,1≤l≤n,41,A2…,n互不相同210EgA=020,B=020,求极小式mA(),m()解:(1)|xI-A|=(x-2)(x-1)极小式为:(x-2)(x-1)或(x-2)(x-1)计算:(4-2/)4-1)=000010k≠000000∴极小式为m4(x)=(x-2)2(x-1)(2)|-B|-=(x-2)2(x-1)00000计算:(B-2)B-1)=000010=000-1八000∴极小式为m(x)=(x-2)(x-1)Eg求下列极小式m(x)4604-60(1)A=-3-50,(2)B=2-303-6100210(3)C,(4)D=000010002000解:(1)特征多项式|x7-A|-(x-1)(x+2)极小式为:(x-1)(x+2)或(x-1)(x+2)共113页矩阵理论A笔记第6页北京航空航天大学张京蕊工程系统工程系验证:(4-D(A+2D=0∴极小式为m(x)-(x-1)(x+2)(3)解法如下引理:A1,A2的极小式为m1(x),m2(x)A10的极小式m(x)等丁m1(x),m2(x)的最小公倍式0A2(此引力可推广到A1,42,43)0100极小式为(x-1)2,0010极小式为(x-1)0取最小公倍式(x-1)2为C的极小式。460(5)F-/40,A1=020|,A00 A0123-6101O引理;设D=,则D的极小式m(x)O验证:先证D的性质(右推公式)设A-(an)xn=(a1,2,…,n)则有AD=(0,01,a2,,.m1)AD2=(0,0,∞1,,x12)AD=(0,….0.,a1,,axn)单位向量技巧:∵AI=A(en,e2…,en)=(el,leAen)=A=(a1, a2,. a,)∴Ae1=01,Ae2=(2,.,A→AD=A(0,e1,e2,…,en-1)=(0,a1,a2…,an-)同理AD2=(AD)D=(0,.01,.12)可知:D-1-(D)Dy2-(0.,0,,e1)≠0D"=(D)D1=0,而特征多项式(x)=|x1-D|=x,极小式为某个x共113页矩阵理论A笔记第7页
- 2020-12-09下载
- 积分:1