声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1198|回复: 1

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

[复制链接]
发表于 2009-3-23 18:06 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
本帖最后由 wdhd 于 2016-9-14 10:49 编辑

  [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 编辑 ]
带阻滤波器.JPG
双线性带阻滤波器1.JPG
加噪声前后频谱.JPG
带阻滤波器.JPG
回复
分享到:

使用道具 举报

 楼主| 发表于 2009-3-24 09:30 | 显示全部楼层
怎么没有人回帖字啊?希望高手多帮忙啊
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-26 12:29 , Processed in 0.058990 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表