Contents

% /home/dimarzio/Documents/working/12303/matlab/prl.m
% Sun Oct 21 21:31:47 2018
%
% Power in an RL Series Circuit
%

Setup

L=100e-3;R=50;f=60;
t=[0:0.1:30]*1e-3;
Vs=120*sqrt(2);

omega=2*pi*f;
vs=real(Vs*exp(1j*omega*t)); % Source

Voltages and currents

Is=Vs/(R+1j*omega*L);
is=real(Is*exp(1j*omega*t));
Vr=Is*R;
Vl=Is*1j*omega*L;
vr=real(Vr*exp(1j*omega*t));
vl=real(Vl*exp(1j*omega*t));

figure;plot(t*1e3,vs,'b-',t*1e3,vr,'g--',t*1e3,vl,'r-.');grid on;
xlabel('t, Time, ms');ylabel('v, Voltage, Volts');
legend('Source','Resistor','Inductor');

figure;plot(t*1e3,vs,'b-',t*1e3,is*100,'m--');grid on;
xlabel('t, Time, ms');ylabel(' ');
legend('Voltage, V','Current*100, A');

Power

ps=-is.*vs;
pr=is.*vr;
pl=is.*vl;

figure;plot(t*1e3,ps,'b-',t*1e3,pr,'g--',t*1e3,pl,'r-.');grid on;
xlabel('t, Time, ms');ylabel('P, Power, W');
legend('Source','Resistor','Inductor');

figure;plot(t,ps+pl+pr);grid on;
xlabel('t, Time, ms');ylabel('P_{total}, Power, W');

Published with MATLAB® R2017a