C = 3e-12; % kapacita
L = 30e-9; % indukcnost
f = logspace([8, 10], 200); % log. frekvencni osa
omega = 2*Pi*f; % uhlova frekvence
% vypocet prenosu
H = (-omega.^2*LC) .../
   [omega.^2*L*C*(3-omega^2*L*C)+1];
H_dB = 20*log_10(abs(H)); % prenos v dB
clear
semilogx(H_dB, f) % logaritmicky graf
grid on % zapnuti mrizky
xlabel('f (Hz)') % popis x-ove osy
Ylabel('|H| (dB)') % popis y-ove osy
