% canceler frequency characteristic
clear
rad=pi/180;
fmhz2=[2.48:.0001:2.52];
fmhz1=[2.4:.001:2.6];
fc=2.5*1000;
fp=50e3;
Tp=1/fp;
arg1=pi*fmhz1*1e6*Tp;
H1=4*sin(arg1).^2;
figure(1)
plot(fmhz1*1000-fc,H1,'-k','LineWidth',2)
grid
axis([min(fmhz1*1000)-fc,max(fmhz1*1000)-fc,0,4])
ylabel('|H(f)|')
xlabel('Doppler Frequency, kHz')
arg2=pi*fmhz2*1e6*Tp;
H2=sin(arg2).^2;
Hdb=20*log10(H2);
figure(2)
plot(fmhz2*1000-fc,Hdb,'k-','LineWidth',2)
grid
axis([min(fmhz2*1000)-fc,max(fmhz2*1000)-fc,-40,0])
ylabel('20*log_1_0(|H(f)|_n_o_r_m)')
xlabel('Doppler Frequency, kHz')
