%C2=(1*10^-10):(1*10^-10):(110*10^-9);
C2=50*10^-9;
C1=32*10^-9;
Xc1=1/(2*pi()*10000*C1);
Xcm=1./(2*pi()*10000*C2);
U1=5;
R=300;
R2=0.5:0.1:100;
R3=0;
Re=U1*((Xcm^2+R*R2+R2.^2)./((R+R2).^2+Xcm^2)-(Xc1^2+R*R3+R3^2)./((R+R3)^2+Xc1^2));
Im=U1*((-Xcm*R)./((R+R2).^2+Xcm^2)-(-Xc1*R)./((R+R3)^2+Xc1^2));
uhel=(180/pi())*angle(Re+Im*1i);
figure(1)
subplot(211)
plot(R2,uhel)
title("Závislost fázového posunu na ESR měřeného kondenzátoru")
xlim([0.5 100])
xlabel("ESR měřeného kondenzátoru [F]")
ylabel("Fázový posun [°]")
subplot(212)
plot(R2,sqrt(Re.^2+Im.^2))
title("Závislost amplitudy na ESR měřeného kondenzátoru")
xlim([0.5 100])
xlabel("ESR měřeného kondenzátoru [F]")
ylabel("Amplituda [V]")
