实验三 卷积.docVIP

实验三 卷积 实验目的： 掌握线性卷积的计算机编程方法，利用卷积的方法观察、分析系统响应的时域特性； 验证卷积定理； 掌握圆周卷积的计算机编程方法，并比较与线性卷积的差别。利用圆周卷积的方法观察、分析系统响应的时域特性。 实验结果： 圆周卷积函数circonv.m function y circonv x1,x2,N if length x1 N error N必须 x1的长度 end if length x2 N error N必须 x2的长度 end x1 [x1 zeros 1,N-length x1 ]; x2 [x1 zeros 1,N-length x2 ]; m [0:1:N-1]; x2 x2 mod -m,N +1 ; H zeros N,N ; for n 1:1:N H n,: cirshift x2,n-1,N ; end y x1*H; 圆周位移函数circshift.m function y cirshift x,m,N if length x N error N必须 x的长度 end x [x zeros 1,N-length x ]; n [0:1:N-1]; n mod n-m,N ; y x n+1 ; M文件： g1 [1 2 3 4 5]; g2 [2 2 0 1 1]; yl conv g1,g2 ; y2 circonv g1,g2,5 ; n1 [0:1:length yl -1]; n2 [0:1:length y2 -1]; subplot 2,1,1 ; stem n1,yl ; ylabel 线性卷积 ; subplot 2,1,2 ; stem n2,y2 ; ylabel 圆周卷积 ; sy3_1a.m M文件： g1 [1 3 5 7 9 11 13 15 17]; g2 [1 -2 3 -2 1]; yl conv g1,g2 ; y2 circonv g1,g2,5 ; n1 [0:1:length yl -1]; n2 [0:1:length y2 -1]; subplot 2,1,1 ; stem n1,yl ; ylabel ‘线性卷积’ ; subplot 2,1,2 ; stem n2,y2 ; ylabel ‘圆周卷积’ ; sy3_1b.m Error using circonv N必须 x1的长度 Error in a at 4 y2 circonv g1,g2,5 ;%用函数circonv计算N1点圆周卷积 分析：N必须大于等于x1的长度。在本程序中，N 5,x1 9，不满足函数的要求，故出现了错误。将N 9带入，便可以消除错误。 更改后的M文件： g1 [1 3 5 7 9 11 13 15 17]; g2 [1 -2 3 -2 1]; yl conv g1,g2 ; y2 circonv g1,g2,9 ; n1 [0:1:length yl -1]; n2 [0:1:length y2 -1]; subplot 2,1,1 ; stem n1,yl ; ylabel 线性卷积 ; subplot 2,1,2 ; stem n2,y2 ; ylabel 圆周卷积 ; sy3_1c.m M文件： N 1:9; g1 [1 2 3 4 5]; g2 [2 2 0 1 1]; g1e [g1 zeros 1,length g2 -1 ]; g2e [g2 zeros 1,length g1 -1 ]; y circonv g1e,g2e,9 sy3_2a.m y 2 6 10 15 21 15 7 9 5 M文件： g1 [1 3 5 7 9 11 13 15 17]; g2 [1 -2 3 -2 1]; g1e [g1 zeros 1,length g2 -1 ]; g2e [g2 zeros 1,length g1 -1 ]; y circonv g1e,g2e,14 sy3_2b.m y 1 1 2 4 5 7 9 11 13 -4 34 -19 17 0 重叠保留法函数overlapsav.m function [y] overlapsav x,h,N Lx length x ;M length h ; M1 M-1;L N-M1;h [h,zeros 1,N-M ]; x [zeros 1,M1 ,x,zeros 1,N-1 ]; K floor Lx+M1-1 / L +1; Y zeros K+1,N ; for k 0:K-1 xk x K*L+1:k*L+N ; Y k+1,: circonv xk,h ;