t=0:0.0000001:0.000005;

y=1.0001*exp(-1999800*t)-0.0001*exp(-200*t);

plot(t,y),grid,xlabel('t'),ylabel('y(t)'),title('charakterystyka czasowa')

aa=-2;

bb=10;

n=1000;

w = logspace(aa,bb,n);

b=(w.^4+399960000*w.^2);

c=((w.^2+1999800^2).*(w.^2+200^2));

P=b./c;

d=(2*10^(6)*w.^3);

e=((w.^2+1999800^2).*(w.^2+200^2));

Q=d./e;

plot(P,Q),grid,xlabel('P(w)'),ylabel('Q(w)'),title('charakterystyka amplitudowo fazowa')