SVPWM算法实现

clc;clear;close;

%%SVPWM算法 %%%判断扇区的 %%%r模拟假设(2/3)Udc r=20; Udc=30;

Ts=62.5%%单位微妙 r=r*cos(pi/6);%%%电压矢量圆环半径 k=720; for i=0:k     angle=(i*pi)/180;     Uta=r*cos(angle);     Utb=r*sin(angle);

    U0=Uta;     U1=Utb;     U2=(sqrt(3)*Uta-Utb)/2;     U3=-((sqrt(3)*Uta+Utb)/2);     %%%I区的判断     if (U1>=0 && U2>=0 && U3<0)        A=1;        B=1;        C=0;     end         %%%二区的判断N=1     if U1>0&&U2<0&&U3<=0        C=0;         B=0;        A=1;      end     

    %%%三区的判断N=5     if U1>=0&&U2<0&&U3>0        C=1;        B=0;        A=1;     end

    %%%四区的判断N=4     if U1<0&&U2<=0&&U3>0         C=1;         B=0;         A=0;     end

    %%%五区的判断N=6     if U1<0&&U2>0&&U3>=0         C=1;         B=1;         A=0;     end          

    %%%六区的判断N=2     if U1<=0&&U2>0&&U3<0        C=0;        B=1;        A=0;     end

    N=4*C+2*B+A;

    if N==3         Section=1;     elseif N==5        Section=3;     elseif N==6          Section=5;     elseif N==2         Section=6;     elseif N==1         Section=2;     elseif N==4         Section=4;     end

    U1=Utb;     U2=(sqrt(3)*Uta-Utb)/2;     U3=-(sqrt(3)*Uta+Utb)/2;     K=(sqrt(3)*Ts)/Udc;

    if Section==1     %%%1区的     T4=K*U2;     T6=K*U1;     T0=((1/2)*(Ts-T4-T6));     T7=T0;     elseif Section==2     %%%2区的     T2=-K*U2;     T6=-K*U3;     T0=((1/2)*(Ts-T2-T6));     T7=T0;     elseif Section==3     %%%3区的     T2=K*U1;     T3=K*U3;     T0=((1/2)*(Ts-T2-T3));     T7=T0;     elseif Section==4     %%%4区的     T1=-K*U1;     T3=-K*U2;     T0=((1/2)*(Ts-T1-T3));     T7=T0;     elseif Section==5     %%%5区     T1=K*U3;     T5=K*U2;     T0=((1/2)*(Ts-T1-T5));     T7=T0;     elseif Section==6     %%%6区     T5=-K*U1;     T4=-K*U3;     T0=((1/2)*(Ts-T4-T5));     T7=T0;     end

    if Section==1        Tx=T4;        Ty=T6;     elseif Section==2        Tx=T2;        Ty=T6;     elseif Section==3        Tx=T2;        Ty=T3;     elseif Section==4        Tx=T1;        Ty=T3;     elseif Section==5        Tx=T1;        Ty=T5;     elseif Section==6        Tx=T4;        Ty=T5;     end

    Ta=(Ts-Tx-Ty)/4;     Tb=Ta+Tx/2;     Tc=Tb+Ty/2;

    if Section==1        Pwm1(i+1)=Ta;        Pwm2(i+1)=Tb;        Pwm3(i+1)=Tc;     elseif Section==2        Pwm1(i+1)=Tb;        Pwm2(i+1)=Ta;        Pwm3(i+1)=Tc;     elseif Section==3        Pwm1(i+1)=Tc;        Pwm2(i+1)=Ta;        Pwm3(i+1)=Tb;     elseif Section==4        Pwm1(i+1)=Tc;        Pwm2(i+1)=Tb;        Pwm3(i+1)=Ta;     elseif Section==5        Pwm1(i+1)=Tb;        Pwm2(i+1)=Tc;        Pwm3(i+1)=Ta;     elseif Section==6        Pwm1(i+1)=Ta;        Pwm2(i+1)=Tc;        Pwm3(i+1)=Tb;     end end hold on plot(0:k,Pwm1,'r'); plot(0:k,Pwm2,'b'); plot(0:k,Pwm3,'y');