Lecture 3 - Power Electronics Applications - Part 2
Lecture 3 - Power Electronics Applications - Part 2
Vs=220; Vm=sqrt(2)*Vs;
R=20;L=20e-3;f=50;
w=2*pi*f;
XL=w*L;
Z=sqrt(R^2+XL^2);
theta=atan(XL/R);
Full-wave uncontrolled RL-loads
ioi=Vm/Z*sin(theta)*((1+exp(-pi/tan(theta)))/(1-exp(-pi/tan(theta))));
io=Vm/Z*(sin(wt-theta))+(ioi+Vm/Z*sin(theta))*exp(-wt/tan(theta));
io=Vm/Z*(sin(wt-pi-theta))+(ioi+Vm/Z*sin(theta))*exp((-wt+pi)/tan(theta));
Full-wave uncontrolled RL-loads
vs2=Vm*sin(wt2);
vo2=-vs2;
%Full-wave uncontrolled RL-loads
clc,clear
Vs=220; Vm=sqrt(2)*Vs;
R=20;L=20e-3;f=50;w=2*pi*f;XL=w*L;
Z=sqrt(R^2+XL^2);theta=atan(XL/R);
ioi=Vm/Z*sin(theta)*((1+exp(-pi/tan(theta)))/(1-exp(-pi/tan(theta))));
wt=0:1e-2:2*pi;
vs=Vm*sin(wt);
for i=1:length(wt)
if wt(i)<=pi
vo(i)=vs(i);
io(i)=Vm/Z*(sin(wt(i)-theta))+(ioi+Vm/Z*sin(theta))*exp(-wt(i)/tan(theta));
vD1(i)=0;
else
vo(i)=-vs(i);
io(i)=Vm/Z*(sin(wt(i)-pi-theta))+(ioi+Vm/Z*sin(theta))*exp((-wt(i)+pi)/tan(theta));
vD1(i)=vs(i);
end
end
Full-wave uncontrolled RL-loads
Vavg=trapz(wt,vo)/(2*pi)
Iavg=Vavg/R
Vrms=sqrt(trapz(wt,vo)/(2*pi))
Irms=sqrt(trapz(wt,io.^2)/(2*pi))
Po=Irms^2*R
PF=Po/(Vs*Irms)
wt=wt*180/pi;
subplot(411),plot(wt,vs,'linewidth',2),ylabel('VS'),grid on
subplot(412),plot(wt,vo,'linewidth',2),ylabel('Vo'),grid on
subplot(413),plot(wt,io,'linewidth',2),ylabel('io'),grid on
subplot(414),plot(wt,vD1,'linewidth',2),ylabel('VD1'),grid on
Full-wave uncontrolled RL-loads (Results)
Full-wave controlled RL-loads
CCM DCM
if theta<alpha
%DCM
if wt(i)<B-pi
vo(i)=-vs(i);
io(i)=Vm/Z*(sin(wt(i)+pi-theta))-
Vm/Z*sin(alpha-theta)*exp((alpha-wt(i)-
pi)/tan(theta));
elseif wt(i)>=B-pi&&wt(i)<=alpha
vo(i)=0;io(i)=0;
elseif wt(i)>alpha&&wt(i)<B
vo(i)=vs(i);
io(i)=Vm/Z*(sin(wt(i)-theta))-Vm/Z*sin(alpha-
theta)*exp((alpha-wt(i))/tan(theta));
elseif wt(i)>B&&wt(i)<alpha+pi
vo(i)=0;io(i)=0;
else
vo(i)=-vs(i);
io(i)=Vm/Z*(sin(wt(i)-pi-theta))-
Vm/Z*sin(alpha-theta)*exp((alpha-
wt(i)+pi)/tan(theta));
end
%CCM
else
if wt(i)<alpha
vo(i)=-vs(i);
io(i)=Vm/Z*(sin(wt(i)+pi-theta))+(ioi-
Vm/Z*sin(alpha-theta))*exp((alpha-wt(i)-
pi)/tan(theta));
elseif wt(i)>alpha&&wt(i)<alpha+pi
vo(i)=vs(i);
io(i)=Vm/Z*(sin(wt(i)-theta))+(ioi-Vm/Z*sin(alpha-
theta))*exp((alpha-wt(i))/tan(theta));
else
vo(i)=-vs(i);
io(i)=Vm/Z*(sin(wt(i)-pi-theta))+(ioi-
Vm/Z*sin(alpha-theta))*exp((alpha-
wt(i)+pi)/tan(theta));
end, end
M-File
clc,clear
global theta alpha
Vs=220; Vm=sqrt(2)*Vs;
R=10;L=10e-3;f=50;
alpha=60*pi/180;
w=2*pi*f;XL=w*L;
Z=sqrt(R^2+XL^2);
theta=atan(XL/R);
ioi=Vm/Z*sin(theta-alpha)*((1+exp(-pi/tan(theta)))...
/(1-exp(-pi/tan(theta))));
wt=0:1e-3:2*pi;
vs=Vm*sin(wt);B=fzero('B_CRL',pi);
M-File
for i=1:length(wt)
if theta<alpha
B=fzero('B_CRL',pi);
if wt(i)<B-pi
vo(i)=-vs(i);
io(i)=Vm/Z*(sin(wt(i)+pi-theta))-Vm/Z*sin(alpha-
theta)*exp((alpha-wt(i)-pi)/tan(theta));
elseif wt(i)>=B-pi&&wt(i)<=alpha
vo(i)=0;io(i)=0;
elseif wt(i)>alpha&&wt(i)<B
vo(i)=vs(i);
io(i)=Vm/Z*(sin(wt(i)-theta))-Vm/Z*sin(alpha-
theta)*exp((alpha-wt(i))/tan(theta));
elseif wt(i)>B&&wt(i)<alpha+pi
vo(i)=0;io(i)=0;
end
M-File
else
if wt(i)<alpha
vo(i)=-vs(i);
io(i)=Vm/Z*(sin(wt(i)+pi-theta))+(ioi-Vm/Z*sin(alpha-theta))*exp((alpha-wt(i)-
pi)/tan(theta));
elseif wt(i)>alpha&&wt(i)<alpha+pi
vo(i)=vs(i);
io(i)=Vm/Z*(sin(wt(i)-theta))+(ioi-Vm/Z*sin(alpha-theta))*exp((alpha-wt(i))/tan(theta));
else
vo(i)=-vs(i);
io(i)=Vm/Z*(sin(wt(i)-pi-theta))+(ioi-Vm/Z*sin(alpha-theta))*exp((alpha-
wt(i)+pi)/tan(theta));
end
end
end
M-File
Vavg=trapz(wt,vo)/(2*pi)
Iavg=Vavg/R
Vrms=sqrt(trapz(wt,vo)/(2*pi))
Irms=sqrt(trapz(wt,io.^2)/(2*pi))
Po=Irms^2*R
PF=Po/(Vs*Irms)
wt=wt*180/pi;
subplot(311),plot(wt,vs,'linewidth',2),ylabel('VS'),grid on
subplot(312),plot(wt,vo,'linewidth',2),ylabel('Vo'),grid on
subplot(313),plot(wt,io,'linewidth',2),ylabel('io'),grid on
Simulink File
Simulink File
Simulink File
Simulink File
------
---
---
Rectifier circuits using Simulink
VSI
VSI
End Lecture (1)