ht7.m

matla窗口制作

function ht7(h_axes1,N,ww1,ww2,ww3,Ws,Wp,Wph,Wsh,b3,c)
n=floor(N);
if b3==1
    f=[0 Ws Wp 1];
    s=[ww1 ww2];
    b=cremez(n,f,'lowpass',s);
    [dp,mag,pha,grd,w]=freqz_m(b,[1]);
elseif b3==2
    f=[0 Ws Wp 1];
    s=[ww1 ww2];
    b=cremez(n,f,'highpass',s);
    [dp,mag,pha,grd,w]=freqz_m(b,[1]);
elseif b3==3
    f=[0 Ws Wp Wph Wsh 1];
    s=[ww1 ww2 ww3];
    b=cremez(n,f,'bandpass',s);
    [dp,mag,pha,grd,w]=freqz_m(b,[1]);
elseif b3==4
    f=[0 Ws Wp Wph Wsh 1];
    s=[ww1 ww2 ww3];
    b=cremez(n,f,'bandstop',s);
    [dp,mag,pha,grd,w]=freqz_m(b,[1]);
end
    axes(h_axes1);
    cla;
    if c==1
    plot(w/pi,dp,'b');
    axis([0,1,-200,10]);
    grid on;
%    set(h_axes1,'Color',[1 1 1]);
    title('Remez-expansion Magnitude Response');
    xlabel('Normalized Frequency(*pi rad/sample)');
    ylabel('Magnitude(dB)');
    elseif c==2
    axes(h_axes1);
    cla;
    plot(w/pi,pha/pi,'b');
    grid on;
    title('Remez-expansion phase Response');
    xlabel('Frequency(*pi rad/sample)');
    ylabel('phase(pi)');
    end
end