www.pudn.com > programradarsystemdisign.zip > fig7_9.m

clear all 
fofr = 0:0.001:1; 
arg = 2.*pi.*fofr; 
nume = 2.*(1.-cos(arg)); 
den11 = (1. + 0.25 * 0.25); 
den12 = (2. * 0.25) .* cos(arg); 
den1 = den11 - den12; 
den21 = 1.0 + 0.7 * 0.7; 
den22 = (2. * 0.7) .* cos(arg); 
den2 = den21 - den22; 
den31 = (1.0 + 0.9 * 0.9); 
den32 = ((2. * 0.9) .* cos(arg)); 
den3 = den31 - den32; 
resp1 = nume ./ den1; 
resp2 = nume ./ den2; 
resp3 = nume ./ den3; 
plot(fofr,resp1,'k',fofr,resp2,'k-.',fofr,resp3,'k--'); 
xlabel('Normalized frequency') 
ylabel('Amplitude response') 
legend('K=0.25','K=0.7','K=0.9') 
grid 
axis tight