ht3.asv

matla窗口制作

function ht3(h_axes1,N,mm,b1,Ws,Wp,Wph,Wsh,c)
N=floor(N);
%delta_w=2*pi/1000;
%Ap=-(min(db(Wp*pi/delta_w+1:1:501)));
%syms c;
if mm==1 
    n=[0:1:N-1];
    wc=((Ws+Wp)/2)*pi;
    Wcl=wc;
    Wch=((Wph+Wsh)/2)*pi;
    if b1==1
        hd=ideal_lp1(wc,N);
    elseif  b1==2
        hd=ideal_hp1(wc,N);
    elseif  b1==3
        hd=ideal_bp1(Wcl,Wch,N);
    elseif  b1==4
        hd=ideal_bs(Wcl,Wch,N);
    end
    h=(bartlett(N))';
    h=hd.*h;
    [dp,mag,pha,grd,w]=freqz_m(h,[1]);
 elseif mm==2
   n=[0:1:N-1];
    wc=((Ws+Wp)/2)*pi;
    Wcl=wc;
    Wch=((Wph+Wsh)/2)*pi;
    if b1==1
        hd=ideal_lp1(wc,N);
    elseif  b1==2
        hd=ideal_hp1(wc,N);
    elseif  b1==3
        hd=ideal_bp1(Wcl,Wch,N);
    elseif  b1==4
        hd=ideal_bs(Wcl,Wch,N);
    end
    h=(blackman(N))';
    h=hd.*h;
   [dp,mag,pha,grd,w]=freqz_m(h,[1]);
elseif mm==3
    n=[0:1:N-1];
    wc=((Ws+Wp)/2)*pi;
    Wcl=wc;
    Wch=((Wph+Wsh)/2)*pi;
    if b1==1
        hd=ideal_lp1(wc,N);
    elseif  b1==2
        hd=ideal_hp1(wc,N);
    elseif  b1==3
        hd=ideal_bp1(Wcl,Wch,N);
    elseif  b1==4
        hd=ideal_bs(Wcl,Wch,N);
    end
    h=(hamming(N))';
    h=hd.*h;
    [dp,mag,pha,grd,w]=freqz_m(h,[1]);
elseif mm==4
    n=[0:1:N-1];
    wc=((Ws+Wp)/2)*pi;
    Wcl=wc;
    Wch=((Wph+Wsh)/2)*pi;
    if b1==1
        hd=ideal_lp1(wc,N);
    elseif  b1==2
        hd=ideal_hp1(wc,N);
    elseif  b1==3
        hd=ideal_bp1(Wcl,Wch,N);
    elseif  b1==4
        hd=ideal_bs(Wcl,Wch,N);
    end
    h=(hanning(N))';
    h=hd.*h;
    [dp,mag,pha,grd,w]=freqz_m(h,[1]);
end 
    axes(h_axes1);
    cla;
    if c==1
    plot(w/pi,dp,'b');
    axis([0,1,-200,10]);
    grid on;
   % set(C,'linewidth',1.5);
    set(h_axes1,'Color',[1 1 1]);
    title(' Window Magnitude Response');
    xlabel('Normalized Frequency(*pi rad/sample)');
    ylabel('Magnitude(dB)');
    elseif c==2
    axes(h_axes1);
    cla;
    %[dp,db,mag,pha,grd,w]=freqz_m2(h,[1]);
    plot(w/pi,pha/pi,'b');
   % axis([0,1,-200,10]);
    grid on;
    set(h_axes1,'Color',[1 1 1]);
    title(' Window phase Response');
    xlabel('Frequency(*pi rad/sample)');
    ylabel('phase(pi)');
    elseif c==3
    axes(h_axes1);
    cla;
    stem(n,h);
   % axis([0,1,-200,10]);
    grid on;
    set(h_axes1,'Color',[1 1 1]);
    title('Impulse Response');
    xlabel('');
    ylabel('phase(pi)');
    end
end