当前位置：文档视界 › 实验七连续系统的复频域分析

实验七连续系统的复频域分析

一、实验目的：

1、了解连续系统的复频域分析的基本实现方法；

2、掌握相关函数调用格式及实现功能。

3、掌握用MATLAB 绘制连续系统零极点图与拉普拉斯变换曲面图

4、掌握用MATLAB 实现连续系统的频率特性及其幅度特性、相位特性。

二、实验仪器

计算机 MATLAB 软件

三、实验原理：

1、系统仿真方法:

根据模型的种类不同分为：

（1）数学仿真－用计算机软件摸拟各种实际系统的数学模型（软件仿真）。（2）物理仿真－用实物模型来摸拟各种实际系统（硬件仿真）。（3）数学--物理仿真—两种的结合

系统仿真分析与设计方法：时域法和频域法

（1）时域法：是以状态方程为基础对系统进行分析设计。系统特性分析包括：李亚谱诺夫（Lyapunov ）稳定性分析，能控能观性分析等。

（2）频域法：主要是借助于传递函数，通过劳斯（Routh ）定理、奈氏图（Nyquist ）、伯德图（

Bode ）

、尼克尔斯图（Nichols ）、根轨迹等概念和方法分析

微分方程

状态空间表达式传递函数

结构图

复域模型

时域模型

频率特性图

频域模型

系统的数学模型

系统的各种特性。如稳定性、动态特性、稳态误差等。

2、实验涉及的MATLAB函数

(1) laplace函数

功能：用符号推理求解拉氏变换。

调用格式：L= laplace（F）F为函数，默认为变量t的函数，返回L为s的函数时，要用syms 命令定义符号变量t。

(2) ilaplace函数

功能：用符号推理求解反拉氏变换。

调用格式：L= ilaplace（F）F为函数，默认为变量t的函数，返回L为s 的函数时，要用syms 命令定义符号变量t。

(3) roots（）函数

功能：求多项式的根

调用格式：r=roots（C）：其中C为多项式的系数向量（从高次到低次），r为根向量。因此可用直接求根来判断系统稳定性。

(4) rlocus函数

功能：求系统的根轨迹

调用格式：rlocus（sys）绘制系统的根轨迹；

rlocus（sys，K）绘制增益为K的闭环极点；

rlocus（sys1，sys2，…）在同一复平面绘制多个系统的根轨迹，为区分各个系统的不同根轨迹，可用不同的颜色来显示，如rlocus（sys1，‘r’，sys2，‘y’…）；

(5) step函数

功能：绘制系统的阶跃曲线。

调用格式：step(num，den，t) ，其num，den分别为传递函数的分子、分母的多项式。t为仿真时间。

（6）impulse函数

功能：绘制系统的冲激响应曲线。

调用格式：impulse(num，den，t)

（7）lsim函数

功能：绘制系统的任意响应曲线。

