% 《數(shù)字信號處理教程——MATLAB釋義與實現(xiàn)》第四章演示程序
%  電子工業(yè)出版社出版  陳懷琛編著 2004年9月
%  數(shù)字濾波器頻帶變換演示
%  call_dlp2lp.m 數(shù)字低通濾波器到數(shù)字低通濾波器

function []=call_dlp2lp(h1,h2)

axes(h1);hold off;
plot([0,0],[-1.2,1.2],'linestyle','-.','color','g');hold on
plot([-1.2,1.2],[0,0],'linestyle','-.','color','g');axis([-1.2,1.2,-1.2,1.2]);
axes(h2);hold off;
plot([0,0],[-1.2,1.2],'linestyle','-.','color','g');hold on
plot([-1.2,1.2],[0,0],'linestyle','-.','color','g');axis([-1.2,1.2,-1.2,1.2]);

w0=0.2*pi;    %低通濾波器的截止頻率(指過渡濾波器，而非最終要設(shè)計的濾波器，它的截止頻率在wl（dtransform.m中)設(shè)置
w=-pi:0.1:pi;  %在這里可以更改其間隔,如果改成0:0.1:2*pi,則下面的判斷要改成abs(w(i))<w0  | abs(w(i))>2*pi-w0
zlp=exp(j*w);

[zhp0,zhp1]=dtransform(w,'lp2lp');


    for i=1:length(w)
        %繪制低通濾波器
        axes(h1);
        if abs(w(i))<w0
            plot(zlp(i),'linestyle','o','color','r');axis([-1.2,1.2,-1.2,1.2]);axis('equal');hold on
        else
            plot(zlp(i),'linestyle','o','color','b');axis([-1.2,1.2,-1.2,1.2]);axis('equal');hold on
        end
        %繪制高通濾波器
        axes(h2);
        if abs(w(i))<w0
            plot(zhp0(i),'linestyle','o','color','r');axis([-1.2,1.2,-1.2,1.2]);axis('equal');hold on
        else
            plot(zhp0(i),'linestyle','o','color','b');axis([-1.2,1.2,-1.2,1.2]);axis('equal');hold on
        end
        pause(0.2);
    end