DSP“数字信号处理”第四章lzh.pptVIP

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
例:画出4点按时间抽取的FFT流图,并用流图计算下列序列的4点DFT X(k)。 §4.10 线性卷积与线性相关的FFT算法 以FIR滤波器为例,它的输出等于有限长单位脉冲响应h(n)与有限长输入信号x(n)的离散线性卷积 y(n)也是有限长序列,其点数为L+M-1点 一、线性卷积的FFT算法 设x(n)为L点,h(n)为M点,输出y(n)为 1) 为了不产生混叠,其必要条件是使x(n),h(n)都补零值点,补到N≥M+L-1, 即: 0≤n≤L-1 L≤n≤N-1 0≤n≤M-1 M≤n≤N-1 2) 然后计算圆周卷积 FFT算法就是用圆周卷积来代替线性卷积 这时,y(n)就代表线性卷积的结果 N N 用FFT计算y(n)的步骤如下 ① 计算 H(k)=DFT[h(n)],N点 ② 计算 X(k)=DFT[x(n)], N点 ③ 计算 Y(k)=X(k)H(k); ④ 计算 y(n)=IDFT[Y(k)],N点 当x(n)的点数很多时,即当LM,需要采用分段卷积或称分段过滤的办法计算卷积。 原因:不能等x(n)全部采集后再进行卷积;否则,使输出相对于输入有太长的延时,需要太多的存储空间; 此外,若N=L+M-1太大,h(n)必须补很多个零值点,很不经济,这时FFT法的优点就表现不出来了。 将x(n)分成点数和h(n)相仿的段,分别求出每段的卷积结果,然后用某种方式把它们合在一起,得到总的输出,其中每一段的卷积均采用FFT方法处理。 有两种分段卷积的办法: 重叠相加法和重叠保留法 iL≤n≤(i+1)L-1 其他n i=0, 1, … 则输入序列可表示成: 用xi(n)表示x(n)的第i 段: 重叠相加法: 设h(n)的点数为M,信号x(n)为很长的序列。将x(n)分为很多段,每段为L点,L选择和M的数量级相同 每一个yi(n)=xi(n)*h(n)都可用上面讨论的快速卷积方法计算。 N L L L 0 L 0 0 L M-1 N 由于xi(n)为L点,而yi(n)为L+M-1点,故相邻两段输出序列必然有(M-1)个点发生重叠,即前一段的后(M-1)个点和后一段的前(M-1)个点相重叠。应该将重叠部分相加再和不重叠的部分共同组成输出y(n)。 线性卷积的特点是:头、尾各有(M-1)长的过渡过程 因此,将x(n)分段后,其每段的卷积结果yi(n)都不能完全和相应的y(n)相等,需要把上一段的后过渡过程和本段的前过渡过程对应相加,才能得到完整的y(n) 重叠相加法图形 h ( n ) 0 N -1 M -1 x ( n ) 0 L 2 L 3 L n n n n n L -1 0 x 0 ( n ) N -1 0 x 1 ( n ) L 2 L -1 L + N -1 3 L -1 0 x 2 ( n ) 2 L 2 L + N -1 重叠相加法图形 ⑤ 将各段yi(n)(包括重叠部分)相加 用FFT法实现重叠相加法的步骤如下 ① 计算N点FFT,H(k)= DFT[h(n)];M→N ② 计算N点FFT,Xi(k)=DFT[xi(n)];L→N ③ 相乘,Yi(k)=Xi(k)H(k); ④ 计算N点IFFT,yi(n)=IDFT[Yi(k)]; 重叠相加是由于各输出段的重叠部分相加而得名的 在重叠相加法中yi(n) 可以完全恢复出y(n)的每一点 由于xi(n)为N=L+M-1点序列,这时用FFT实现h(n)和xi(n)的N点圆周卷积,圆周卷积结果的前(M-1)个值发生混叠,不等于线性卷积值,必须舍去。 即用每段yi(n)的后L个值,首尾相接构成y(n)。 此时得到的y(n)与x(n)的长度(点数)相同。 重叠保留法: 与重叠相加法比较: 相同之处:先将x(n)分段,每段L个点。 不同之处:序列中补零处不补零,而是在每一段的前边补上前一段最后(M-1)个值,组成L+M-1点序列xi(n) 0 M-1 M-1 M-1 N M-1 L L L L N 重叠保留法示意图 重叠保留法示意图 重叠保留法的名称是因为每一组相继的输入段均由L个新数据和前一段保留下来的(M-1)个旧数据所组成而得名的。 在重叠保留法中yi(n) 不能恢复出y(n)的后过渡部分 例 x (n)={3,4,2,5,6,7,3,2,8,3,4,5,4,2,3,7,6,9,2,2,4,5,3,1} x 0(n)={0,0,0,0,3, 4, 2, 5, 6, 7, 3, 2} x 1(n)={6,7,3,2,8, 3, 4, 5, 4, 2, 3, 7} x 2(n)={4,2,3,7,6, 9,

文档评论(0)

1243595614 + 关注
实名认证
文档贡献者

文档有任何问题,请私信留言,会第一时间解决。

版权声明书
用户编号:7043023136000000

1亿VIP精品文档

相关文档