调用格式：lsim(num ， den ， u ，t ）， u 为任意输入信号。（8）gensing （）函数功能：用来产生一些信号。

调用格式：[u ，t]= gensing （type ，tau ，Tf ，Ts ）。Type 为产生信号的类型（’sin ’正弦波、’square ’为方波、’pulse ’为脉冲序列）。Tau 为信号周期，Tf 信号持续时间，Ts 采样周期。u 为所产生的信号。

gensing （）函数和lsim(）联合起来使用。四、实验内容

经典的拉普拉斯变换分析方法即先从时域变换到复频域，在复频域经过处理后再利用拉普拉斯反变换从复频域变换到时域，完成对时间问题的求解，涉及的函数有laplace 函数和ilaplace 函数等

1、求解时域信号对应的拉普拉斯变换、已知象函数求拉普拉斯反变换

例1、求信号e -3t sin2tu(t)拉氏变换 MATLAB 程序： >> clear; >> syms t s

>> f=exp(-3*t)*sin(2*t); >> F=laplace(f) 运行结果： F =

1/2/(1/4*(s+3)^2+1)

例2、已知象函数，求原函数f(t)

MATLAB 程序： clear; syms s t

F=(s-2)/(s^4+ 3*s^3+3* s^2+s) f_t=ilaplace(F) 运行结果：

f_t=-2+3/2*t^2*exp(-t)+2*t*exp(-t)+2*exp(-t)

s s s s s F +++-=2

332)(

2、利用拉普拉斯变换求解系统全响应例

3、已知)(3)(t u t e =，6

()()(2

++=

=s s s s E s R s H ,初始状态y(0)=1 y ’(0)=1;

求系统零状态响应)(t r zs 。

%% 零输入响应 num=[1 0 ]; den=[1 5 6 ]; sys=tf(num,den); t=0:0.01:3; sys1=ss(sys); y=[1 1 ];

u=zeros(1,length(t)); rzi=lsim(sys1,u,t,y);% subplot(311); plot(t,rzi);

title('零输入响应yzi(t)'); ylabel('rzi(t)') %% 零状态响应 syms s f

f=ilaplace(3/((s+2)*(s+3)));t=0:0.01:3; rzs=3*exp(-2*t)-3*exp(-3*t); subplot(312); plot(t,rzs) title('零状态响应'); ylabel('rzs(t)')

%% 全响应=零状态响应+零输入响应 r=rzi+rzs'; subplot(313); plot(t,r);

title('全响应');

ylabel('r(t)');xlabel('时间（秒）');ylabel('r(t)')

例4、已知求解y(t)

MATLAB 程序参考： clear;

b=［1］;a=［1 3 2］; ［A B C D ］=tf2ss(b,a); sys=ss(A,B,C,D); t=0:0.01:10; f=-2*exp(-t); zi=［2 1］; y=lsim(sys,f,t,zi); plot(t,y);grid; xlabel(′时间(t)′); ylabel(′y(t)′); title(′系统的全响应′); 运行结果：

、使用部分分式法求解拉普拉斯反变换

使用部分分式法求解拉普拉斯反变换，方法简单，易于实现，但是对于分母

,2_)0(',1_)0(,0e 2)(2)(d d 3)(d d 22

===+++-y y t y t y t

t y t

多项式复杂难于分解的、或者分解后有多重极点的，手动计算的运算量就很大了，这里介绍使用matlab 方便实现部分分式的分解，进而易于就拉斯反变换。例5、已知象函数求F(s)使用部分分式展开的表示形式

分析：关键是使用部分分式法裂项分开 Matlab 语句：

b=［0 0 0 1 -2］; %分子系数 a=［1 3 3 1 0］; %分母系数

［r,p,k ］=residue(b,a) %留数运算，其中r 是系数， p 是极点， k 是直接项运行实验结果： r =

2.0000 2.0000

3.0000 -2.0000 p =

-1.0000 -1.0000 -1.0000 0 k =

[ ] 因此s

s s s s F 21

1(2)

1(3)(2

例6、已知象函数 ,求原函数f(t) 分析：关键是使用部分分式法裂项分开 Matlab 语句：

b=［1 5 9 7］; %分子系数 a=［0 1 3 2］; %分母系数

［r,p,k ］=residue(b,a) %留数运算，其中r 是系数， p 是极点， k 是直接项运行结果： r = -1 2

237

95)(22

3+++++=s s s s s s F 237

95)(22

+++++=s s s s s s F

p = -2 -1 k =

1 2

因此用部分分式法裂项后的表示形式为： 2

22)(+-+

+++=s s s s F

4、利用MATLAB 绘制系统传递函数H(s)的零极点图（1）传递函数：

用多项式形式表示的开环传递函数，

在Matlab 中函数tf （）可用来建立传递函数模型，其调用格式为：G （S ）=tf （num ，den ）。（2）用多项式形式表示的闭环传递函数：

2111

211

21)

(1)()(+-+-+-++++++++++=

n n n

m m m

a s a s a

b s

b b s

b s

b s G s G s φ

（3）用零极点增益形式表示的开环传递函数

)

())(()())(()(2121n m p s p s p s z s z s z s k

s G ------=

在Matlab 中函数zpk （）可用来建立传递函数零、极、增益模型，其调用格式为：G （S ）=zpk （z ，p ，k ）。例7、已知某系统的传递函数10

111184)

()()(23

+++++=

=s s s s s s E s R s H ，绘制零极点图

MATLAB 参考程序：

num=[1,4,8] %传递函数分子多项式系数行向量 den=[1,11,11,10] %传递函数分母多项式系数行向量 G=tf(num,den) %建立传递函数模型

[z,p,k]=tf2zp(num,den) %提取传递函数的零极点和增益 pzmap(G) %绘制其零极点图 grid on %打开绘图网格

211121)

()()(+-+-++++++=

=n n n

m m m a s

a s a

b s

b s den s num s G

运行结果：零极点图

5、用MATLAB 绘制拉普拉斯变换的曲面图

从三维几何空间的角度来看，)(s F 和)(s ?对应着复平面的两个曲面，如果我们能绘出它们的三维曲面图，我们就可以直观的分析连续信号的拉普拉斯变换F （s ）随复变量s 的变化。

上述过程我们可以利用MATLAB 的三维绘图功能来实现。现在考虑如何用MATLAB 来绘制s 平面的有限区域上连续时间信号f （t ）的拉普拉斯变换F(s)的曲面图，我么以单位阶跃信号u （t ）为例来说明实现过程。

我们知道，对单位阶跃信号f （t ）= u （t ），其拉普拉斯变换为s

s F 1)(=

。

首先，我们用两个向量来确定绘制曲面图的s 平面的横、纵坐标的范围。然后再调用meshgrid()函数来产生矩阵S ，并用矩阵来表示绘制曲面图的附平面区域，对应的MATLAB 命令如下：

[x,y]=meshgrid(x1,y1)

S=x+i*y

上述命令产生的矩阵S 包含了复平面2.02.02.02.0<<-<<-ωδj 、范围以内间隔0.03取样的所有样点。

最后我们在计算出信号拉普拉斯变换在复平面上的这些样点上的值，即可用函数mesh （）来绘出其曲面图。

例8、已知连续时间信号f(t)=sin(t)u(t),求出该信号的拉普拉斯变换，并用MATLAB 绘制拉普拉斯变换的曲面图

参考程序：

a=-0.5:0.08:0.5;

b=-1.99:0.08:1.99;

[a,b]=meshgrid(a,b);

d=ones(size(a));

c=a+i*b;

c=c.*c;

c=c+d;

c=1./c;

c=abs(c);

mesh(a,b,c);

surf(a,b,c);

axis([-0.5,0.5,-2,2,0,15]);

title('单边正弦信号拉斯变换曲面');

colormap(hsv);

执行上述命令后，绘制的拉普拉斯变换曲面图为

五、实验报告要求：

1、认真学习实验原理和实验内容

2、仔细学习实验例题，要求独立上机练习

3、熟练使用MATLAB求解信号的复频域象函数、绘制频率函数零极点图

4、规范化的书写实验报告

大作业1(机电控制系统时域频域分析)

《机电系统控制基础》大作业一基于MATLAB的机电控制系统响应分析哈尔滨工业大学 2013年11月4日

1 作业题目 1. 用MATLAB 绘制系统2 ()25()() 425 C s s R s s s Φ== ++的单位阶跃响应曲线、单位斜坡响应曲线。 2. 用MATLAB 求系统2 ()25 ()()425 C s s R s s s Φ==++的单位阶跃响应性能指标：上升时间、峰值时间、调节时间和超调量。 3. 数控直线运动工作平台位置控制示意图如下： X i 伺服电机原理图如下： L R (1)假定电动机转子轴上的转动惯量为J 1，减速器输出轴上的转动惯量为J 2，减速器减速比为i ，滚珠丝杠的螺距为P ，试计算折算到电机主轴上的总的转动惯量J ； (2)假定工作台质量m ，给定环节的传递函数为K a ，放大环节的传递函数为K b ，包括检测装置在内的反馈环节传递函数为K c ，电动机的反电势常数为K d ，电动机的电磁力矩常数为K m ，试建立该数控直线工作平台的数学模型，画出其控制系统框图； (3)忽略电感L 时，令参数K a =K c =K d =R=J=1，K m =10，P/i =4π，利用MATLAB 分析kb 的取值对于系统的性能的影响。

2 题目1 单位脉冲响应曲线单位阶跃响应曲线

源代码 t=[0:0.01:1.6]; %仿真时间区段和输入 nC=[25]; dR=[1,4,25]; fi=tf(nC,dR); %求系统模型 [y1,T]=impulse(fi,t); [y2,T]=step(fi,t); %系统响应 plot(T,y1); xlabel('t(sec)'),ylabel('x(t)'); grid on; plot(T,y2); xlabel('t(sec)'),ylabel('x(t)'); grid on; %生成图形 3 题目2 借助Matlab，可得： ans = 0.4330 0.6860 25.3826 1.0000 即

电力系统分析实验报告四(理工类)

西华大学实验报告（理工类）开课学院及实验室：实验时间：年月日一、实验目的 1)初步掌握电力系统物理模拟实验的基本方法。 2)加深理解功率极限的概念，在实验中体会各种提高功率极限措施的作用。 3)通过对实验中各种现象的观察，结合所学的理论知识，培养理论结合实际及分析问题的能力。二、实验原理所谓简单电力系统，一般是指发电机通过变压器、输电线路与无限大容量母线联接而且不计各元件的电阻和导纳的输电系统。对于简单系统，如发电机至系统d 轴和g 轴总电抗分别为d X ∑和q X ∑，则发电机的功率特性为当发电机装有励磁调节器时，发电机电势q E 随运行情况而变化，根据一般励磁调节器的性能，可认为保持发电机'q E （或' E ）恒定。这时发电机的功率特性可表示成或这时功率极限为随着电力系统的发展和扩大，电力系统的稳定性问题更加突出，而提高电力系统稳定性和输送能力的最重要手段之一，就是尽可能提高电力系统的功率极限。从简单电力系统功率极限的表达式看，要提高功率极限，可以通过发电机装设性能良好的励磁调节器，以提高发电机电势、增加并联运行线路回路数；或通过串联电容补偿等手段，以减少系统电抗，使受端系统维持较高的运行电压水平；或输电线采用中继同步调相机、中继电力系统等手段以稳定系统中继点电压。 (3)实验内容 1)无调节励磁时，功率特性和功率极隈的测定 ①网络结构变化对系统静态稳定的影响（改变戈）：在相同的运行条件下（即系统电压U-、发电机电势E 。保持不变．罚芳赆裁Ll=E 。），分别测定输电线单回线和双回线运行时，发电机的功一角特性曲线，＆豆甍辜授冁蝮和达到功率极限时的功角值。同时观察并记录系统中其他运行参数（如发电极端毫玉萼蔫交化。将两种情况下的结果加以比较和分析。实验步骤如下： a)输电线路为单回线； b)发电机与系统并列后，调节发电机，使其输出的有功和无ZZ 蔓专零： c)功率角指示器调零； d)逐步增加发电机输出的有功功率，而发电机不调节震磁： e)观察并记录系统中运行参数的变化，填入表1.3中： f)输电线路为双回线，重复上述步骤，将运行参数填入表l 。毒=：

