《信号与系统》实验指导书

金陵科技学院实验一MATLAB语言工作环境和基本操作一、实验目的1．初步了解MATLAB开发环境和常用菜单的使用方法；2．熟悉MATLAB常用窗口，包括命令窗口、历史窗口、当前工作窗口、工作空间浏览器窗口、数组编辑器窗口和M文件编辑/调试窗口等；3．了解MATLAB的命令格式；4．熟悉MATLAB的帮助系统。二、实验原理1.简介MATLAB，MatrixLaboratory的缩写，是由Mathworks公司开发的一套用于科学工程计算的可视化高性能语言，具有强大的矩阵运算能力。与大家常用的Fortran和C等高级语言相比，MATLAB的语法规则更简单，更贴近人的思维方式，被称之为“草稿纸式的语言”。MATLAB安装完成后，会自动在WINDOWS桌面上生成一个快捷方式，它是指向安装目录下\bin\win32\matlab.exe的链接，双击它即可来到MATLAB集成环境的基本窗口，通常称之为命令窗口。MATLAB的退出与普通WIN32的程序一样，值得一提的是它有一个自身专有的快捷键Ctrl+Q。MATLAB启动后界面如图所示。命令窗口(CommandWindow)：(1)用于执行MATLAB命令，正常情况下提示符为“”，表示MATLAB进入工作状态。(2)在提示符后输入运算指令和函数调用等命令（不带“；”），MATLAB将迅速显示出结果并再次进入准备工作状态。(3)若命令后带有“；”，MATLAB执行命令后不显示结果。(4)在准备工作状态下，如果按上下键，MATLAB会按顺序依次显示以前输入的命令，若要执行它，则直接回车即可。工作空间(Workspace)：(1)显示计算机内存中现有变量的名称、类型、结构及其占用子节数等。(2)如果直接双击某变量，则弹出ArrayEditor窗口供用户查看及修改变量内容。(3)该窗口上有工具条支持用户将某变量存储到文件中或者从文件中载入某变量。命令历史记录(CommandHistory)：(1)保存并显示用户在命令窗口中输入过的命令，以及每次启动MATLAB的时间等信息(2)若双击某条命令记录，则MATLAB会再次执行该命令。当前路径窗口(CurrentDirectory)：(1)先是当前路径内的所有文件。(2)用户可以在这里新建或删除一个文件，也可以双击一个文件，在编辑/调试窗口中打开。设置当前路径(CurrentDirectory)：(1)用于选择当前工作路径。(2)可以在命令窗口中输入文件名来直接调用工作路径下的文件。使用MATLAB时获取帮助的两种方法：(1)直接在命令窗口中输入help函数名或命令；(2)在帮助窗口中浏览或搜索相应信息。还可以参考MATLAB的Demo程序来学习MATLAB编程。2.基本语法标识符标识符是标识变量名、常量名、函数名和文件名的字符串的总称。在MATLAB中，变量和常量的标志符最长允许19个字符；函数和文件名通常不超过8个字符（受操作系统文件管理器的限制）。这些字符包括全部的英文字母（区分大小写）、阿拉伯数字和下划线等符号。标识符的第一个字符必须是英文字母。（1）变量及其赋值变量代表一个活若干个内存单元。变量名必须是以字母开头，后接字母、数字或下划线的字符序列，最多63个字符。而且区分大小写。赋值：变量名=表达式（2）矩阵及其赋值最简单的矩阵建立方法就是直接输入矩阵的元素。具体方法是：将矩阵的所有元素用方括号括起来，按矩阵行的顺序输入各元素，同一行的各元素之间用空格或逗号分隔，不同行的元素之间用分号分隔。（3）数组及其赋值数组通常是指单行或单列的矩阵。MATLAB提供两种为等间隔数组赋值的简易方法。①用两个冒号组成等增量语句，其格式为：t=[初值：增量：终值]或者t=初值：增量：终值。②用linspace函数。调用格式为：linspace(初值，终值，点数)。如：（4）元素群的四则运算和幂次运算元素群运算也就是矩阵中所有元素按照单个元素进行运算。在运算符*、/、\、^前加上一个点符号“.”，就表示做元素群运算。参与元素群运算的两个矩阵必须是同阶的（只有标量除外，它会自动扩展为同阶矩阵参与运算）。运算式输出结果Z=X.*YZ=41018Z=X.\YZ=4.00002.50002.0000Z=X.^YZ=132729Z=X.^2Z=149Z=2.^[XY]Z=248163264运算符MATLAB中的数学运算符名称说明名称说明+-矩阵加，矩阵减*矩阵乘/\矩阵右除，矩阵左除^矩阵求幂.*.^群乘，群求幂./.\数组右除，数组左除'.'共轭转置，转置=赋值MATLAB中的比较和逻辑运算符名称说明名称说明==等于~=不等于=大于，大于等于=小于，小于等于&与|或~非xor(a,b)a和b异或any(a)a中有元素非零则为真all(a)a中所有元素都非零则为真3.M文件及程序调试三、主要实验仪器及材料微型计算机、Matlab6.5以上版本的编程环境。四、实验内容1.熟悉MATLAB开发环境和常用菜单的使用方法；了解help功能使用方法。2.赋值一个4阶矩阵，使其主对角线上预算都为0，与主对角线相邻元素都为1，其余为2；3.编写一个m文件，求自然数前N项和。N=100。五、思考题1.MATLAB的文件名和变量名有何要求？2.如何绘制双坐标系？六、实验报告要求1.简述实验原理及目的。主程序文件(scriptfile)（主程序）：用户为解决特定问题而编制函数文件(functionfile)（子程序）：必须由其他M文件调用；往往具有一定的通用性M文件由MATLAB语句构成的文本文件，以m作为文件扩展名；是ASCII码文本文件，可以直接阅读并用任何编辑器来建立。行命令模式：键入一行命令后，让系统立即执行该命令。缺点：程序可读性差，难以存储。运行模式程序文件模式：可存储，尤其适于复杂问题2.记录调试运行情况及所遇问题的解决方法。3.简要回答思考题。«信号与系统»课程实验指导书实验一常用连续信号的表示适用专业：通信工程,电子信息工程学时：2学时一、实验目的熟悉MATLAB软件；利用MATLAB软件，绘制出常用的连续时间信号。二、实验内容1.绘制正弦信号)tAsint(f0（），其中A=1，2，6/；例：A=1;w0=2*pi;phi=pi/6;t=0:0.001:8;ft=A*sin(w0*t+phi);plot(t,ft);gridon;2.绘制指数信号atAet(f），其中A=1，0.4a；调用形式：A*exp(a*t)3.绘制矩形脉冲信号，脉冲宽度为2；调用函数：rectpuls(t-2*T,2*T)4.绘制三角波脉冲信号，脉冲宽度为4；斜度为0.5；调用函数：tripuls(t,4,0.5)5.对上题三角波脉冲信号进行尺度变换，分别得出)2t(f，)2t2(f；6.绘制抽样函数Sa（t）,t取值在-3到+3之间；（调用函数sinc）7.绘制周期矩形脉冲信号，参数自定；调用函数：square(2*pi*30*t,75)8.绘制周期三角脉冲信号，参数自定。调用函数：sawtooth(t,0.5)三、实验设备586以上计算机，装有MATLAB7.0软件四、实验原理和方法上机操作五、实验步骤1、打开MATLAB界面，建立新文件。2、根据实验要求，编写程序。3、调试运行程序，并修改完善。4、保存设计成果，撰写实验报告。六、实验注意事项爱护机器，遵守规章制度金陵科技学院«信号与系统»课程实验指导书实验二连续时间系统分析适用专业：通信工程,电子信息工程学时：2学时一、实验目的利用MATLAB软件对系统进行冲激响应、零状态响应、频率响应的求解并领会其物理意义；掌握周期信号的傅立叶级数系数以及功率谱的计算。二、实验内容1.一力学系统，其系统微分方程为(t)100100y(t))t('2y)t(''y，求该系统冲激响应；sys=tf([10],[12100]);t=0:0.01:5;y=impulse(sys,t);plot(t,y);2.一系统微分方程为f(t)100y(t))t('2y)t(''y，若t10sin2)t(f，求该系统零状态响应；sys=tf([1],[12100]);f=10*sin(2*pi*t);y=lsim(sys,f,t);3.一系统微分方程为3x(t)3y(t)t)('y，求该系统的频率响应（绘出其幅值和相位特性曲线）；a=[13];b=3;freqs(b,a);4.一系统微分方程为5x(t)(t)'x'y(t))t('4yt)(''3y，求该系统的频率响应（绘出其幅值和相位特性曲线）；b=[105];5.计算教材例3-1-2周期三角波信号的傅立叶级数系数并绘出其频谱。n1=-N:-1;c1=-4*j*sin(n1*pi/2)/pi^2./n1.^2;%计算n=-N到－1时的傅立叶级数系数c0=0;%计算n=0到N时的傅立叶级数系数n2=1:N;c2=-4*j*sin(n2*pi/2)/pi^2./n2.^2;%计算n=1到N时的傅立叶级数系数cn=[c1c0c2];stem(n,abs(cn));stem(n,angle(cn));6.计算信号t)fcos(2t)f2(cos)t(x21的功率并绘制其功率谱，其中f1=37Hz,f2=219Hz。ts=0.001;fs=1/ts;t=[0:0.001:10];x=cos(2*pi*47*t)+cos(2*pi*219*t);p=(norm(x)^2)/length(x);psd=spectrum(x,1024);specplot(psd,fs)三、实验设备586以上计算机，装有MATLAB7.0软件四、实验原理和方法上机操作五、实验步骤1、打开MATLAB界面，建立新文件。2、根据实验要求，编写程序。3、调试运行程序，并修改完善。4、保存设计成果，撰写实验报告。六、实验注意事项爱护机器，遵守规章制度金陵科技学院«信号与系统»课程实验指导书实验三连续时间系统的复频域分析适用专业：通信工程,电子信息工程学时：2学时一、实验目的利用MATLAB软件掌握Laplace正变换与反变换，对F（s）进行部分分式展开，求系统函数的零极点并画出零极点分布图，并判断系统是否稳定。二、实验内容1.求函数tf(t)esin(at)(t)u的Laplace变换；f=sym('exp(-t)*sin(a*t)');F=laplace(f);pauseF2.求1)/(ssF(s)22的Laplace反变换；F=sym('s^2/(s^2+1)');ft=ilaplace(F);3.对3s)4s2)/(ss(F(s)23进行部分分式展开；formatrat;%将分数以近似的整数之比的形式显示num=[12];%num为分子多项式的系数向量den=[1430];%den为分母多项式的系数向量[r,p,k]=residue(num,den);%r为所得部分分式展开式的系数向量，p为极点，k为分式的直流分量4.求出系统函数)22ss/()1s()s(H2的零极点并画出零极点分布图；b=[1-1];a=[122];zs=roots(b);ps=roots(a);plot(real(zs),imag(zs),'o',real(ps),imag(ps),'kx','markersize',12);b=[1-1];a=[122];sys=tf(b,a);pzmap(sys);5.画出系统函数)12s2ss/(1)s(H23的零极点分布图，求出系统的单位冲激响应h(t)和幅频响应)j(H，并判断系统是否稳定。num=[1];den=[1221];sys=tf(num,den);poles=roots(den);figure(1);%系统函数的零极点分布图pzmap(sys);t=0:0.02:10;h=impulse(num,den,t);figure(2);%系统的单位冲激响应plot(t,h);title('ImpulseRes