%Punto 5
syms AH AL t1 t2 T n
eq1=AH-AL-1;
eq2=AH*0.33+AL*0.67;
sol2=solve(eq1,eq2,'AH,AL');
AH=sol2.AH; AL=sol2.AL; f=100; T=1/f; wq=2*pi*f; t1=0.33*T; t2=0.67*T;
% %Calcolo funzione di rete
R23=(R2*R3)/(R2+R3)
H=(R3*1/(j*n*wq*C2))/((1/(j*n*wq*C2)+R23)*(R2+R3));
%Dall'analisi della forma della funzione di rete vediamo che è un filtro
%passa basso quindi ci vanno bene 20 armoniche
N=20;
An_sym=(2/T)*(int(AH*cos(n*wq*t),t,0,t1)+int(AL*cos(n*wq*t),t,t1,T));
Bn_sym=(2/T)*(int(AH*sin(n*wq*t),t,0,t1)+int(AL*sin(n*wq*t),t,t1,T));
j_t=0;
iR2=0;
for n=1:N
An=subs(An_sym);
Bn=subs(Bn_sym);
FAN=An;
FBN=Bn*exp(-j*pi/2);
FCN=FAN+FBN;
j_t=j_t+abs(FCN)*cos(n*wq*t+angle(FCN));
FIRn=subs(H,'n',n)*FCN;
iR2=iR2+abs(FIRn)*cos(n*wq*t + angle(FIRn));
end
figure(3), ezplot(j_t,[0,3*T]),hold on, grid on, axis auto, ezplot(iR2,[0,3*T]), axis auto, grid on