实验五信号与系统的复频域分析

实验五信号与系统的复频域分析王靖 08通信 12号实验目的（1）掌握利用MA TLAB 进行连续时间信号与系统的复频域分析。（2）掌握利用MA TLAB 进行离散系统的复频域分析。实验环境安装MATLAB7.0以上版本的计算机实验内容 1. 利用help 命令了解以下命令的基本用法 residue ，roots ，pzmap ，cart2pol ，residuez ，tf2zp ，zplane 2. 部分分式展开的MATLAB 实现用部分分式展开法求X(s)的反变换。 2321 ()452s X s s s s +=+++ 步骤一：建立新的m 文件，保存并命名为program1.m 。步骤二：输入以下命令，理解每条命令的含义。 %program1，部分分式展开法求反变换 [10 1];[1452];[,,](,) n u m d en r p k resid u e n u m d en === 步骤三：保存程序并运行，记录得到的结果。如右图所示步骤四：由得到的结果可以直接获得X(s)展开表示式 25 4 2 ()21(1)X s s s s =-++++：步骤五：由此可得到X(s)反变换的原函数，记录。 X(t)=(5exp(-2*t)-4exp(-t)+2texp(-t)) 思考：将其转换成极坐标形式，应该如何使用cart2pol 命令？离散系统的部分分式展开，如何使用命令residuez ，得到的结果如何利用？将笛卡尔坐标转化为极坐标用 [angle,mag]=cart2pol(real(r),imag(r)) [r,p,k] = residuez(nun,,den)

(实验三)连续时间LTI系统的频域分析汇总

