clear; close all;

%% 
K = 0.1;
T1 = 1;
T2 = 0.001;
T3 = 0.1

% (K(T1p + 1)(T2p + 1))/(T2p + 1)

N = 200;
W10min = -1;
W10max = 4;

%% À×Õ è Ô×Õ ïî ÏÔ
w = logspace(W10min, W10max, N);

LMR = (K.*(w).*sqrt( (1 - (T1*T2).*(w.^2)).^2 + ((T1 + T2).*w).^2 ))./sqrt( (1 - (T1.*T2)));
LMR = 20*log10(LMR);

LPR =pi./2 + atan(w) + atan(0.001.*w) - atan(0.1.*w); 
LPR( LPR > 0 ) = LPR( LPR > 0) - pi;
LPR = LPR * 180/pi;