《带通滤波器设计》word版_第1页
《带通滤波器设计》word版_第2页
《带通滤波器设计》word版_第3页
《带通滤波器设计》word版_第4页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1、.以前我做的一个完整的信号滤波前后的程序，可能数据有些不一样，但用MATLAB做滤波器的大体思路都差不多，希望对你有用。（完全自创，要给分哦）f=8*103; %信号频率f1=30*103; %噪声频率l=100; %信号的长度fs=100*103; %采样频率%生成信号t=(0:l-1)/fs;s=0.5*sin(2*pi*f*t)+cos(2*pi*f1*t);figure(1);subplot(1,2,1);plot(fs*t(1:60),s(1:60); title('时域波形');xlabel('时间');ylabel('幅度');gr

2、id on;nf=2nextpow2(l);Y=fftshift(fft(s,nf);%k=-fs/2+(0:nf-1)*fs/nf; k=fs*linspace(0,1,nf)-fs/2;subplot(1,2,2);plot(k,abs(Y); %信号的频谱图title('信号频谱');xlabel('频率');ylabel('幅度');grid on;%滤波器设定ws=10*103; %通带截止频率ws1=20*103; %阻带起始频率ds=50; %阻带衰减wt=2*pi*ws/fs; %经采样后的通带截止角频率wz=2*pi*ws1/f

3、s; %阻带的wc=(wt+wz)/2; %归一化后的滤波器截止频率n=ceil(6.6*pi/(wz-wt); %t=(n-1)/2; q=fir1(n-1,wc/pi,hanning(n);%滤波器时域函数w=linspace(0,pi,512);h1=freqz(q,1,512); %进行512个点的傅里叶变换figure(2);plot(w/pi,20*log10(abs(h1);title('滤波器频谱图');xlabel('频率');ylabel('幅度');grid ;%经滤波器后的信号signal=conv(s,q);nx=(0:n+l-2)/fs;figure(3);subplot(1,2,1);plot(nx(n:60)*fs,signal(n:60);title('经滤波后的信号时域图');xlabel('时间');ylabel('幅度');grid on;lx=2nextpow2(nx);Y=fftshift(fft(signal,lx);k=-fs/2+(0:lx-1)*fs/lx; %subplot(1,2,2);plot(k,abs(Y); %信号的频谱图title('滤波后

人人文库> 全部分类> 应用文书 > 工作计划

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论