利用FFT实现快速卷积(杭电)_第1页
利用FFT实现快速卷积(杭电)_第2页
利用FFT实现快速卷积(杭电)_第3页
利用FFT实现快速卷积(杭电)_第4页
利用FFT实现快速卷积(杭电)_第5页
已阅读5页,还剩1页未读 继续免费阅读

下载本文档

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

文档简介

1、信号、系统与信号处理实验实验报告实验名称:利用FFT实现快速卷积一、 实验目的1.加深理解FFT在实现数字滤波(或快速卷积)中的重要作用,更好地利用FFT进行数字信号处理。2.进一步掌握圆周卷积和线性卷积两者之间的关系。二、 实验内容与要求1、给定两个序列x(n)=2,1,1,2,h(n)=1,-1,-1,1。首先直接在时域计算两者的线性卷积;其次用FFT快速计算二者的线性卷积,验证结果。2、数字滤波器的脉冲响应为 ,N可自定,本实验取N=17,输入序列x(n)可选下列几种情况:,N可自取16,N=16, N=163、实验前,预先编制一个应用FFT实现数字滤波器的通用程序。4、上机独立调试,并

2、打印或记录实验结果。5、将实验结果与预先笔算的结果比较,验证其正确性。三、 实验程序与结果1、 给定两个序列x(n)=2,1,1,2,h(n)=1,-1,-1,1。首先直接在时域计算两者的线性卷积;其次用FFT快速计算二者的线性卷积,验证结果。clear all;n=2,1,1,2;m=1,-1,-1,1;y1n=conv(n,m);ny1=0:1:length(y1n)-1;stem(ny1,y1n);title(利用conv函数求线性卷积)figure; N=length(n)+length(m)-1;NK=fft(n,N);MK=fft(m,N);YK=NK.*MK;y2n=ifft(Y

3、K,N);x=0:N-1;stem(x,y2n);title(利用fft方法求线性卷积)2、数字滤波器的脉冲响应为 ,N可自定,本实验取N=17,输入序列x(n)可选下列几种情况:,N可自取16 ,N=16 , N=16主程序:clear all;n=0:1:15;m=0:1:16; xn1=ones(1,16);xn2=cos(2*n*pi/16);xn3=(1/3).n;hn=(-1/2).n;y1=ccc(xn1,hn);y2=ccc(xn2,hn);y3=ccc(xn3,hn);定义函数实现利用FFT实现数字滤波器:function yf = ccc(xn,hn)N1=length(xn);N2=length(hn);N=N1+N2-1;XK=fft(xn,N);HK=fft(hn,N);YK=XK.*HK;yf=ifft(YK,N);if all (imag(xn)=0)&(all(imag(hn)=0) yf=real(y

温馨提示

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

评论

0/150

提交评论