实验三连续时间LTI 系统的频域分析一、实验目的 1、掌握系统频率响应特性的概念及其物理意义； 2、掌握系统频率响应特性的计算方法和特性曲线的绘制方法，理解具有不同频率响应特性的滤波器对信号的滤波作用； 3、学习和掌握幅度特性、相位特性以及群延时的物理意义； 4、掌握用MA TLAB 语言进行系统频响特性分析的方法。基本要求：掌握LTI 连续和离散时间系统的频域数学模型和频域数学模型的MATLAB 描述方法，深刻理解LTI 系统的频率响应特性的物理意义，理解滤波和滤波器的概念，掌握利用MATLAB 计算和绘制LTI 系统频率响应特性曲线中的编程。二、实验原理及方法 1 连续时间LTI 系统的频率响应所谓频率特性，也称为频率响应特性，简称频率响应（Frequency response ），是指系统在正弦信号激励下的稳态响应随频率变化的情况，包括响应的幅度随频率的变化情况和响应的相位随频率的变化情况两个方面。上图中x(t)、y(t)分别为系统的时域激励信号和响应信号，h(t)是系统的单位冲激响应，它们三者之间的关系为：)(*)()(t h t x t y =，由傅里叶变换的时域卷积定理可得到： )()()(ωωωj H j X j Y = 3.1 或者： ) () ()(ωωωj X j Y j H = 3.2 )(ωj H 为系统的频域数学模型，它实际上就是系统的单位冲激响应h(t)的傅里叶变换。即 ? ∞ ∞ --= dt e t h j H t j ωω)()( 3.3 由于H(j ω)实际上是系统单位冲激响应h(t)的傅里叶变换，如果h(t)是收敛的，或者说是绝对可积（Absolutly integrabel ）的话，那么H(j ω)一定存在，而且H(j ω)通常是复数，

3-系统分析实验报告

管理信息系统实验报告实验3 系统分析课程名称：管理信息系统指导教师：王玮班级：信管1401 学号：姓名：唐赛赛时间： 2016.04.06 地点： 3 号机房

一、实验目的 1．了解开发Visio解决方案的基本概念和关于Visio工具的一些基本的操作和应用； 2．掌握系统分析阶段数据流程图的画法；二、实验步骤和实验结果：使用Visio中提供的“组织结构图”模具，绘制下面例题的组织结构图，附在图后。 2、使用Visio绘制“业务流程图模具”和“数据流程图模具”（1）创建“业务流程图模具” 先在“框图”－〉“基本形状”中找到圆角矩形，右击选择“添加到我的形状”－〉“添加到新模具”。之后出现“另存为”对话框，把新模具命名为“业务流程图”，把圆角矩形形添加到了新模具“业务程图”中。用同样的思路，先在“框图”－〉“基本形状”中找到圆形，右击选择“添加到我的形状”－〉“添加到模具“业务程图”中；在“框图”－〉“基本形状”找到矩形，在“流程图”中的“IDEFO图表形状”找到动态连接线，在“流程图”中的“SDL图表形状”中找到文档，多文档，添加到模具“业务程图”中。可以通过设置“动态连接线”属性来改变其形状。如下图：

添加完成后，我们就可以在画业务流程图时打开该模具，业务流程图所有的元素都会在一个模具中显示出来。（2）创建“数据流程图模具”先在“框图”－〉“基本形状”中找到圆形（或是“流程图”中的“混合流程图形状”中找到外部实体2 ），右击选择“添加到我的形状”－〉“添加到新模具”（注，使用外部实体2来表示外部实体的时候，请将之旋转180度使用）。之后出现“另存为”对话框，把新模具命名为“数据流程图”，这样我们就把圆形形添加

北京理工大学信号与系统实验报告5-连续时间系统的复频域分析

实验5连续时间系统的复频域分析（综合型实验）一、实验目的 1）掌握拉普拉斯变换及其反变换的定义并掌握MATLAB 实现方法。 2）学习和掌握连续时间系统函数的定义及复频域分析方法。 3）掌握系统零极点的定义，加深理解系统零极点分布与系统特性的关系。二、实验原理与方法 1.拉普拉斯变换连续时间信号x(t)的拉普拉斯变换定义为 (s)(t)e st X x dt +∞ --∞ = ? （1）拉普拉斯反变换为1 (t)(s)e 2j st j x X ds j σσπ+∞ -∞ =? （2） MATLAB 中相应函数如下： (F) L laplace = 符号表达式F 拉氏变换，F 中时间变量为t ，返回变量为s 的结果表达式。 (F,t)L laplace =用t 替换结果中的变量s 。 () F ilaplace L =以s 为变量的符号表达式L 的拉氏反变换，返回时间变量为t 的结果表达式。 (,) F ilaplace L x =用x 替换结果中的变量t 。

的连续时间系统，其系统函数为s 的有理函数 110 110 ...(s)...M M M M N N N N b s b s b H a s a s a ----+++= +++ （7） 3.连续时间系统的零极点分析系统的零点指使式（7）的分子多项式为零的点，极点指使分母多项式为零的点，零点使系统的值为零，极点使系统的值为无穷大。通常将系统函数的零极点绘在s 平面上，零点用O 表示，极点用?表示，这样得到的图形为零极点分布图。可以通过利用MATLAB 中的求多项式根的roots 函数来实现对（7）分子分母根的求解，调用格式如下： r=roots(c)，c 为多项式的系数向量，返回值r 为多项式的根向量。求取零极点以及绘制系统函数的零极点分布图可以采用pzmap 函数，调用格式如下： pzmap(sys)绘出由系统模型sys 描述的系统的零极点分布图。 [p,z]=pzmap(sys)这种调用方式返回极点与零点，不绘出零极点分布图。还有两个专用函数tf2zp 和zp2tf 可实现系统的传递函数模型和零极点增益模型的转换。调用格

连续时间LTI系统的频率特性及频域分析

实验报告实验项目名称：运用Matlab进行连续时间信号卷积运算（所属课程：信号与系统）学院：电子信息与电气工程学院专业： 10电气工程及其自动化姓名： xx 学号： 201002040077 指导老师： xxx

一、实验目的 1、学会运用MATLAB 分析连续系统的频率特性。 2、掌握相关函数的调用。二、实验原理 1、一个连续LTI 系统的数学模型通常用常系数线性微分方程描述，即 )()()()()()(01 )(01)(t e b t e b t e b t r a t r a t r a m m n n +'++=+'++ (1) 对上式两边取傅里叶变换，并根据FT 的时域微分性质可得： )(])([)(])([0101ωωωωωωE b j b j b R a j a j a m m n n +++=+++ 101)()()()()(a j a j a b j b j b j E j R j H n n m m ++++++==ωωωωωωω H ( j ω )称为系统的频率响应特性，简称系统频率响应或频率特性。一般H ( j ω )是复函数，可表示为： )()()(ω?ωωj e j H j H = 其中， )(ωj H 称为系统的幅频响应特性，简称为幅频响应或幅频特性；)(ω?称为系统的相频响应特性，简称相频响应或相频特性。H ( j ω )描述了系统响应的傅里叶变换与激励的傅里叶变换间的关系。H ( j ω )只与系统本身的特性有关，与激励无关，因此它是表征系统特性的一个重要参数。 MATLAB 信号处理工具箱提供的freqs 函数可直接计算系统的频率响应的数值解，其语句格式为：H=freqs(b,a,w)其中，b 和a 表示H ( j ω )的分子和分母多项式的系数向量；w 为系统频率响应的频率范围，其一般形式为w1:p:w2，w1 为频率起始值，w2 为频率终止值，p 为频率取值间隔。 H 返回w 所定义的频率点上系统频率响应的样值。注意，H 返回的样值可能为包含实部和虚部的复数。因此，如果想得到系统的幅频特性和相频特性，还需要利用abs 和angle 函数来分别求得。

控制系统的频域分析实验报告

实验名称：控制系统的频域分析实验类型：________________同组学生姓名：__________ 一、实验目的和要求用计算机辅助分析的方法，掌握频率分析法的三种方法，即Bode 图、Nyquist 曲线、Nichols 图。二、实验内容和原理（一）实验原理 1．Bode(波特)图设已知系统的传递函数模型： 1 1211121)(+-+-+???+++???++=n n n m m m a s a s a b s b s b s H 则系统的频率响应可直接求出： 1 1211121)()()()()(+-+-+???+++???++=n n n m m m a j a j a b j b j b j H ωωωωω MATLAB 中，可利用bode 和dbode 绘制连续和离散系统的Bode 图。 2．Nyquist(奈奎斯特)曲线 Nyquist 曲线是根据开环频率特性在复平面上绘制幅相轨迹，根据开环的Nyquist 线，可判断闭环系统的稳定性。反馈控制系统稳定的充要条件是，Nyquist 曲线按逆时针包围临界点(-1，j0)p 圈，为开环传递函数位于右半s 一平面的极点数。在MATLAB 中，可利用函数nyquist 和dnyquist 绘出连续和离散系统的乃氏曲线。 3．Nicho1s(尼柯尔斯)图根据闭环频率特性的幅值和相位可作出Nichols 图，从而可直接得到闭环系统的频率特性。在 MATLAB 中，可利用函数nichols 和dnichols 绘出连续和离散系统的Nichols 图。（二）实验内容 1．一系统开环传递函数为 ) 2)(5)(1(50)(-++=s s s s H 绘制系统的bode 图，判断闭环系统的稳定性，并画出闭环系统的单位冲击响应。 2．一多环系统 ) 10625.0)(125.0)(185.0(7.16)(+++=s s s s s G 其结构如图所示试绘制Nyquist 频率曲线和Nichols 图，并判断稳定性。（三）实验要求

北京理工大学信号与系统实验实验5连续时间系统地复频域分析报告报告材料

实验5 连续时间系统的复频域分析一、实验目的 1.掌握拉普拉斯变换及其反变换的定义，并掌握MATLAB 实现方法。 2.学习和掌握连续时间系统系统函数的定义及复频域分析方法。 3.掌握系统零极点的定义，加深理解系统零极点分布与系统特性的关系。二、实验原理与方法 1.拉普拉斯变换连续时间信号)(t x 的拉普拉斯变换定义为 )1.....(..........)()(dt e t x s X st ? +∞ ∞ --= 拉普拉斯反变换定义为 )2....(..........)(21)(ds e s X j t x j j st ?∞ +∞ -=σσπ 在MATLAB 中，可以采用符号数学工具箱的laplace 函数和ilaplace 函数进行拉氏变换和反拉氏变换。 L=laplace(F)符号表达式F 的拉氏变换，F 中时间变量为t ，返回变量为s 的结果表达式。 L=laplace(F,t)用t 替换结果中的变量s 。 F=ilaplace(L)以s 为变量的符号表达式L 的拉氏反变换，返回时间变量为t 的结果表达式。 F=ilaplace(L,x)用x 替换结果中的变量t 。除了上述ilaplace 函数，还可以采用部分分式法，求解拉普拉斯逆变换，具体原理如下：当 X (s )为有理分式时，它可以表示为两个多项式之比： )3.(..........)()()(0 110 11a s a s a b s b s b s D s N s X N N N N M M M M +?+++?++==---- 式（3）可以用部分分式法展成一下形式 )4.....(.............)(2211N N p s r p s r p s r s X -++-+-= 通过查常用拉普拉斯变换对，可以由式（1-2）求得拉普拉斯逆变换。利用 MATLAB 的residue 函数可以将 X (s )展成式（1-2）所示的部分分式展开式，该函数的调用格式为：[r,p,k] = residue(b,a) 其中b 、a 为分子和分母多项式系数向量，r 、p 、k 分别为上述展开式中的部分分式系数、极点和直项多项式系数。 2.连续时间系统的系统函数

系统分析实验报告

天津职业技术师范大学课程设计大学学籍管理系统的设计与开发专业：软件工程班级学号：软件1002-17 学生姓名：靳利强指导教师：龚良波老师二〇一三年七月

一．需求分析 1.课程名称：大学教务信息系统的设计与开发 2.设计目的：为方便学校做好学生学籍管理工作，设计一个学生学籍管理系统，在设计过程中作了系统分析和总体设计，软件设计采取模块化的设计思路。 3.需求概述该学生学籍管理系统主要对学生学籍信息、成绩信息进行管理，提供一个平台，供学籍管理人员增删改查学生信息、学生成绩信息。系统分为学生信息管理、学生成绩管理、信息查询等几个模块。学籍管理人员登录成功后可以对学生信息管理、学生成绩管理、信息查询等模块进行操作，如学生信息添加、修改、删除和查询；学生成绩登记、修改、删除和查询；查询信息等。 4功能需求： 1)功能齐全：界面操作灵活方便，设计包括以下基本功能： 2)学生信息管理、教师信息管理、财务信息管理、班级信息管理、课程信息管理、成绩信息管理、打印信息管理、教室信息管理、综合信息查询、系统管理等，至少实现其中的三个功能，且每个功能至少包括两个子功能。 3)按照软件工程的要求进行分析、设计和开发。 4)界面友好：界面友好、输入有提示、尽量展示人性化。 5)可读性强：源程序代码清晰、有层次、主要程序段有注释。

