Solution 10.8.9.7a The MATLAB program w=logspace(-2,3,200);; s=j*w;; K =78 z1 = 3 p1 = 0 p2 = 1 p3 = 18 mag = 20*log10( abs((K.*( +z1))./((s+p1).*(s+p2).*(s+p3))));; phase=(angle(s+z1)-angle(s+p1)-angle(s+p2)-angle(s+p3))*180/pi;; % % 0.5 dB % M=10^(0.5/20);; c=-M^2/(M^2 - 1);; r =abs(M/(M^2-1));; theta = linspace(0,2*pi,360);; x=c+r*cos(theta);; y=r*sin(theta);; s=x+j*y;; gamma = angle(s)*180/pi;; b=zeros(1,360);; c=gamma > b;; gamma = gamma -360*c;; rho = 20*log10(abs(s));; G=[gamma'];; R=[rho'];; % % 1dB % M=10^(1/20);; c=-M^2/(M^2 - 1);; r =abs(M/(M^2-1));; theta = linspace(0,2*pi,360);; x=c+r*cos(theta);; y=r*sin(theta);; s=x+j*y;; gamma = angle(s)*180/pi;; b=zeros(1,360);; c=gamma > b;; gamma = gamma -360*c;; rho = 20*log10(abs(s));; 1 G=[G gamma'];; R=[R rho'];; M=1.0593;; c=-M^2/(M^2 - 1);; r =abs(M/(M^2-1));; theta = linspace(0.00001,2*pi,360);; x=c+r*cos(theta);; y=r*sin(theta);; s=x+j*y;; gamma = angle(s)*180/pi;; b=zeros(1,360);; c=gamma > b;; gamma = gamma -360*c;; rho = 20*log10(abs(s));; G=[G gamma'];; R=[R rho'];; % % 2dB % M=10^(2/20);; c=-M^2/(M^2 - 1);; r =abs(M/(M^2-1));; theta = linspace(0.00001,2*pi,360);; x=c+r*cos(theta);; y=r*sin(theta);; s=x+j*y;; gamma = angle(s)*180/pi;; b=zeros(1,360);; c=gamma > b;; gamma = gamma -360*c;; rho = 20*log10(abs(s));; G=[G gamma'];; R=[R rho'];; % % 3dB % M=10^(3/20);; c=-M^2/(M^2 - 1);; r =abs(M/(M^2-1));; theta = linspace(0.00001,2*pi,360);; x=c+r*cos(theta);; y=r*sin(theta);; s=x+j*y;; gamma = angle(s)*180/pi;; b=zeros(1,360);; 2 c=gamma > b;; gamma = gamma -360*c;; rho = 20*log10(abs(s));; G=[G gamma'];; R=[R rho'];; % % 4dB % M=10^(4/20);; c=-M^2/(M^2 - 1);; r =abs(M/(M^2-1));; theta = linspace(0.00001,2*pi,360);; x=c+r*cos(theta);; y=r*sin(theta);; s=x+j*y;; gamma = angle(s)*180/pi;; b=zeros(1,360);; c=gamma > b;; gamma = gamma -360*c;; rho = 20*log10(abs(s));; G=[G gamma'];; R=[R rho'];; % % 6dB % M=10^(6/20);; c=-M^2/(M^2 - 1);; r =abs(M/(M^2-1));; theta = linspace(0.00001,2*pi,360);; x=c+r*cos(theta);; y=r*sin(theta);; s=x+j*y;; gamma = angle(s)*180/pi;; b=zeros(1,360);; c=gamma > b;; gamma = gamma -360*c;; rho = 20*log10(abs(s));; G=[G gamma'];; R=[R rho'];; % % 9dB % M=10^(9/20);; c=-M^2/(M^2 - 1);; r =abs(M/(M^2-1));; 3 theta = linspace(0.00001,2*pi,360);; x=c+r*cos(theta);; y=r*sin(theta);; s=x+j*y;; gamma = angle(s)*180/pi;; b=zeros(1,360);; c=gamma > b;; gamma = gamma -360*c;; rho = 20*log10(abs(s));; G=[G gamma'];; R=[R rho'];; % % 12 dB % M=10^(12/20);; c=-M^2/(M^2 - 1);; r =abs(M/(M^2-1));; theta = linspace(0.00001,2*pi,360);; x=c+r*cos(theta);; y=r*sin(theta);; s=x+j*y;; gamma = angle(s)*180/pi;; b=zeros(1,360);; c=gamma > b;; gamma = gamma -360*c;; rho = 20*log10(abs(s));; G=[G gamma'];; R=[R rho'];; plot(phase,mag,'k-',G,R,'k--') grid on axis([-360 0 -10 30]) print -deps 10897a.eps Draws the plot shown in Figure 1. As can be seen the logmagnitude plot just touches the 2 dB circle. Then 10 2=20 =1:2589: Based on this calculation wewould expect the maximum overshoot to be about 25%, whichiscloseto what it actually is. 4 -350 -300 -250 -200 -150 -100 -50 0 -10 -5 0 5 10 15 20 25 30 0.5 dB 1 dB 2 dB 3 dB 4 dB 6 dB 9 dB 12 dB Figure 1: Large scale polar plot 5