%zobrazení frewkvenční charakteristiky RLC obvodu
%       C1 = 100 nF   L = 0,01 H   R = 1000 Ohm
%     ,------||---------vvvvv-----|OOOO|-------,
%  o--|           C2 = 2 nF                    |---o
%     `-------------||-------------------------´
%  U = 1000 V
%  f = proměnná pro vyšetření charakteristiky

f = logspace(1,5,30)  %generuj logaritmiockou řadu hodnot frekvence (Hz)

%výpočet obvodu - násobení a dělení po složkách ( ./  .*)
om = 2*pi*f
Z1 = 1./(i*om*100e-9) + i.*om*0.01 + 1000
Z2 = 1./(i*om*2e-9)
Z = Z1.*Z2./(Z1 + Z2)
I = 1000./Z

%zobrazení grafu
figure(1)   %otevře okno
subplot(2,1,1)  %v okně budou 2 grafy + přepne vykreslování do prvního
semilogx(f, abs(I),'r*-') %graf s logarimickou osou X
xlabel('f (Hz)')  %popisky os
ylabel('|I| (A)')
grid on  %zapnutí mřížky
zoom on

subplot(2,1,2)
semilogx(f, angle(I)*180/pi,'b*-')
xlabel('f (Hz)')
ylabel('fáze (°)')
grid on
zoom on