实验二 利用DFT计算线性卷积_第1页
实验二 利用DFT计算线性卷积_第2页
实验二 利用DFT计算线性卷积_第3页
实验二 利用DFT计算线性卷积_第4页
实验二 利用DFT计算线性卷积_第5页
已阅读5页,还剩8页未读 继续免费阅读

下载本文档

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

文档简介

1y(n)=x(n)*h(n)=Σx(m)h(n-m)my(n)x(n)和h(n)的N点圆周卷积为:yN(n)=x(m)h((n-m))RN(n)完成卷积运算,要求N=2y。采用补零的办法是x(n)和h(n)的长度均为N。(2)计算x(n)和h(n)的N点FFTFFTx(n)喻X(k)FFTh(n)喻H(k)Y(k)=X(k)H(k)Y(k)—我们考察单位取样响应为h(n)的线性系统,输入为x(n),输出为y(n),则y(n)=x(n)*h(n)2当输入序列时再开始进行卷积,会使输出相对输入有较大的延时,再者如果序列太长,需要大量的存储单元。为此,我们把,分别求出每段的卷积,合在一起其到最后的总输出。重叠保留法:设x(n)的长度为NX,h(n)的长度为M。我们把序列x(n)分成多段N点序列xi(n),每段与前一段重叠M-1个样本。由于第一段没有前一段保留信号,为了修正,我们在第一个输入段前面填充M-1个零。计算每一段h(n)的圆周卷积,则其每段卷积结果(1)在x(n)前面填充M-1个零,一一一—M-1个^FFTxi(n)喻Xi(k)FFTh(n)喻H(k)Y(k)=Xi(k)H(k)IFFTY(k)喻yi(n)(4)舍去每一段卷积结果的前M-1个样本,连接剩下样本,得到卷积结果y(n)。iM数量级相同。以x(n)表示每段信号,则i伪3伪x(n)*y(n)=xi(n)*h(n)每一段卷积yi(n)的长度为L+M-1,所以在做求和时,相邻两段序列有M-1个样本不重叠部分共同组成输出y(n)。i(1)将x(n)分为若干L点子段x(n)。i(2)计算每一段与h(n)的卷积yi(n),根据快速卷积算法利用FFT计算卷积。伪y(n)=yi(n-iL)4h=[cos(0.2*pi*n2),zeros(1h=cos(0.2*pi.*ones(15h=cos(0.2*pi.*ones(1x_temp=x(((count-1)*subLen+1):(count*6y(((count-1)*subLen+1):(count*subLen+M-1))=y(((count-1)*subL结论:与第二问的结果进行比较,发现运用重叠相加x_temp=x_bar(((count-1)*useLen+1):(count-1)*useLeny(((count-1)*useLen+1):(count*useLen))=y_tempx_temp=x_bar((count*us

温馨提示

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

评论

0/150

提交评论