数字信号 实验四.docVIP

数字信号处理 第 第 PAGE 1 页 共 NUMPAGES 1 页 实 验 四 实验项目名称： 线性卷积与循环卷积 实验日期： 2015.05.11 实验学时： 2 实验地点： S1-326 姓名： 秦显富 学号： 201201100118 实验成绩： 指导教师签字： 实验目的和要求： （１）进一步掌握线性卷积的计算机编程方法，利用卷积的方法观察系统响应的时域特性。 （２）掌握循环卷积的计算机编程方法，并比较与线性卷积的差别，验证二者之间的关系。利用循环卷积的方法观察、分析系统响应的时域特性。 实验内容与原理： 实验原理： （１）线性卷积： 线性时不变系统（Linear Time-Invariant System, or L. T. I系统）输入、输出间的关系为：当系统输入序列为，系统的单位脉冲响应为，输出序列为，则系统输出为：；上式称为线性卷积。 （２）循环卷积 设两个有限长序列和，长度分别为和，　　。如果，则 上式称为和的循环卷积。 （３）两个有限长序列的线性卷积 序列和，长度分别为点和点，为这两个序列的线性卷积，则为且线性卷积的非零值长度为点。 （４）循环卷积与线性卷积的关系 序列为点长，序列为点长，若序列和进行点的循环卷积，其结果是否等于该两序列的线性卷积，完全取决于循环卷积的长度。 由教材相关推导，得，也就是说，循环卷积是线性卷积的周期延拓序列再取主值区间。 当时循环卷积等于线性卷积，即； 当时，两者不等。 2.实验内容： 已知长为4的两个有限长序列， （1）利用MATLAB的conv（）函数求线性卷积，并绘出图形。 （２）利用MATLAB构建的循环卷积函数计算下述4种情况下循环卷积，并绘出图形。 ⑤　⑥　⑦　⑧ 调用fft()函数利用循环卷积定理计算下述4种情况下循环卷积，并绘出图形。 ⑤ ⑥ ⑦ ⑧ 比较线性卷积和循环卷积的结果，总结结论。 实验软硬件环境： 计算机一台、MATLAB2014a仿真软件 实验编程及调试： . x=[1,2,3,4];nx=0:3; h=[4,3,2,1];nh=0:3; y=conv(x,h);ny=0:6; figure;stem(ny,y); title(线性卷积); ylabel(Y); xlabel(n) . 先构建circonv.m文件 function y=circonv(x1,x2) xn2=[x2(1),fliplr(x2)]; xn2(length(xn2))=[]; C=xn2; R=x2; M=toeplitz(C,R); y=x1*(M); ?clear; x=[1,2,3,4];nx=0:3; h=[4,3,2,1];nh=0:3; x5=[x,0];nx5=0:4; h5=[h,0];nh5=0:4; yn=circonv(x5,h5);nyn=0:4; figure;stem(nyn,yn); title(五点循环卷积); ylabel(yn); xlabel(n); ?clear; x=[1,2,3,4];nx=0:3; h=[4,3,2,1];nh=0:3; x6=[x,0,0];nx6=0:5; h6=[h,0,0];nh6=0:5; yc6=circonv(x6,h6);nyc6=0:5; figure;stem(nyc6,yc6); title(六点循环卷积); xlabel(n); ylabel(yc6); ?x=[1,2,3,4];nx=0:3; h=[4,3,2,1];nh=0:3; x7=[x,0,0,0];nx7=0:6; h7=[h,0,0,0];nh7=0:6; yc7=circonv(x7,h7);nyc7=0:6; figure;stem(nyc7,yc7); title(七点循环卷积); xlabel(n); ylabel(yc7); ④x=[1,2,3,4];nx=0:3; h=[4,3,2,1];nh=0:3; x8=[x,0,0,0,0];nx8=0:7; h8=[h,0,0,0,0];nh8=0:7; yc8=circonv(x8,h8);nyc8=0:7; figure;stem(nyc8,yc8); title(八点循环卷积); xlabel(n); ylabel(yc8); . ?x=[1,2,3,4,0];nx=0:4; h=[4,3,2,1,0];nh=0:4; X=fft(x,5); H=fft(h,5); Y=X.*H; yc5=ifft(Y,5);nyc5=0:4; figure;stem(nyc5,yc5); title(快速五点循环卷积); xlabel(yc5); ylabel(n); ?x=[1,2,3,4,0,0];nx=