% Prova d'esame del 2018-03-19
function q

clear
format short g
format compact

[ST,I] = dbstack;
%fprintf(2,['$$$$ ' num2str(100*rand) ' $$$$\n'])
fprintf(2,'$$$$ %0.5g $$$$ %s\n', 100*rand,ST.file)

%% Quesiti A1 e B1
clear
E = 20;    J = 3;
R1 = 10;   R2 = 15;
R3 = 10;   R4 = 15;
%
dispah(' Quesito A1, Metodo 1: Millmann ')
EJ = R1*J; % trasformazione della coppia R1 J da Norton a Thevenin
V3 = (E/R4+EJ/(R1+R2))/(1/R3+1/R4+1/(R1+R2)); % Millmann
I3 = V3/R3;
V2 = (V3-EJ)*R2/(R1+R2);
VJ = V3-V2;
PJ = VJ*J;
disp(table(V3,I3,V2,VJ,PJ))

disp(' ')
dispah(' Quesito B1, Metodo 1: Millmann ')
V4 = (E/R3-EJ/(R1+R2))/(1/R3+1/R4+1/(R1+R2));
I3 = (V4-E)/R3;
V2 = (V4+EJ)*R2/(R1+R2);
VJ = V2-V4;
PJ = VJ*J;
disp(table(V4,I3,V2,VJ,PJ))

%% Quesiti A2 e B2
clear
disp(' ')
dispah(' Quesito A2 ')
aEm = 300;   omega = 250;
R1 = 1;   R2 = .5;
L1 = 4e-3;
C  = 4e-3;
%
dispah(' t < 0 ')
Em = aEm; % * exp(1i*pi/4);
ZL1 =  1i*omega*L1;
ZC  = -1i/(omega*C);
Z   = R1 + R2 + ZL1;
I   = Em/Z;
disp(table(ZL1,ZC,Z,I))
absI  = abs(I);
beta = angle(I);
betadeg = beta/pi*180;
I0 = absI*cos(beta); % t=0
disp(table(absI,beta,betadeg,I0))
str  = sprintf('i(t) = %0.4g*cos(%0.4g*t + %0.4g)', absI,omega,beta);
disp(str)
%
dispah(' t > 0 ')
Zr = R1 + ZL1;
Ir = Em/Zr;
absIr = abs(Ir);
betar = angle(Ir);
betardeg = betar/pi*180;
Ir0 = absIr*cos(betar);
disp(table(Zr,Ir,absIr,betar,betardeg,Ir0))
str  = sprintf('ir(t) = %0.4g*cos(%0.4g*t + %0.4g)', absIr,omega,betar);
disp(str)
tau = L1/R1;
lambda = -1/tau;
disp(table(omega,tau,lambda))
disp('i(t) = (I0-Ir0)*exp(-t/tau) + ir(t)')
str = sprintf('     = %0.4g*exp(-t/tau) + ir(t)', I0-Ir0);
disp(str)

disp(' ')
dispah(' Quesito B2 ')
clear
omega = 200;
R1 = 1;   R2 = .5;
L1 = 5e-3;
C  = 5e-3;
disp('rispetto al quesito A2 cambia solo:')
tau = L1/R1;
lambda = -1/tau;
disp(table(omega,tau,lambda))

%% Quesiti A3 e B3
clear
disp(' ')
dispah(' Quesito A3 ')
E = 20;
R1 = 10;   R2 = 5;   g = .5;
%
G1 = 1/R1;   G2 = 1/R2;
G = [G1   -G1;
     g-G1 G1+G2];
Jcc = [1; -1]*G1*E;
disp(table(G,Jcc))

disp(' ')
dispah(' Quesito B3 ')
G = [G1+G2 g-G2;
     -G2   G2];
Jcc = [1; -1]*G2*E;
disp(table(G,Jcc))

%% Quesito C1
clear
disp(' ')
dispah(' Quesiti C1 e D1 ')
E = 230;   f = 50;
P = 1000;  Q = 1000;
R = 10;   XL = 10;   XC = -10;
%
A = sqrt(P^2+Q^2);
absI = A/3/E;
phi = atan(Q/P);
alpha1 = 0;                beta1 = alpha1-phi;
alpha2 = alpha1-2*pi/3;    beta2 = alpha2-phi;
alpha3 = alpha2-2*pi/3;    beta3 = alpha3-phi;
E1 = E*exp(1i*alpha1);  I1 = absI*exp(1i*beta1);
E2 = E*exp(1i*alpha2);  I2 = absI*exp(1i*beta2);
E3 = E*exp(1i*alpha3);  I3 = absI*exp(1i*beta3);
ZL = 1i*XL;
ZC = 1i*XC;
VOOp = -(E1/R+E2/ZL+E3/ZC)/(1/R+1/ZL+1/ZC);
VR = VOOp+E1;
VL = VOOp+E2;
VC = VOOp+E3;
IR = VR/R;
IL = VL/ZL;
IC = VC/ZC;
I123 = [I1; I2; I3];  absI123 = abs(I123);  beta123  = angle(I123);  beta123deg  = beta123*180/pi;
VRLC = [VR; VL; VC];  absVRLC = abs(VRLC);  alphaRLC = angle(VRLC);  alphaRLCdeg = alphaRLC*180/pi;
IRLC = [IR; IL; IC];  absIRLC = abs(IRLC);  betaRLC  = angle(IRLC);  betaRLCdeg  = betaRLC*180/pi;
I1tot = I1+IR;
I2tot = I2+IL;
I3tot = I3+IC;
I123tot = [I1tot; I2tot; I3tot];  absI123tot = abs(I123tot);  beta123tot = angle(I123tot);  beta123totdeg = beta123tot*180/pi;
V23 = E2-E3;
disp(table(I123,absI123,beta123,beta123deg))
disp(table(VOOp,V23))
disp(table(VRLC,absVRLC,alphaRLC,alphaRLCdeg))
disp(table(IRLC,absIRLC,betaRLC,betaRLCdeg))
disp(table(I123tot,absI123tot,beta123tot,beta123totdeg))
dispah(' C1 ')
W = real(V23*conj(I2tot));
disp(table(V23,W))
str  = sprintf('i3(t) = sqrt(2)*%0.4g*cos(2*pi*f*t + %0.4g)', abs(I3tot),angle(I3tot));
disp(str)
dispah(' D1 ')
W = real(-V23*conj(I3tot));
disp(table(W))
str  = sprintf('i1(t) = sqrt(2)*%0.4g*cos(2*pi*f*t + %0.4g)', abs(I1tot),angle(I1tot));
disp(str)


return

publish('q_180918_AD','pdf')

%break

function c=paral(a,b)
c=a.*b./(a+b);

function dispah(str)
disp(sprintf('<a href="">%s</a>', str))

function dispc(str,k)
%
% aqdispc('[variablename1 variablename2 ...]') displays the following:
% [variablename1 variablename2]
% values1 values2

if nargin==1 || k==1
  %disp(str)
  disp(['<a href="">' str '</a>'])
end
disp(evalin('caller',str))