clc
clear all
close all
alpha = 0.1;
syms s w
%w = logspace(0,5,501);
%s = i*w;

%prvi blok normirano
br1=[1 0];
naz1=s+2.1056;
naz1=sym2poly(naz1);
h1=tf(br1,naz1)
%drugi blok normirano
br2=[1 0 0];
naz2=s^2+1.5559*s+2.02476;
naz2=sym2poly(naz2);
h2=tf(br2,naz2)
%treci blok normirano
br3=[1 0 0]
naz3=s^2+0.316282*s+1.07753;
naz3=sym2poly(naz3)
h3=tf(br3,naz3)

h=h1*h2*h3

%bode normirano
w1n=logspace(-1,1,1001); %od dekade ispod wg do dekade iznad wg
bode(h,h1,h2,h3,w1n)
grid on
legend('H(s)','H1(s)','H2(s)','H3(s)')
%w2n=logspace(0,1,1001); %od wg do dekade iznad
%bode(h,h1,h2,h3,w2n)
%grid on
%legend('H(s)','H1(s)','H2(s)','H3(s)')

hwabs=vpa((sqrt((w^10 + 7.536*w^8 + 7.06*w^6)^2 + (3.978*w^9 - 9.885*w^7 + 4.594*w^5)^2))/(w^10 + 30.897*w^8 - 7.734*w^6 + 240.67*w^4 - 40.98*w^2 + 21.105));
pretty(hwabs)

%phiw = atan((3.978*w^9 )/())
[br,naz] = tfdata(h,'v')
br=subs(poly2sym(br),'x',s)
naz=subs(poly2sym(naz),'x',s)
hsym=vpa(br/naz);
pretty(hsym)
pretty(simplify(ilaplace(hsym)))

%denormirana

denh=subs(hsym,s,(s/(2*pi*300)));
denh=simplify(denh);
fprintf('\nDENORMIRANA\n')
pretty(vpa(denh))
%za bode potrebno posebno denormirati brojnik i nazivnik zbog tf-a
br=subs(br,s,(s/(2*pi*300)))
naz=subs(naz,s,(s/(2*pi*300)))
denbode=tf(sym2poly(br),sym2poly(naz))
bode(denbode)
xlim([10^2.27 10^4.27])