sd4 zad4

% Paweł Pytlowski

% czyszczenie workspace'u i otwartych wykresów
clc;
close all;
clear;


f_gr = 1000; % podana przez prowadzacego częstotliwość graniczna (1 kHz)

mnoznik = 60; % znaleziony doświadczalnie mnożnik RC
              % dla którego f graniczna odpowiada magnitutde = -20 dB

%f_gr = 1/(2*pi*R*C)
RC = mnoznik * 1/(2*pi*f_gr);

% L(s)/M(s) = G(s)
Licznik = [1];
Mianownik = [RC, 1];

G = tf(Licznik, Mianownik);

figure;
hold on;

% zaznaczam punkt dla częstotliwości granicznej na czerwono
bode(G,f_gr,'ro');
legend('f graniczna');

% wyznaczam diagram Bodego dla G
bode(G,{1e1,1e4});
title(['Diagram Bodego dla RC = ' num2str(RC)]);
hold off;