6)健壮性好：用户输入非法数据时，系统应及时给出警告信息。二.概要设计 1.功能模块： 2数据流图：（1）学生端

（2）管理员端

学生端功能： A 登录,学生登录后，验证成功，进入其信息展示页。管理员端功能： B 登录，管理员登录后，验证成功，进入学生信息列表，可以对学生信息进行修改，删除，按班级查询，按学号查询，按名字查询。上传图片，更新图片等操作。三．详细设计及实现数据库设计：学生表：教师表：

实验六-信号与系统复频域分析

实验六信号与系统复频域分析一、实验目的 1.学会用MATLAB进行部分分式展开； 2.学会用MATLAB分析LTI系统的特性； 3.学会用MATLAB进行Laplace正、反变换。 4.学会用MATLAB画离散系统零极点图； 5.学会用MATLAB分析离散系统的频率特性；二、实验原理及内容 1．用MATLAB进行部分分式展开用MATLAB函数residue可以得到复杂有理分式F(s)的部分分式展开式，其调用格式为其中，num,den分别为F(s)的分子和分母多项式的系数向量，r为部分分式的系数，p为极点，k为F(s)中整式部分的系数，若F(s)为有理真分式，则k为零。例6-1 用部分分式展开法求F(s)的反变换解:其MATLAB程序为 format rat; num=[1,2]; den=[1,4,3,0]; [r,p]=residue(num,den) 程序中format rat是将结果数据以分数形式显示

F(s)可展开为 210.536()13 F s s s s --=++++ 所以，F(s)的反变换为 3211()()326t t f t e e u t --??=--???? 2．用MATLAB 分析LTI 系统的特性系统函数H （s ）通常是一个有理分式，其分子和分母均为多项式。计算H （s ）的零极点可以应用MATLAB 中的roots 函数，求出分子和分母多项式的根，然后用plot 命令画图。在MATLAB 中还有一种更简便的方法画系统函数H （s ）的零极点分布图，即用pzmap 函数画图。其调用格式为 pzmap(sys) sys 表示LTI 系统的模型，要借助tf 函数获得，其调用格式为 sys=tf(b,a) 式中，b 和a 分别为系统函数H （s ）的分子和分母多项式的系数向量。如果已知系统函数H （s ），求系统的单位冲激响应h(t)和频率响应H ω（j ）可以用以前介绍过的impulse 和freqs 函数。例6-2 已知系统函数为 321221 s s s +++H(s)= 试画出其零极点分布图，求系统的单位冲激响应h(t)和频率响应H ω（j ），并判断系统是否稳定。解：其MATLAB 程序如下： num=[1];

