基于Matlab的信号与系统实验指导

1、基于Matlab的信号与系统实验指导实验一 连续时间信号在Matlab中的表示一、实验目的1、学会运用Matlab表示常用连续时间信号的方法2、观察并熟悉这些信号的波形和特性二、实验原理及实例分析1、信号的定义与分类2、如何表示连续信号？连续信号的表示方法有两种；符号推理法和数值法。从严格意义上讲，Matlab数值计算的方法不能处理连续时间信号。然而，可利用连续信号在等时间间隔点的取样值来近似表示连续信号，即当取样时间间隔足够小时，这些离散样值能被Matlab处理，并且能较好地近似表示连续信号。3、Matlab提供了大量生成基本信号的函数。如：（1）指数信号：K*exp(a*t)（2）正弦信号

2、：K*sin(w*t+phi)和K*cos(w*t+phi)（3）复指数信号：K*exp(a+i*b)*t)（4）抽样信号：sin(t*pi)注意：在Matlab中用与Sa(t)类似的sinc(t)函数表示，定义为： （5）矩形脉冲信号：rectpuls(t,width)（6）周期矩形脉冲信号：square(t,DUTY)，其中DUTY参数表示信号的占空比DUTY%，即在一个周期脉冲宽度（正值部分）与脉冲周期的比值。占空比默认为0.5。（7）三角波脉冲信号：tripuls(t, width, skew)，其中skew取值范围在-1+1之间。（8）周期三角波信号：sawtooth(t, widt

3、h)（9）单位阶跃信号：y=(t=0)三、实验内容1、验证实验内容直流及上述9个信号2、程序设计实验内容（1）利用Matlab命令画出下列连续信号的波形图。（a）（b）（c）（2）利用Matlab命令画出复信号的实部、虚部、模和辐角。四、实验报告要求1、格式：实验名称、实验目的、实验原理、实验环境、实验内容、实验思考等2、实验内容：程序设计实验部分源代码及运行结果图示。实验二 连续时间信号在Matlab中的运算一、实验目的1、学会运用Matlab进行连续时间信号的时移、反褶和尺度变换；2、学会运用Matlab进行连续时间信号微分、积分运算；3、学会运用Matlab进行连续时间信号相加、相乘运算

4、；4、学会运用Matlab进行连续时间信号卷积运算。二、实验原理及实例分析1、信号的时移、反褶和尺度变换信号的平移、反转和尺度变换是针对自变量时间而言的，其数学表达式和波形变换中存在着一定的变化规律。从数学表达式上来看，信号的上述所有计算都是自变量的替换过程。所以在使用Matlab进行连续时间信号的运算时，只需要进行相应的变量代换即可完成相关工作。2、连续时间信号的微分和积分符号运算工具箱有强大的积分运算和求导功能。连续时间信号的微分运算，可使用diff命令函数来完成，其语句格式为：diff(function, variable,n)其中，function表示需要进行求导运算的函数，或者被赋值

5、的符号表达式；variable为求导运算的独立变量；n为求导阶数，默认值为一阶导数。连续时间信号积分运算可以使用int命令函数来完成，其语句格式为：int(function, variable, a, b)其中，function表示被积函数，或者被赋值的符号表达式；variable为积分变量；a为积分下限，b为积分上限，a和b默认时则求不定积分。3、信号的相加和相乘运算信号的相加和相乘是信号在同一时刻取值的相加和相乘。因此Matlab对于时间信号的相加和相乘都是基于向量的点运算。4、连续信号的卷积运算卷积积分是信号与系统时域分析的重要方法之一。定义为：Matlab进行卷积计算可通过符号运算方法

6、和数值计算方法实现。（1）Matlab符号运算法求连续信号卷积从卷积定义出发，可以利用Matlab符号运算法求卷积积分，但要注意积分变量和积分限的选取。例：试用Matlab符号运算法求卷积y(t)=u(t)-u(t-1)*u(t)-u(t-1)。（2）Matlab数值计算法求连续信号的卷积例：试用Matlab数值计算法求信号和的卷积。三、实验内容1、已知信号的波形（课本P11例题），画出的波形图。2、使用微分命令求关于变量x的一阶导数；使用积分命令计算不定积分 ，定积分。3、已知，使用命令画出两信号和及两信号乘积的波形图。其中，4、四、实验报告要求1、格式：实验名称、实验目的、实验原理、实验环

