ch3example2a.m

清华大学出版社 邵玉斌编写的《通信系统建模与仿真实例分析》一书的所有MATLAB和SIMULINK代码

% ch3example2A.m
f_N=8000;                          % 采样率
f_p=2100; f_s=2500; R_p=3; R_s=25; % 设计要求指标
Ws=f_s/(f_N/2); Wp=f_p/(f_N/2);    % 计算归一化频率
[n, Wn]=buttord(Wp,Ws,R_p,R_s);    % 计算阶数和截止频率
[b,a]=butter(n, Wn);               % 计算H(z)
figure(1);
freqz(b,a, 1000, 8000)             % 作出H(z)的幅频相频图, freqz(b,a, 计算点数, 采样率)
subplot(2,1,1); axis([0 4000 -30 3])
figure(2);                              % 第二种作图方法
f=0:40:4000;                            % 计算频率点和频率范围
z=exp(j*2*pi*f./(f_N));                 % 
H_z=polyval(b,z)./polyval(a,z);         % 计算相应频率点处H(s)的值
subplot(2,1,1); plot(f, 20*log10(abs(H_z))); % 幅频特性
axis([0 4000 -40 1]);
xlabel('频率 Hz');ylabel('幅度 dB');
subplot(2,1,2); plot(f, angle(H_z));         % 相频特性
xlabel('频率 Hz');ylabel('相角 rad');