www.gusucode.com > matlab通信工程仿真源码（张德丰等编著）程序书籍 > matlab_code/matlab通信工程仿真源码（张德丰等编著）/第9章/li9_2.m

    clear all;
f_N=8000;         %采样率
f_p=1000;f_s=700;  %设计要求指标
R_p=3;R_s=20;
Ws=f_s/(f_N/2);Wp=f_p/(f_N/2);   % 计算归一化频率
[n,Wn]=cheb1ord(Wp,Ws,R_p,R_s)   %计算阶数和截止频率
[b,a]=cheby1(n,R_p,Wn,'high');       %计算阶数和截止频率
freqz(b,a,1000,8000);    %作出H(z)的幅频相频图
subplot(2,1,1);axis([0 4000 -30 3]);
subplot(2,1,1);xlabel('频率Hz');ylabel('幅度dB');
subplot(2,1,2);xlabel('频率Hz');ylabel('相解rad');