序列的运算——卷积.ppt

序列的运算——卷积 一、实验目的 1．掌握信号的线性卷积运算。 2．掌握信号的循环卷积运算。 3. 掌握信号循环卷积计算线性卷积的条件 二、实验仪器 1．PC机； 2．MATLAB软件。 三、实验原理 学习使用MATLAB软件进行建模；学习信号卷积运算的MATLAB实现；实现信号的线性卷积运算、应用DFT实现线性卷积运算、验证循环卷积计算线性卷积的条件。 如果x(n)=x1(n)*x2(n),利用循环卷积计算，可以用circonvt函数实现如下： 四、实验内容 设卷积下面信号 选定循环卷积的长度为N=21。确定这两个序列的循环卷积与线性卷积中哪些数值是相同的。编写程序代码，并输出图形，怎么样使两者数据相同。 五、思考题 （1）如果x(n)=x1(n)*x2(n),利用线性卷积计算，用conv函数实现如下：y=conv(x1,x2)。 （2）设卷积下面信号 选定循环卷积的长度为N=21。确定这两个序列的循环卷积与线性卷积中哪些数值是相同的。编写程序代码，并输出图形，分析错误数据的原因。怎么样使两者数据相同。 * * unction?fn=circonvt(x1,x2,N)? %circonvt 函数实现输入序列 x1 和 x2 的循环卷积 ,fn 为输出序列 ? %N? 为循环卷积长度 ? % 实现方法 if?(length(x1)N|length(x2)N)% 判断输入信号的长度 ? error(N 的长度必须大于输入数据的长度 );? end? x1=[x1,zeros(1,N-length(x1))];? x2=[x2,zeros(1,N-length(x2))];? m=0:N-1;? x=zeros(N,N);? for?n=0:N-1? ???? x(:,n+1)=x2(mod((n-m),N)+1);? end;? fn=x1*x;?% 循环计算卷积 ? function y=myconv(x1,x2)　　　x1=input(x1=);　　　x2=input(x2=);　　　N=input(N=);　　　x1=[x1,zeros(1,N-length(x1))];　　　x2=[x2,zeros(1,N-length(x2))];　　　 V=circlel(x2)　　　Z=x1*V;　　　stem(Z);xlabel(n);ylabel(Z);grid on;　　　title(循环卷积结果Z) y1=circonvt(x1,h,21); y2=circonvt(x1,h,N); x1=[x1 zeros(1,N-length(x1))]; h=[h zeros(1,N-length(h))]; X1=fft(x1,N); H=fft(h,N); X=X1.*H; x=ifft(X); x=real(x); Subplot(2,2,1);stem(n,x1);title(x1(n));axis([0,33,0,1]); Subplot(2,2,2);stem(n,h);title(h(n));axis([0,33,0,1]); Subplot(2,2,3);stem(ny,y1,‘fill’);title(’21点循环卷积’);axis([0,33,0,8]); hold on; Subplot(2,2,4);stem(n,x);title(‘线性卷积); axis([0,33,0,8]); Subplot(2,2,3);stem(n,y2,r,--); axis([0,33,0,8]); hold off