%  Some simple transfer function analysis
%
%
%
c = 1.2;
a = 0.1;
rp = 1.0;
dp = 0.1;
rv = 5.0;
dv = 4.0;
c = 1.0;
a = 0.1;
d = 0.1;
w = 0.1*(1:10000);
%cutter = ones(1,10000);
lw = log(w);
s = i*w;
G = a*c./(s - (rp-dp));
F = c./(s - (rv - dv));
LD = exp(-i*w*d);
xfer = G./(1 + (F.*LD).*G);
subplot(2,1,1);
plot(lw,abs(xfer));
phase = atan(imag(xfer)./real(xfer));
subplot(2,1,2);
plot(lw,phase);