7、境、实验内容（上述4部分代码及结果图形）、实验思考等。实验三 连续时间LTI系统的时域分析一、实验目的1、学会使用符号法求解连续系统的零输入响应和零状态响应2、学会使用数值法求解连续系统的零状态响应3、学会求解连续系统的冲激响应和阶跃响应二、实验原理及实例分析1、连续时间系统零输入响应和零状态响应的符号求解连续时间系统可以使用常系数微分方程来描述，其完全响应由零输入响应和零状态响应组成。MATLAB符号工具箱提供了dsolve函数，可以实现对常系数微分方程的符号求解，其调用格式为：dsolve(eq1,eq2,cond1,cond2,v)其中参数eq表示各个微分方程，它与MATLAB符号表达式

8、的输入基本相同，微分和导数的输入是使用Dy，D2y，D3y来表示y的一价导数，二阶导数，三阶导数；参数cond表示初始条件或者起始条件；参数v表示自变量，默认是变量t。通过使用dsolve函数可以求出系统微分方程的零输入响应和零状态响应，进而求出完全响应。 实例1试用Matlab命令求齐次微分方程的零输入响应，已知起始条件为。注意，程序中绘图的时间区间一定要t0，本程序中取0, 8，程序运行后结果如下。2、连续时间系统零状态响应的数值求解在实际工程中使用较多的是数值求解微分方程。对于零输入响应来说，其数值解可以通过函数initial来实现，而该函数中的参量必须是状态变量所描述的系统模型，由于现

9、在还没有学习状态变量相关内容，所以此处不做说明。对于零状态响应，MATLAB控制系统工具箱提供了对LTI系统的零状态响应进行数值仿真的函数lsim，利用该函数可以求解零初始条件下的微分方程的数值解。其调用格式为：y=lsim(sys,f,t)，其中t表示系统响应的时间抽样点向量，f是系统的输入向量；sys表示LTI系统模型，用来表示微分方程、差分方程或状态方程。在求解微分方程时，sys是有tf函数根据微分方程系数生成的系统函数对象，其语句格式为：sys=tf(a,b)。其中，a和b分别为微分方程右端和左端的系数向量。例如，对于微分方程可以使用获得其LTI模型。注意，如果微分方程的左端或者右端表

10、达式有缺项，则其向量a或者b中对应元素应该为零，不能省略不写。3、连续时间系统冲激响应和阶跃响应的求解在连续时间LTI系统中，冲激响应和阶跃响应是系统特性的描述。在MATLAB中，对于冲激响应和阶跃响应的数值求解，可以使用控制工具箱中提供的函数impulse和step来求解。其中t表示系统响应的时间抽样点向量，sys表示LTI系统模型。三、实验内容1、已知系统的微分方程和激励信号，使用MATLAB命令画出系统的零状态响应和零输入响应（零状态响应分别使用符号法和数值法求解，零输入响应只使用符号法求解）。要求题目2必做，题目1选做。2、已知系统的微分方程，使用MATLAB命令画出系统的冲激响应和阶

11、跃响应（数值法）。要求题目2必做，题目1选做。四、实验报告要求1、格式：实验名称、实验目的、实验原理、实验环境、实验内容（上述几部分代码及结果图形）、实验思考等。实验四 傅里叶变换(FT)及其性质一、实验目的1、学会运用Matlab求连续时间信号的傅里叶2、学会运用Matlab求连续时间信号的频谱图3、学会运用Matlab分析连续时间信号的傅里叶变换的性质二、实验原理及实例分析（一）傅里叶变换的实现 例1：用Matlab符号运算求解法求单边指数信号的FT。例2：用Matlab符号运算求解法求的IFT。例3：用Matlab命令绘出例1中单边指数数信号的频谱图。例4：用Matlab命令求图示三角脉

