性');
subplot(2,2,3),plot(angle(YY3));xlabel('时间/s');ylabel('相位');title('blackman窗函数滤波相位特性');
subplot(2,2,4),plot(angle(YY4));xlabel('时间/s');ylabel('相位');title('kaiser窗函数滤波相位特性');
3.分别设计hanning hamming blackman kaiser窗函数highpass_FIR function highpassfilter clc;
clear all;
Fs=100;%采样频率
fs=35;%高通阻带模拟截止频率 fp=40;%高通通带模拟起始频率 ws=2*pi*fs/Fs; wp=2*pi*fp/Fs; wn=(wp+ws)/2/pi; Bt=wp-ws;
N0=ceil(55*pi/Bt); N=N0+mod(N0+1,2); %调用窗函数
window1=hanning(N); window2=hamming(N); window3=blackman(N);
[n,Wn,beta,ftype]=kaiserord([35,40],[0 1],[0.01 0.01],100); window4=kaiser(n+1,beta); %设计加窗函数fir1
b1=fir1(N-1,wn,'high',window1); b2=fir1(N-1,wn,'high',window2); b3=fir1(N-1,wn,'high',window3);
b4=fir1(n,Wn,'high',window4 ,'noscale'); %求取频率响应
[H1,W1]=freqz(b1,1,512,2); [H2,W2]=freqz(b2,1,512,2); [H3,W3]=freqz(b3,1,512,2); [H4,W4]=freqz(b4,1,512,2);
figure(1);
subplot(2,2,1),plot(W1,20*log10(abs(H1)));%绘制频率响应图形 axis([0,1,-100,100]);
title('高通hanning窗的频率响应图形'); xlabel('频率（Hz)'); ylabel('幅值');
subplot(2,2,2),plot(W2,20*log10(abs(H2)));%绘制频率响应图形 axis([0,1,-100,100]);
title('高通hamming窗的频率响应图形'); xlabel('频率（Hz)'); ylabel('幅值');
subplot(2,2,3),plot(W3,20*log10(abs(H3)));%绘制频率响应图形 axis([0,1,-100,100]);
title('高通blackman窗的频率响应图形'); xlabel('频率（Hz)'); ylabel('幅值');
subplot(2,2,4),plot(W4,20*log10(abs(H4)));%绘制频率响应图形 axis([0,1,-100,100]);
title(' 高通kaiser窗的频率响应图形'); xlabel('频率（Hz)'); ylabel('幅值');
T=1/Fs;
L=100;%信号长度
t=(0:L-1)*T;%定义时间范围和步长
y=sin(2*pi*5*t)+5*sin(2*pi*15*t)+8*sin(2*pi*40*t);%滤波前的图形 NFFT = 2^nextpow2(L); % Next power of 2 from length of y Y = fft(y,NFFT)/L;%将时域信号变换到频域 f = Fs/2*linspace(0,1,NFFT/2+1);%频域采样 figure