实验报告-数字信号处理实验1.docVIP

大连理工大学 本科实验报告 实验名称 ： 学院（系）： 信息与通信学院 专 业： 电子信息工程 班 级： 电子1301 学 号： 201383022 学生姓名： 同 组 人： 年 月 日 实验项目列表 序号 实验项目名称 学时 成 绩 指导教师 预习 操作 结果 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 总计 学分：  大连理工大学实验报告 学院（系）： 信息与通信工程 专业： 电子信息工程 班级： 电子1103 姓 名： 于梦月 学号： 201158021 组： ___ 实验时间： 实验室： 实验台： 指导教师签字： 成绩： 实验一 线性卷积与圆周卷积 实验目的和要求 见预习报告 实验原理和内容 见预习报告 主要仪器设备 Matlab软件 四、实验步骤与操作方法 给出序列x=[3,11,7,0,-1,4,2]，h=[2,3,0,-5,2,1]；用两种方法求两者的线性卷积y，对比结果。 直接调用matlab内部函数conv来计算。 解： clear clc nx=0:1:6; nh=0:1:5; ny=0:11; x=[3,11,7,0,-1,4,2]; h=[2,3,0,-5,2,1]; y=conv(x,h); 结果：y = 6 31 47 6 -51 -5 41 18 -22 -3 8 2 根据线性卷积的步骤计算。 将函数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 ) 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) 解： function yc=circonv(x1,x2,N); if length(x1)N error(N should higher than or equal to the length of x1!); end if length(x2)N error(N should higher than or equal to the length of x2!); end x1=[x1,zeros(1,N-length(x1))]; x2=[x2,zeros(1,N-length(x2))]; n=[0:N-1]; x2=x2(mod(-n,N)+1); H=zeros(N,N); for n=1:1:N H(n,:)=cirshifted(x2,n-1,N); end yc=x1*H; 其中cirshifted函数如下： function y = cirshifted( x,m,N ) if length(x)N error(x的长度必须小于N); end x=[x,zeros(1,N-length(x))]; n=[0:1:N-1]; y=x(mod(n-m,N)+1); end 求序列和，（）的线性卷积 解： Clear clc nx=-5:50; nh=-5:50; x=uct(nx)-uct(nx-10); h=0.8