shiyan31.m

来自「数字信号处理课程的有关实验」· M 代码 · 共 25 行

    wp=0.2*pi;
    ws=0.3*pi;
    Rp=1;
    Rs=15;
    Ts=0.02*pi;
    Fs=1/Ts;
%选择滤波器的最小阶数
    [N,Wn]=buttord(wp,ws,Rp,Rs,'s');
    %创建butterworth模拟滤波器
    [Z,P,K]=buttap(N);
    %把滤波器零极点模型转化为传递函数模型
    [Bap,Aap]=zp2tf(Z,P,K);
    %把模拟滤波器原型转换成截至频率为Wn的低通滤波器
    [b,a]=lp2lp(Bap,Aap,Wn);
    %用双线性变换法实现模拟滤波器到数字滤波器的转换
    [bz,az]=bilinear(b,a,Fs);
    %绘制频率响应曲线
    [H,W]=freqz(bz,az);
    plot(W*Fs/(2*pi),abs(H));
    grid
    xlabel('频率／Hz')
    ylabel('频率响应幅度')