求助模拟滤波器到数字滤波器

lianmei1984

　　[x1,fs,Nbits]=wavread('121.wav');
　　fs=8000;
　　sound(x1,fs); %播放语音信号
　　y1=fft(x1,1024); %对信号做1024点FFT变换
　　f=fs*(0:511)/1024;
　　plot(abs(y1(1:512))) %做原始语音信号的FFT频谱图
　　加噪声t=0:1/fs:(size(x1)-1)/fs;
　　Au=0.03;
　　d=[Au*cos(2*pi*2000*t)]';
　　x2=x1+d;
　　sound(x2,fs)%播放加噪声后的语音信号
　　y2=fft(x2,1024);
　　figure(1)
　　plot(t,x2)
　　title('加噪后的信号');
　　xlabel('time n');
　　ylabel('fuzhi n');
　　figure(2)
　　subplot(2,1,1);
　　plot(f,abs(y1(1:512)));
　　title('原始语音信号频谱');
　　xlabel('Hz');
　　ylabel('fuzhi');
　　subplot(2,1,2);
　　plot(f,abs(y2(1:512)));
　　title('加噪后的信号频谱');波形见附件3
　　请问高手如何设计滤波器滤掉噪声阿。应该设计带阻滤波器，我这有个归一化设计数字带阻滤波器的程序
　　f0=2000/4000;
　　>> fc=100/4000;
　　>> ws=[f0-0.5*fc f0+0.5*fc];
　　>> wp=[f0-fc f0+fc];
　　>> Ap=1;
　　>> As=30;
　　>> [N,wc]=buttord(wp,ws,Ap,As,'s');
　　>> [num,den]=butter(N,wc,'stop');
　　>> [H,W]=freqz(num,den,fs); plot(W*fs/(2*pi),20*log10(abs(H)))波形见附件4
　　请问高手如何设计模拟滤波器滤掉噪声阿?需要设计模拟带阻滤波器的程序，希望高手帮忙啊
　　[ 本帖最后由 lianmei1984 于 2009-3-23 18:31 编辑 ]

lianmei1984

怎么没有人回帖字啊？希望高手多帮忙啊