终稿信号与系统MATLAB实验报告YAO.doc

信号与系统 MATLAB 仿真实验报告 姓名：姚俊涛 学号：2009302540313 班级：0910 学院：电气工程学院 专业：电气工程与自动化 目录 实验1： 连续时间信号的表示及可视化 4 1.1、实验内容 4 1.2、实验程序与实验图形 4 (1)冲激信号 4 (2)阶跃信号 5 (3)指数信号 5 (4)门函数 6 (5)f(t)=Sa(nwt) 7 (6)f(t)=Sin(2*pi*f*t) 8 1.3．实验体会： 9 实验2： 离散时间信号的表示及可视化 10 2.1、实验内容 10 2.2、实验程序与实验图形 10 (1)离散抽样序列 10 (2) 阶跃序列 11 (3)指数序列 11 (4) f(n)=R(n) (分别取不同的N值) 12 (5) f(n)=Sa(nw) 14 (6) f(n)=Sin(nw) (分别取不同的w值) 15 2.3．实验心得： 17 实验3： 系统的时域求解 17 3.1、实验内容 17 3.2、实验程序与实验图形 17 3.2.1 画出x(n)、h(n)、y(n)的波形 17 3.2.2 绘出y(n)的单位抽样响应，并绘出H（）的幅频及相频特性曲线 18 3.3．实验心得 20 实验4： 信号的DFT分析 20 4.1、实验内容 20 4.2、实验程序与实验图形 20 4.2.1 N=10时 21 4.2.2 N=16时 22 4.2.3 N=22时 23 4.3 扩展实验 N=32()和N=64()时 24 4.3.1 N=32时 24 4.3.2 N=64时 25 4.4 结果对比与差别原因分析 26 4.4.1进行理论值与计算值的比较 26 4.4.2查阅相关资料讨论信号频谱分析过程中误差原因及改善方法 26 4.3．实验心得 27 实验5： 系统时域解的快速卷积求法 27 5.1、实验内容 27 5.2、问题分析 27 5.3、实验程序与实验图形 28 5.3.1 序列、的图形 28 5.3.2 N=25时 29 5.3.3 N=34时 30 5.3.4 N=45时 31 5.4正确卷积结果 32 5.5 结果对比与差别原因分析 33 5.4、实验心得体会 34 六、信号与系统MATLAB仿真实验心得体会 34 实验1： 连续时间信号的表示及可视化 1.1、实验内容 ；；（分别取）； ；；（分别画出不同周期个数的波形）； 1.2、实验程序与实验图形 (1)冲激信号 syms t f=dirac(t) ezplot(f,[-5:5]); xlabel((t)); ylabel((f)); 冲激函数示意图如下 (2)阶跃信号 创建Heaviside 的M 文件，该文件如下： function f=Heaviside(t) 正确定义出该函数并保存运行后，就可调用该函数了。如先定义向量： t=-5:0.01:5 然后调用 Heaviside 函数表示出该信号并绘出波形 f=heaviside(t) plot(t,f) axis([-5,5,-0.2,1.2]) (3)指数信号 将a=1和a=2时的图形分别画于同一界面中，便于形成对比， syms t subplot 121 f=exp(t) % a=1 ezplot(f,[-5:5]); xlabel((t)); ylabel((f)); title(f(t)=exp(t)) subplot 122 f=exp(-t) %a=-1 ezplot(f,[-5:5]); xlabel((t)); ylabel((f)); title(f(t)=exp(t)) (4)门函数 t=-5:0.01:5; %设定时间变量t的范围及步长 y=rectpuls(t,2); %用rectpuls（t a)命令表示门函数，默认以零点为中心，宽度为a plot(t,y); %用plot函数绘制连续函数 title(门函数); %用title函数设置图形的名称 axis([-5 5 -0.5 1.5]); (5)f(t)=Sa(nwt) syms t f=Sinc(pi/8*t) subplot 121 ezplot(f,[-15:15]); xlabel((t)); ylabel((f)); axis([-15 15 -0.3 1.2]) title(f(t)=Sa(pi/8*t)); f=sinc(pi/18*t) subplot 122 ezplot(f,[-15:15]); xlabel((t)); ylabel((f)); title(f(t)=Sa(pi/18*t)); (6)f(t)=Sin(2*pi*f*t) 周期信号(分别画出不同周期个数的波形) sym t a %频率记为a f=sin(2*pi*