fc=500;
fc1=fc/2000;%即截止頻率為0.25，最大頻率為1000hz，采樣頻率為2000hz
BW=0.05;%即sinc濾波器的降落帶?50hz
M=4/BW;%即濾波器內核長度
w=zeros（1M+1）;
h=zeros（1M+1）;
X1=zeros（1M+1）;
X2=zeros（1M+1）;
%a=zeros（M1）;
y=zeros（M+11）;
?t=0:0.0005:0.0005*M;
?a=1*sin（2*pi*500*t）+1*sin（2*pi*1000*t）+1*sin（2*pi*1100*t）;
?b=sin（2*pi*500*t）;
?%plot（a）
for?i=1:M+1
????w（i）=0.42-0.5*cos（2*pi*i/M）+0.08*cos（4*pi*i/M）;%布萊克曼窗
??