ht4.m

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

function ht4(h_axes1,N,w1,w2,up1,up2,up3,b2)
n=floor(N);
if b2==1
    f=[0 w1 1];
    amp=[1 0];
    up=[(up1+1) up2];                    %假设通带波纹和阻带波纹的上沿和下沿对称
    lo=[((up1+1)-2*up1) -up2];
    b=fircls(n,f,amp,up,lo);
    [h,w]=freqz(b);
elseif b2==2
    f=[0 w1 1];
    amp=[0 1];
    up=[up2 (up3+1)];
    lo=[-up2 ((up3+1)-2*up3)];
    b=fircls(n,f,amp,up,lo);
    [h,w]=freqz(b);
elseif b2==3
    f=[0 w1 w2 1];
    amp=[0 1 0];
    up=[up1 (up2+1) up3];
    lo=[-up1 ((up2+1)-2*up2) -up3];
    b=fircls(n,f,amp,up,lo);
    [h,w]=freqz(b);
elseif b2==4
    f=[0 w1 w2 1];
    amp=[1 0 1];
    up=[(1+up1) up2 (1+up3)];
    lo=[((1+up1)-2*up1) -up2 ((1+up3)-2*up3)];
    b=fircls(n,f,amp,up,lo);
    [h,w]=freqz(b);                 %要改一下h,w变量的参数名字   记住！
end