信号与线性系统实验指导书

《信号与系统》课内实验指导

--MATLAB仿真

李小平 编写

南京工程学院通信工程学院

2007年9月 南京

05级1/2班（本科）使用对象：康尼学院通信工程专业

实验一 信号的表示、可视化及时域运算

一、实验目的：熟悉MATLAB软件的使用，学会信号的表示及用MATLAB来产生信号，

实现信号的可视化。

二、实验时数：2学时 三、实验内容：

信号按照自变量的取值是否连续可分为连续时间信号和离散时间信号。信号与系统的仿真运算对信号与系统课程的学习有着重要意义。MATLAB以其科学与技术计算通用的交互系统和编程语言，具有良好的交互数学计算和易于使用的集成图形，编程简单，功能连贯，成为信号与系统的仿真运算首选。

本次实验完成以下内容：

1． 熟悉MATLAB语言。A）双击快捷图标，启动MATLAB。B）了解界面功能。

主界面包括标题栏、菜单栏、工具栏和输入输出区等。C）简单的信号处理可以通过在Command Window中输入相关指令按Enter键执行来达到处理目的。复杂的信号处理可以通过编写M文件，执行M文件达到处理目的。 M文件的编写类似于Word文档编写。选中命令窗口File菜单中M-File打开M文件编辑器或单击命令窗口工具栏中New M-File图标编辑M-文件。

D）M-文件可以在命令窗口输入已保存的M-文件名回车执行，也可使用M-文件编辑器窗口中选择Debug菜单中的RUN运行M-文件。

2．连续信号的表示与可视化。MATLAB提供了大量的生成基本信号的函数，一种常用的方法是调用相关函数及运算通过向量形式来表示信号，用适当的MATLAB语句表示信号后，可以利用MATLAB的绘图命令绘制出直观的信号波形。 本次实验完成指数信号Aeat，正弦信号Asin(?0t??)，Acos(?0t??)，抽样

函数Sa(t)?sin(?t)，矩形脉冲函数（产生以t=o为中心的宽度为width的矩形?t波）rectpuls(t,width)，三角波脉冲函数tripuls(t,width,skew)。

对于连续时间信号f(t)，可以用两个行向量f和t来表示，其中向量t是形如t=t1:p:t2的MATLAB命令定义的时间范围向量，t1为信号起始时间，t2为信号终止时间，p为时间间隔。向量f为连续信号f(t)在向量t所定义的时间点上的样值。用向量对连续信号进行表示后，就可以用plot命令来绘制出信号的时域波形。plot命令可将点与点间用直线连接，当点与点间的距离很小时，绘出的波形就成了光滑的曲线。 MATLAB命令如下：

plot(t,f);

title(‘f(t)=Sa(t)’); xlabel(‘t’);

axis([-10,10,-0.4,1.1]);

以下是有关波形产生程序： %program1_1 exponential signal t=0:0.001:10;A=1;a=-0.6; ft=A*exp(a*t); plot(t,ft)

%program1_2 sinusoidal signal t=0:0.001:8;A=1;w0=2*pi;phi=pi/6; ft=A*sin(w0*t+phi); plot(t,ft)

%program1_3 sample function t=-3:.001:3; ft=sinc(t); plot(t,ft)

%program1_4 rectpuls function t=-10:0.001:10;A=1; ft=rectpuls(t,2); ft1=rectpuls(t-2-1,4); plot(t,ft,t,ft1)

%program1_5 tripuls function t=-3:.001:3; ft= tripuls (t,4,0.5); plot(t,ft)

3．离散信号的表示与可视化。离散信号的表示与可视化和连续信号类似，画图略有不同。这里仅处理单位脉冲序列和单位阶跃序列。 %program1_6 delta function n=-30:30;

delta=[zero(1,30),1, zero(1,30)];

stem(n,delta)

%program1_6 un function n=-30:30;

un=[zero(1,30), one(1,31)]; stem(n,un)

4. 信号的运算。以下通过实例编程实现信号的平移、翻转、尺度变换及微分、积分，离散序列的差分与求和。

例 通过MATLAB编程对program1_5中的三角波f(t)实现如下变换：f(t?2)，

f(2t)，f(?t)，f(2?2t)幵画出变换后的曲线。

%program1_6 t=-3:.001:3; ft= tripuls (t,4,0.5); ft1= tripuls (t-2,4,0.5); ft2= tripuls (2t,4,0.5); ft3= tripuls (-t,4,0.5); ft4= tripuls (2-2t,4,0.5); subplot(2,3,1) plot(t,ft) title(‘f(t)') subplot(2,3,2) plot(t,ft1)