ht6.m

来自「matla窗口制作」· M 代码 · 共 26 行

function ht6(h_axes1,N,ww1,ww2,ww3,Ws,Wp,Wph,Wsh,b3)
n=floor(N);         %奇数要报错（高通和带阻）
if b3==1
    f=[0 Ws Wp 1];
    s=[ww1 ww2];
    b=cremez(n,f,'lowpass',s);
    [h,w]=freqz(b);
elseif b3==2
    f=[0 Ws Wp 1];
    s=[ww1 ww2];
    b=cremez(n,f,'highpass',s);
    [h,w]=freqz(b);
elseif b3==3
    f=[0 Ws Wp Wph Wsh 1];
    s=[ww1 ww2 ww3];
    b=cremez(n,f,'bandpass',s);
    [h,w]=freqz(b);
elseif b3==4
    f=[0 Ws Wp Wph Wsh 1];
    s=[ww1 ww2 ww3];
    b=cremez(n,f,'bandstop',s);
    [h,w]=freqz(b);
end