光信1501.docVIP

光信1501 陈志强 实验一 MATLAB方法实现信号波形和微分方程的建立与求解 一、实验目的 1、熟悉和掌握常用的用于信号与系统时域仿真分析的MATLAB函数； 2、掌握连续时间和离散时间信号的MATLAB产生，掌握用周期延拓的方法将一个非周期信号进行周期信号延拓形成一个周期信号的MATLAB编程； 3、牢固掌握系统的单位冲激响应的概念，掌握LTI系统的卷积表达式及其物理意义，掌握卷积的计算方法、卷积的基本性质； 4、掌握利用MATLAB计算卷积的编程方法，并利用所编写的MATLAB程序验证卷积的常用基本性质； 掌握MATLAB描述LTI系统的常用方法及有关函数，并学会利用MATLAB求解LTI系统响应，绘制相应曲线。 基本要求：掌握用MATLAB描述连续时间信号和离散时间信号的方法，能够编写MATLAB程序，实现各种信号的时域变换和运算，并且以图形的方式再现各种信号的波形。掌握线性时不变连续系统的时域数学模型用MATLAB描述的方法，掌握卷积运算、线性常系数微分方程的求解编程。 三、实验内容及步骤 实验前，必须首先阅读本实验原理，读懂所给出的全部范例程序。实验开始时，先在计算机上运行这些范例程序，观察所得到的信号的波形图。并结合范例程序应该完成的工作，进一步分析程序中各个语句的作用，从而真正理解这些程序。 实验前，一定要针对下面的实验项目做好相应的实验准备工作，包括事先编写好相应的实验程序等事项。 Q1-1：修改程序Program1_1，将dt改为0.2，再执行该程序，保存图形，看看所得图形的效果如何？ dt = 0.01时的信号波形 dt = 0.2时的信号波形 这两幅图形有什么区别，哪一幅图形看起来与实际信号波形更像？ 答：图如上，两者变化幅度不一样，第一幅与实际波形更像。 Q1-2：以Q1_2为文件名存盘，产生实门信号和信号。 要求在图形中加上网格线，并使用函数axis()控制图形的时间范围在-2~2秒之间。然后执行该程序，保存所的图形。 代码: t=-2:0.01:2; width=1; g2t=rectpuls(t,width); ft=cos(10*pi*t).*g2t; plot(t,g2t); hold on; plot(t,ft); axis([-2 2 -2 2]); grid on; Q1-3：将实验原理中所给的单位冲激信号和单位阶跃信号的函数文件在MATLAB文件编辑器中编写好，并分别以以文件名delta和Heaviside存入work文件夹中以便于使用。 抄写函数文件delta如下： 抄写函数文件u如下： function y=delta(t) dt=0.01; y=(Heaviside(t)-Heaviside(t-dt))/dt; Q1-4：仿照前面的示例程序的编写方法，编写一个MATLAB程序，以Q1_4为文件名存盘，使之能够在同一个图形窗口中的两个子图中分别绘制信号x[n]=0.5|n| 和x(t)=cos(2πt)[u(t)-u(t-3)]。要求选择的时间窗能够表现出信号的主要部分（或特征）。 编写的程序Q1_4如下： function y = Q1-4 t=0.001; n = -2:t:2; x = 0.5.^(abs(n)); subplot(1,2,1); plot(n,x,-) title(x(n)) xlabel(n (s)) ylabel(y(n)); grid on； axis([-2 2 -2 2]) t = 0.001; t = -2:t:5; x = cos(2*pi*t).*((t=0)-(t=3)); subplot(1,2,2); plot(t,x,-) title(x(t)) xlabel(t (s)) ylabel(y(n)); grid on axis([-2 5 -2 2]) end 信号x[n]=0.5|n| 的波形图和信号x(t)=cos(2πt)[u(t)-u(t-3)]的波形图 Q1-5：根据示例程序的编程方法，编写一个MATLAB程序，以Q1_5为文件名存盘，由给定信号 x(t) = e-0.5tu(t) 求信号y(t) = x(1.5t+3)，并绘制出x(t) 和y(t)的图形。 编写的程序Q1_5如下： function y = Q1_5 dt