实验八系统的复频域分析

一、实验目的 1、掌握系统的复频域分析方法。 2、掌握测试系统的频率响应的方法。二、预习内容 1、系统频响的方法。（见第四章波特图的介绍）三、实验原理 1. N 阶系统系统的传递函数用微分方程描述的N 阶系统为：根据零状态响应（起始状态为零），则对其进行拉氏变换有：则系统传递函数可表达为：用差分方程描述的N 阶系统为：根据零状态响应（起始状态为零），则对其进行拉氏变换有：则系统传递函数可表达为： 2.根据系统传递函数的零极点图分析系统零点：传递函数分子多项式的根。极点：传递函数分母多项式的根。根据零极点图的不同分布分析系统。 3.涉及到的Matlab 函数 (1)freqz 函数：实验六中出现过，可用来求单位圆上的有理z 变换的值。调用格式：同实验六 (2)zplane 函数：得到有理z 变换的零极点图。调用格式：zplane(num,den)

其中，num和 den是按z ?1 的升幂排列的、z 变换分子分母多项式系数的行向量。 (3)roots 函数：求多项式的根。调用格式：r=roots(c), c 为多项式系数向量；r 为根向量。四、实验内容 1.系统零极点的求解 (1)求解系统和的零极点，验证下面程序的运行结果，根据系统零极点图分析系统性质。 b=[1,0,-1]; a=[1,2,3,2]; zr=roots(b); pr=roots(a); plot(real(zr),imag(zr),'go',real(pr),imag(pr),'mx','markersize',12,'linewidth',2); grid; legend('零点','极点'); figure; zplane(b,a); (2)参考上述程序，绘制系统和的零极点图，并分析系统性质。与用zplane 函数直接绘制系统零极点图（注:圆心的圆圈并非系统的零点）做比较。

系统分析实验报告2016

本科实验报告课程名称：系统分析与设计实验项目：《》实验实验地点：专业班级：学号：学生姓名：指导教师： 2016年11月日

一、实验目的通过《系统分析与设计》实验，使学生在实际的案例中完成系统分析与系统设计中的主要步骤，并熟悉信息系统开发的有关应用软件，加深对信息系统分析与设计课程基础理论、基本知识的理解，提高分析和解决实际问题的能力，使学生在实践中熟悉信息系统分析与设计的规范，为后继的学习打下良好的基础。二、实验要求学生以个人为单位完成，自选题目，班内题目不重复，使用UML进行系统分析与设计，并完成实验报告。实验报告（A4纸+电子版)在最后一次上课时提交（10周）。三、实验主要设备：台式或笔记本计算机四、实验内容 1 选题及项目背景学生填写自选题目 2 定义学生填写（对自选项目系统进行描述200-400字） 3 参考资料学生填写 4 系统分析与设计 4.1需求分析 4.1.1识别参与者学生填写 4.1.2 对需求进行捕获与描述学生填写时删除以下括号内容（内容要求1：对每个用例进行概要说明，参考以下格式：用例名称：删除借阅者信息执行者：管理员目的：完成一次删除借阅者信息的完整过程。）（内容要求2：选择其中一个用例（如下订单）给出其用例描述。格式参考下表

） 4.1.3 用例图通过已掌握的需求，初步了解系统所要完成的功能。下面给出用例图。 4.1.4 分析与讨论 1)建模用例图的步骤、方法？ 2)如何识别系统的参与者？应该如何划分用例，应注意哪些问题？ 3)心得 4.2 建立对象模型 4.2.1 候选类的数据字典学生填写 4.2.2定义类（内容以“书籍信息”类为例列出该类的属性和操作如下： “书籍信息”类 ?属性国际标准书号(ISBN)：文本(String) 书名(name)：文本

理工大学信号与系统实验报告连续时间系统的复频域分析

理工大学信号与系统实验报告连续时间系统的复频域分析 Document serial number【UU89WT-UU98YT-UU8CB-UUUT-UUT108】

实验5连续时间系统的复频域分析（综合型实验）一、实验目的 1）掌握拉普拉斯变换及其反变换的定义并掌握MATLAB 实现方法。 2）学习和掌握连续时间系统函数的定义及复频域分析方法。 3）掌握系统零极点的定义，加深理解系统零极点分布与系统特性的关系。二、实验原理与方法 1.拉普拉斯变换连续时间信号x(t)的拉普拉斯变换定义为(s)(t)e st X x dt +∞ --∞ =? （1）拉普拉斯反变换为1 (t)(s)e 2j st j x X ds j σσπ+∞ - ∞ = ? （2） MATLAB 中相应函数如下： (F)L laplace = 符号表达式F 拉氏变换，F 中时间变量为t ，返回变量为s 的结果表达式。 (F,t)L laplace =用t 替换结果中的变量s 。 ()F ilaplace L =以s 为变量的符号表达式L 的拉氏反变换，返回时间变量为t 的结果表达式。 (,)F ilaplace L x =用x 替换结果中的变量t 。拉氏变换还可采用部分分式法，当(s)X 为有理分式时，它可以表示为两个多项式之比： 110 1 10 ...(s)(s)(s)...M M M M N N N N b s b s b N X D a s a s a ----+++==+++ （3）

上式可以采用部分分式法展成以下形式 1212(s)...N N r r r X s p s p s p = +++--- （4）再通过查找常用拉氏变换对易得反变换。利用residue 函数可将X(s)展成（4）式形式，调用格式为： [r,p,k]residue(b,a)=其中b 、a 为分子和分母多项式系数向量，r 、p 、k 分别为上述展开式中的部分分式系数、极点和直项多项式系数。 2.连续时间系统的系统函数连续时间系统的系统函数是指系统单位冲激响应的拉氏变换 (s)(t)e st H h dt +∞ --∞ = ? （5）连续时间系统的系统函数还可以由系统输入与输出信号的拉氏变换之比得到。 (s)(s)/X(s)H Y = （6）单位冲激响应(t)h 反映了系统的固有性质，而(s)H 从复频域反映了系统的固有性质。由（6）描述的连续时间系统，其系统函数为s 的有理函数 110 1 10 ...(s)...M M M M N N N N b s b s b H a s a s a ----+++=+++ （7） 3.连续时间系统的零极点分析系统的零点指使式（7）的分子多项式为零的点，极点指使分母多项式为零的点，零点使系统的值为零，极点使系统的值为无穷大。通常将系统函数的零极点绘在s 平面上，零点用O 表示，极点用?表示，这样得到的图形为零极点分布图。可以通过利用MATLAB 中的求多项式根的roots 函数来实现对（7）分子分母根的求解，调用格式如下：

