Scribd
Upload
7 views2 pages

Power_Lab_2

The document outlines the calculation of instantaneous current and power for an inductive load with a given supply voltage. It provides the expressions for current i(t) and power p(t), along with MATLAB code to plot these functions over an interval. The plots include instantaneous voltage, current, power, real power, and reactive power.

Uploaded by

uyelikicin2013
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
7 views2 pages

Power_Lab_2

The document outlines the calculation of instantaneous current and power for an inductive load with a given supply voltage. It provides the expressions for current i(t) and power p(t), along with MATLAB code to plot these functions over an interval. The plots include instantaneous voltage, current, power, real power, and reactive power.

Uploaded by

uyelikicin2013
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 2

Figure 1

The supply voltage in Figure 1 is given by v(t) = 100 cos (ωt) and the load is
inductive with impedance Z = 1.25∠60° Ω. Determine the expression for the
instantaneous current i(t) and the instantaneous power p(t). Use MATLAB to plot
i(t), v(t), p(t), pR(t), and pX(t) over an interval of 0 to 2π.

⏟
𝜃 = 𝜃⏟𝑣 − 𝜃⏟𝑖
𝑁𝑒𝑡 𝑃ℎ𝑎𝑠𝑒 𝐴𝑛𝑔𝑙𝑒 𝑉𝑜𝑙𝑡𝑎𝑔𝑒 𝑃ℎ𝑎𝑠𝑒 𝐴𝑛𝑔𝑙𝑒 𝐶𝑢𝑟𝑟𝑒𝑛𝑡 𝑃ℎ𝑎𝑠𝑒 𝐴𝑛𝑔𝑙𝑒

v(t) = 100 cos (ωt)


i(t) = 80 cos (ωt - 60°) A
p(t) = v(t) i(t) = 8000 cos ωt cos (ωt - 60°) W
clear all
clc
Vm = 100; thetav = 0; % Voltage amplitude and phase angle
Z = 1.25; gama = 60; % Impedance magnitude and phase angle
thetai = thetav - gama; % Current phase angle in degree
theta = (thetav - thetai)*pi/180; % Degree to radian
Im = Vm/Z; % Current amplitude
wt = 0:.05:4*pi; % wt from 0 to 2*pi
v = Vm*cos(wt); % Instantaneous voltage
i = Im*cos(wt + thetai*pi/180); % Instantaneous current
p = v.*i; % Instantaneous power
V = Vm/sqrt(2); I=Im/sqrt(2); % rms voltage and current
P = V*I*cos(theta); % Average power
Q = V*I*sin(theta); % Reactive power
S = P + j*Q; % Complex power
pr = P*(1 + cos(2*(wt + thetav)));
px = Q*sin(2*(wt + thetav));
PP = P*ones(1, length(wt)); % Average power of length w for plot
xline = zeros(1, length(wt)); % generates a zero vector
wt=180/pi*wt; % converting radian to degree

subplot(2,2,1), plot(wt, v, wt,i,wt, xline), grid


title(['v(t)=Vm coswt, i(t)=Im cos(wt+',num2str(thetai), ')'])
xlabel('wt, degree'), legend('v(t)','i(t)')

subplot(2,2,2), plot(wt, p, wt, xline)


title('p(t)=v(t) i(t)'),xlabel('wt, degree'), grid, legend('p(t)')

subplot(2,2,3), plot(wt, pr, wt, PP,wt,xline), grid


title('pr(t)'), xlabel('wt, degree'), legend('pr(t)')

subplot(2,2,4), plot(wt, px, wt, xline), grid


title('px(t)'), xlabel('wt, degree'), legend ('px(t)')

You might also like