p308.m

天气雷达的地物杂波的模拟及其仿真

%program p308
wp=2000*2*pi;
ws=3000*2*pi;
Rp=3;
Rs=15;
Fs=10000;
Nn=128;
[N,Wn]=buttord(wp,ws,Rp,Rs,'s')%compute order and cutoff freq.
[z,p,k]=buttap(N);%compute the analog filter
[Bp,Ap]=zp2tf(z,p,k);
[b,a]=lp2lp(Bp,Ap,Wn);
[bz,az]=impinvar(b,a,Fs)
freqz(bz,az,Nn,Fs)