example01.m

matla窗口制作

%___________________________________________________利用fircls函数，自行设计窗函数，对信号进
%___________________________________________________行滤波(带通滤波器）____________________
t=0:1/600:1;                            % 信号的采样频率为600hz,(w=2*pi/fs)
sig=sin(2*pi*15*t)+0.5*sin(2*pi*90*t)+0.2*sin(2*pi*200*t);
n=50;                                  %窗的宽度决定所设计的滤波器的过度带的宽度
f=[0 0.2 0.4 1];                        %设置窗函数的四个临界频率
a=[0 1 0];                              %确定滤波器的类型
up=[0.02 1.02 0.01];                    %确定滤波器的上限参数
lo=[-0.02 0.98 -0.01];                  %确定滤波器的下限参数
b=fircls(n,f,a,up,lo);                  %返回其单位脉冲响应
newsig=fftfilt(b,sig);                  %通过b滤波器对向量sig进行滤波
ft=t(300:350);                          %后面四行表示对滤除后的信息进行处理，让其能在规定的采样频率中显示出来
ns=newsig(300:350); 
zns=interp(ns,12);                      %在300到350之间每一中插入12个数据，以至于最终有600个数据，使得向量匹配
znt=interp(ft,12);
%_________________________________绘制图象_____________________________________
subplot(321);
plot(t,sig);
title('原始信号');
%——————————————————————————————————————————
subplot(322);
[H,W]=freqz(sig);
plot(W/pi,abs(H));
title('原始信号频谱');
grid
%——————————————————————————————————————————
subplot(323);
stem(b);
title('窗函数的脉冲响应');
%——————————————————————————————————————————
subplot(324);
[H,W]=freqz(b);
plot(W/pi,abs(H));
title('窗函数的频谱');
grid
%—————————————————————————————————————————
subplot(325);
plot(znt,zns);
title('滤波后的信号');
subplot(326);
%__________________________________________________________________________
[H,W]=freqz(newsig);
plot(W/pi,abs(H));
title('滤除后的频谱');
grid