连续系统的复频域分析

实验四：连续系统的复频域分析一、实验目的： 1、掌握连续与离散时间系统的正反复频域与Z域变换 2、掌握利用MATLAB进行零极点分析，进一步了解零极点对整个系统的影响 3、掌握simulink环境下系统建模与仿真以及系统求解。二、实验内容： 1、已知某连续系统的系统函数为：（1）利用[r, p, k]=residue(num, den),求H（s）的极零点以及多项式系数; （2）画出系统的零极点分布图，判断系统得稳定性。（3）求h(t)，判断系统得稳定性。 2、已知某离散系统的系统函数为：，（1）利用[r, p, k]=residuez(num, den)求H（z）的极零点以及多项式系数; （2）画出零极点分布图，判断系统得稳定性。（3）求单位函数响应用impz(b, a)，判断系统是否稳定； 3、已知线性时不变微分方程在Simulink环境下搭建起系统的仿真模型，并查看仿真结果曲线。（1）写出传递函数H（s），绘出系统模拟框图；（2）当f(t)分别为，，的零状态响应；且当与课本P81的结果进行比较（3）方程的初值为, ,求全响应； 4、已知某信号，n(t)为正态噪声干扰且服从N(0,0.22)分布，对此信号进行采样，采样间隔为0.001s，之后对此信号进行Botterworth低通滤波，从信号中过滤10HZ的输出信号，试对系统进行建模与仿真。三、实验数据处理与结果分析：第一题：题1_1：

>> num=[2,5]; den=[1,1,3,2]; [r,p,k]=residue(num,den) r = -0.5750 - 0.7979i -0.5750 + 0.7979i 1.1499 p =-0.1424 + 1.6661i -0.1424 - 1.6661i -0.7152 k =[]

管理信息系统分析实验报告

《管理信息系统》实验二题目：系统分析专业：信息管理与信息系统班级：1106班姓名 ************************* 指导教师：贺玉珍老师完成日期：2014.4.28

运城学院超市管理系统设计分析说明书一、系统目标：随着小超市规模的发展不断扩大，商品数量急剧增加，有关商品的各种信息量也成倍增长。超市时时刻刻都需要对商品各种信息进行统计分析。而大型的超市管理系统功能过于强大而造成操作繁琐降低了小超市的工作效率。超市管理系统是市场上最流行的超市上常用的系统之一，它主要包含以下几个模块：系统权限的设定、原始数据录入、数据的汇总及查询等。从而，实现对进货、销售及员工信息等实现全面、动态、及时的管理。本文系统的分析了软件开发的背景以过程；首先介绍了软件的开发环境，其次介绍了本软件的详细设计过程：数据库的设计、各个模块的设计和实现，以及具体界面的设计和功能。二、系统的初步调查通过实地参观和学习，对超市的整体情况进行调研。了解超市的组织机构划分，充分了解超市进销存的流程的整体情况，对开发新系统的态度等。通过召开座谈会和个人访谈方法了解各个部门的主要职能及具体运作方式、过程等。进行初步调研的具体内容为：（1）员工的规模：大约有多少员工，有多少是稳定的，有多少是浮动的；（2）员工管理人员的数量；（3）超市的商品销售状况（4）客户编码方式；三、可行性分析： 1.技术可行性研究，在IT行业中从业的工作人员一般都要求掌握计算机技术，具有一定的软硬件基础，会使用各种管理软件，熟悉IT产品。因为，有的超市对员工的素质要求比较高，从管理层到下面的销售人员，都要求具有一定的计算机基础，所以在新系统投入使用时，只要对员工进行少量的培训，系统的功能和使用方法就基本上能够是系统顺利运行。 2经济可行性研究，因为通过网络传递销售信息可以不受距离的限制，因此可以借阅许多的人力和物力，方便管理，由此可以减少不必要的开支，同时该系统可以提高超市的销售效率，即提高了超市的经济效益，所以从经济上完全是可行的，(1)超市有能力承担系统开发费用，(2)新系统将为企业带来经济效益3操作可行性研究，本系统采用基于Windows的图形用户界面，而该系统是大家熟悉的操作系统，对于那些有一般的计算机知识的人员就可以轻松上手。而整个超市管理系统采用最友好的交互界面，简介明了，不需要对数据库进行深入的

连续系统的时域、频域分析

学生实验报告实验课程:信号与系统E D A 实验地点:东1教 414 学院: 专业: 学号 : 姓名 :

