MATLAB-FILTER.rar_butter_filter_mobilcommuniction

% 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');