实验线性卷积与圆周卷积的计算
实验二 DFT(FFT)的应用—利用FFT实现快速卷积
姓名：高铭遥 班级：16131701 学号：1120171450 成绩：实验二 DFT/FFT 的应用-利用FFT 实现快速卷积[实验目的]1．深刻理解DFT/FFT 的概念和性质，进一步掌握圆周卷积和线性卷积两者之间的关系。
2．掌握DFT/FFT 的应用。
理解FFT 在实现数字滤波（或快速卷积）中的重要作用，更好地利用FFT 进行数字信号处理。
[实验内容及要求]1．给定两个序列()[]2,1,1,2x n =，()[]1,1,1,1h n =--。
首先直接在时域计算两者的线性卷积；然后用FFT 快速计算二者的线性卷积，验证结果。
（1）线性卷积 程序代码：figure(1);N1=4; N2=4; xn=[2,1,1,2]; hn=[1,-1,-1,1];N=N1+N2-1;%卷积后的序列长度 yn=conv(xn,hn);%线性卷积 x=0:N-1;stem(x,yn);title('线性卷积'); 运行结果：（2）FFT 卷积快速卷积 程序代码： figure(1); n=0:1:3; m=0:1:3;N1=length(n);%xn 的序列长度 N2=length(m);%hn 的序列长度 xn=[2,1,1,2]; hn=[1,-1,-1,1];姓名：高铭遥 班级：16131701 学号：1120171450 成绩：N=N1+N2-1;%卷积后的序列长度XK=fft(xn,N);%xn 的离散傅里叶变换 HK=fft(hn,N);%hn 的离散傅里叶变换 YK=XK.*HK;yn=ifft(YK,N);%逆变换if all(imag(xn)==0)&&(all(imag(hn)==0))%实序列的循环卷积仍为实序列 yn=real(yn); endx=0:N-1;stem(x,yn);title('FFT 卷积'); 运行结果：结果分析：对比（1）和（2）直接线性卷积和FFT 快速卷积的结果可以验证，用FFT 线性卷积的结果是与直接卷积的结果相同的，FFT 可以实现快速卷积，提高运算速度。
线性卷积与圆周卷积的计算
数字信号处理实验报告
实验名称：线性卷积与圆周卷积的计算（实验四）
专业班 级： 学生姓名： 学 号： 指导教师：
二○一 年 月 日
1
一、 实验原理 1.线性卷积
当系统输入序列为 x(n),系统的单位冲激响应为 h(n),输出序列为 y(n),则线性时不变系统 y(n) = x(n)*h(n). 2.圆周卷积
圆周移位代替线性移位的好处： 时域圆周卷积在频域上相当于两序列的 DFT 的相乘，而计算 DFT 可 以采用它的快速算法——快速傅立叶变换（FFT），因此圆周卷积和线 性卷积相比，计算速度可以大大加快。
7
调试中所遇到的问题： 图形输出与算得的不符，经过结果逆推回去，发现计算圆周卷积的一 个矩阵没有转置。 圆周卷积与线性卷积的关系： 根据实验结果，可以发现，当 N>=L+P-1 时，圆周卷积等于线性卷积； 而当 N<L+P-1 时，圆周卷积等于两个序列的线性卷积加上时间的混 叠。
6
线性卷积的运算步骤： 求 x1(n)与 x2（n） 的线性卷积：对 x1(m)或 x2（m）先进行镜像移 位 x1（-m），对移位后的序列再进行从左至右的依次平移 x(n-m),当 n=0,1,2.…N-1 时，分别将 x(n-m)与 x2（m）相乘，并在 m=0,1,2.… N-1 的区间求和，便得到 y（n）。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
圆周卷积通用程序
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function yc = circonv(x1,x2,N) % 定义函数 if length(x1)>N
线性卷积和圆周卷积的关系
1、线性卷积和圆周卷积的关系（参考书P122）dsp31：ppt91设)(1n x 是1N 点的有限长序列)10(1-≤≤N n ，设)(2n x 是2N 点的有限长序列)10(2-≤≤N n ，（A ）)(1n x 和)(2n x 的线性卷积为∑∑-=∞-∞=-=-=1021211)()()()()(N m m l m n x m x m n x m x n y 则线性卷积)(n y l 的长度为121-+N N 。
（B ）两个有限长序列)(1n x 和)(2n x 做L 点的圆周卷积：首先将两个序列补零，扩成长度为L 的序列： ⎩⎨⎧-≤≤-≤≤=1,010),()(1111L n N N n n x n x ⎩⎨⎧-≤≤-≤≤=1,010),()(2222L n N N n n x n x 圆周卷积为：)(]))(()([)(1021n R m n x m x n y L L m L ⋅-=∑-=这里必须将一个序列变成L 点周期延拓序列，这里采用)(2n x 序列：∑∞-∞=+==r L rL n x n x n x )())(()(~222 把它带入到)(n y 中并考虑到前面的线性卷积公式，可得到：)(])()([)(]))(()([)(10211021n R m rL n x m x n R m n x m x n y L L m r L L m L ⋅-+=⋅-=∑∑∑-=∞-∞=-=)(])()([2101n R m rL n x m x Lr L m ⋅-+=∑∑∞-∞=-= )(])([n R rL n y L r l ⋅+=∑∞-∞=所以L 点圆周卷积)(n y 是线性卷积)(n y l 以L 为周期的周期延拓序列的主值序列。
因为)(n y l 有121-+N N 个非零值，所以延拓周期L 必须满足：121-+≥N N L 。
这时各延拓周期才不会交叠，而)(n y 的前121-+N N 个值正好是)(n y 的全部非零序列值，也正是线性卷积)(n y l 。
实验四 有限长序列的线性卷积、圆周卷积及分段卷积(数字信号处理)
电子信息与自动化学院《数字信号处理》实验报告学号： 姓名：实验名称： 实验四 有限长序列的线性卷积、圆周卷积及分段卷积一、 实验目的(1) 在理论学习的基础上，通过本实验，加深对线性卷积、圆周卷积、分段卷积的理解；(2) 掌握计算线性卷积、圆周卷积、分段卷积的方法；(3) 体会有限长序列卷积运算的关系；二、 实验原理1、有限长序列卷积有两种形式：线性卷积和圆周卷积然而现实中要解决的实际问题是要计算两个有限长序列的线性卷积，如信号通过线性系统，系统的输出 y(n)是输入信号 x(n)与系统抽样响应 h(n)的线性卷积：y(n)=x(n)*h(n)。
设n x 1和n x 2是两个长度分别为 M 和 N 的有限长序列，则其线性卷积为)(*)()(211n x n x n y =。
)(1n y 是一个长度为 L1=N+M-1 点的有限长序列．将n x 1和n x 2均补零成 L 点的有限长序列，其中 L ≥max(M,N),则其 L 点的圆周卷积为)(]))(()([)()()(1021212n R m n x m x n x n x n y L L m L ∑-=-=⊗=，现在讨论)(1n y 和)(2n y 的关系。
显然]∑∑∑∑∑∑∑∑∞-∞=∞-∞=∞∞=-=-=∞-∞=-=-=+=+=+-=+-=-=-=r r L r M m L M m L r L M m L L L m L rL n y n R rL n x n xn R rL m n x n R rL m n x m x n R m n x m x n R m n x m x n y )([)()](*)([)()()()()()(]))(()([)(]))(()([)(121102102112110212由此可见，L 点的圆周卷积)(2n y 是线性卷积)(2n y 以 L 为周期，进行周期延拓后在区间 0 到 L-1 范围内所取的主值序列。
数字信号处理实验线性卷积圆周卷积
数字信号处理实验线性卷积圆周卷积⼤连理⼯⼤学实验报告学院（系）：电信专业：⽣物医学⼯程班级：**1101姓名：**** 学号：201181*** 组：___实验时间：实验室：实验台：指导教师签字：成绩：实验⼀线性卷积和圆周卷积⼀、实验程序1.给出序列x=[3,11,7,0,-1,4,2]，h=[2,3,0,-5,2,1]；⽤两种⽅法求两者的线性卷积y，对⽐结果。
a)直接调⽤matlab内部函数conv来计算。
b)根据线性卷积的步骤计算。
clear;clc;x=[3 11 7 0 -1 4 2];n1=0:1:length(x)-1;h=[2 3 0 -5 2 1];n2=0:1:length(h)-1;y=conv(x,h);n3=0:1:length(x)+length(h)-2;figure(1);subplot(121);stem(n1,x,'.');axis([0 6 -15 15]);title('x(n)序列');grid;subplot(122);stem(n2,h,'.');axis([0 5 -10 10]);title('h(n)序列');grid;figure(2);subplot(121);stem(n3,y,'.');axis([0 12 -60 60]);title('调⽤conv函数的线性卷积后序列');grid;N=length(x);M=length(h);L=N+M-1;for(n=1:L)y1(n)=0;for(m=1:M)k=n-m+1; if(k>=1&k<=N)y1(n)=y1(n)+h(m)*x(k); end; end; end;subplot(122);stem(n3,y1,'*');axis([0 12 -60 60]);title('按步骤计算的线性卷积后序列');grid; 结果2.卷积后结果y=[ 6 , 31 , 47 , 6 , -51 , -5 , 41 , 18 , -22 , -3 , 8 , 2]。
实验二 利用DFT计算线性卷积
实验二 利用DFT 计算线性卷积一、实验目的1.掌握利用FFT 计算线性卷积的原理及具体实现方法。
2.加深理解重叠相加法和重叠保留法。
3.考察利用FFT 计算线性卷积各种方法的适用范围。
二、实验设备与环境计算机、MATLAB 软件环境三、实验基础理论1、线性卷积与圆周卷积设)(n x 为L 点序列，)(n h 为M 点序列，)(n x 和)(n h 的线性卷积为：∑∞-∞=-==m m n h m x n h n x n y )()()(*)()()(n y 的长度为：L+M-1，)(n x 和)(n h 的N 点圆周卷积为：)())(()()(10n R m n h m x n y N N m N -=∑-= 其中：1-+≥L M N此时圆周卷积等于线性卷积，而圆周卷积可利用FFT 计算。
2、快速卷积利用FFT 计算线性卷积步骤如下：（1）为了使线性卷积可以利用圆周卷积来计算，必须选择；同时为了能使用基2-FFT 完成卷积运算，要求γ2=N 。
采用补零的办法是)(n x 和)(n h 的长度均为N 。
（2）计算)(n x 和)(n h 的N 点FFT)()()()(k H n h k X n x FFTFFT →→（3）组成卷积 )()()(k H k X k Y =（4）利用IFFT 计算IDFT ，得到线性卷积y(n)(k)()IFFT Y y n −−−→3、分段卷积我们考察单位取样响应为)(n h 的线性系统，输入为)(n x ，输出为)(n y ，则)(*)()(n h n x n y =当输入序列时再开始进行卷积，会使输出相对输入有较大的延时，再者如果序列太长，需要大量的存储单元。
为此，我们把，分别求出每段的卷积，合在一起其到最后的总输出。
这种方法称为分段卷积。
分段卷积可细分为重叠相加法和重叠保留法。
重叠保留法：设)(n x 的长度为X N ，)(n h 的长度为M 。
我们把序列)(n x 分成多段N 点序列)(n x i ，每段与前一段重叠M-1个样本。
线性卷积与圆周卷积的计算(杭电)
信号、系统与信号处理实验Ⅱ实验报告实验名称：线性卷积与圆周卷积的计算一、实验目的1、通过编程，上机调试程序，进一步增强使用计算机解决问题的能力。
2、掌握线性卷积与圆周卷积软件实现的方法，并验证两者之间的关系。
二、实验内容与要求已知两个有限长序列：x(n)= δ(n)+2δ(n-1)+3δ(n-2)+4δ(n-3)+5δ(n-4)；h(n)= δ(n)+2δ(n-1)+δ(n-2)+2δ(n-3)1.编制一个计算两个线性卷积的通用程序，计算x(n)*h(n)。
2.编制一个计算圆周卷积的通用程序，计算上述4种情况下两个序列x(n)与h(n)的圆周卷积。
3.上机调试并打印或记录实验结果。
4.将实验结果与预先笔算的结果比较，验证真确性。
三、实验程序与结果1、计算两个线性卷积的通用程序，计算x(n)*h(n)。
xn=[1 2 3 4 5]hn=[1 2 1 2]N=length(xn);M=length(hn);L=N+M-1;for(n=1:L)y(n)=0;for(m=1:M)k=n-m+1;if(k>=1&k<=N)y(n)=y(n)+hn(m)*xn(k);endendendy=conv(xn,hn);ny=0:L-1;stem(ny,y) ;xlabel('n ');ylabel('y(n) ');figurestem(ny,yn) ;xlabel('n ');ylabel('y ');根据定义编写循环实现线性卷积结果：01234567n y (n )Conv 函数实现线性卷积结果：01234567n y2. 计算圆周卷积的通用程序，计算上述4种情况下两个序列x(n)与h(n)的圆周卷积。
主程序：clear allN=[5 6 9 10];xn=[1 2 3 4 5];hn=[1 2 1 2];yc1=circonv(xn,hn,N(1))yc2=circonv(xn,hn,N(2))yc3=circonv(xn,hn,N(3))yc4=circonv(xn,hn,N(4))figurestem(0:N(1)-1,yc1);xlabel('时间序号n');ylabel('信号幅度');title('5点圆周卷积');figurestem(0:N(2)-1,yc2);xlabel('时间序号n');ylabel('信号幅度');title('6点圆周卷积');figurestem(0:N(3)-1,yc3);xlabel('时间序号n');ylabel('信号幅度');title('9点圆周卷积');figurestem(0:N(4)-1,yc4);xlabel('时间序号n');ylabel('信号幅度');title('10点圆周卷积');定义函数：function yc=circonv(x1,x2,N)if length(x1)>Nerror('N必须大于等于x1的长度'); endif length(x2)>Nerror('N必须大于等于x2的长度'); endx1=[x1,zeros(1,N-length(x1))];x2=[x2,zeros(1,N-length(x2))];n=[0:N-1];x2=x2(mod(-n,N)+1);H=zeros(N,N);for n=1:1:NH(n,:)=cirshiftd(x2,n-1,N);yc=x1*H';function y=cirshiftd(x,m,N)if length(x)>Nerror('x 的长度必须小于N');endx=[x,zeros(1,N-length(x))];n=[0:1:N-1];y=x(mod(n-m,N)+1);时间序号n 信号幅度5点圆周卷积00.51 1.52 2.533.54 4.55时间序号n 信号幅度时间序号n 信号幅度时间序号n 信号幅度四、仿真结果分析编写的线性卷积程序和conv 函数的结果相同，也与笔算结果相同。
关于线性卷积及圆周卷积的简便竖式法计算
关于线性卷积及圆周卷积的简便竖式法计算
线性卷积和圆周卷积是数字信号处理中常见的两种卷积操作。
简单来说，线性卷积可以把两个信号之间的关系映射到输出上，而圆周卷积是一种更为复杂的运算，它可以寻找两个旋转的信号之间的关系。
下面就描述一下这两种卷积的简便竖式法计算。
线性卷积：
输入：
f(n)=x(n)*h(n)
f：输入信号；
x：样本函数；
h：滤波器。
步骤：
（1）将输入信号f分段；
（2）用滤波器在f的每一段输入取值上乘以x；
（3）对f的每一段结果求和，最终得到f的线性卷积输出。
圆周卷积：
输入：
F(n)=X(n)*H(n)
F：输入信号；
X：变换函数；
H：滤波器。
步骤：
（1）将输入信号F分段，每一段变换为正弦、余弦等函数；
（2）对每一段变换后的函数，用滤波器H乘以X；
（3）对每一段变换后函数结果求叠加和，以得到F的圆周卷积输出。
总结：
上述简便竖式法计算描述了两种卷积的计算步骤，即线性卷积和圆周卷积，在结果求叠加和时，用来表示信号实际上与自身的旋转有关的圆周卷积结果是不同的。
因此，这两种卷积的计算采用的步骤也有所不同。
以上就是线性卷积及圆周卷积的简便竖式法计算的长文描述。
用matlab验证卷积定理
用matlab验证卷积定理
卷积定理
一、实验目的
通过本实验，验证卷积定理，掌握利用DFT和FFT计算线性卷积的方法。
二、实验原理
时域圆周卷积在频域上相当于两序列DFT的相乘，因而可以采用FFT的算
法来计算圆周卷积，当满足
121
L N N
≥+-时，线性卷积等于圆周卷积，因此可利用FFT计算线性卷积。
三、实验内容和步骤
1．给定离散信号()
x n和()
h n，用图解法求出两者的线性卷积和圆周卷积；2．编写程序计算线性卷积和圆周卷积；
3．比较不同列长时的圆周卷积与线性卷积的结果，分析原因。
四、实验设备
计算机、Matlab软件
五、实验报告要求
1．整理好经过运行并证明是正确的程序，并且加上详细的注释。
2．给出笔算和机算结果对照表，比较不同列长时的圆周卷积与线性卷积的结果对照，作出原因分析报告。
3．结出用DFT计算线性卷积的方法。
圆周卷积与周期卷积、线性卷积的关系与计算
圆周卷积与周期卷积、线性卷积的关系及计算一、三者关系设：1122()01()01x n n N x n n N ≤≤-≤≤-N ：圆周卷积的点数⏹ 圆周卷积是周期卷积的主值序列。
周期卷积：1120()()()N m y n x m x n m -==-∑ （1）圆周卷积：1120()()()[()(())]()N c N N N m y n y n R n x m x n m R n -===-∑1210[()(())]()N N N m x m x n m R n -==-∑ （2）注意：（2）式直接使用的前提是圆周卷积的点数N 应满足：12max[,]N N N ≥（一般题目均符合此种情况）⏹ 周期卷积是线性卷积的周期延拓。
线性卷积：1112120()()*()()()N l m y n x n x n x m x n m -===-∑212121()()()*()N m x m x n m x n x n -==-=∑ （4）圆周卷积与线性卷积的关系：()[()]()c l N r y n y n rN R n ∞=-∞=+∑ （5）注意：上述关系式对任意长度的圆周卷积均适合。
二、举例说明1、对于12max[,]N N N ≥的情况，各教材例题很多，不再举例。
2、12N N N N <<或的情况。
习题8.已知序列()()2(1)(4)3(5)x n n n n n δδδδ=+-+-+-,4()()y n R n =，求：（1）()()*()z n x n y n =（2）()()f n x n =○5()y n （5点圆周卷积）。
解：(){1,2,0,0,4,3},(){1,1,1,1}x n y n ==（1）()()(){1,3,3,3,3,4,4,4,3}z n x n y n =*=（过程略） （2）()()f n x n =○5()y n （5点圆周卷积），N =5。
线性卷积与圆周卷积演示程序的设计
1.设计基本原理1.1课题研究的背景卷积运算广泛的应用于通讯、电子、自动化等领域的线性系统的仿真、分析及数字信号处理等方面。
在MATLAB中可以使用线性卷积和圆周卷积实现离散卷积。
线性卷积是工程应用的基础，但圆周卷积实现线性离散卷积具有速度快等优势。
圆周卷积采用循环移位，在MATLAB中没有专用函数，需要根据圆周卷积的运算过程编制程序代码。
本实验主要围绕线性卷积和圆周卷积的演示程序设计来展开，给出了线性卷积和圆周卷积演示的程序及动态实现。
在线性时不变连续系统中，利用系统的冲激响应和叠加原理来求系统对任意激励信号作用时的零状态响应，这就是卷积方法的原理。
因此，在时域内，卷积运算是求解线性非时变系统零状态响应的重要方法，特别是激励信号为时限信号时尤其如此。
卷积运算的计算比较复杂，是信号与系统分析中的重点和难点，特别适合用于计算机来计算。
以往的卷积积分多用fortran、c、VB等语言编程，不仅编程繁琐，而且可视性差。
用MATLAB来计算卷积积分问题要比用C、FORTRAN 等语言完成相同的事情简洁的多。
在MATLAB中，有很多现成的函数可以直接调用，而且在计算机方面，可以直接用相应的计算机符号即可。
在编写程序语言方面，它与其他语言相比更为简单。
正因为上述原因，使他深受工程技术人员及科学专家的欢迎，并很快成为应用学科计算机辅助分析、设计、仿真、教学等领域不仅可缺少的基础软件。
1.2课题研究意义本课程为电子信息工程专业的独立实践课，是建立在信号与系统、数字信号处理等课程的基础上，加强实践环节而开设的。
其目的在于通过本课程设计使学生进一步巩固数字信号处理的基本概念、理论、分析方法和实现方法；使学生能有效地将理论和实际紧密结合；增强学生软件编程实现能力和解决实际问题的能力。
通过课程设计，主要达到以下的目的：（1）使学生增进对MATLAB的认识，加深对信号处理理论的理解。
（2）使学生掌握数字信号处理中频谱分析的概念和方法。
线性卷积和圆周卷积
x1 ( n)
DSP
1
0
x2 ( n)
N 1
n
1
0
x2 ((0 m)) N RN ( m)
N 1
n
1
x2 ((1 m)) N RN ( m)
1
N 1
m
0
y ( n) x1 ( n) N x2 ( n)
N 1
m
3 2 1 0
N 1
m
结论
DSP
设两序列为x(n)和h(n)，长度分别为N、M；则其 线形卷积的长度为N+M-1，而圆周卷积的长度为 K=Max（N，M）。实际上还可以计算K+1、 K+2、。。。、N+M-1、N+M、N+M+1….等点 的圆周卷积，只有L≥N1＋N2－1，则L点圆周卷 积和线性卷积相等。
DSP
圆周卷积过程
DSP
补零：若 x(n) 的长度是 N，h(n) 的长度是 M， 取，对序列补零至H点。 周期延拓：先在哑变量坐标 m 上作出 x(m) 和 h(m)，将h(m)周期延拓。 翻褶，取主值序列：对 h(m)以 m=0的垂直轴 为对称轴翻褶成h (-m)，然后取主值序列。 圆周移位：对得到的序列进行圆周移位。 相乘相加：与 x(m) 对应项相乘，并累加，得到 y(n)
DSP
线性卷积和圆周卷积
付海燕
实验目的
DSP
理解离散序列的线性卷积与圆周卷积的原理，比 较其相同和不同点 掌握线性卷积与圆周卷积的计算步骤和计算方法 能熟练使用Matlab的相关命令。
卷积和的运算步骤
DSP
翻 褶 ： 先 作 出 x(m) 和 h(m)， 将 h(m) 以 m=0的垂直轴为对称轴翻褶成h(-m)。 移位：将h(-m)移位，即得h(n-m)。当n为正 整数时，右移n位。当为负整数时，左移n位。 相乘：再将 h(n-m) 和 x(m) 的相同 m 值的对应 点值相乘 相加：把以上所有对应点的乘积叠加起来，即得 y(n)值。
线性卷积与圆周卷积
线性卷积与圆周卷积伦敦奥运会也将接近尾声，中国体育代表团获得的⾦牌屈居第⼆，貌似后⾯也没有⼏个夺⾦点了，不过这已经很不错了。
然⽽奥运期间⼏场⽐赛却让⼈难忘。
⼥双打“假球”事件惹来国内⼀⽚哗然，有骂她们的，也有为她们不值的，其实错也不应该让她们承担，奥运作为竟技⽐赛，不可能完全不考虑胜负⽽每场⽐赛却拼个你死我活，当然得考虑战略战术，谁愿意⼀上来就和本国另⼀队拼个你死我活，想想当年中国队与中国⾹港队⾜球在⼩组中为让中国队出线打成7:0，就觉得这确实没什么；林丹决赛和李宗伟那场巅峰对决确实扣⼈⼼弦，不愧是超级丹！中国⼥排进四强⽐赛中在2:3情况下不敌⽇本队，让⼈觉得挽惜，中国⼥排打得也不错，貌似第⼀局有点轻敌，更没有想到的是⽇本⼥排异常的顽强；刘翔摔倒那⼀刻估计很多⼈像会不禁惊叫出来，当时为了看刘翔⽐赛，提前离开教研室去吃饭，回来在正好在⼤屏幕上看到刘翔马上要⽐赛，⼤屏幕前已经站了很多同学，当时⼼想进决赛应该没问题，但摔倒那⼀刻许多⼈都表⽰失望，其实也不是怪刘翔，⽽是当希望变成失望时⼀时⽆法接受，因此刘翔承载了太多⼈的希望，也可能是由于⼤家让刘翔承载太多，让他折断飞翔的两翼.....继续关注奥运⽐赛。
话⼊正题，前两天看了下线性卷积和圆周卷积，并对他们之间的关系作出验证。
线性卷积与圆周卷积离散线性卷积的定义：设长度为N1的序列x(n)和长度为N2的序列h(n)进⾏线性卷积，得到长度为N1+N2-1的y(n):离散圆周卷积的定义：圆周卷积是定义在有限长序列之间的。
设有限长序列x(n)和h(n)的长度分别为N1和N2，取N>=max(N1,N2)，定义它们的N点圆周卷积为：圆周卷积定理：设有限长序列x(n)和h(n)的长度分别为N1和N2，取N>=max(N1,N2)，分别对x(n)和h(n)取N点的DFT，将结果取N点的IDFT得到y(n)，且y(n)=.圆周卷积定理建⽴起圆周卷积与DFT之间的关系，因此求圆周卷积只须⽤DFT进⾏计算即可，⽽DFT可⽤FFT实现。
线性卷积与圆周卷积实验报告
大连理工大学实验报告学院（系）：电子信息与电气工程学部专业：电子信息工程(英语强化) 班级：电英1001班姓名：刘志旋学号：201081510 组：___ 实验时间：实验室：实验台：指导教师签字：成绩：实验二电话拨号音的合成与识别五、实验数据记录和处理在函数num1_OpeningFcn中定义全局变量handles.NUM实验程序如下:1.数字键1的代码wl1=2*pi*697/8192; %数字1 的行频wh1=2*pi*1209/8192; %数字1 的列频n=[1:410]; % 每个数字用410 个采样点表示d1=sin(wl1*n)+sin(wh1*n); % 对应行频列频叠加space=zeros(1,410); %410 个0 模拟静音信号phone=[handles.NUM,d1];handles.NUM=[phone, space]; %将得到的820点存入全局变量handles.NUM中n1=strcat(get(handles.numshow,'string'),'1'); % 获取数字号码set(handles.numshow,'string',n1); % 显示号码guidata(hObject, handles); %刷新wavplay(d1,8192); %按键音其他的数字键同1，只需要修改其中的行列频和显示数字2.删除键（*键）的代码wl1=2*pi*941/8192; %*键的行频wh1=2*pi*1209/8192; %*键的列频n=[1:410];d1=sin(wl1*n)+sin(wh1*n);num=get(handles.numshow,'string'); %得到屏幕上的内容l=length(num); %得到长度n11=strrep(num,num,num(1:l-1)); %删除最后的一个字符set(handles.numshow,'string',n11); %显示现在的结果L=length(handles.NUM);handles.NUM=handles.NUM(1:L-820); %删除820点，该为要删除点的采样值guidata(hObject, handles);wavplay(d1,8192);3.全部清除键（#键）的代码wl1=2*pi*941/8192;wh1=2*pi*1477/8192;n=[1:410]; % 每个数字用 410 个采样点表示d1=sin(wl1*n)+sin(wh1*n); % 对应行频列频叠加handles.NUM=[]; %将全局变量handles.NUM清空set(handles.numshow,'string',[]); % 将显示清空guidata(hObject, handles);wavplay(d1,8192);4.识别键（receive 键）的代码L=length(handles.NUM);n=L/820;number=zeros(1,n);for i=1:nj=(i-1)*820+1;d=handles.NUM(j:(j+410-1)); % 截取出每个数字f=fft(d,8192); % 以N=2048 作FFT 变换a=abs(f);p=a.*a/8192; % 计算功率谱num(1)=find(p(1:1000)==max(p(1:1000))); % 找行频%找列频num(2)=1000+find(p(1000:1700)==max(p(1000:1700)));if (num(1) < 730)row=1; % 确定行数elseif (num(1) < 810)row=2;elseif (num(1) < 900)row=3;elserow=4;endif (num(2) < 1260)column=1; % 确定列数elseif (num(2) < 1400)column=2;elsecolumn=3;end z=[row,column]; % 确定数字if z==[4,2]phone=0;elseif z==[1,1]phone=1;elseif z==[1,2]phone=2;elseif z==[1,3]phone=3;elseif z==[2,1]phone=4;elseif z==[2,2]phone=5;elseif z==[2,3]phone=6;elseif z==[3,1]phone=7;elseif z==[3,2]phone=8;elseif z==[3,3]phone=9;end%将整型数转化成字符number=strcat(number,int2str(phone)); end% 显示号码set(handles.shownum1,'string',number);5.数字拨号音的频谱图产生代码a=handles.NUM(1:420);a=fft(a,8192);a=abs(a);figure(1)stem(a)title('第一个输入数字拨号音的频谱图')xlabel('频率Hz'),ylabel('幅度')figure(2)a=a(1:1700);stem(a)title('第一个输入数字拨号音的频谱图')xlabel('频率Hz'),ylabel('幅度')六、实验结果与分析1．GUI下的实验结果图上面的六张图从左到右，从上到下分别是：运行程序后的界面、输入1008611后的界面、按1次*键（删除键）后的界面、按1次Receive（识别键）后的界面、按1次#键（清除键）后的界面、再按1次Receive（识别键）后的界面。
20点圆周卷积,并与线性卷积结果比较。
1、计算下列信号的圆周卷积，分别求出15点，20点圆周卷积，并与线性卷积结果比较。
)()(),5(7.0)(),(*)()()()(),(7.0)(),(*)()(1021552222101151111n R n h n R n x n h n x n y n R n h n R n x n h n x n y n n =-=====-分析线性卷积和圆周卷积之间的关系。
2.利用MALAB 分别求120()0.9()n x n R n =，225()0.9()n x n R n =与10()()h n R n =卷积和，并对两个线性卷积的结果进行比较，并给出其中一个用笔算时卷积和conv 计算线性卷积对照图；。
3. 已知系统的差分方程和输入信号分别为用递推法计算系统的零状态响应和零输入相应。
4. 有四个分系统T 1、 T 2、 T 3和T 4， 四个分系统分别用下面的单位脉冲响应或者差分方程描述：编写程序计算整个系统的单位脉冲响应h (n ), 0≤n ≤99。
5.假设系统函数如下式：(1) 画出极、 零点分布图， 并判断系统是否稳定；(2) 用输入单位阶跃序列u (n )检查系统是否稳定。
6. 下面四个二阶网络的系统函数具有一样的极点分布：)2(2)()1(21)(-+=-+n x n x n y n y {} 1 ,2 ,4 ,3 ,2 ,1 )(=n x ⎪⎩⎪⎨⎧==其它05,4,3,2,1,021)(:11n n h T n⎩⎨⎧==其它05,4,3,2,1,01)(:22n n h T )2(41)1(21)(41)(:33-+-+=n x n x n x n y T )1()()2(81.0)1(9.0)(:4-++---=n v n v n y n y n y T 5147.13418.217.098.22505)(2342-++--+=z z z z z z z H试用MA TLAB 语言研究零点分布对于单位脉冲响应的影响。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
题目：已知两个有限长序列
x(n>=δ(n>+2δ(n-1>+3δ(n-2>+4δ(n-3>+5δ(n-4>
h(n>=δ(n>+2δ(n-1>+δ(n-2>+2δ(n-3>
计算以下两个序列的线性卷积和圆周卷积
<1）x(n>⑤y(n> (2>x(n>⑥y(n> (3>x(n>⑨y(n>
(4>x(n>⑩y(n>b5E2RGbCAP
●调用函数circonv
function yc=circonv(x1,x2,N>
%用直接法实现圆周卷积
%y=circonv(x1,x2,N>
%y:输出序列
%x1,x2:输入序列
%N:圆周卷积的长度
if length(x1>>N
error。
end
if length(x2>>N
error。
end
%以上语句判断两个序列的长度是否小于N
x1=[x1,zeros(1,N-length(x1>>]。
%填充序列x1(n>使其长度为N，序列h(n>的长度为N1,序列x(n>的长度为N2p1EanqFDPw
x2=[x2,zeros(1,N-length(x2>>]。
%填充序列x2(n>使其长度为N
n=[0:1:N-1]。
x2=x2(mod(-n,N>+1>。
%生成序列x2((-n>>N，镜像，可实现对x(n>以N为周期的周期延拓，加1是因为MATLAB向量下标只能从1开始。
DXDiTa9E3d
H=zeros(N,N>。
%生成N行N列的零矩阵
for n=1:1:N
H(n,:>=cirshiftd(x2,n-1,N>。
%该矩阵的k行为x2((k-1-n>>N end
yc=x1*H'。
%计算圆周卷积
●调用函数cirshiftd
function y=cirshiftd(x,m,N>
%直接实现序列x的圆周移位
%y=cirshiftd(x,m,N>
%x:输入序列，且它的长度小于N
%m:移位位数
%N：圆周卷积的长度
%y:输出的移位序列
if length(x>>N
error('x的长度必须小于N'>。
end
x=[x,zeros(1,N-length(x>>]。
n=[0:1:N-1]。
y=x(mod(n-m,N>+1>。
•函数<1）x(n>⑤y(n>
clear all。
N1=5。
N2=4。
xn=[1 2 3 4 5]。
%生成x(n>
hn=[1 2 1 2]。
%生成h(n>
yln=conv(xn,hn>。
%直接用函数conv计算线性卷积
ycn=circonv(xn,hn,5>。
%用函数circonv计算N1点圆周卷积ny1=[0:1:length(yln>-1]。
ny2=[0:1:length(ycn>-1]。
subplot(2,1,1>。
%画图
stem(ny1,yln>。
ylabel('线性卷积'>。
subplot(2,1,2>。
stem(ny2,ycn>。
ylabel('圆周卷积'>。
•函数<2）x(n>⑥y(n>
clear all。
N1=5。
N2=4。
xn=[1 2 3 4 5]。
%生成x(n>
hn=[1 2 1 2]。
%生成h(n>
yln=conv(xn,hn>。
%直接用函数conv计算线性卷积
ycn=circonv(xn,hn,6>。
%用函数circonv计算N1点圆周卷积ny1=[0:1:length(yln>-1]。
ny2=[0:1:length(ycn>-1]。
subplot(2,1,1>。
stem(ny1,yln>。
ylabel('线性卷积'>。
subplot(2,1,2>。
stem(ny2,ycn>。
ylabel('圆周卷积'>。
•函数<3）x(n>⑨y(n>
clear all。
N1=5。
N2=4。
xn=[1 2 3 4 5]。
%生成x(n>
hn=[1 2 1 2]。
%生成h(n>
yln=conv(xn,hn>。
%直接用函数conv计算线性卷积
ycn=circonv(xn,hn,9>。
%用函数circonv计算N1点圆周卷积
ny1=[0:1:length(yln>-1]。
ny2=[0:1:length(ycn>-1]。
subplot(2,1,1>。
stem(ny1,yln>。
ylabel('线性卷积'>。
subplot(2,1,2>。
stem(ny2,ycn>。
ylabel('圆周卷积'>。
•函数<4）x(n>⑩y(n>
clear all。
N1=5。
N2=4。
xn=[1 2 3 4 5]。
%生成x(n>
hn=[1 2 1 2]。
%生成h(n>
yln=conv(xn,hn>。
%直接用函数conv计算线性卷积
ycn=circonv(xn,hn,10>。
%用函数circonv计算N1点圆周卷积ny1=[0:1:length(yln>-1]。
ny2=[0:1:length(ycn>-1]。
subplot(2,1,1>。
stem(ny1,yln>。
ylabel('线性卷积'>。
subplot(2,1,2>。
stem(ny2,ycn>。
ylabel('圆周卷积'>。
思考题：
①圆周卷积与线性卷积的关系：
若有x1(n>与x2<n）两个分别为N1与N2的有限长序列，则它们的线性卷积y1<n）为N1+N2-1的有限长序列，而它们的N点圆周卷积y2<n）则有以下两种情况：1，当N<N1+N2-1时，y2<n）是由y1<n）的前N点和后<N1+N2-1-N）点圆周移位后的叠加而成；N> N1+N2-1时，y2<n）的前N1+N2-1的点刚好是y1<n）的全部非零序列，而剩下的N-(N1+N2-1>个点上的序列则是补充的零。
RTCrpUDGiT
②线性卷积运算步骤：
求x1(n>与x2<n）的线性卷积：对x1(m>或x2<m）先进行镜像移位x1<-m），对移位后的序列再进行从左至右的依次平移x(n-m>,当n=0,1,2.…N-1时，分别将x(n-m>与x2<m）相乘，并在m=0,1,2.…N-1的区间求和，便得到y<n）5PCzVD7HxA
③圆周卷积运算步骤：
圆周卷积过程中，求和变量为m，n为参变量，先将x2(m>周期化，形成x2((m>>N,再反转形成x2((-m>>N,取主值序列则得到x2((-m>>NRN(m>,通常称之为x2(m>的圆周反转。
对x2(m>圆周反转序列圆周右移n,形成x2((n-m>>NRN(m>,当n=0,1,2,…,N-1时，分别将x1(m>与x2((n-m>>NRN(m>相乘，并在m=0到N-1区间内求和，便得到圆周卷积y(n>。
jLBHrnAILg
④用圆周移位代替线性移位的好处：
时域圆周卷积在频域上相当于两序列的DFT的相乘，而计算DFT可以采用它的快速算法——快速傅立叶变换<FFT），因此圆周卷积和线性卷积相比，计算速度可以大大加快。
xHAQX74J0X
实验总结：
通过本次实验，我掌握了线性卷积与圆周卷积软件实现的方法，并验证了两者之间的关系，同时，通过上机调试程序，进一步增强了我使用计算机解决问题的能力。
在编程过程中照葫芦画瓢是没有用的，只有真正理解程序中的精华与奥妙，才能真正掌握，把书本上的知识转换为自己的知识。
LDAYtRyKfE
申明：
所有资料为本人收集整理，仅限个人学习使用，勿做商业用途。