2．信号卷积,根据PPT 中的实验2、2与2、3内容完成课堂练习,写出程序及运行结果。用Matlab 实现卷积运算)(*)(t h t f ,其中 )()()],2()([2)(t e t h t t t f t εεε-=--=,)2 ()(2t h t h =;对比说明信号)( t f 分别输入系统)(和)(2t h t h 时的输出有什么区别并分析原因。 >> p=0、01; nf=0:p:4; f=2*(heaviside(nf)-heaviside(nf-2)); nh=0:p:6; h=exp(-nh)、*(nh>0); y=conv(f,h);

t=0:length(y)-1; subplot(3,1,1),stairs(nf,f);title('f(t)');axis([0 6 0 2、1]); subplot(3,1,2),plot(nh,h);title('h(t)');axis([0 6 0 1、1]); subplot(3,1,3),plot(0、01*t,y); title('y(t)=f(t)*h(t)'); >> p=0、01; nf=0:p:4; f=2*(heaviside(nf)-heaviside(nf-2)); nh=0:p:6; h=exp(-2*nh)、*(2*nh>0); y=conv(f,h); t=0:length(y)-1; subplot(3,1,1),stairs(nf,f);title('f(t)');axis([0 6 0 2、1]);

管理信息系统实验报告分析

实验报告课程：管理信息系统一、实验目的验证有关概念和理论，加深对概念和知识的理解和认识；熟悉和掌握Visual Basic 6.0 软件的使用方法；初步具备信息管理知识和制作数据字典、系统数据流程图的能力。运用课程讲授的管理信息系统的系统分析方法、模块化系统设计方法以及系统的调试方法进行人事档案管理信息系统的分析、设计、开发、实现与调试。二、实验方法面向对象法三、实验环境及开发工具 1.硬件环境在最低配置的情况下，系统的性能往往不尽如人意，但现在的硬件性能已经相当的出色，而且价格便宜，因此通常给服务器的配置高性能的硬件。处理器：Interl Pentium II 266 MX 或更高内存：64M 硬盘空间：2 GB 显卡：SVGA 显示适配器显示器：液晶17寸 2.软件环境操作系统：Windows/98/ME/2000/XP或更高版本数据库：Microsoft Access 2000 3.实验开发工具：Visual Bisic 6.0程序系统四、实验内容

（一）、系统分析 1、系统数据流程图 2、数据字典 3、系统中所有实体（包括实体的属性）以及实体之间的联系类型分析人员的个人资料经过专业的处理部门的处理形成个人档案。档案包括自然情况，工作情况，简历，政治情况等各方面信息，内容比较庞大复杂。将档案信息传送到人员信息库。同时还综合考虑档案管理工作的性质，总结归纳出所需实现

的功能。为人事档案进行服务，对人事的变动、人事资料、以及人事资料的查询，统计等功能。总体上说具有编辑，查询，用户管理，图表统计等功能。然后将最终结果提交到人力资源管理部门，由人力资源管理人员进行审查，以便于对职工的调配。 4、典型处理的表达档案完整添加用户档案到档案库个人信息成功添加到档案库修改用户档案信息失败退回用户档案退回用户档案（二）、系统设计 1、子系统划分（或功能划分或模块划分）功能划分 1、用户管理功能：设置使用人事管理系统的用户及其使用权限。整个人事管理系统由多个功能模块组成，不同的模块完成不同的功能，所以可以为不同的职工分配不同的功能，使其具有不同的权限，完成其权限所对应的功能，从而很好地管理好整个系统。 2、辅助表管理功能：通过它的这个功能可以有效的对本单位人事部门的扩充进行及时的计算机管理。只要管理员进行简单的数据字段添加即可。辅助表管理功能是高级管理员及中级管理员拥有的权限，它的功能是对数据库进行新表的添加。 3、档案编辑功能：档案编辑模块中有4个子模块。他们是档案卡片、个人简历、家庭成员、历史档案等功能。这些功能因管理员的权限不同所表示出的功能使用也不同，普通管理员没有数据修改及删除的权利。在这些功能里详细的记录了所有单位员工的资料。 4、档案查询功能：对档案卡片的查询功能，在这里可以查到符合程序要求的任何信息。

第5章_用MATLAB进行控制系统频域分析

第5章用MATLAB 进行控制系统频域分析一、基于MATLAB 的线性系统的频域分析基本知识（１）频率特性函数)(ωj G 。设线性系统传递函数为： n n n n m m m m a s a s a s a b s b s b s b s G ++???++++???++=---1101110)( 则频率特性函数为： n n n n m m m m a j a j a j a b j b j b j b jw G ++???++++???++=---)()()()()()()(1101110ωωωωωω 由下面的MATLAB 语句可直接求出G(jw)。 i=sqrt(-1) % 求取-1的平方根 GW=polyval(num ，i*w)./polyval(den ，i*w) 其中（num ，den ）为系统的传递函数模型。而w 为频率点构成的向量，点右除（./）运算符表示操作元素点对点的运算。从数值运算的角度来看，上述算法在系统的极点附近精度不会很理想，甚至出现无穷大值，运算结果是一系列复数返回到变量GW 中。（２）用MATLAB 作奈魁斯特图。控制系统工具箱中提供了一个MATLAB 函数nyquist( )，该函数可以用来直接求解Nyquist 阵列或绘制奈氏图。当命令中不包含左端返回变量时，nyquist （）函数仅在屏幕上产生奈氏图，命令调用格式为： nyquist(num,den) nyquist(num,den,w) 或者 nyquist(G) nyquist(G,w) 该命令将画出下列开环系统传递函数的奈氏曲线： ) () ()(s den s num s G = 如果用户给出频率向量w,则w 包含了要分析的以弧度/秒表示的诸频率点。在这些频率点上，将对系统的频率响应进行计算，若没有指定的w 向量，则该函数自动选择频率向量进行计算。 w 包含了用户要分析的以弧度/秒表示的诸频率点,MATLAB 会自动计算这些点的频率响应。当命令中包含了左端的返回变量时，即: [re,im,w]=nyquist(G) 或

信号与系统报告实验5 连续系统的复频域分析实验

信号与系统实验报告实验五连续系统的复频域分析实验五连续系统的复频域分析一、实验目的 1. 深刻理解拉普拉斯变换、逆变换的定义，掌握用MATLAB实现拉普拉斯变换、逆变换的方法。 2会求几种基本信号的拉氏变换。 3 掌握用MATLAB绘制连续系统零、极点的方法。 4 求解系统函数H(s)。二

1已知连续时间信号f（t）=sin（t）u（t）、求出该信号的拉普拉斯变换，并用MATLAB 绘制拉普拉斯变换的曲面图。 syms t; ft=sin(t)*heaviside(t); Fs=Laplace(ft); a=-0.5:0.08:0.5; b=-2:0.08:2; [a,b]=meshgrid(a,b); c=a+i*b; d=ones(size(a)); c=c.*c; c=c+d; c=1./c; c=abs(c); mesh(a,b,c); surf(a,b,c) axis([-0.5,0.5,-2,2,0,10]) colormap(hsv

) 2求[(1-e^(-at))]/t的拉氏变换。 syms t s a f1=(1-exp(-a*t))/t; F=laplace(f1,t,s) F = log(s+a)-log(s) 3求F（s）=-log（s）+ log（s+a）的拉氏逆变换syms t s a F =log(s+a)-log(s); f1=ilaplace(F,s,t) f1 = (1-exp(-a*t))/t

4已知某连续系统的系统函数为： H（s）=(s^2+3s+2)/(8s^4+2s^3+3s^2+5)试用MATLAB求出该系统的零极点，画出零极点分布图。 b=[1 3 2]; a=[8 2 3 0 5]; zs=roots(b); ps=roots(a); hold on plot(real(zs),imag(zs),'o'); plot(real(ps),imag(ps),'x'); grid axis([-2.5,1,-1,1]) 5已知H(s)=（s+1）/（s^2+s+1）,绘制阶跃响应图形，冲激响应图形，频率激响应图形。 syms t s H=(s+1)/(s^2+s+1); f1=ilaplace(H,s,t); f2=heaviside(t);

实验七 连续系统的复频域分析