a2d_filter.m

很多matlab的源代码

h0=figure('position',[10 10 800 500],...
    'name','IIR滤波器的数模转换');
 
wp = 0.2*pi;
ws = 0.3*pi;
Rp = 1;
As = 15;
T = 1;
OmegaP = wp/T;
OmegaS = ws/T;
[cs,ds] = afd_butt(OmegaP,OmegaS,Rp,As);
[b,a] = imp_invr(cs,ds,T);
[C,B,A] = dir2par(b,a)
[db,mage,pha,grd,w] = freqz_m(b,a);
subplot(2,3,4);
plot(w/pi,db);grid
xlabel('w');
ylabel('dB');
axis([0 1 -45 5])
subplot(2,3,1);
plot(w/pi,mage);grid
xlabel('w');
ylabel('|H|');
axis([0 1 0 1.1])
subplot(2,3,2);
plot(w/pi,pha);grid
xlabel('w');
ylabel('相位');
axis([0 1 min(pha) max(pha)]);
subplot(2,3,5);
plot(w/pi,grd);grid
xlabel('w');
ylabel('样本');


Ts_text=uicontrol('style','text',...
    'position',[600 460 100 20],...
    'string','Method');

Ts_1=uicontrol('style','radiobutton',...
    'position',[600 440 100 20],...
    'string','脉冲响应法',...
    'value',1,...
    'callback',[...
        'set(Ts_1,''value'',1),',...
        'set(Ts_2,''value'',0),',...
        'group_id=10,',...
        'position=group_id+type,',...
        'mfile_name=strcat(''a8'',int2str(position)),',...
        'eval(mfile_name);']);


Ts_2=uicontrol('style','radiobutton',...
    'position',[600 420 110 20],...
    'string','双线性变换法',...
    'value',0,...
    'callback',[...
        'set(Ts_1,''value'',0),',...
        'set(Ts_2,''value'',1),',...
        'group_id=16,',...
        'position=group_id+type,',...
        'mfile_name=strcat(''a8'',int2str(position)),',...
        'eval(mfile_name);']);

type_text=uicontrol('style','text',...
    'position',[600 380 100 20],...
    'string','Prototype Filter');

type_1=uicontrol('style','radiobutton',...
    'position',[600 360 100 20],...
    'string','Butterworth',...
    'value',1,...
    'callback',[...
        'set(type_1,''value'',1),',...
        'set(type_2,''value'',0),',...
        'set(type_3,''value'',0),',...
        'set(type_4,''value'',0),',...
        'type=1;',...
        'position=group_id+type,',...
        'mfile_name=strcat(''a8'',int2str(position)),',...
        'eval(mfile_name);']);
type_2=uicontrol('style','radiobutton',...
    'position',[600 340 100 20],...
    'string','Chebyshev-I',...
    'value',0,...
    'callback',...
        ['set(type_1,''value'',0),',...
        'set(type_2,''value'',1),',...
        'set(type_3,''value'',0),',...
        'set(type_4,''value'',0),',...
        'type=2;',...
        'position=group_id+type,',...
        'mfile_name=strcat(''a8'',int2str(position)),',...
        'eval(mfile_name);']);

type_3=uicontrol('style','radiobutton',...
    'position',[600 320 100 20],...
    'string','Chebyshev-II',...
    'value',0,...
    'callback',[...
        'set(type_1,''value'',0),',...
        'set(type_2,''value'',0),',...
        'set(type_3,''value'',1),',...
        'set(type_4,''value'',0),',...
        'type=3;',...
        'position=group_id+type,',...
        'mfile_name=strcat(''a8'',int2str(position)),',...
        'eval(mfile_name);']);

type_4=uicontrol('style','radiobutton',...
    'position',[600 300 100 20],...
    'string','Elliptic',...
    'value',0,...
    'callback',[...
        'set(type_1,''value'',0),',...
        'set(type_2,''value'',0),',...
        'set(type_3,''value'',0),',...
        'set(type_4,''value'',1),',...
        'type=4;',...
        'position=group_id+type,',...
        'mfile_name=strcat(''a8'',int2str(position)),',...
        'eval(mfile_name);']);

Hc_close=uicontrol(gcf,'style','push',...
   'position',[650 80 80 25],...
   'string','Close',...
   'Callback','close');
Hc_info=uicontrol(gcf,'style','push',...
'position',[650 150 80 25],...         
'string','info',...                 
'Callback','web http://www.bupt.edu.cn -browser');                 


%A2D_filter_sub;

group_id=10;
type=1;