12、冲的FT，并画出其幅度谱。例5：用Matlab数值计算法求例3的三角脉冲幅度频谱图。（二）FT的性质1、尺度变换例6：设矩形信号，利用Matlab命令绘出该信号及其频谱图。同时绘出的频谱图，并加以比较。下面利用Matlab将常规矩形脉冲信号的频谱和其调制信号（课本例3-4信号）频谱进行比较。Matlab源程序如下：傅里叶变换的其它性质可用类似的方法验证，希望大家课下练习。三、实验内容 注意：(1)写代码时ji3、分别利用Matlab符号运算求解法和数值计算法求下图所示信号的FT，并画出其频谱图。4、已知门函数自身卷积为三角波信号，试用Matlab命令验证FT的时域卷积定理。四、实验报告要求实验

13、名称、实验目的、实验原理、实验环境、实验内容（上述几部分代码及结果图形）、实验思考等。实验五 信号抽样及抽样定理一、实验目的学会运用MATLAB完成信号抽样以及对抽样信号的频谱进行分析学会运用MATLAB改变抽样时间间隔，观察抽样后信号的频谱变化学会运用MATLAB对抽样后的信号进行重建二、实验原理及实例分析（一）信号抽样 信号抽样是利用抽样脉冲序列从连续信号中抽取一系列的离散值，通过抽样过程得到的离散值信号称为抽样信号，记为。从数学上讲，抽样过程就是信号相乘的过程，即因此，可以使用傅里叶变换的频域卷积性质来求抽样信号的频谱。常用的抽样脉冲序列有周期矩形脉冲序列和周期冲激脉冲序列。上式表明，信

14、号在时域被抽样后，它的频谱是原连续信号频谱以抽样角频率为间隔周期的延拓，即信号在时域抽样或离散化，相当于频域周期化。在频谱的周期重复过程中，其频谱幅度受抽样脉冲序列的傅里叶系数加权，即被加权。可以看出，是以为周期等幅地重复。程序运行结果，如下页图示。 很明显，升余弦脉冲信号的频谱抽样后发生了周期延拓，频域上该周期为。（二）抽样定理如果是带限信号，带宽为，则信号可以用等间隔的抽样值来唯一表示。经过抽样后的频谱就是将的频谱在频率轴上以抽样频率为间隔进行周期延拓。因此，当时，周期延拓后频谱不会产生频率混叠；当时，周期延拓后频谱将产生频率混叠。通常把满足抽样定理要求的最低抽样频率称为奈奎斯特频率，把最

15、大允许的抽样间隔称为奈奎斯特间隔。（三）信号重建抽样定理表明，当抽样定理小于奈奎斯特间隔时，可以使用抽样信号唯一表示原信号，即信号的重建。为了从频谱中无失真的恢复原信号，可以采用截止频率为的理想低通滤波器。上式表明连续信号可展开为抽样函数的无穷级数，该级数的系数为抽样值。利用MATLAB中的函数来表示，所以可获得由重建的表达式，即实例9-3 对实例9-1中的升余弦脉冲信号，假设其截止频率为，抽样间隔，采用截止频率的低通滤波器对抽样信号滤波后重建信号，并计算重建信号与原升余弦脉冲信号的绝对误差。程序运行结果如下图所示。从图中可以看出，重建后的信号与原升余弦脉冲信号的误差在以内，因为当选取升余弦脉

16、冲信号带宽为时，实际上已经将很少的高频分量忽略了。程序运行结果如下图所示。结果表明信号不满足抽样定理时，会产生较大的失真，并且绝对误差十分明显。三、实验内容1、设有三个不同频率的正弦信号，频率分别为，；现在使用抽样频率对这三个信号进行抽样，使用MATLAB命令画出各抽样信号的波形和频谱，并分析其频率混叠现象。2、结合抽样定理，利用MATLAB编程实现信号经过冲激脉冲抽样后得到的抽样信号及其频谱，并利用构建信号。四、实验报告要求实验名称、实验目的、实验原理、实验环境、实验内容（上述几部分代码及结果图形）、实验思考等。实验六 连续时间LTI系统的频率特性及频域分析一、实验目的1、学会运用MATLA

