用matlab实现序列的圆周卷积.docx

数字信号处理实验报告 实验项目名称： 用MATLAB 实验项目名称： 用MATLAB实现序列的| 周卷积 实验日期： 20124J28 实验成绩： 实验评定标准: 1)实验结果是否正确 A ( ) B ( ) C () 2)实验结果分析 A ( ) B ( ) C () 3)实验报告是否按照规定格式 A ( ) B ( ) C () 一、 实验目的 通过本实验，掌握一些基本而且重要的离散时间信号，熟悉基本离散时间 信号的MATLAB实现方法。 二、 实验器材 PC 机，MATLAB 软件。 三、 实验内容 计算两序列 xl (n) ={1,23,4,5}, x2 (n) ={1,2,3,4,5,4,3,2,1}的圆周卷积。 四、 实验结果 实验代码： clear all n =O:N-1 close all n l=0:N-2; clc n 2=0:N-3; xl司1,2,3,4,5,6,7,8]; yl=circonvt(xl,x2,N); x2=[l,2,3,4,5,6,7,8,7,6,5,4,3,2,1]; y2=circonvt(xl,x2,N-l); N=length(xl)+length(x2); y3=circonvt(xl,x2,N-2); xl=[xl zeros(l,N-length(xl))]; x2=[x2 zeros(l,N-length(x2))]; Xfl=dft(xl,N); Xf2=dft(x2?N); Xf=Xfl.*Xf2; x=idft(Xf,N); x=real(x); subplot(2,3,l) stem(n’xl); title(xl(n)); subplot(2,3,2) stem(n,x2); title(x2(n)) function y=circonvt(xl/x2/N) if length(xl)N error(N必须 xxl的长度) end if length(x2)N error(N必须=x2的长度) end xl=[xl zeros(l,N-length(xl))]; x2=[x2 zeros(l,N-length(x2))]; m=[O:l:N-l]; x2=x2(mod(-m,N)+l); H=zeros(N,N); for n=l:l:N H(n,:)=cirshift(x2,n-1,N); end y=xl*H; function y=cirshift(x/m/N) subplot(2/3,3); stem(n,x); title(x(n)=IDFT(X(k))); subplot(23/4); stem(n,yl); title(N点圆周卷积J; subplot(2,3,5); stem(nl,y2); title(N-l点圆周卷积)； subplot(2,3,6); stem(n2,y3); title(N-2点圆周卷积)； if length(x)N error(N必须=x的长度) end x=[x zeros(l,N-length(x))]; n=[O:l:N-l]; n=mod (n-m,N); y=x(n+l); function (Xk]=dft(xn,N) n=[O:l:N-l];k=[O:l:N-l]; WN=exp(-j*2*pi/N); nk=n**k; WNnk=WN.Ank; Xk= xn * WNnk; function [xn]=idft(Xk,N) %计算逆离散傅里叶变换 %[xn]=idft(Xk,N) n=[O:l:N-l]; k=[O:l:N-l]; WN=exp(-j*2*pi/N); nk=n*k; 实验结果：  WNnk=WN.A(-nk); xn=(Xk*WNnk)/N; xi(n) x2(n xjni-DFTO^n 冷卷円 点 IB 崗◎帜 实验结果分析-