function[]=Filter1();
[x,fs,bits]=wavread('bluetooth_input.wav');%播放原始信號
N=length(x);%返回采樣點數
t=(1:N)/fs;
df=fs/N;%采樣間隔
n1=1:N/2;
f=(n1-1)*df;%頻帶寬度
y1=awgn(x,45);%加入高斯白噪聲，信噪比為45
%*****************************布拉克曼濾波器設計*************************
Window=blackman(35); %長度為35的布拉克曼窗Window
b=fir1(34,0.1,Window);%產生低通布拉克曼濾波器
H=freqz(b,1,512);
figure(6);
plot(abs(H));
title('布拉克曼濾波器的頻率響應');
%**********************************************************************
ys=fftfilt(b,y1);%信號送入濾波器濾波，ys為輸出
fftwave=fft(ys);%將濾波后的語音信號進行快速傅立葉變換
figure(7);
plot(f,20*log10(abs(fftwave(n1)))),grid on;
title('布拉克曼濾波器濾波后信號的頻譜圖');
xlabel('頻率/Hz');
ylabel('幅值/db');
grid;
wavwrite(ys,fs,'bluetooth_input.wav');%把濾波后的信號寫入d:\cc1.wav中
[g,fs,bits]=wavread('bluetooth_input.wav');%語音回放
wavplay(g,fs);