17、B分析连续系统地频率特性2、学会运用MATLAB进行连续系统的频域分析二、实验原理及实例分析（一）连续时间LTI系统的频率特性 一个连续LTI系统的数学模型通常用常系数线性微分方程描述，即 （*）对上式两边取傅里叶变换，并根据FT的时域微分性质可得：定义为：可见为两个的多项式之比。其中，分母、分子多项式的系数分别为（*）式左边与右边相应项的系数，称为系统的系统函数，也称为系统的频率响应特性，简称系统频率响应或频率特性。一般是复函数，可表示为：其中，称为系统的幅频响应特性，简称为幅频响应或幅频特性；称为系统的相频响应特性，简称相频响应或相频特性。描述了系统响应的傅里叶变换与激励的傅里叶变换间的关

18、系。只与系统本身的特性有关，与激励无关，因此它是表征系统特性的一个重要参数。MATLAB信号处理工具箱提供的freqs函数可直接计算系统的频率响应的数值解，其语句格式为：H=freqs(b,a,w)其中，b和a表示的分子和分母多项式的系数向量；w为系统频率响应的频率范围，其一般形式为w1:p:w2，w1为频率起始值，w2为频率终止值，p为频率取值间隔。H返回w所定义的频率点上系统频率响应的样值。注意，H返回的样值可能为包含实部和虚部的复数。因此，如果想得到系统的幅频特性和相频特性，还需要利用abs和angle函数来分别求得。例1：已知某连续LTI系统的微分方程为：求该系统的频率响应，并用MAT

19、LAB绘出其幅频特性和相频特性图。解：对上式两端取FT，得：因此，频率响应为：利用MATLAB中的freqs函数可求出其数值解，并绘出其幅频特性和相频特性图。MATLAB源程序和程序运行结果如下例2：下图是实用带通滤波器的一种最简单形式。试求当时该滤波器的幅频特性和相频特性。解：带通滤波器的频率响应为：代入参数，带通滤波器的谐振频率为：带通滤波器的幅频特性和相频特性的MATLAB源程序如下：程序运行结果如上右图所示，可以看到，该带通滤波器的特性是让接近谐振频率的信号通过而阻止其它频率的信号。（二）连续时间LTI系统的频率特性 连续LTI系统的频域分析法，也称为傅里叶变换分析法。该方法是基于信号

20、频谱分析的概念，讨论信号作用于线性系统时在频域中求解响应的方法。傅里叶分析法的关键是求取系统的频率响应。傅里叶分析法主要用于分析系统的频率响应特性，或分析输出信号的频谱，也可用来求解正弦信号作用下的正弦稳态响应。下面通过实例来研究非周期信号激励下利用频率响应求零状态响应。例3：下图(a)为RC低通滤波器，在输入端加入矩形脉冲如图(b)所示，利用傅里叶分析法求输出端电压。解：RC低通滤波器的频率响应为：，其中激励信号的FT为：因此，响应的FT为：MATLAB源程序和程序运行结果如下所示：由上图可看出，时域中输出信号与输入信号的波形产生了失真，表现在波形的上升和下降部分，输出信号的波形上升和下降部

21、分比输入波形要平缓许多。而在频域，激励信号频谱的高频分量与低频分量相比受到较严重的衰减。这正是低通滤波器所起的作用。对于周期信号激励而言，可首先将周期信号进行傅里叶级数展开，然后求系统在各傅里叶级数分解的频率分量作用下系统的稳态响应分量，再由系统的线性性质将这些稳态响应分量叠加，从而得到系统总的响应。该方法的理论基础是基于正弦信号作用下系统的正弦稳态响应。对于正弦激励信号，当经过系统，其稳态响应为：例4：设系统的频率响应为，若外加激励信号为，用MATLAB命令求其稳态响应。解：MATLAB源程序和程序运行结果如下：从图形可看出，信号通过该系统后，其高频分量衰减较大，说明该系统是低通滤波器。三、

22、实验内容1、试用MATLAB命令求下图所示电路系统的幅频特性和相频特性。已知。2、已知系统微分方程和激励信号如下，试用MATLAB命令求系统的稳态响应。（1）（2）四、实验报告要求实验名称、实验目的、实验原理、实验环境、实验内容（上述几部分代码及结果图形）、实验思考等。实验七 拉普拉斯变换（LT）一、实验目的1、学会运用MATLAB求拉普拉斯变换（LT）2、学会运用MATLAB求拉普拉斯反变换（ILT）二、实验原理及实例分析 LT是分析连续信号与系统的重要方法。运用LT可以将连续LTI系统的时域模型简便地进行变换，经求解在还原为时域解。从数学角度看，LT是求解常系数线性微分方程的工具。由LT导

