% Prova d'esame Albanese, de Magistris del 14-03-2016

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)

%% Problema 1
R1 = 2;
R2 = 8;
R3 = 4;
R4 = 4;
R5 = 6;
J  = 2;
E  = 20;
%
% Metodo 1: Thevenin
disp('<a href=""> Traccia A, Esercizio 1, Metodo 1: Thevenin </a>')
R12 = paral(R1,R2);
Vac0 = -R12*J;
Vbc0 = -E*R4/(R3+R4);
Vab0 = Vac0 - Vbc0;
E0 = Vab0
R0 = R12 + paral(R3,R4)
Vab = E0*R5/(R0+R5)
%
% Metodo 2: equivalenza e Millmann
disp('<a href=""> Traccia A, Esercizio 1, Metodo 2: equivalenza e Millmann </a>')
EJ = R12*J;
Vcb = (EJ/(R12+R5) + E/R3) / (1/(R12+R5) + 1/R3 + 1/R4)
I5 = (Vcb-EJ) / (R12+R5)
Vab = R5*I5

%% Problema 2
Em = 100;  omega = 500;
R1 = 20;
R2 = 2;
C  = 20e-6;
L  = 100e-3;
n  = 5;
%
% Metodo 1: riporto al primario
disp('<a href=""> Traccia A, Esercizio 2, Metodo 1: riporto al primario </a>')
E_   = 100/sqrt(2);    % fasore in termini del valore efficace
ZL_  =  1i*omega*L
ZC_  = -1i/(omega*C)
R2p  = n^2*R2;         % R2 riportata al primario
ZLC_ = paral(ZL_,ZC_)
Z_   = paral(R2p,ZLC_)
V1_  = E_*Z_/(R1+Z_);  % tensione al primario
V1   = abs(V1_)        % valore efficace
PR2p = V1^2/R2p;       % potenza assorbita da R2p
PR2  = PR2p
%
% Metodo 2: Thevenin ai capi di R2
disp('<a href=""> Traccia A, Esercizio 2, Metodo 2: Thevenin ai capi di R2 </a>')
V10_ = E_*ZLC_/(R1+ZLC_); % tensione al primario quando il secondario è a vuoto
V20_ = V10_/n             % tensione al secondario a vuoto
E0_  = V20_               % tensione a vuoto di Thevenin
Zp_  = paral(R1,ZLC_)
Z0_  = Zp_/n^2;
V2_  = E0_*R2/(R2+Z0_);   % tensione ai capi di R2
V2   = abs(V2_)           % valore efficace
PR2_chk  = V2^2/R2