%?SPWM電壓型逆變器的數(shù)學模型
??clear?all
??clc
??%?----------------------------------------------------
%?參數(shù)輸入
??Ud=390;%直流輸入電壓
??F1=50;%調(diào)制波頻率
??h=0.00001;%采樣時間間隔
??Tend=1.2/F1;%仿真時間
??Kf=21;%載波比K
??Ka=4/3;%調(diào)制比（Ka>1）
??F2=Kf*F1;%載波頻率
??T1=1/F1;
??T2=1/F2;
??P（1）=0;%相位
??P（2）=2*pi/3;
??P（3）=4*pi/3;
??%?----------------------------------------------------
%?主程序
??t=0;%采樣時間
??i=1;%采樣次數(shù)
??while?t<=Tend
	tt2=mod（tT2）;
	if?tt2<=T2/2%確定三角波信號表達式
???????A2=Ka*（4*tt2/T2-1）;
????else
	???A2=Ka*（-4*tt2/T2+3）;
????end
%確定逆變器三個開關變量S（k）
	for?k=1:3
??	??A1=sin（2.*pi*t/T1-P（k））;%確定正弦波信號表達式
??	??if?A1>A2
?????????S（k）=1;
??????else
?????????S（k）=-1;
??????end
????end
%?逆變器輸出線電壓（負載電機線電壓）
????Uab=（S（1）-S（2））*Ud/2;
	Ubc=（S（2）-S（3））*Ud/2;
	Uca=（S（3）-S（1））*Ud/2;
%?逆變器輸出端至直流電源中性點的電壓
????Uag=S（1）*Ud/2;
????Ubg=S（2）*Ud/2;
????Ucg=S（3）*Ud/2;
%?負載電機相電壓
????Ua0=（2*Uag-Ubg-Ucg）/3;
????Ub0=（-Uag+2*Ubg-Ucg）/3;
????Uc0=（-Uag-Ubg+2*Ucg）/3;
%?----------------------------------------------------
%建立時間、線電壓、相電壓、載波信號與調(diào)制信號的數(shù)組????
????tx（i）=t;
????Uaby（i）=Uab;
????Ubcy（i）=Ubc;
????Ucay（i）=Uca;
????Uagy（i）=Uag;
????Ubgy（i）=Ubg;
????Ucgy（i）=Ucg;
????Ua0y（i）=Ua0;
????Ub0y（i）=Ub0;
????Uc0y（i）=Uc0;
????A1y（i）=A1;
????A2y（i）=A2;
????t=t+h;
????i=i+1;
??end??
%?---------------------------------------------------------------------------------
??figure（1）;??plot（txUaby）?title（‘AB線電壓‘）;?xlabel（‘t（s）‘）;?ylabel（‘Uab‘）;?grid;
??figure（2）;??plot（txUbcy）?title（‘BC線電壓‘）;?xlabel（‘t（s）‘）;?ylabel（‘Ubc‘）;?grid;
??figure（3）;??plot（txUcay）?title（‘CA線電壓‘）;?xlabel（‘t（s）‘）;?ylabel（‘Uca‘）;?grid;
??figure（4）;??plot（txUagy）?title（‘AG電壓‘）;???xlabel（‘t（s）‘）;?ylabel（‘Uag‘）;?grid;
??figure（5）;??plot（txUbgy）?title（‘BG電壓‘）;???xlabel（‘t（s）‘）;?ylabel（‘Ubg‘）;?grid;
??figure（6）;??plot（txUcgy）?title（‘CG電壓‘）;???xlabel（‘t（s）‘）;?ylabel（‘Ucg‘）;?grid;
??figure（7）;??plot（txUa0y）?title（‘A相電壓‘）;??xlabel（‘t（s）‘）;?ylabel（‘Ua0‘）;?grid;
??figure（8）;??plot（txUb0y）?title（‘B相電壓‘）;??xlabel（‘t（s）‘）;?ylabel（‘Ub0‘）;?grid;
??figure（9）;??plot（txUc0y）?title（‘C相電壓‘）;??xlabel（‘t（s）‘）;?ylabel（‘Uc0‘）;?grid;
??figure（10）;?plot（txA1y‘c+:‘txA2y‘c+:‘）;??title（‘正弦波‘）;???xlabel（‘t（s）‘）;?ylabel（‘A‘）;??grid;
%?***********************************************************************************