大连理工大学数字信号处理实验一实验报告.docVIP

大连理工大学实验报告 学院（系）： 信息与通信工程学院专业：电子信息工程 班级： 姓 名： 学号： 组： 实验时间： 实验室： 实验台： 指导教师签字： 成绩： 实验一 线性卷积和圆周卷积 实验目的和要求 （与预习报告相同，省略） 实验原理和内容 （与预习报告相同，省略） 三、实验程序和结果 给出序列x=[3,11,7,0,-1,1,2],h=[2,3,0,-5,2,1]，用两种方法求两者的线性卷积y，对比结果。 1.1、直接调用matlab内部函数conv函数，程序为： x=[3,11,7,0,-1,1,2]; h=[2,3,0,-5,2,1]; N=length(x); %求x序列长度 M=length(h); %求h序列长度 L=N+M-1; %y序列的长度 y=conv(x,h); %调用conv函数 nx=0:N-1;nh=0:M-1;ny=0:L-1; %各序列向量显示范围 subplot(311);stem(nx,x,*r); %画出序列波形 xlabel(n);ylabel(x(n));grid on; subplot(312);stem(nh,h,*r); xlabel(n);ylabel(h(n));grid on; subplot(313);stem(ny,y,*r); xlabel(n);ylabel(y(n));grid on; y=conv(x,h); 程序运行结果为： y = 6 31 47 6 -51 -11 32 18 -7 -9 5 2 1.2、根据线性卷积步骤计算，程序为： x=[3,11,7,0,-1,1,2]; h=[2,3,0,-5,2,1]; N=length(x); M=length(h); L=N+M-1; ny=0:L-1; y=zeros(1,L); %先将y序列置0 for n=1:1:L %外重循环，依次得到y(n) for m=1:1:N %內部循环，叠加求和 k=n-m+1; %线性移位得到对应h序号 if ((k=1)(k=M)) %保证k的范围 y(n)=y(n)+x(m)*h(k); %叠加求和求得y(n) end end end stem(ny,y,*r); xlabel(n);ylabel(y(n)); grid on; 程序结果为： y = 6 31 47 6 -51 -11 32 18 -7 -9 5 2 2、将函数conv扩展成为函数conv_m,可以对任意基底的序列求卷积。 function [y,ny]=conv_m(x,nx,h,nh) ny1=nx(1)+nh(1); %输出序列的起始位置 ny2=nx(length(x))+nh(length(h)); %输出序列的截止位置 ny=[ny1:ny2]; y=conv(x,h) end 创建函数circonv，实现序列的圆周卷积，格式： function y=circonv(x1,x2,N) 3.1、采用圆周卷积步骤法 function y=circonv(x1,x2,N) if length(x1)N %x1的长度必须小于N error(N must be bigger than the length); end if length(x2)N %x2的长度必须小于N error(N must be bigger than the length); end x1=[x1,zeros(1,N-length(x1))]; %将x1补零补到长度N x2=[x2,zeros(1,N-length(x2))]; %将x2补零补到长度N y=zeros(1,N); %先将y序列置0 for n=1:1:N %外重循环，依次得到y(n) for m=1:1:N %内重循环，叠加求和 k=mod(n-m,N)+1; %圆周移位得到对应h序号 y(n)=y(n)+x1(m)*x2(k); %叠加求和得到y(n) end end 3.2采用FFT变换法 function y=circonv2(x1,x2,N) i