大连理工大学数字信号处理实验实验一线性卷积与圆周卷积..docVIP

大连理工大学实验报告 学院（系）： 电信 专业： 电子信息工程 班级： 电子1204 姓 名： 梁宇 学号： 201281001 组： ___ 实验时间： 实验室： 实验台： 指导教师签字： 成绩： 实验一、线性卷积与圆周卷积 实验题目和代码 给出序列x=[3,11,7,0,-1,4,2]，h=[2,3,0,-5,2,1]；用两种方法求两者的线性卷积y，对比结果。 直接调用matlab内部函数conv来计算。 根据线性卷积的步骤计算。 实验源程序： （a） clc clear all x=[3,11,7,0,-1,4,2]; nx=[0：length(x)-1]; h=[2,3,0,-5,2,1]; nh=[0:length(h)-1]; y=conv(x,h); ny=[0:1:length(y)-1]; subplot(1,3,1) stem(nx,x); xlabel(n) ylabel(幅度) title(X) subplot(1,3,2) stem(nh,h) xlabel(n) ylabel(幅度) title(h) subplot(1,3,3) stem(ny,y); xlabel(n) ylabel(幅度) title(线性卷积) 运行结果： （b） clc clear all x=[3,11,7,0,-1,4,2]; h=[2,3,0,-5,2,1]; x1=zeros(1,length(h)-1); x2=zeros(1,length(h)-1); x3=[x1 x x2]; %补零 nx=[0:length(x)-1]; h1=zeros(1,length(x)+(length(h)-2)); h2=[h h1]; nh=[0:length(h)-1]; %补零 y=[];i=1; for j=0:1:16 M=j+2; a=0; for t=1:1:M-1 a=a+x3(t)*h2(M-t); end y(i)=a; i=i+1; end ny=[-(length(h)-1):1:length(y)-6]; subplot(1,3,1) stem(nx,x); xlabel(n) ylabel(幅度) title(X) subplot(1,3,2) stem(nh,h) xlabel(n) ylabel(幅度) title(h) subplot(1,3,3) stem(ny,y); xlabel(n) ylabel(幅度) title(线性卷积) axis([0 length(h)+length(x)-2 0 60]) %去掉补位0 运行结果： 将函数conv稍加扩展为函数conv_m，它可以对任意基底的序列求卷积。格式如下： function [y,ny]=conv_m(x,nx,h,nh) % 信号处理的改进卷积程序 % [y,ny]=conv_m(x,nx,h,nh) % [y,ny]=卷积结果 % [x,nx]=第一个信号 % [h,nh]=第二个信号 源程序： function [y,ny]=conv_m(x,nx,h,nh) x=input(x=); h=input(h=); nx=length(x); nh=length(h); ny=nx+nh-1; for(n=1:ny) y(n)=0; for(m=1:nh) k=n-m+1; if(k=1k=nx) y(n)=y(n)+h(m)*x(k); end; end; end; n=0:1:ny-1; xx=0:nx-1; hh=0:nh-1; subplot(311); stem(xx,x,.); title(x序列); grid; subplot(312); stem(hh,h,.); title(h序列); grid; subplot(313); stem(n,y,.); title(调用co