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

信号与系统与技术学院信息科学与技术201313020220 任课教师 杨斯涵 设计（论文）成绩 教务处 制 2015年 1 月 6 日 摘要 此次的信号与系统课程设计的任务是在MATLAB软件下进行离散系统的分析仿真及卷积运算。 关键字：课程设计函数程序，分析仿真。、 离散系统的响应及卷积运算 1.1、已知某LTI离散系统的单位序列响应，试用MATLAB中的函数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)的时域宽度是k)与h(k)之和减一。1.2、已知某LTI离散系统的差分方程如下： ， 试用MATLAB中的impz（）函数绘出该系统0～60时间范围内的单位序列响应波形。 如下：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） 2.1、使用Matlab中的imread命令读取图像，并用imshow显示图像。 代码：imshow(f); title(原图) 2.2、将图像数据变换到二维频域（采用命令fft），判断该图像的截止频率。 在频域进行低通滤波和高通滤波，恢复空域结果（采用命令ifft），比较滤波前后的图像差异。 f=imread(F:\邓泉铃\灰度图.jpg); imshow(f); C=fft2(double(f)); D=fftshift(fft2(double(f))); E=ifft2(D); figure,imshow(log(abs(D)+1),[]); title(变换频谱图); 从图像中看出x与的交点就是灰度图的。 ）低通滤波处理F=fft2(M); F=fftshift(F); figure;imshow(log(abs(F)),[]); title(快速傅里叶变换（FFT）图像) [M,N]=size(F); nn=2; d0=100; m=fix(M/2);n=fix(N/2); for i=1:M for j=1:N d=sqrt((i-m)^2+(j-n)^2