Ex 2.16

L = 0.1e-3;
C = 0.25e-9;
R = 50;
V0 = 10e-3;
f = 0:10:1e7;
omega = 2*pi*f;
VR = R./sqrt(R^2 + (omega*L - 1./(omega*C)).^2);
plot(f,VR);
xlabel('Frekvens f');
ylabel('V_R');
[maxVR,i] = max(VR);
f0 = f(i);
disp(['Resonansfrekvens f0 = ',num2str(f0),' Hz']);
V1 = VR(round(i*1.1));
disp(['Spänning vid 10% högre frekvens = ',num2str(V1),' V']);
f1 = f(0.49995*maxVR<VR & VR<=0.5*maxVR);
disp(['Frekvens vid 50% av max spänning = ',num2str(f1),' Hz']);

Resonansfrekvens f0 = 1006580 Hz
Spänning vid 10% högre frekvens = 0.38261 V
Frekvens vid 50% av max spänning = 940020  1077860 Hz