成都理工大学matlab信号与系统课程设计实验报告.docVIP

信号与系统 课程设计(论文) 设计（论文）题目 离散系统的响应及卷积运算和灰度图像频域处理 学院名称 信息科学与技术学院 专业名称 电子信息科学与技术 学生姓名 邓泉铃 学生学号 201313020220 任课教师 杨斯涵 设计（论文）成绩 教务处 制 2015年 1 月 6 日 摘要 此次的信号与系统课程设计的任务是在MATLAB软件下进行离散系统的分析仿真及卷积运算。。技术内容是：根据时域分析原理，利用MATLAB软件求解系统零状态响应，冲激响应。对MATLAB软件进行程序操作，同时利用MATLAB软件也能对书本上的知识进行验证，在MATLAB软件下编写函数程序，然后运行程序，与书本上的信号的求解进行对照分析和比较。对MATLAB软件进行一定的了解和运用之后，开始做此次课程设计——LTI连续系统的分析仿真，用MATLAB软件对此次课程设计的系统零状态响应、冲激响应进行绘图求解，并且记录其分析过程。 关键字：课程设计，函数程序，分析仿真。 一、基础设计 离散系统的响应及卷积运算 1.1、已知某LTI离散系统的单位序列响应，试用MATLAB中的函数conv( )求解该系统在激励时的零状态响应，绘出其时域波形图，并说明序列和的时域宽度与的时域宽度的关系。 调用conv函数代码： function[f,k]=dconv(f1,f2,k1,k2) %The function of computer f=f1*f2 f=conv(f1,f2) k0=k1(1)+k2(1); k3=length(f1)+length(f2)-2; k=k0:k0+k3 subplot(2,2,1) stem(k1,f1) title(f1(k)) xlabel(k) ylabel(f1(k)) subplot(2,2,2) stem(k2,f2) title(f2(k)) xlabel(k) ylabel(f2(k)) subplot(2,2,3) stem(k,f); title(f(k)=f1(k)*f2(k)) xlabel(k) ylabel(f(k)) MATLAB代码如下： f1=ones(1,4); k1=0:3; f2=ones(1,3); k2=0:2; [f,k]=dconv(f1,f2,k1,k2) 运行结果如右： f = 1 2 3 3 2 1 k = 0 1 2 3 4 5 f = 1 2 3 3 2 1 k =0 1 2 3 4 5 y(k)的时域宽度是f(k)与h(k)之和减一。 1.2、已知某LTI离散系统的差分方程如下： ， 试用MATLAB中的impz（）函数绘出该系统0～60时间范围内的单位序列响应波形。 MATLAB代码如下： a=[2,-1,3] b=[1,0,2]; subplot(1,1,1) impz(b,a,60) 输出图像如下： 1.3、已知某离散系统的差分方程为 ， 若输入信号为，试用MATLAB中的filter( )和stem( )函数绘出该系统在时间间隔内的零状态响应。 代码如下： a=[1,0.5,1] b=[1,-1,0]; k=0:20; x=exp(-2*k); dk=zeros(1,length(k)); dk(1)=1; uk=ones(1,length(k)); subplot(2,1,1) stem(k,x) title(x(k)) y=filter(b,a,x) subplot(2,1,2) stem(k,y) title(零状态响应) 二、综合设计 灰度图像频域处理 实验目的：正确理解二维傅里叶变换及滤波的基本概念，掌握低通、高通滤波器。 内容：提供一幅灰度图像； （1） 使用Matlab中的imread命令读取图像，并用imshow显示图像； （2） 将图像数据变换到二维频域（采用命令fft），判断该图像的截止频率。 （3） 在频域进行低通滤波和高通滤波，恢复空域结果（采用命令ifft），比较滤波前后的图像