信号与系统上机实验完整版
信号与系统上机实验指导书
《信号与系统》上机实验指导书山东建筑大学信息与电气工程学院电子信息教研室目录实验一、连续信号的时域描述与运算 (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文件函数。
电子科技大学信号与系统上机实验报告
信号与系统上机实验报告班级：2010023010学号:2010020030016姓名：李思豪时间：2011年12月27日第一题一、实验目的1．熟悉MATLAB5.3的软件操作环境和编程方法。
2.学习画出信号波形并从中确定信号性质。
3.计算卷积，频率响应以及利用函数：conv, freqz, freqs and filter所得到的输出信号。
二、实验内容（一）实验程序n=[0:31];x1=sin(pi*n/4).*cos(pi*n/4);x2=(cos(pi*n/4)).^2;x3=sin(pi*n/4).*cos(pi*n/8);subplot(3,1,1);stem(n,x1)title('figures of signals')xlabel('n')ylabel('x1[n]')subplot(3,1,2);stem(n,x2)xlabel('n')ylabel('x2[n]')subplot(3,1,3);stem(n,x3)xlabel('n')ylabel('x3[n]')（二）实验图像：三、实验结论x1的信号周期是4，x2信号周期是4，x3周期是16第二题一、实验目的：1.确定离散时间信号的性质2.证明不满足线性性质 3.证明不是因果的 二、实验内容 （一） 1.试验程序]1[][][++=n x n x n y ()()][2/sin ][n x n y π=x1=[1 zeros(1,10)];x2=[2 zeros(1,10)];y1=sin((pi/2).*x1);y2=sin((pi/2).*x2);x3=2*x1+6*x2;y3=sin((pi/2).*x3);y4=2*y1+6*y2; subplot(2,1,1);stem(n,y3) title('figure of y3') xlabel('n')ylabel('y3')subplot(2,1,2);stem(n,y4) title('figure of y4') xlabel('n')ylabel('y4')2.实验图像（二）1.实验程序：nx=[-5:9];x1=[zeros(1,5) 1 ones(1,9)]; x2=[zeros(1,4) 1 ones(1,10)]; y=x1+x2; subplot(3,1,1); stem(nx,x1); xlabel('n'); ylabel('x1[n]'); subplot(3,1,2); stem(nx,x2); xlabel('n'); ylabel('x2[n]'); subplot(3,1,3) ; stem(nx,y); xlabel('n'); ylabel('y[n]');2.实验图象：三、实验结论1.系统不是线性的。
武大电气学院信号系统实验报告
《信号与系统》上机实验实验一连续时间信号的表示及可视化一．实验目的熟练掌握连续时间信号的表示及可视化处理。
二．实验源程序δf(t))=)(tf=@(t)dirac(t) %定义函数ezplot(f,[-5:5]); %利用ezplot( )命令绘制连续图形xlabel('(t)'); %横坐标ylabel('(f)'); %纵坐标f(t)= ε(t)（f=Heaviside(n)）f=@(t)heaviside(t) %定义函数ezplot(f,[-5:5]); %利用ezplot( )命令绘制连续图形xlabel('(t)'); %横坐标ylabel('(f)'); %纵坐标f(t)=at e(分别取a>0及a<0)a=1时f=@(t)exp(t) %定义函数ezplot(f,[-5:5]); %利用ezplot( )命令绘制连续图形xlabel('(t)'); %横坐标ylabel('(f)'); %纵坐标a=-1时f=@(t)exp(-t) %定义函数ezplot(f,[-5:5]); %利用ezplot( )命令绘制连续图形xlabel('(t)'); %横坐标ylabel('(f)'); %纵坐标f(t)=R(t)t=-5:0.01:5; %设定时间变量t的范围及步长y=rectpuls(t,2); %用rectpuls（t a)命令表示门函数，默认以零点为中心，宽度为aplot(t,y); %用plot函数绘制连续函数grid on; %显示网格命令title('门函数'); %用title函数设置图形的名称axis([-5 5 -0.5 1.5]);f(t)=Sa(wt)w=5时，f=Sa(5*t)f=@(t)Sinc(5*t) %定义函数ezplot(f,[-5:5]); %利用ezplot( )命令绘制连续图形xlabel('(t)'); %横坐标ylabel('(f)'); %纵坐标axis([-5 5 -1.2 1.2])w=8时，f=Sa(8*t)f=@(t)sinc(8*t) %定义函数ezplot(f,[-4:4]); %利用eaplot( )命令绘制连续图形xlabel('(t)'); %横坐标ylabel('(f)'); %纵坐标f(t)=Sin(2πft)(分别画出不同周期个数的波形)f(t)=Sin(t)f=@(t)sin(t) %定义函数ezplot(f,[-15:15]); %利用eaplot( )命令绘制连续图形xlabel('(t)'); %横坐标ylabel('(f)'); %纵坐标axis([-15 15 -1.2 1.2])三．程序运行结果(1)(2)（3）-5-4-3-2-1012345-1-0.8-0.6-0.4-0.200.20.40.60.81(t)dirac(t)(f )-5-4-3-2-101234500.20.40.60.81(t)heav iside(t)(f )（4）-5-4-3-2-1012345010********607080(t)exp(t)(f )-5-4-3-2-1012345010********607080(t)exp(-t)(f )（5）-5-4-3-2-1012345-1-0.8-0.6-0.4-0.200.20.40.60.81(t)Sinc(5 t)(f )(6)-4-3-2-101234 -1-0.8-0.6-0.4-0.20.20.40.60.81(t)sinc(8 t)(f)-15-10-5051015 -1-0.8-0.6-0.4-0.20.20.40.60.81(t)sin(t)(f)实验二离散时间信号的表示及可视化一．实验目的学会对离散时间信号进行标识和可视化处理。
信号与系统上机实验报告
姓名：学号：学院：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)结果：由图形可知周期T1=4由图形可知周期T2=4由图形可知周期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;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)并画图。
信号与系统上机实验1
实验四傅里叶变换(FT)及其性质一、实验目的1、学会运用Matlab求连续时间信号的傅里叶2、学会运用Matlab求连续时间信号的频谱图3、学会运用Matlab分析连续时间信号的傅里叶变换的性质二、实验原理及实例分析（一）傅里叶变换的实现例1：用Matlab 符号运算求解法求单边指数信号)()(2t u e t f t-=的FT 。
例2：用Matlab 符号运算求解法求211)(ωω+=j F 的IFT 。
例3：用Matlab命令绘出例1中单边指数数信号的频谱图。
例4：用Matlab命令求图示三角脉冲的FT，并画出其幅度谱。
例5：用Matlab数值计算法求例3的三角脉冲幅度频谱图。
（二）FT的性质1、尺度变换例6：设矩形信号)5.0tuf，利用Matlab命令=ttu)5.0+()-(-(绘出该信号及其频谱图。
同时绘出)2()2/f和的频谱图，并加以t(tf比较。
下面利用Matlab将常规矩形脉冲信号的频谱和其调制信号（课本例3-4信号）频谱进行比较。
Matlab源程序如下：傅里叶变换的其它性质可用类似的方法验证，希望大家课下练习。
三、实验内容[注意：(1)写代码时j i]3、分别利用Matlab符号运算求解法和数值计算法求下图所示信号的FT，并画出其频谱图。
4、已知门函数自身卷积为三角波信号，试用Matlab命令验证FT的时域卷积定理。
四、实验报告要求实验名称、实验目的、实验原理、实验环境、实验内容（上述几部分代码及结果图形）、实验思考等。
五、开发环境Matlab 7.0六、实验内容：1．（1）ft = sym('sin(2*pi*(t-1))/(pi*(t-1))');Fw = fourier(ft);subplot(2,2,1);ezplot(abs(Fw));grid on;title('幅度谱');phase = atan(imag(Fw)/real(Fw)); subplot(2,1,2);ezplot(phase);grid on;title('相位谱');幅度谱-505w相位谱w 2．（1）syms t;Fw = sym('10/(3+i*w) - 4/(5+i*w)'); ft = ifourier(Fw,t);ezplot(ft,t);-10123456x y x = 2 heaviside(t) (5 exp(-3 t)-2 exp(-5 t)), y = t3.f=sym('(t-2)*(Heaviside(t+2)-Heaviside(t+1))+Heaviside(t+1)-He aviside(t-1)+(-t+2)*(Heaviside(t-1)-Heaviside(t-2))');Fw = simplify(fourier(ft));subplot(2,1,1);ezplot(abs(Fw));grid on;title('幅度谱');phase = atan(imag(Fw)/real(Fw));subplot(2,1,2);ezplot(phase);grid on;title('相位谱');-6-4-20246w幅度谱-6-4-20246w 相位谱。
信号与系统-上机实验
信号与系统上机实验报告
2012029170009
丁启程
实验一：表示信号、系统的matlab函数、工具箱
用matlab表示基本单位冲激函数、单位阶跃函数、正弦函数及其相加、相乘和奇偶分离
单位冲激信号相加：
单位冲激信号相乘：
单位阶跃信号相加：
单位阶跃信号相乘：
正弦信号相乘：
X(n)=sin(pi*n/4)*cos(pi*n/4)的图像
根据图像可知，周期为4
实验二：离散时间的冲激响应、卷积和
利用conv函数计算卷积和，利用filter函数计算差分方程的系统输出并研究使用filter函数计算卷积的方法。
用conv、filter函数计算卷积
用filter函数对两个离散信号进行卷积的方法：
》x=[m:n] %m,n表示需要卷积的信号范围
》y=filter(y1,y2,x) %将y1，y2两个信号卷积，得到y
实验三：离散函数的转移函数、零极点图及模拟由离散信号的转移函数绘制零极点图并讨论其频率响应。
绘制系统零极点图。
信号与系统 实验报告2
信号与系统实验报告实验二连续时间LTI系统的时域分析一、实验目的1．学会用MATLAB 求解连续系统的零状态响应；2. 学会用MATLAB 求解冲激响应及阶跃响应；3．学会用MATLAB 实现连续信号卷积的方法；二、实验原理1．连续时间系统零状态响应的数值计算我们知道，LTI连续系统可用如下所示的线性常系数微分方程来描述，在MATLAB 中，控制系统工具箱提供了一个用于求解零初始条件微分方程数值解的函数lsim。
其调用格式y=lsim(sys,f,t) 式中，t表示计算系统响应的抽样点向量，f是系统输入信号向量，sys是LTI系统模型，用来表示微分方程，差分方程或状态方程。
其调用格式sys=tf(b,a) 式中，b 和a分别是微分方程的右端和左端系数向量。
例如，对于以下方程可用获得其LTI模型。
注意，如果微分方程的左端或右端表达式中有缺项，则其向量a或b 中的对应元素应为零，不能省略不写，否则出错。
2．连续时间系统冲激响应和阶跃响应的求解在MATLAB 中，对于连续LTI系统的冲激响应和阶跃响应，可分别用控制系统工具箱提供的函数impluse和step来求解。
其调用格式为y=impluse(sys,t)y=step(sys,t)式中，t表示计算系统响应的抽样点向量，sys是LTI系统模型。
3. 用MATLAB实现连续时间信号的卷积信号的卷积运算有符号算法和数值算法，此处采用数值计算法，需调用MATLAB 的conv( )函数近似计算信号的卷积积分。
连续信号的卷积积分定义是如果对连续信号进行等时间间隔∆均匀抽样，则分别变为离散时间信号其中，m为整数。
当∆足够小时，既为连续时间信号。
因此连续时间信号卷积积分可表示为采用数值计算时，只求当∆= nt 时卷积积分)(tf 的值)( ∆ nf ，其中，n为整数，既其中，实际就是离散序列的卷积和。
当∆足够小时，序列就是连续信号的数值近似，既上式表明，连续信号的卷积，可用各自抽样后的离散时间序列的卷积再乘以抽样间隔∆。
信号与系统上机实验(12级)
《信号与系统》上机实验实验指导（12级试用）实验1 信号变换与系统非时变性质的波形绘制1、实验目的掌握MATLAB画连续信号波形的方法，并能绘制信号变换和系统线性非时变性质的输出波形。
2、预习阅读内容《信号与系统实验教程》实验1、实验2。
3、参考程序《信号与系统实验教程》第2章例2.1-2程序（如何用MA TLAB画波形），例2.2-3程序（如何进行波形变换），例2.2-4程序（如何进行波形分解）。
4、MATLAB学习要点●用MA TLAB画波形的三要素（1）横轴赋值，两种方法t=0:0.1:30; t = 起点：间隔：终点，强调横轴的间隔。
t=linspace(0,30,300);t有300点，从起点0到终点30,等间隔自动计算，强调横轴的总点数。
（2）纵轴赋值f=sin(pi*t.*sign(t)); 根据所给函数写出MA TLAB表达式（3）波形绘制figure(1)，plot(t,f,'linewidth',2)●MATLAB的内部函数和外部函数（1）内部函数：MA TLAB软件内部就存在，调用就可以了，本次介绍几个常用函数门函数rectpuls（）调用格式 y=rectpuls(t) 产生高度为1,宽度为1的门函数调用格式 y=rectpuls(t,W) 产生高度为1,宽度为W的门函数三角脉冲函数tripuls（）调用格式 y=tripuls(t) 产生高度为1,宽度为1的三角脉冲函数调用格式 y=tripuls(t,w) 产生高度为1,宽度为w的三角脉冲函数调用格式 y=tripuls(t,w,s) 产生高度为1,宽度为w的三角脉冲函数,-1<s<1， s=0时为对称三角形; s=-1, 三角形顶点左边。
符号函数sign（）, 指数函数exp（）, 正弦函数sin（），余弦函数cos（）（2）外部函数：根据编程的需要，自己编写的函数。
介绍几个自编函数阶跃函数u( ), 文件名u.m，经常要用到，调用：y=u(t)function y=u(t)y=(t>=0);画图函数myplot( ),函数名为myplot.m ，可以自动将波形放在图的中间，使其美观。
华工电信学院信号与系统实验一报告参考模板
华工电信学院信号与系统实验信号与系统实验报告（一）实验项目名称：MATLAB 编程基础及典型实例 上机实验题目：信号的时域运算及MA TLAB 实现 一、实验目的学习并掌握使用MATLAB 产生基本信号、绘制信号波形、实现信号的可视化表示，为信号分析和系统设计奠定基础。
二、实验内容1. 利用Matlab 产生下列连续信号并作图。
(1) 51),1(2)(<<---=t t u t x(2) 2000,)8.0cos()1.0cos()(<<=t t t t x ππ 2. 利用Matlab 产生下列离散序列并作图。
(1) ⎩⎨⎧≤≤-=其他,055,1][k k x , 设1515-≤<k 。
(2) )]25.0cos()25.0[sin()9.0(][k k k x k ππ+=,设2020-≤<k 。
3. 已知序列]3,2,1,0,1,2;2,3,1,0,2,1[][--=-=k k x ， ]21,0,1,1,1[][=-=k k h 。
(1) 计算离散序列的卷积和][][][k h k x k y *=，并绘出其波形。
(2) 计算离散序列的相关函数][][][n k y k x k R k xy +=∑∞-∞=，并绘出其波形。
(3) 序列相关与序列卷积有何关系？三、实验细节1. 利用Matlab 产生下列连续信号并作图。
(1) 51),1(2)(<<---=t t u t xt=-1:0.01:5;x=-2.*((t-1)>=0); plot(t,x);axis([-1,5,-2.2,0.2])-112345-2-1.5-1-0.5(2) 2000,)8.0cos()1.0cos()(<<=t t t t x ππ t=0:2:200;x=cos(0.1*pi*t).*cos(0.8*pi*t); plot(t,x);20406080100120140160180200-1-0.8-0.6-0.4-0.200.20.40.60.812. 利用Matlab 产生下列离散序列并作图。
信号与系统(matlab上机实验报告)
通信与信息工程学院信号与系统分析MATLAB上机实验报告班级学号姓名电信1603班通信与信息工程学院二〇一八年实验题目一：信号的表示及可视化一、实验目的1.掌握连续信号的 MATLAB 表示方法（表达式及图形描述）；2. 掌握离散序列的 MATLAB 表示方法（表达式及图形描述）；二、实验原理向量法表示信号和符号法表示信号后直观的绘出图形1.连续信号的表示及可视化在 MATLAB 中，是用连续信号在等时间间隔点的样值来近似地表示连续信号的，当取样时间间隔足够小时，这些离散的样值就能较好的近似出连续信号。
（1）：向量法表示以用两个行向量 f 和 t 来表示，其中 t 向量是形如t=t1:p:t2 的 MATLAB 命令定义的时间范围向量，t1为信号起始时间，t2为终止时间，p 为时间间隔。
向量 f 为连续信号 f(t)在向量 t 所定义的时间点上的样值。
然后用plot画出图。
（2）：符号运算法表示：用 ezplot 命令绘制出信号的波形。
（3）：常用信号的表示：单位阶跃信号方法一：一种得到单位阶跃信号的方法是在 MATLAB 的 Symbolic MathToolbox 中调用单位阶跃函数 Heaviside，这样可方便地表示出单位阶跃信号。
但是，在用函数 ezplot 实现其可视化时，就出现一个问题：函数 ezplot 只能画出既存在于 Symbolic Math 工具箱中，又存在于总 MATLAB 工具箱中的函数，而 Heaviside 函数仅存在 Symbolic MathToolbox 中，因此，就需要在自己的工作目录 work 下创建 Heaviside 的M 文件。
方法二：用符号函数 sgn(t)的表示可调用 MATLAB 中的符号函数 sign 来实现。
单位冲激信号： MATLAB 是不能表示单位冲激信号的，但我们可用时间宽度为 dt ,高度为 1/dt的矩形脉冲来近似地表示冲激信号。
信号与系统-Matlab上机实验报告
江南大学物联网工程学院信号与系统—Matlab
上机实验报告
班级：物联网工程1302班
姓名：明坤
学号：1030613221
实验任务：
1.画出f(t) = 2sin200πt·cos20πt·e-2t + cos50πt·sin70πt的频谱。
2.画出f(t) = 2sin200πt·cos20πt·e-2t * cos50πt·sin70πt的波形和频谱。
源码：
%信号与系统Matlab作业
t=0:1e-3:100;%选择分析的时域范围离散化ft
%为方便计算把ft分成两部分
f1 = 2*sin( 200*pi*t ) .* cos( 20*pi*t ).*exp( -2*t ) ;
f2 = cos( 50*pi*t ) .* sin( 70*pi*t ) ;
ft_1=f1+f2;%第一题ft
y=fft(ft_1);%快速傅里叶变换
subplot(3,1,1);%分割绘图区域
plot(abs(y));%求出Fourier变换后的振幅
grid on;title('第一题频谱图');
ft_2=conv(f1,f2);%第二题ft 利用卷积函数
subplot(3,1,2);plot(ft_2);grid on;title('第二题时域波形');
axis([0 600 -5 5]);%为时域图选择一个合适的绘图范围
y=fft(ft_2);
subplot(3,1,3);plot(abs(y));grid on;title('第二题频谱图');
图形结果：。
信号与系统实验报告(实验11)
长江大学电工电子实验中心电路与系统(2)实验报告姓名班级序号指导教师成绩实验名称：非周期信号的频谱一、实验目的1.掌握用matlab绘制离散信号波形图的方法，学会常见波形的绘制。
2.掌握用matlab观测离散信号的周期性的方法。
3.掌握离散信号能量的计算方法。
二、实验内容实验11-1的程序 r11-1.mn=-5:15;f1=(1-0.8.^(n-3)).*u(n-3);figure(1),mystem(n,f1)n=-20:5;f2=(0.9.^(-n).*cos(n.*pi/8)).*u(-n+3);figure(2),mystem(n,f2);n=-12:12;f3=0.5*n.*rectpuls(n,16);figure(3),mystem(n,f3);n=-18:18;f4=sawtooth(2*n.*pi/14);figure(4),mystem(n,f4);-551015-20-15-10-55-10-50510-15-10-5051015 %实验11-2的程序r11-2.mclear k0=-7;k1=7; k=k0:k1; f=fd2(k+2);subplot(1,4,1),mystem(k,f), title('f(k+2)') f=fd2(k+2).*u(-k-2); subplot(1,4,2),mystem(k,f), title('f(k+2)*u(-k+2)') f=fd2(-k+2);subplot(1,4,3),mystem(k,f), title('f(-k+2)') f=fd2(-k+2).*u(k-1); subplot(1,4,4),mystem(k,f), title('f(-k+2)*u(k-1)')function f=fd2(k)f=delta(k+2)+delta(k+1)-1/3*(k-3).*(u(k)-u(k-3));-55f(k+2)-55f(k+2)*u(-k+2)-55f(-k+2)-55f(-k+2)*u(k-1)%实验11-3的程序r11-3.mfigure(1)k=-25:25;f=sin(k*pi/4)-2*cos(k*pi/6);mystem(k,f)[x,y]=ginput(2)gtext(['\bf 周期:N=',num2str(round(x(2)-x(1)))])figure(2)k=-25:25;f=cos(3*k/7-pi/8);mystem(k,f)-25-20-15-10-50510152025-25-20-15-10-50510152025%实验11-4的程序r11-4.mclearfigure(1)k=-6:6;f1=[0 0 0 0 0 0 6 4 2 2 0 0 0];mystem(k,f1)figure(2)f2=[0 0 0 0 -3 -2 -1 0 1 0 0 0 0];mystem(k,f2)figure(3)k=-15:15;f3=cos(pi*k/2);mystem(k,f3);figure(4)k=-2:15;f4=8*(0.5).^k.*u(k);mystem(k,f4)disp('离散信号的功率和能量’)E1=sum(abs(f1).^2)E2=sum(abs(f2).^2)k=0:3;f3=cos(pi*k/2);E3=sum(abs(f3).^2);P3=E3/4 E4=sum(abs(f4).^2)离散信号的功率和能量E1 =60E2 =15P3 =0.5000E4 =85.3333-6-4-20246-6-4-20246-15-10-5051015-202468101214三、上机实验的收获与体会通过进行matlab上机实验，学会了用matlab绘制离散信号波形图的方法，学会常见波形的绘制为进一步学习信号与系统分析打下了良好的基础。
- 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