基于TI公司DSP例程算法的SVPWM技术的研究及其实现.doc

SVPWM算法技术的研究及其实现 -By 左文全2011年3月06日星期日SVPWM的输入为两相旋转坐标系下的变量U，U，对于为什么选择两相旋转坐标系下的变量U，U作为输入变量，其原因在于后面的所有中间变量的计算都是基于两相旋转坐标系下的变量U，U计算的，所有的算法都是这样的，也是最简单的，在矢量控制应用中，整个系统的输入为ID,IQ，经过PID调节器调节和IPARK逆变换后得到U，U，这就是SVPWM的输入，SVPWM输出为6路PWM波。从输入到输出的算法实现即为本文的重点。PMSM矢量控制系统如下图0所示：(1) 由，扇区确定电压矢量所在扇区：先计算三个中间变量B0，B1，B2：在使用符号函数计算扇区中间变量P值：P值与扇区号之间的对应关系为：P123456扇区号261435 图1 扇区分布图2扇区变量P与扇区SECTOR之间对应关系图3 扇区变量P与扇区SECTOR之间关系(2) 确定单位矢量切换时间t1，t2：先计算三个中间变量X，Y，Z：根据电压矢量所在扇区确定切换时间t1，t2：扇区号123456t1-ZZX-X-YYt2XY-YZ-Z-X如果t1+t21，则；如果t1+t2UALFA,UBETA获得给定信号/= Ualfa = _IQmpy(Ud,Cosine) - _IQmpy(Uq,Sine);Ubeta = _IQmpy(Uq,Cosine) + _IQmpy(Ud,Sine); /=/通过UALFA,UBETA确定扇区变量P/= B0=Ubeta;B1=_IQmpy(_IQ(0.8660254),Ualfa)- _IQmpy(_IQ(0.5),Ubeta);/ 0.8660254 = sqrt(3)/2 B2=_IQmpy(_IQ(-0.8660254),Ualfa)- _IQmpy(_IQ(0.5),Ubeta); / 0.8660254 = sqrt(3)/2if(B0_IQ(0) P0=1;else P0=0;if(B1_IQ(0) P1=1;else P1=0;if(B2_IQ(0) P2=1;else P2=0;P=P0+2*P1+4*P2;/=/通过扇区变量P确定扇区Sector/= if(P=1) Sector=2; else if(P=2) Sector=6; else if(P=3) Sector=1; else if(P=4) Sector=4; else if(P=5) Sector=3; else if(P=6) Sector=5; if(i=256) i=0;/利用CCS Graph观察扇区变量P与扇区/=/根据扇区Sector计算Ta,Tb,Tc/=X=Ubeta;Y=_IQmpy(_IQ(0.8660254),Ualfa)+ _IQmpy(_IQ(0.5),Ubeta);/ 0.8660254 = sqrt(3)/2 Z=_IQmpy(_IQ(-0.8660254),Ualfa)+ _IQmpy(_IQ(0.5),Ubeta); / 0.8660254 = sqrt(3)/2 if(Sector=1)t1=-Z;t2=X;Taon=_IQmpy(_IQ(0.5),(_IQ(1)-t1-t2);Tbon=Taon+t1;Tcon=Tbon+t2;Ta=Taon;Tb=Tbon;Tc=Tcon;else if(Sector=2)t1=Z;t2=Y;Taon=_IQmpy(_IQ(0.5),(_IQ(1)-t1-t2);Tbon=Taon+t1;Tcon=Tbon+t2; Ta=Tbon;Tb=Taon;Tc=Tcon; else if(Sector=3)t1=X;t2=-Y;Taon=_IQmpy(_IQ(0.5),(_IQ(1)-t1-t2);Tbon=Taon+t1;Tcon=Tbon+t2; Ta=Tcon;Tb=Taon;Tc=Tbon; else if(Sector=4)t1=-X;t2=Z;Taon=_IQmpy(_IQ(0.5),(_IQ(1)-t1-t2);Tbon=Taon+t1;Tcon=Tbon+t2;Ta=Tcon;Tb=Tbon;Tc=Taon; else if(Sector=5)t1=-Y;t2=-Z;Taon=_IQmpy(_IQ(0.5),(_IQ(1)-t1-t2);Tbon=Taon+t1;Tcon=Tbon+t2; Ta=Tbon;Tb=Tcon;Tc=Taon; else if(Sector=6)t1=Y;t2=-X;Taon=_IQmpy(_IQ(0.5),(_IQ(1)-t1-t2);Tbon=Taon+t1;Tcon=Tbon+t2; Ta=Taon;Tb=Tcon;Tc=Tbon; /=/根据占空比Ta,Tb,Tc计算全比较器持?/= Ta=_IQmpy(_IQ(2),(Ta-_IQ(0.5);Tb=_IQmpy(_IQ(2),(Tb-_IQ(0.5);Tc=_IQmpy(_IQ(2),(Tc-_IQ(0.5);Tmp = (int32)T1Period*(int32)Modulation; / Q15 = Q0*Q15，占空比设置MPeriod = (int16)(Tmp15); / Q0 = (Q15-Q0)Tmp = (int32)MPeriod*(int32)Ta; / Q15 = Q0*Q15，计算全比较器CMPR1赋值EvaRegs.CMPR1 = (int16)(Tmp16) + (int16)(MPeriod1); / Q0 = (Q15-Q0)/2 + (Q0/2)Tmp = (int32)MPeriod*(int32)Tb; / Q15 = Q0*Q15，计算全比较器CMPR2赋值EvaRegs.CMPR2 = (int16)(Tmp16) + (int16)(MPeriod1); / Q0 = (Q15-Q0)/2 + (Q0/2)Tmp = (int32)MPeriod*(int32)Tc; / Q15 = Q0*Q15，计算全比较器CMPR3赋值EvaRegs.CMPR3 = (int16)(Tmp16) + (int16)(MPeriod