Circuit Theory/1Source Excitement/Example 7/code

From Wikibooks, open books for an open world
Jump to navigation Jump to search

In Matlab:

R = 10;
L = .01;
w = 377;
Vm = 120*2^.5;
omega = 2*pi/3;
Im = Vm/(R^2+(L*w)^2)^.5;
alpha = omega - atan(L*w/R);
syms t A C;
I = Im * cos(w*t+alpha) + A*exp(-t/(L/R)) + C;
leftSide = Vm*cos(w*t + omega) - I*R;
rightSide = L * diff(I,t);
vpa(leftSide, 3);
vpa(rightSide,3);
t=0;
teq = subs(leftSide) - subs(rightSide)
vpa(teq,3)
Ieval = subs(I);
vpa(Ieval,3)
S = solve(Ieval,teq,A,C);
S = [S.A S.C];
vpa(S(1),3)
vpa(S(2),3)

In MuPad and MatLab (have to cut and paste)

R := 10;
L := .01;
w := 377;
Vm := 120*2^.5;
omega := 2*pi/3;
Imagnitude := Vm/(R^2+(L*w)^2)^.5;
alpha := omega - atan(L*w/R);
Itime := Imagnitude * cos(w*t+alpha) + A*exp(-t/(L/R)) + C;
dItime := diff(Itime,t);
t:=0;
numeric::solve([Itime=0,Vm*cos(w*t + omega) - Itime*R -  L * dItime = 0],[A,C],FixedPrecision)

sig1 = 5.986596875*sin(2*pi/3 - atan(0.377));
sig2 = cos(2*pi/3 - atan(0.377));
A = -16.97056275*cos(2*pi/3) - 2.775557562*10^-17*sig2 - sig1
C =   16.97056275*cos(2*pi/3) - 15.87956731*sig2+ sig1