信号与系统上机实验完整版
信号与系统上机实验指导书
《信号与系统》上机实验指导书山东建筑大学信息与电气工程学院电子信息教研室目录实验一、连续信号的时域描述与运算 (3)实验二、离散信号的时域描述与运算 (14)实验三、连续信号的频域分析 (19)实验四、离散信号的频域分析…………………………………………………….实验五、连续线性时不变系统分析 (22)实验一 连续信号的时域描述与运算一、 实验目的1，通过绘制典型信号的波形，了解这些信号的基本特征；2，通过绘制信号运算结果的波形，了解这些信号运算对信号所起的作用。
二、 实验原理及方法1． 基于matlab 的信号描述方法如果一个信号在连续时间范围内（除有限个间断点外）有定义，则称该信号为连续时间信号，简称为连续信号。
从严格意义上讲，matlab 数值计算方法并不能处理连续信号，但是可以利用连续信号在等时间间隔点的采样值来近似表示连续信号，即当采样间隔足够小时，这些离散采样值能够被matlab 处理，并且能够较好地近似表示连续信号。
（1）向量表示法对于连续时间信号f(t),可以定义两个行向量f 和t 来表示，其中向量t 是形如t=t1:△t:t2的matlab 命令定义的时间范围向量，t1为信号的起始时间，t2为终止时间，△t 为时间间隔；向量f 为连续时间信号f(t)在向量t 所定义的时间点上的采样值。
例如对于连续正弦信号)4sin()(t t f π=，可以用向量表达式表示为：t=-8:1:8;y=sin(pi*t/4); plot(t,y); grid on;axis([-8 8 -1.1 1.1]);绘制的信号波形如实验图1－1所示，当把时间间隔△t 取得更小（如△t ＝0。
01）时，就可以得到f(t)较好的近似波形，如实验图1－2所示。
（2）符号运算表示法如果信号可以用一个符号表达式来表示，则可用ezplot 命令绘制出信号的波形。
例如对于连续信号)4sin()(t t f π=，可以用符号表达式表示为：f=sym(‘sin(pi/4*t)’);ezplot(f,[-8 8]);该命令绘制的信号波形如实验图1－3所示。
信号与系统上机实验报告-信号与系统实验报告实验一
姓名：学号：学院：Q1:0<=n<=31x1(n)=sin(pi*n/4)*cos(pi*n/4);x2(n)=cos(pi*n/4)*cos(pi*n/4)x3(n)=sin(pi*n/4)*cos(pi*n/8)分别画出图形,求出其周期。
代码：n=0:31;x1=sin(pi*n/4).*cos(pi*n/4);x2=cos(pi*n/4).*cos(pi*n/4);x3=sin(pi*n/4).*cos(pi*n/8);stem(n,x1)stem(n,x2)stem(n,x3)结果：05101520253035由图形可知周期T1=4由图形可知周期T2=405101520253035由图形可知周期T3=16Q2:当0<=n<=5时，h(n)=n;其他h(n)=0;x(n)=h(n);求y(n)=x(n)*h(n);用stem函数画出y(n).代码：n=0:5;h=n;x=h;y=conv(x,h);stem(y)结果：Q3：（a）.定义用向量a1和b1描述差分方程y(n)-0.8y(n-1)=2x(n)-x(n-2)表征的因果LTI系统，(b).用在（a）中的系数向量，利用freqz定义H1是在0和pi之间4个等份频率上频率响应的值，omega1是这些频率值。
(c).用在（a）中的系数向量，利用freqz定义H2是在0和2*pi之间4个等份频率上频率响应的值，omega2是这些频率值。
代码：n=4;a1=[5,0,-4]b1=[10,0,-5][H1,W1] = freqz(b1,a1,n)[H2,W2] = freqz(b1,a1,n,'whole')结果：a1 =5 0 -4b1 =10 0 -5H1 =5.00001.7073 - 0.3659i1.66671.7073 + 0.3659iW1 =0.78541.57082.3562H2 =5.00001.66675.00001.6667W2 =1.57083.14164.7124Q4: X1(n)=u(n)-u(n-8); 其周期N1=8,X2(n)=u(n)-u(n-8); 其周期N2=16,X3(n)=u(n)-u(n-8); 其周期N3=32,(1) 画出这些周期信号在0<=n<=63的图形(2) 求其对应的付氏级数，(分别为a1,a2,a3)并画图。
北航信号与系统上机实验报告
信号与系统上机实验报告我是 buaa 快乐的小2B目录实验一、连续时间系统卷积的数值计算 (3)一、实验目的 (3)二、实验原理 (3)三、实验程序源代码、流图实验程序源代码 (4)4.1源代码与程序框图： (4)4.2数据与结果 (5)4.3数据图形 (6)实验二、信号的矩形脉冲抽样与恢复 (7)一、实验目的： (7)二、实验原理： (7)三、实验内容 (9)四、实验程序流程图和相关图像 (9)4.1、画出f(t)的频谱图即F(W)的图像 (9)4.2、对此频域信号进行傅里叶逆变换，得到相应的时域信号，画出此信号的时域波形f(t) (11)4.3、三种不同频率的抽样 (14)4.4、将恢复信号的频谱图与原信号的频谱图进行比较 (17)实验五、离散时间系统特性分析 (21)一、实验目的： (21)二、实验原理： (21)三、实验内容 (21)四、程序流程图和代码 (22)五、实验数据： (23)5.1单位样值响应 (23)5.2幅频特性 (24)六、幅频特性和相频特性曲线并对系统进行分析。
(25)6.1幅频特性曲线 (25)6.2相频特性曲线 (26)实验一、连续时间系统卷积的数值计算一、实验目的1 加深对卷积概念及原理的理解；2 掌握借助计算机计算任意信号卷积的方法。
二、实验原理1 卷积的定义卷积积分可以表示为2 卷积计算的几何算法卷积积分的计算从几何上可以分为四个步骤：翻转→平移→相乘→叠加。
3 卷积积分的应用卷积积分是信号与系统时域分析的基本手段，主要用于求系统零状态响应，它避开了经典分析方法中求解微分方程时需要求系统初始值的问题。
设一个线性零状态系统，已知系统的单位冲激响应为h(t)，当系统的激励信号为e(t)时，系统的零状态响应为由于计算机技术的发展，通过编程的方法来计算卷积积分已经不再是冗繁的工作，并可以获得足够的精度。
因此，信号的时域卷积分析法在系统分析中得到了广泛的应用。
卷积积分的数值运算实际上可以用信号的分段求和来实现，即：如果我们只求当t ）时r(t)的值，则由上式可以得到：1 1 2t = nΔt （n为正整数, nΔt 记为当 1 Δt 足够小时，( ) 2 r t 就是e(t)和h(t)卷积积分的数值近似，由上面的公式可以得到卷积数值计算的方法如下：1、将信号取值离散化，即以Ts 为周期，对信号取值，得到一系列宽度间隔为Ts 的矩形脉冲原信号的离散取值点，用所得离散取值点矩形脉冲来表示原来的连续时间信号；2、将进行卷积的两个信号序列之一反转，与另一信号相乘，并求积分，所得为t=0 时的卷积积分的值。
《信号与系统》实验报告（完整版）
《信号与系统》实验报告（完整版）长江大学电工电子实验中心电路与系统(2)实验报告姓名高文昌班级电信10909班序号06指导教师黄金平老师成绩实验名称：连续信号的绘制一、实验目的1．掌握用Matlab 绘制波形图的方法，学会常见波形的绘制。
2．掌握用Matlab 编写函数的方法。
3．周期信号与非周期信号的观察。
加深对周期信号的理解。
二、实验内容1、用MATLAB 画出下列信号的波形。
(a) ][cos )(1t t f ε=; (b) )]2()2([2||)(2--+=t t t t f εε; (c) )]2()([sin )(3t t t t f ---=εεπ; (d) )sgn()()(24t t G t f =; (e) )2()(265-=t Q t G f ; (f) )sin(|)|2()(6t t t f πε-= （a ）t=linspace(-10,10,400);f1=u(cos(t));figure(1),myplot(t,f1)xlabel('Time(sec)'),ylabel('f1(t)')(b)t=linspace(-4,4,400);f2=abs(t)/2.*(u(t+2)-u(t-2)); figure(2),myplot(t,f2)xlabel('Time(sec)'),ylabel('f2(t)');(c)t=linspace(-1,3,400);f3=sin(pi*t).*(u(-t)-u(2-t)); figure(3),myplot(t,f3)xlabel('Time(sec)'),ylabel('f3(t)')(d)t=linspace(-2,2,400); f4=sign(t).*rectpuls(t,2); figure(4),myplot(t,f4)xlabel('Time(sec)'),ylabel('f3(t)')(e)t=linspace(-1,4,400);f5=rectpuls(t,6).*tripuls(t-2,4); figure(5),myplot(t,f5)xlabel('Time(sec)'),ylabel('f5(t)')(f)t=linspace(-4,4,400); f6=u(2-abs(t)).*sin(pi*t) figure(6),myplot(t,f6)xlabel('Time(sec)'),ylabel('f6(t)')2、用基本信号画出图2.1-10中的信号。
信号与系统上机课程实验1和2报告
实验一1.连续时间信号的运算实验程序如下syms tf1=('(-t+4)*(u(t)-u(t-4))'); f2=('sin(2*pi*t)');y1=subs(f1,t,-t);f3=y1+f1;y2=subs(f1,t,-2*t);f4=-(y2+f1);f5=f2*f3;f6=f1*sin(2*pi*t);y3=subs(f6,t,t-2);f7=y3+f2;subplot(2,2,1);ezplot(f4,[-10,10]);subplot(2,2,2);ezplot(f5,[-6,6]);subplot(2,2,3);ezplot(f6,[-1,5]);subplot(2,2,4);ezplot(f7,[0,8]);其中阶跃函数定义子程序为function f=u(t)f=(t>0);程序运行结果如下图2. 信号的直流/交流分量实验程序如下t=0:0.1:500;f=100*abs(sin(2*pi*t/50));f1=fDC(f);f2=fAC(f,f1);plot(t,f,'r',t,f1,'b',t,f2,'g') 其中函数定义的子程序如下function fD=fDC(f)fD=mean(f);function fA=fAC(f,fD)fA=f-fD;程序运行结果如下图实验二1．f1 与f2 的卷积实验程序如下t0=-1;t1=2;dt=0.01;t=t0:dt:t1;f1=u(t+0.5)-u(t-1);f2=(0.5*t).*(u(t)-u(t-2)); f=dt*conv(f1,f2);subplot(2,2,1);plot(t,f1)title('f1');axis([t0,t1,-0.1,1.2]); subplot(2,2,2);plot(t,f2)title('f2');axis([t0,t1,-0.1,1.2]); subplot(2,1,2);t=2*t0:dt:2*t1;plot(t,f)title('f1*f2');axis([2*t0,2*t1,-0.1,1.2]);其中阶跃函数定义子程序如下function y=u(t)y=(t>0);程序运行结果如下2.求冲激响应h(t)实验程序a=[2 1 8];b=[0 0 1];impulse(b,a)程序运行结果如下3.求阶跃响应g(t) 实验程序a=[2 1 8];b=[0 0 1];step(b,a)程序运行结果如下4.求零状态响应实验程序a=[1 4 4];b=[0 1 3];p=0.01;t=0:p:8;x=exp(-t);lsim(b,a,x,t)实验程序运行结果如下。
西工大信号与系统上机实验3
上机实验3 连续LTI 系统的频域分析 一、实验目的 （1）掌握连续时间信号傅立叶变换和傅立叶逆变换的实现方法，以及傅立叶变换的时移特性，傅立叶变换的频移特性的实现方法； （2）了解傅立叶变换的特点及其应用； （3）掌握函数fourier 和函数ifourier 的调用格式及作用； （4）掌握傅立叶变换的数值计算方法，以及绘制信号频谱图的方法。
二、实验内容与方法1.验证性实验（1）傅立叶变换。
已知连续时间信号()2t f t e -=，通过程序完成()f t 的傅立叶变换。
MATLAB 程序：syms t;f=fourier(exp(-2*abs(t)));ezplot(f);运行结果如下：试画出()()323t f t e U t -=的波形及其幅频特性曲线。
MATLAB 程序：Syms t v w ff=2/3*exp(-3*t)*sym(‘Heaviside(t)’);F=fourier(f);subplot(2,1,1);ezplot(f);subplot(2,1,2);ezplot(abs(F));信号()()323t f t e U t -=的波形及其幅频特性曲线如图所示：（2）傅立叶逆变换。
已知()211f t ω=+，求信号()F j ω的逆傅立叶变换。
MATLAB 程序：syms t wifourier(1/(1+w^2),t)结果如下：()()()()11*exp **exp *22ans t U t t U t =-+ (3)傅立叶变换数值计算。
已知门函数()()()()211f t g t U t U t ==+--，试采用数值计算方法确定信号的傅立叶变换()F j ω。
MATLAB 程序：R=0.02;t=-2:R:2;f=stepfun(t,-1)-stepfun(t,1);W1=2*pi*5;N=500;k=0:N;W=k*W1/N;F=f*exp(-j*t'*W)*R;F=real(F);W=[-fliplr(W),W(2:501)];F=[fliplr(F),F(2:501)];subplot(2,1,1);plot(t,f); axis([-2,2,-0.5,2]);xlabel('t');ylabel('f(t)'); title('f(t)=U(t+1)-U(t-1)');subplot(2,1,2);plot(W,F); axis([-40,40,-1,2]);title('f(t)的傅立叶变换');ylabel('F(w)');xlabel('w');信号的傅立叶变换如图：（4）连续函数的傅立叶变换。
信号与系统实验报告最终版
实验报告2015年 6 月实验1 常见信号观测实验一、实验目的1.观察和测量各种典型信号；2.掌握有关信号的重要性，了解其在信号与系统分析中的应用。
二、实验原理说明 1．正弦函数信号； 2．指数函数信号； 3．指数衰减震荡函数信号； 4．抽样函数信号； 5.钟形函数信号； 三、实验原理波形产生原理框图如下图所示四、实验步骤1．打开实验箱，调节SW101（程序选择）按钮，使程序指示灯显示D3D2D1D0=0001，对应信号观测；（实验箱上电时默认D3D2D1D0=0001，因此不用调节）2．将跳线开关K801，K802，K803和K804连续到左侧；3. 用示波器分别测量TP801，TP802，TP803，TP804，TP805的波形，并记录下来。
测试点说明如下：（1）TP801：测试正弦函数信号波形（2）TP802：测试指数函数信号波形（3）TP803：测试指数衰减震荡函数信号波形（4）TP804：测试抽样函数信号波形（5）TP805：测试种形函数信号波形五、实验设备1.双踪示波器2.信号系统实验箱六、实验结果实验2 冲激响应与阶跃响应一、实验目的1.观察和测量RLC串联电路的阶跃响应与冲激响应的波形和有关参数，并研究其电路元件参数变化对响应状态的影响；2.掌握有关信号时域的测量方法。
二、实验原理说明实验如图1-1所示为RLC串联电路的阶跃响应与冲激响应的电路连接图，图2-1（a）为阶跃响应电路连接示意图；图2-1（b）为冲激响应电路连接示意图。
三、实验内容1.阶跃响应波形观察与参数测量设激励信号为方波，其幅度为1.5V，频率为500Hz。
实验电路连接图如图2-1（a）所示。
①连接P04与P914。
②调节信号源，使P04输出f=500Hz，占空比为50%的脉冲信号，幅度调节为1.5V；（注意：实验中，在调整信号源的输出信号的参数时，需连接上负载后调节）③示波器CH1接于TP906，调整W902，使电路分别工作于欠阻尼、临界和过阻尼三种状态，并将实验数据填入表格2-1中。
西工大信号系统上机实验一实验二
西工大信号系统上机实验一实验二上机实验1 连续时间信号的时域分析一、实验目的（1）掌握连续时间信号的时域运算的基本方法；（2）掌握相关函数的调用格式及作用；（3）掌握连续信号的基本运算；（4）掌握利用计算机进行卷积运算的原理和方法；（5）熟悉连续信号卷积运算函数conv的应用。
二、实验原理信号的基本运算包括信号的相加（减）和相乘（除）。
信号的时域变换包括信号的平移、翻转、倒相、尺度变换等，由以下公式所描述。
（1）加（减）：f(t)=f1(t)±f2(t)（2）乘：f(t)=f1(t)×f2(t)（3）延时或平移：f(t)→f(t-t0) t0>0时右移；t0<0时左移（4）翻转：f(t)→f(-t)（5）尺度变换：f(t)→f(at)（6）标量相乘：f(t)→af(t)（7）倒相：f(t)→-f(t)（8）微分：f(t)→df(t)/d(t)（9）积分：f(t)（10）卷积：f(t)=f1(t)*f2(t)三、涉及的MATLAB函数及实现1.stepfun函数功能：产生一个阶跃信号调用格式：stepfun(t,t0)，其中t是时间区间，在该区间内阶跃信号一定会产生；t0是信号发生从0到1的条约的时刻。
2.diff函数调用格式：diff(f)：求函数f对预设的独立变数的一次微分值。
diff(f，’t’)：求函数f对独立变数t的一次微分。
3.int函数调用格式：Int(f)：函数F对预设独立变数的积分值。
Int(f,’t’)：函数f对独立和变数t的积分值。
4.conv函数功能：实现信号的卷积运算。
调用格式：w=conv(u,v)：计算两个有限长度序列的卷积。
说明：该函数假定两个序列都从零开始。
四、试验内容与结果（1）相加:f(t)=f1(t)+f2(t)MATLAB程序：clear all;t=0:0.001:3;b=3;t0=1;u=stepfun(t,t0);n=length(t);for i=1:nu(i)=b*u(i)*(t(i)-t0);endy=sin(2*pi*t);f=y+u;plot(t,f);xlable('时间(t)');ylabel('幅值f(t)');title('连续信号的相加'); 实验结果截屏如下：（2）相乘：f(t)=f1(t)*f2(t)MATLAB程序如下：clear all;t=0:0.001:3;b=3;t0=1;u=stepfun(t,t0);n=length(t);for i=1:nu(i)=b*u(i)*(t(i)-t0);endy=sin(2*pi*t);f=y.*u;plot(t,f);xlabel('时间(t)');ylabel('幅值f(t)');title('连续信号的相乘'); 实验结果截屏如下：(3)倒相：由f(t)得到-f(t)MATLAB程序如下：clear all;t=-1:0.02:1;g1=3.*t.*t;g2=-3.*t.*t;grid on;plot(t,g1,'r-',t,g2,'b--');xlabel('t');ylabel('g(t)');title('倒相');实验结果截屏如下：(4)综合：f(t)=(1+t/2)*[U(t+2)-U(t-2)];绘制f(t+2);f(t-2);f(-t);f(2t);-f(t)MATLAB程序：syms tf=sym('(t/2+1)*(heaviside(t+2)-heaviside(t-2))');subplot(2,3,1);ezplot(f,[-3,3]);y1=subs(f,t,t+2);subplot(2,3,2);ezplot(y1,[-5,1])y2=subs(f,t,t-2);subplot(2,3,3);ezplot(y2,[-1,5]);y3=subs(f,t,-t);subplot(2,3,4);ezplot(y3,[-3,3]);y4=subs(f,t,2*t);subplot(2,3,5);ezplot(y4,[-2,2]);y5=-f;subplot(2,3,6);ezplot(y5,[-3,3]);运行程序前建立一个Heaviside的M文件函数。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
0.2
0
-0.2
-0.4
-0.6
-0.8
-1
-10 -8 -6 -4 -2
0
2
4
6
8 10
练习：
e （1）画信号 f (t) = −0.4t sin(2πt) 在 0≤t≤8 范围内的连续波形。
0.2
程序三：产生周期锯齿波
t=(0:0.001:2.5); y=sawtooth(2*pi*30*t); plot(t,y),
axis([0 0.2 -1 1])
注解：x = sawtooth(t)
功能：产生一个周期为 2π、幅度在-1 到+1 之间的 周期性三角波信号。若x = sawtooth(Ωt)，则对应 的周期为 2π/Ω。
2、掌握用 matlab 软件产生基本信号
信号按照自变量的取值是否连续可分为连续时间信号和离散时间信号。对信
号进行时域分析，首先需要将信号随时间变化的规律用二维曲线表示出来。对于
简单信号可以通过手工绘制其波形，但对于复杂的信号，手工绘制信号波形显得
十分困难，且难以绘制精确的曲线。
利用 Matlab 软件的信号处理工具箱(Signal Processing Toolbox)中的专用
ctrl+c
中断一个 Malab 任务
操作符
分号操作符 分号用来取消打印。如果语句的最后一个符号是分号，则打印 被取消，但是命令仍在执行，而结果不再显示。这是一个有益的特性，因为打印 中间结果可能不必要。此外，在输入矩阵时，除非最后一行，分号用来指示一行 的结束。
冒号操作符 冒号操作符在 MATLAB 中起着重要作用。该操作符用来建立向 量，赋予矩阵下标和规定叠代。例如，j:k 表示[j j+1 … k]，A(:，j)表示矩阵 A 的第 j 列，A(i，:)表示矩阵 A 的第 i 行。
列出已定义的变量 与 who 相同，但仅列出全局变量 给出所有以 a 开头的变量的一个列表 给出比命令 who 更详细的列表，如显示矩阵的维数 与 whos 相同，但仅列出全局变量 现实某个 MATLAB 文件所在的目录 查询指定目录下各类文件的性质 查询某个 MATLAB 语法、变量及工具箱函数
第 1 页 共 20 页
函数产生信号并绘出波形。
1
请画出以下实验的实验结果图，
0.8
并注生正弦波
0.4
0.2
t=0:0.01:3*pi; y=sin(2*t); plot(t,y)
注解：先在（0~3π）范围内按等步长0.01 生成向量t，然后计算向量y的值，最后调 用plot（）函数绘制二维曲线图。
第 2 页 共 20 页
4. 进一步加深对信号与系统的基本原理、方法及应用的理解。
二、实验内容和步骤
1、使用 Matlab 进行基本的矩阵运算
依此完成以下实验内容，通过实验学习掌握 Matlab 软件的基本使用方法以及 使用 Matlab 进行基本的矩阵运算的方法和步骤。 程序一：
在命令窗口中输入： >>A=[1,2,3;4,5,6]，B=A’，C=sqrt(A) 窗口中显示的 A、B、C 的值为： >> A=[1,2,3;4,5,6], A=
>>r=rank(A),e=eig(A)
输出矩阵 A 的秩和特征值。R= ,e’=
注解： eig(a) 求矩阵 a 的特征值
第 4 页 共 20 页
inv(a)或 a ^ (-1) 求矩阵 a 的逆矩阵 rank(a) 求矩阵 a 的秩 trace(a) 求矩阵 a 的迹(对角线元素之和)
>> A=magic(4) A=
0
-0.2
-0.4
-0.6
-0.8
-1
0
1
2
3
4
5
6
7
8
9 10
对于连续时间信号f(t)，可以用两个行向量f和t来表示，其中向量t是形如t=t :p:t 的
1
2
MATLAB命令定义的时间范围向量，t 为信号起始时间，t 为信号终止时间，p为时间间隔。向
1
2
量f为连续信号f(t)在向量t所定义的时间点上的样值。
Inf
定义为 1/0 （正无穷大）。当出现被零除时，Matlab 就返回 inf，
并不中断执行而继续计算
NaN
定义为“Not a Number”（表示不定值），这个非数值要么是％类
型，要么是 inf/inf
变量列表 Who Who global who a* whos whos global which what lookfor
请分别写出在命令窗口中输入以下内容，MatLab 的执行结果中 A~G 变量的值: >>A=fix(6*rand(3,3))； >>B=fix(9*rand(3,3))； >>C=A+B； >>D=A-B； >>E=A*B； >>F=A.*B； >>G=A.^(1/2)； 提示：fix（）函数为取整函数，如 fix（2.5）=2，fix（-2.5）=-2
Matlab 的程序设计基础概述
用户可在 Matlab 的命令窗口键入一个命令，也可以由它定义的语言在编辑器 中编写应用程序（如：记事本、MatLab 源程序编辑器），Matlab 软件对此进行解 释后，在 Matlab 环境下对它进行处理，最后返回结果。
Matlab 的基本语句结构 Matlab 的实质就是对矩阵进行运算处理，以复数矩阵作为最基本的运算单元，
0.2
f(t)=Sa(t)
0
如果把 t 改为 t=-10:0.11:10; 则波形变为： -0.2
f(t)=Sa(t)
-0.4
-10 -8 -6 -4 -2
0
2
4
6
8 10
t
1
0.8
0.6
0.4
0.2
0
-0.2
-0.4
-10 -8 -6 -4 -2
0
2
4
6
8 10
t
如果需要在一幅图上画出多条曲线，可采用具有多个自变量的 plot 命令。多
x = sawtooth(t, width) 功能：产生一个周期为 2π、幅度在-1 到+1 之间的 周期性三角波信号。其中width表示最大幅度出现的
0 -0.2 -0.4 -0.6 -0.8
-1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
1 0.8 0.6 0.4 0.2
rand（）函数为随机数生成函数。 分析“*”与“.*”的区别；“；”及“who”等变量的使用方法。 .^为矩阵中对所有的单个元素进行乘方运算符，程序中为开平方，如果为 A.^2，则表示为 2 次方。 >> A
第 3 页 共 20 页
A= 522 150 344
>> B B=
331 871 001 >> C C= 853 9 12 1 345 >> D D= 2 -1 1 -7 -2 -1 343 >> E E= 31 29 9 43 38 6 41 37 11 >> F F= 15 6 2 8 35 0 004 >> G G= 2.2361 1.4142 1.0000 2.2361 1.7321 2.0000
序设计流程、高质量的图形可视化界面设计、丰富的交互式仿真集成环境，以及
与其他程序和语言便捷的接口，是通用的科学计算、数值仿真及数据可视化的重
要工具。其广泛应用于数学计算；新算法研究与开发；建模、仿真及样机开发；
数据分析、探索及可视化、科技与工程的图形功能；友好图形界面的应用程序开
发，其主要特点包括：
《信号与系统》实验上机指导
MATLAB 的基本介绍
MATLAB (MATrix LABoratory)由美国 The MathsWorks 公司开发，是当今国际
上在科学界和教育界中最具有影响力、也最具有活力的软件，它起源于矩阵运算，
现已发展成为一种高度集成的计算机语言；它提供了强大的科学运算、灵活的程
0 -0.2 -0.4 -0.6 -0.8
-1 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2
位置：即在一个周期内，信号从t=0 到width×2π时
函数值从-1 到+1 线性增加，而从width×2π到 2π又是从+1 到-1 线性下降。width取值在 0 ~
用上述向量对连续信号进行表示后，就可以用plot命令来绘制出信号的时域波形。plot 命令可将点与点间用直线连接，当点与点间的距离很小时，绘出的波形就成了光滑的曲线。
第 5 页 共 20 页
程序二: 产生周期方波
1
0.8
t=0:0.01:1; 0.6
y=square(4*pi*t); 0.4
plot(t,y)
16 2 3 13 5 11 10 8 9 7 6 12 4 14 15 1 >> [sum(A),sum(A'),trace(A),trace(rot90(A))] ans = 34 34 34 34 34 34 34 34 34 34 >> r=rank(A),e=eig(A) r= 3 e= 34.0000 8.9443 -8.9443 -0.0000
0.1
0
0
1
2
3
4
5
6
7
8
9 10
t
第 6 页 共 20 页
程序五：产生连续信号 f(t)=Sa(t)=sin(t)/t
1