% 《數字信號處理教程——MATLAB釋義與實現》第八章模擬濾波器頻帶變換演示程序
% 電子工業出版社出版  陳懷琛編著 2004年9月
% lvboqi.m    從模擬低通濾波器到模擬高通，帶通，帶阻濾波器的演示! 

% 在本演示程序中，用的為巴特沃斯濾波器，階數為 N=5（應該根據實際數計算出來，具體的計算公式和方法見上面的參考書）
% 如果需要根據實際的要求來設計相應的濾波器的話，所設計的低通濾波器的參數在每一個子函數中進行修改
% 因為對于不同的條件,所要顯示的時間長度--T,  時間間隔---dt;  通帶截止頻率--fc;  阻帶截止頻率fs 
% 及在設計帶通,帶阻濾波器時中心頻率--f0    歸一化后都是不一樣的

handles=figure('name','模擬濾波器設計演示');
set(handles,'unit','normalized','position',[0.1,0.2,0.5,0.5]);
set(handles,'defaultuicontrolunits','normalized');

%坐標軸y一用來繪制低通濾波器
haxes1=axes('position',[0.05,0.55,0.7,0.4],...
             'ylim',[0,2]);               
                  
%在這里繪制出一般的模擬低通濾波器的圖形，及其一些設計參數
p1=[0:0.1:3]*j; fp=1;  fs=2;
for i=1:length(p1)
    hl1(i)=1/abs((p1(i)^2 + 0.6180*p1(i) + 1)*(p1(i)^2 + 1.6180*p1(i) + 1)*(p1(i)+1));
end
h=[fliplr(hl1),hl1];
p=[-fliplr(p1),p1];
plot(p/j,h);axis([-3,3,0,1.2]);hold on
plot([0,0],[0,1.2],'linestyle','-','color','r');
plot([fp,fp],[0,hl1(fp/0.1+1)],'linestyle','-.','color','g');
plot([-fp,-fp],[0,hl1(fp/0.1+1)],'linestyle','-.','color','g');
plot([fs,fs],[0,hl1(fs/0.1+1)],'color','r');
plot([-fs,-fs],[0,hl1(fs/0.1+1)],'color','r');
plot([-3,fp],[0.707,0.707],'linestyle','-.','color','g');


haxes2=axes('position',[0.05,0.1,0.7,0.4]);

hpushbutton1=uicontrol(handles,'style','pushbutton',...
                       'position',[0.78,0.6,0.2,0.1],...
                       'string','模擬低通到高通',...
                      'callback',['call_l2hp4(haxes1,haxes2);']);
hpushbutton2=uicontrol(handles,'style','pushbutton',...
                       'position',[0.78,0.4,0.2,0.1],...
                       'string','模擬低通到帶通',...
                      'callback',['call_l2bp5(haxes1,haxes2);']);
hpushbutton3=uicontrol(handles,'style','pushbutton',...
                       'position',[0.78,0.2,0.2,0.1],...
                       'string','模擬低通到帶阻',...
                      'callback',['call_l2bs6(haxes1,haxes2);']);
hpushbutton4=uicontrol(handles,'style','pushbutton',...
                       'position',[0.78,0.8,0.2,0.1],...
                       'string','模擬低通到低通',...
                      'callback',['call_l2lp(haxes1,haxes2);']);