23、出的系统函数对系统特性分析也具有重要意义。（一）拉普拉斯变换（LT） 对于一些不满足绝对可积条件的时域信号，是不存在傅里叶变换的。为了使更多的函数存在变换，并简化某些变换形式或运算过程，引入衰减因子，其中，为任意实数，使得满足绝对可积条件，从而求的傅里叶变换，即把频域扩展为复频域。连续时间信号的LT定义为： （*）ILT定义为： （*）式（*）和（*）构成了拉普拉斯变换对，称为的像函数，而称为的原函数。可以将拉普拉斯变换理解为广义的傅里叶变换。考虑到实际问题，人们用物理手段和实验方法所能记录和处理的一切信号都是有起始时刻的，对于这类单边信号或因果信号，我们引入单边LT，定义为：如果连续信号可用

24、符号表达式表示，则可用MATLAB的符号数学工具箱中的laplace函数来实现其单边LT，其语句格式为：。式中L返回的是默认符号为自变量s的符号表达式，f则为时域符号表达式，可通过sym函数来定义。例1：用MATLAB的laplace函数求的FT。解：MATLAB的源程序为：f=sym(exp(-t)*sin(a*t);L=laplace(f)或syms a tL=laplace(exp(-t)*sin(a*t);laplace函数另一种语句格式为：。它返回的函数L是关于符号对象v的函数，而不是默认的s。对上例中如果要求FT后的表达式自变量为v，则MATLAB源程序为：syms a t vf=

25、exp(-t)*sin(a*t);L=laplace(f,v)注：请自行验证结果正确与否（二）拉普拉斯反变换（ILT） 1、基于MATLAB符号数学工具箱实现ILT如果连续信号可用符号表达式表示，则可用MATLAB的符号数学工具箱中的ilaplace函数来实现其ILT，其语句格式为：。式中f返回的是默认符号为自变量t的符号表达式，L则为s域符号表达式，也可通过sym函数来定义。例2：试用MATLAB的ilaplace函数求的ILT。解：MATLAB源程序为：F=sym(s2/(s2+1);ft=ilaplace(F)或syms sft=ilaplace(s2/(s2+1)注：请自行验证结果正确

26、与否2、基于MATLAB部分分式展开法实现ILT用MATLAB函数residue可得到复杂有理式F(s)的部分分式展开式，其语句格式为：其中B、A分别表示F(s)的分子和分母多项式的系数向量；r为部分分式的系数；p为极点；k为F(s)中整式部分的系数。若F(s)为有理真分式，则k为0。例3：利用MATLAB部分分式展开法求的ILT。解：MATLAB源程序为：format rat;B=1,2;A=1,4,3,0;r,p=residue(B,A)程序中的format rat是将结果数据以分数的形式表示，其运行结果为：r= -1/6-1/22/3p=-3-10从上述结果可知，F(s)有3个单实极点，

27、即，其对应部分分式展开系数为：-1/6、-1/2、2/3。因此，F(s)可展开为：。所以，F(s)的反变换为：例4：利用MATLAB部分分式展开法求的ILT。解：F(s)的分母不是标准的多项式形式，可利用MATLAB的conv函数将因子相乘的形式转换为多项式的形式，其MATLAB源程序为：B=1,-2;A=conv(conv(1,0,1,1),conv(1,1,1,1);r,p=residue(B,A)程序运行结果（略）根据程序运行结果，F(s )可展开为：所以，F(s)的ILT为：（三）拉普拉斯变换法求解微分方程拉普拉斯变换法是分析连续LTI系统的重要手段。LT将时域中的常系数线性微分方程，

28、变换为复频域中的线性代数方程，而且系统的起始条件同时体现在该代数方程中，因而大大简化了微分方程的求解。借助MATLAB符号数学工具箱实现拉普拉斯正反变换的方法可以求解微分方程，即求得系统的完全响应。例5：已知某连续LTI系统的微分方程为：，且已知激励信号，起始条件为，求系统的零输入响应、零状态响应和全响应。解：对原方程两边进行拉普拉斯变换，并利用起始条件，得：将起始条件及激励变换代入整理可得：其中，第一项为零输入响应的拉普拉斯变换，第二项为零状态响应的拉普拉斯变换。利用MATLAB求其时域解，源程序如下：syms t sYzis=(3*s+13)/(s2+3*s+2);yzi=ilaplace

29、(Yzis)yzi= -7*exp(-2*t)+10*exp(t)xt=4*exp(-2*t)*Heaviside(t);Xs=laplace(xt);Yzss=Xs/(s2+3*s+2);yzs=ilaplace(Yzss)yzs= 4*(-1-t)*exp(-2*t)+4*exp(-t)yt=simplify(yzi+yzs)yt= -11*exp(-2*t)+14*exp(-t)-4*t*exp(-2*t)系统的零输入响应为：系统的零状态响应为：系统的完全响应为：三、实验内容1、试用MATLAB命令求课本习题4-1各函数的LT。注：要求至少任选2个小题。2、试用MATLAB命令求课本习题

30、4-4各函数的ILT。注：要求至少任选2个小题，且分别用两种方法求解。3、试用MATLAB命令和拉普拉斯变换法求课本习题2-6。注：要求至少任选一种情况求解。四、实验报告要求实验八 Z变换及离散时间系统的Z域分析一、实验目的1、学会运用MATLAB求离散时间信号的z变换和z反变换2、学会运用MATLAB分析离散时间系统的系统函数的零极点3、学会运用MATLAB分析系统函数的零极点分布与其时域特性的关系4、学会运用MATLAB进行离散时间系统的频率特性分析二、实验原理及实例分析（一）Z正反变换序列的双边z变换定义为：序列的单边z变换定义为：MATLAB符号数学工具箱提供了计算离散时间信号的单边z

31、变换的函数ztrans和z反变换函数iztrans，其语句格式分别为：Z=ztrans(x)X=iztrans(Z)上式中的x和Z分别为时域表达式和z域表达式的符号表示，可以通过sym函数来定义。x=sym(an*cos(pi*n);z=ztrans(x);simplify(z)ans=z/(z+a)(2) z变换的MATLAB程序为x=sym(2(n-1)-(-2)(n-1);Z=ztrans(x);simplify(z)ans=z2/(z-2)/(z+2)如果信号的z域表示式是有理数，则进行z变换的另外一个办法就是对X(z)进行部分分式展开，然后求各简单分式的z变换。设X(z)的有理分式表

32、示为MATLAB信号工具箱提供了一个对X(z)进行部分分式展开的函数residuez，其语句格式为R,P,K=residuez(B,A)其中，B，A分别表示X(z)的分子与分母多项式的系数向量，R为部分分式的系数向量，P为极点向量，K为多项式的系数。若X(z)为有理真分式，则K为0。（二）系统函数的零极点分析离散时间系统的系统函数定义为如果系统函数的有理函数表达式为在MATLAB中系统函数的零极点就可以通过函数roots得到，也可以借助函数tf2zp得到，tf2zp的语句格式为:Z,P,K=tf2zp(B,A)其中，B与A分别表示为H(z)的分子与分母多项式的系数向量。它的作用是将H(z)的有

33、理分式表示为转换为零极点增益形式：（三）系统函数的零极点分布与其时域特性的关系在离散系统中，z变换建立了时域函数h(n)与z域函数H(z)之间的关系。因此，H(z)从形式上可以反映h(n)的部分内在性质。下面通过讨论H(z)的一阶极点情况，来说明系统函数的零极点分布与系统时域特性的关系。从图14-2可以看出，当极点位于单位圆内时，h(n)为衰减序列；当极点位于单位圆上时，h(n)为等幅序列；当极点位于单位圆外时，h(n)为增幅序列。若H(z)有一价实数极点，则h(n)为指数序列；若H(z)为一价共轭极点，则h(n)为指数震荡序列；则H(z)的极点位于虚轴左边，则h(z)序列按一正一负的规律交替

34、变化。（四）离散时间LTI系统的频率特性分析对于因果稳定的离散时间系统，如果激励为正弦序列，则系统的稳态响应为。其中，通常为复数。离散时间系统的频率响应定义为，其中，称为离散时间系统的幅频响应特性；称为离散时间系统的相频响应特性；是以为周期的周期函数。因此，只要分析在范围内的情况，便可知道整个系统在频域的特性。MATLAB提供了求离散时间系统频响特性的函数freqz，调用freqz的格式主要有两种形式。一种形式为：H,w=freqz(B,A,N), 其中B与A分别表示H(z)的分子与分母多项式的系数向量；N为正整数，默认值为512；返回值包含范围内的N个频率等分点；返回值H则是离散时间系统频率

35、响应在范围内N的频率处对应的值。另外一种形式为H,w=freqz(B,A,N，whole)。与第一种方式的不同之处在于角频率的范围由扩展到。三、实验内容1、使用MATLAB的residuez函数，求出的部分分式展开和。2、使用MATLAB画出下列因果系统的系统函数的零极点图，并判断系统的稳定性。3、使用MATLAB绘制出的频率响应曲线。四、实验报告要求实验名称、实验目的、实验原理、实验环境、实验内容（上述几部分代码及结果）、实验思考等。附录资料：matlab绘图指令大全绘图指令1 二维曲线图1.1 绘制折线图plot指令图例Y=1,3,6,5,9,0,2;plot(Y);X=0: pi/10:

36、 pi*2;Y=sin(X);plot(X,Y);X=0: pi/10: pi*2;Y1=sin(X);Y2=cos(X);Plot(X,Y1,X,Y2);调整坐标范围：axisaxis(0,300,0,2)1.2 绘制自定义函数DrawCircle.mfunction DrawCircle(Point,Radius) Hold on t=0: pi/10: 2*pi; x=Point(1)+ Radius*cos(t); y=Point(2)+ Radius*sin(t); plot(x,y);DrawCircle(10 10,1)DrawCircle(20 10,2)DrawCircle(

37、10 20,3)1.3 绘制符号函数显函数ezplot(sin(x),0,2*pi)隐函数ezplot(x2+y2-10,-5,5,-6,6)参数方程ezplot(cos(t)3,sin(t)3,0,2*pi)1.4 绘制自定义函数function y=myf1(x) y=sqrt(100-x2);fplot(myf1,-15 15)fplot(sin(x) cos(x) myf1(x),-15 15)1.5 图形修饰 设置颜色 y m c r g b w k 设置线型 - : -. - 设置标记 . o x + * 指令图例Y=1,3,6,5,9,0,2;plot(Y, r-+);X=0:

38、pi/10: pi*2;Y=sin(X);plot(X,Y, b-.);X=0: pi/10: pi*2;Y1=sin(X); Y2=cos(X);plot(X,Y1,r+-, X,Y2,b-*); 在指定坐标处，书写文字：text(3.5, 0.6, 曲线比较);x=1.6*pi, 1.6*pi; y=-0.3, 0.8;s=曲线cos; 曲线sin; text(x,y,s);1.6 更多类型的二维图指令图例bar直方图X=0:pi/10:2*pi;Y=sin(X);bar(X,Y);polar极坐标图T=0: pi/10: 4*pi;R=T;polar(T, R);误差棒棒图X=0:pi/

39、10:2*pi;Y=sin(X);e=0.2*rand(size(X);errorbar(X,Y,e);火柴杆图X=0:pi/10:2*pi; Y=sin(X);stem(X,Y);stairs楼梯图X=0:pi/10:2*pi; Y=sin(X);stairs(X,Y);多边形填色图X=1,2,3,4,5; Y=3,5,2,1,6;fill(X,Y,r);hold on; % 保持图形plot(X,Y,o)1.7 数值函数的二维图 可用于绘图，更可用于采样取点。 fplot(0.5*cos(x),-pi,pi) % 绘图X,Y = fplot(0.5*cos(x),-pi,pi); % 返回

40、点坐标fplot(cos(x),-pi,pi,r-+); % 观察点的位置控制采样点的密度fplot(cos(x),-pi,pi,r-+,0.05);fplot(cos(x),-pi,pi,r-+,0.1); 可绘制系统函数，也可绘制自定义函数的图形。2 三维曲线图2.1 三维曲线plot3指令图例X=0: 0.1: 8*pi;Y=sin(X);Z=cos(X); plot3(X,Y,Z,r);X=0: 0.1: 8*pi;Y=sin(X);Z1=cos(X);Z2=2*cos(X); plot3(X,Y,Z1,r, X,Y,Z2,b);2.2 三维面填色fill3指令图例X1=2,2,1;Y

41、1=0,2,1; Z1=0,0,1;fill3(X1,Y1,Z1,r);hold on;X2=1,0,0;Y2=1,2,0;Z2=1,0,0;fill3(X2,Y2,Z2,r);X3=0,2,1;Y3=2,2,1;Z3=0,0,1;fill3(X3,Y3,Z3,b);text(1,1,1,1,1,1);3 曲面图形3.1 网格点坐标的表示x=1:2:7y=2:2:6X,Y = meshgrid(x,y)X = 1 3 5 7 1 3 5 7 1 3 5 7Y = 2 2 2 2 4 4 4 4 6 6 6 63.2 三维网格mesh、meshc、meshz 用途：数据场的观察分析命令图例随机数

42、据的网格Z=rand(5,5);mesh(Z);% 设置颜色colormap(1,0,0);自定义函数的网格x=-4: 1: 4;y=-5: 1: 5;X,Y=meshgrid(x,y);Z=X.2+Y.2;mesh(X,Y,Z); 消影开关：hidden on / hidden off 利用peaks(50)作为模拟数据矩阵;命令图例 带等高线的网格Z=peaks(50);meshc(Z);Z=peaks(50);meshc(Z);colormap(1,0,0);带基准面的网格Z=peaks(50);meshz(Z);剪孔Z=peaks(50);Z(30:45,15:30)=NaN*ones

43、(16,16);meshc(Z);3.3 着色表面图surf、surfc命令图例表面着色的网格Z=peaks(50);surf(Z);自定义函数的着色网格x=-2: 0.1: 2;y=-2: 0.1: 2;X,Y=meshgrid(x,y);Z=sqrt(X.2+Y.2);surfc(X,Y,Z);3.4 二元函数的伪彩色图pcolor 用途：污染浓度场的观察分析。命令图例Z=peaks(50); pcolor(Z);colorbar(hor);colorbar(vec);3.5 等高线contour不仅可用于绘图，更可以用以求截面数据。命令图例以矩阵下标为x、y分量的等高线Z=peaks(5

44、0);C=contour(Z);colormap(1,0,0);C：保存了全部等高线上的点坐标。均分n条等高线，并标注之Z=peaks(50);n=5;C=contour(Z,n);colormap(1,0,0);clabel(C);在指定高度绘制等高线Z=peaks(50);V=-10: 2: 10;C=contour(Z, V);colormap(1,0,0);clabel(C);完整图形数据的等高线x=-2: 0.1: 2;y=-2: 0.1: 2;X,Y=meshgrid(x,y);Z=sqrt(X.2+Y.2);n=5;C=contour(X,Y,Z,n);三维等高线x=-2:0.1

45、:2;y=-2:0.1:2;X,Y=meshgrid(x,y);Z=sqrt(X.2+Y.2);n=10;C=contour3(X,Y,Z,n); 3.6 矢量场图quiver用于挖掘数据变化趋势。命令图例构造起伏跌宕的曲面x=-2:0.2:2;y=-1:0.2:1;X,Y=meshgrid(x,y);Z=X.*exp(-X.2-Y.2);mesh(X,Y,Z);xlabel(X轴);ylabel(Y轴);colormap(1,0,0);计算曲面的梯度px,py= gradient(Z,0.2,0.2);绘制矢量场图quiver(x,y,px,py);3.7 视角控制view视点控制方式及效果

46、：view(1 1 1)view(2 1 1)view(3 1 1)view(1 1 1)view(1 2 1)view(1 3 1)view(1 1 1)view(1 1 2)view(1 1 3)方位角、仰角控制方式及效果：缺省为(-37.5,30)。view(-37.5,30)view(-17.5,30)view(-5.5,30)view(-37.5,30)view(-37.5,45)view(-37.5,60)3.8 多视区控制subplotsubplot(2,1,1); mesh(X,Y,Z);subplot(2,1,2); quiver(x,y,px,py);3.9 制作、播放动画x,y,z=peaks(30); s