自动控制原理实验三用MATLAB实现线性系统的频域分析

自动控制原理实验三用MATLAB实现线性系统的频域分析
自动控制原理实验三用MATLAB实现线性系统的频域分析

实验三 用MATLAB 实现线性系统的频域分析

实验目的:

1.掌握MATLAB 平台下绘制典型环节及系统开环传递函数的Bode 图和Nyquis 图(极坐标图)绘制方法;

2.掌握利用Bode 图和Nyquis 图对系统性能进行分析的理论和方法。 实验要求:

1.根据内容要求,写出调试好的MATLAB 语言程序,及对应的结果。 2. 记录显示的图形,根据实验结果与各典型环节的频率曲线对比分析。 3. 记录并分析ζ对二阶系统bode 图的影响。

4.根据频域分析方法分析系统,说明频域法分析系统的优点。 5.写出实验的心得与体会。 实验内容:

1.典型二阶系统

2

2

22)(n

n n

s s s G ωζωω++= 绘制出6=n ω,1.0=ζ,0.3,0.5,0.8,2的bode 图,记录并分析ζ对系统bode 图的影响。 Bode 图程序:

wn=6;znb=[0.1 0.3 0.5 0.8 2];

w=logspace(0,2,10000);figure(1);n=[wn^2];

for k=znb d=[1 2*k*wn wn^2];sys=tf(n,d);bode(sys,w);hold on; end

运行结果:

结果分析:

从图中可看出ζ越小,中频段振荡越剧烈。该二阶系统是典型的振荡环节,

谐振频率 )2

20(2122

2≤<*-*=ζζωωn r ,谐振峰值

)220(121222≤<-**=ζζζr M 当22

02

<<ζ时, r ω,r M 均为ζ的减函数,ζ越小,r M ,r ω越大,振荡幅度越大,超调量越大,过程越不平稳且系统响应速

度越慢,当

12

2

2

<<ζ时。)(ωA 单调减小,此时无谐振峰值和谐振频率,过程较平稳。

2.系统的开环传递函数为

)5)(15(10

)(2+-=

s s s s G

)106)(15()

1(8)(22

++++=s s s s s s G

)11.0)(105.0)(102.0()

13/(4)(++++=

s s s s s s G

绘制系统的Nyquist 曲线、Bode 图,说明系统的稳定性。

)5)(15(10

)(2+-=

s s s s G

Bode 图: 程序代码:

num=[10];

den=[5 24 -5 0 0]; w=logspace(-2,3,100); bode(num,den,w) grid 实验结果:

Nyquist 图: 程序代码:

num=[10];

den=[5 24 -5 0 0]; nyquist(num,den) 实验结果:

结果分析:

因为开环传递函数在S 右半平面有一个极点,即P=1,从Nyquist 曲线可看出,奈氏曲线没有包围(-1,0),即R=0,根据奈氏稳定判据,Z=P-R=1,不等于0,所以该系统不稳定,从阶跃响应曲线上也可以看出,系统不稳定。

2

)106)(15()

1(8)(22++++=

s s s s s s G

Bode 图: 程序代码:

num=[8 8];

den=[1 21 100 150 0 0]; w=logspace(-2,3,100); bode(num,den,w) grid 实验结果:

Nyquist图:

程序代码:

num=[8 8];

den=[1 21 100 150 0 0];

nyquist(num,den)

实验结果:

结果分析:

因为开环传递函数在S 右半平面没有极点,即P=0,从Nyquist 曲线可看出,奈氏曲线逆时针包围(-1,0)一圈,即R=1,根据奈氏稳定判据,Z=P-R=-1,不等于0,所以该系统不稳定,从阶跃响应曲线上也可以看出,系统不稳定。

3

)11.0)(105.0)(102.0()

13/(4)(++++=

s s s s s s G

Bode 图: 程序代码:

num=[4/3 4];

den=[0.0001 0.008 0.152 1 0]; w=logspace(-2,3,100); bode(num,den,w) grid 实验结果:

Nyquist图:

程序代码:

num=[4/3 4];

den=[0.0001 0.008 0.152 1 0];

nyquist(num,den)

实验结果:

结果分析:

因为开环传递函数在S右半平面没有极点,即P=0,从Nyquist曲线可看出,奈氏曲线没有包围(-1,0),即R=0,根据奈氏稳定判据,Z=P-R=0,所以该系统不稳定,从阶跃响应曲线上也可以看出,系统阶跃响应最终趋于稳定,所以系统稳定。

体会:本次实验内容划算简单,就是实验结果的分析有点复杂,本次实验总结有:频域分析法分析系统具有很多优点,控制系统及其元部件的频率特性可以用分析法和实验法获得,并可用多种形式的曲线表示,因而系统分析和控制器的设计可以应用图解法进行。频域分析法不仅适用于线性定常系统,还可以推广应用于某些非线性控制系统。通过这次实验,我学会了用MATLAB来分析系统的频域特性,频域特性的图解法主要有,Nyquist曲线、Bode图,Nyquist曲线和Bode图主要用来分析系统的开环频率特性,手工绘制Nyquist曲线、Bode图很麻烦,而高阶系统只能大概地绘出,这给我们分析系统带来了很大的不便,使用MATLAB 软件可以方便而精确地绘制出Nyquist曲线、Bode图和,使得我们分析和设计系统更加方便。

建议:希望老师在以后的实验中多增加几个课时

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

北京理工大学信号与系统实验报告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 可实现系统的传递函数模型和零极点增益模型的转换。调用格

用MATLAB解决线性代数问题实验报告

实验三使用MATLAB解决线性代数问题学院:数计学院班级:1003班姓名:黄晓丹学号:1051020144 实验目的: 学习MATLAB有关线性代数运算的指令,主要学习运用MATLAB解决矩阵除法,线性方程组的通解,矩阵相似 对角化问题,以及解决投入产出分析等应用问题。 实验内容: 矩阵转置:A=[1 2;3 4];B=[4 3;2 1]; >> A',B' ans = 1 3 2 4 ans = 4 3 3 1 矩阵加减:A-B ans= -3 -1 1 3 矩阵乘法:A*B,A.*B(数组乘法)||比较矩阵乘法与数组乘法的区别ans= 8 5 20 13 ans= 4 6 6 4 矩阵除法:A\B,B./A ans=

-6 -5 5 4 ans= 4 1.5 0.6667 0.25 特殊矩阵生成:zeros(m,n)||生成m行n列的矩阵 ones(m,n)||生成m行n列的元素全为一的矩阵 eye(n)||生成n阶单位矩阵 rand(m,n)||生成m行n列[0 ,1]上均匀分布随 机数矩阵 zeros(2,3) ans = 0 0 0 0 0 0 >> ones(3,3) ans = 1 1 1 1 1 1 1 1 1 >> eye(3)

ans = 1 0 0 0 1 0 0 0 1 >> rand(2,4) ans = Columns 1 through 3 0.9501 0.6068 0.8913 0.2311 0.4860 0.7621 Column 4 0.4565 0.0185 矩阵处理:trace(A)||返回矩阵的迹 diag(A)||返回矩阵对角线元素构成的向量 tril(A)||提取矩阵的下三角部分 triu(A)||提取矩阵的上三角部分 flipud(A)||矩阵上下翻转 fliplr(A)||矩阵左右翻转 reshape(A,m,n)||将矩阵的元素重排成m行n列矩阵A=[1 2 3;4 5 6;7 8 9]; >> t=trace(A),d=diag(A),u=triu(A)

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

西华大学实验报告(理工类) 开课学院及实验室: 实验时间 : 年 月 日 一、实验目的 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 。毒=:

(实验三)连续时间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度使用)。之后出现“另存为”对话框,把新模具命名为“数据流程图”,这样我们就把圆形形添加

连续时间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 函数来分别求得。

系统分析实验报告

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

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

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

(2)管理员端

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

2018年自动控制原理期末考试题[附答案解析]

. 2017 年自动控制原理期末考试

卷与答案 一、填空题(每空1分,共20分) 1、对自动控制系统的基本要求可以概括为三个方面,即:稳定性、快速性和准确性。 2、控制系统的输出拉氏变换与输入拉氏变换在零初始条件下的比值称为传递函数。 3、在经典控制理论中,可采用劳斯判据(或:时域分析法)、根轨迹法或奈奎斯特判据( 或:频域分析 法) 等方法判断线性控制系统稳定性。 4、控制系统的数学模型,取决于系统参数, 与外作用及初始条件无关。和结构 A( )L( )lg) 或:,横坐标为( 、线性系统的对数幅频特性,纵坐标取值为5。20lg

,其中P 是指开环传函中具有正实部的极点的个数,6、奈奎斯特稳定判据中,Z = P - RZ是指闭环 传函中具有正实部的极点的个数,R 指奈氏曲线逆时针方向包围(-1, j0 )整圈数。 定义为调整时间。%是超调量。、在二阶系统的单位阶跃响应图中,7t s K)A(22 (T)1K),则其开环幅频特性为(T8、设系统的开环传递函数为1,相12 s(Ts 1)(T s 1)21110) (T ) tg (Ttg 。频特性为()9021 9、反馈控制又称偏差控制,其控制作用是通过给定值与反馈量的差值进行的。 0.5t0.2 t,则该系统的传递函数G(s) 为、若某系统的单位脉冲响应为10。5e g (t) 10e 510 0.2 sss 0.5s

11、自动控制系统有两种基本控制方式,当控制装置与受控对象之间只有顺向作用而无反向联系时,称 为开环控制系统;当控制装置与受控对象之间不但有顺向作用而且还有反向联系时,称为闭环控制系 统;含有测速发电机的电动机速度控制系统,属于闭环控制系统。 、根轨迹起始于开环极点,终止于开环零点。12 、稳定是对控制系统最基本的要求,若一个控制系统的响应曲线为衰减振荡,则该系统13稳定。判断 页脚

Matlab 使用之线性代数综合实例讲解

一、上机目的 1、培养学生运用线性代数的知识解决实际问题的意识、兴趣和能力; 2、掌握常用计算方法和处理问题的方法; 二、上机内容 1、求向量组的最大无关组; 2、解线性方程组; 三、上机作业 1、设A=[2 1 2 4; 1 2 0 2; 4 5 2 0; 0 1 1 7]; 求矩阵A列向量组的一个最大无关组. >> A=[2 1 2 4;1 2 0 2;4 5 2 0;0 1 1 7] A = 2 1 2 4 1 2 0 2 4 5 2 0 0 1 1 7 >> rref(A) ans = 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 所以矩阵A的列向量组的一个最大无关组就是它本身; 2、用Matlab解线性方程组 (1) >> A=[2 4 -6;1 5 3;1 3 2] A = 2 4 -6 1 5 3 1 3 2 >> b=[-4;10;5]

b = -4 10 5 >> x=inv(A)*b x = -3.0000 2.0000 1.0000 >> B=[3 41 -62;4 50 3;11 38 25] B = 3 41 -62 4 50 3 11 38 25 >> c=[-41;100;50] c = -41 100 50 >> x=inv(B)*c x = -8.8221 2.5890 1.9465 3、(选作)减肥配方的实现 设三种食物每100克中蛋白质、碳水化合物和脂肪的含量如下表,表中还给出了20世纪80年代美国流行的剑桥大学医学院的简捷营养处方。现在的问题是:如果用这三种食物作为每天的主要食物,那么它们的用量应各取多少才能全面准确地实现这个营养要求? 四、上机心得体会

系统分析实验报告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):文本

自动控制原理-线性系统的频域分析实验报告

自动调节系统频域分析 班级11081801 学号1108180135 姓名王佳炜 日期2014.1.5

线性系统的频域分析 一、实验目的 1.掌握用MATLAB 语句绘制各种频域曲线。 2.掌握控制系统的频域分析方法。 二、实验内容 1.典型二阶系统 2 2 22)(n n n s s s G ωζωω++= 绘制出6=n ω,1.0=ζ,0.3,0.5,0.8,2的bode 图,记录并分析ζ对系统bode 图的影响。 解: 程序如下: num=[0 0 36];den1=[1 1.2 36];den2=[1 3.6 36]; den3=[1 6 36];den4=[1 9.6 36];den5=[1 24 36]; w=logspace(-2,3,100); bode(num,den1,w) grid hold bode(num,den2,w) bode(num,den3,w) bode(num,den4,w) bode(num,den5,w)

-100-80-60-40-200 20M a g n i t u d e (d B )10 -2 10 -1 10 10 1 10 2 10 3 P h a s e (d e g ) Bode Diagram Frequency (rad/sec) 分析:随着.0=ζ的增大 ,伯德图在穿越频率处的尖峰越明显,此处用渐近线代替时误差越大. 2.系统的开环传递函数为 ) 5)(15(10 )(2 +-= s s s s G ) 106)(15() 1(8)(22++++= s s s s s s G ) 11.0)(105.0)(102.0() 13/(4)(++++= s s s s s s G 绘制系统的Nyquist 曲线、Bode 图,说明系统的稳定性,并通过绘制阶跃响应曲线验证。 解: 程序如下 奈氏曲线: (1) num1=[0,0,10];den1=conv([1,0],conv([1,0],conv([5,-1],[1,5]))); w=logspace(-1,1,100); nyquist(num1,den1,w)

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

理工大学信号与系统实验报告连续时间系统的 复频域分析 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)分子分母根的求解,调用格式如下:

Matlab线性代数实验指导书

Matlab线性代数实验指导书 理学院线性代数课程组 二零零七年十月

目录 一、基础知识 (1) 1.1、常见数学函数 (1) 1.2、系统在线帮助 (1) 1.3、常量与变量 (2) 1.4、数组(矩阵)的点运算 (3) 1.5、矩阵的运算 (3) 二、编程 (4) 2.1、无条件循环 (4) 2.2、条件循环 (5) 2.3、分支结构 (5) 2.4、建立M文件 (6) 2.5、建立函数文件 (6) 三、矩阵及其运算 (7) 3.1、矩阵的创建 (7) 3.2、符号矩阵的运算 (11) 四、秩与线性相关性 (14) 4.1、矩阵和向量组的秩以及向量组的线性相关性 (14) 4.2、向量组的最大无关组 (14) 五、线性方程的组的求解 (16) 5.1、求线性方程组的唯一解或特解(第一类问题) (16) 5.2、求线性齐次方程组的通解 (18) 5.3、求非齐次线性方程组的通解 (19) 六、特征值与二次型 (22) 6.1、方阵的特征值特征向量 (22) 6.2、正交矩阵及二次型 (23)

一、基础知识 1.1常见数学函数 函数数学计算功能函数数学计算功能 abs(x) 实数的绝对值或复数的幅值floor(x) 对x朝-∞方向取整acos(x) 反余弦arcsinx gcd(m,n) 求正整数m和n的最大公约数acosh(x) 反双曲余弦arccoshx imag(x) 求复数x的虚部angle(x) 在四象限内求复数x的相角lcm(m,n)求正整数m和n的最小公倍 自然对数(以e为底数) asin(x) 反正弦arcsinx log(x) 常用对数(以 10 为底数) asinh(x) 反双曲正弦arcsinhx log10(x) atan(x) 反正切arctanx real(x) 求复数 x 的实部atan2(x,y) 在四象限内求反正切rem(m,n) 求正整数m和n的m/n之余数atanh(x) 反双曲正切arctanhx round(x) 对x四舍五入到最接近的整数 符号函数:求出 x 的符号ceil(x) 对x朝+∞方向取整 sign(x) conj(x) 求复数x的共轭复数 sin(x) 正弦sinx 反双曲正弦sinhx cos(x) 余弦cosx sinh(x) cosh(x) 双曲余弦coshx sqrt(x) 求实数x的平方根exp(x) 指数函数e x tan(x) 正切tanx fix(x) 对 x 朝原点方向取整 tanh(x) 双曲正切tanhx 如:输入 x=[-4.85 -2.3 -0.2 1.3 4.56 6.75],则: ceil(x)= -4 -2 0 2 5 7 fix(x) = -4 -2 0 1 4 6 floor(x) =-5 -3 -1 1 4 6 round(x) = -5 -2 0 1 5 7 1.2 系统的在线帮助 1.2.1 help 命令: 1.当不知系统有何帮助内容时,可直接输入 help以寻求帮助: >> help(回车) 2.当想了解某一主题的内容时,如输入: >> help syntax (了解Matlab的语法规定) 3.当想了解某一具体的函数或命令的帮助信息时,如输入: >> help sqrt (了解函数sqrt的相关信息) 1.2.2 lookfor 命令 现需要完成某一具体操作,不知有何命令或函数可以完成,如输入: >> lookfor line (查找与直线、线性问题有关的函数) 1.3 常量与变量

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

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

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

实验4:连续系统的频域分析

实验4:连续系统的频域分析 一、实验目的 (1)掌握连续时间信号的傅里叶变换和傅里叶逆变换的实现方法。 (2)掌握傅里叶变换的数值计算方法和绘制信号频谱的方法。 二、实验原理 1.周期信号的分解 根据傅里叶级数的原理,任何周期信号都可以分解为三角级数的组合——称为 ()f t 的傅里叶级数。在误差确定的前提下,可以由一组三角函数的有限项叠加而得到。 例如一个方波信号可以分解为: 11114111 ()sin sin 3sin 5sin 7357E f t t t t t ωωωωπ?? = ++++ ??? 合成波形所包含的谐波分量越多,除间断点附近外,它越接近于原波形,在间断点附近,即使合成的波形所含谐波次数足够多,也任存在约9%的偏差,这就是吉布 斯现象(Gibbs )。 2.连续时间信号傅里叶变换的数值计算 由傅里叶变换的公式: ()()lim ()j t j n n F j f t e dt f n e ωωττωττ∞ ∞ ---∞ →=-∞ ==∑ ? 当 ()f t 为时限信号时,上式中的n 取值可以认为是有限项N ,则有: ()(),0k N j n n F k f n e k N ωτττ-==≤≤∑,其中2k k N π ωτ = 3.系统的频率特性 连续LTI 系统的频率特性称为频率响应特性,是指在正弦信号激励作用下稳态响应随激励信号频率的变化而变化的情况,表示为 () ()() Y H X ωωω= 三、实验内容与方法 1.周期信号的分解 【例1】用正弦信号的叠加近似合成一个频率为50Hz 的方波。 MATLAB 程序如下: clear all; fs=10000; t=[0:1/fs:0.1]; f0=50;sum=0; subplot(211) for n=1:2:9 plot(t,4/pi*1/n*sin(2*pi*n*f0*t),’k ’); hold on; end title(‘信号叠加前’); subplot(212) for n=1:2:9;

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

学生实验报告实验课程:信号与 系统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、档案查询 功能:对档案卡片的查询功能,在这里可以查到符合程序要求的任何信息。

自动控制原理线性系统的频域分析实验四

武汉工程大学实验报告专业电气自动化班号指导教师 姓名同组者无

M a g n i t u d e (d B ) 10 1010101010P h a s e (d e g )Frequency (rad/sec) 当3.0=ζ时,程序如下: num=[0 0 36];den=[1 3.6 36];w=logspace(-2,3,100);bode(num,den,w) grid M a g n i t u d e (d B ) 10 1010101010P h a s e (d e g )Bode Diagram Frequency (rad/sec) 当5.0=ζ时,程序如下: num=[0 0 36];den=[1 6 36];w=logspace(-2,3,100);bode(num,den,w) grid

M a g n i t u d e (d B ) 10 1010101010P h a s e (d e g )Frequency (rad/sec) 当8.0=ζ时,程序如下: num=[0 0 36];den=[1 9.6 36];w=logspace(-2,3,100);bode(num,den,w) grid M a g n i t u d e (d B ) 10 1010101010P h a s e (d e g ) Bode Diagram Frequency (rad/sec) 当2=ζ时,程序如下: num=[0 0 36];den=[1 24 36];w=logspace(-2,3,100);bode(num,den,w) grid

电力系统分析实验报告金科

学生实验报告 (理工类) 课程名称:专业班级: 学生学号:学生: 所属院部:指导教师: 20 13 ——20 14 学年第二学期 金陵科技学院教务处制 实验一电力系统分析计算 实验项目名称:电力系统分析计算实验学时: 2

同组学生:实验地点: C208 实验日期: 2014 6 23 实验成绩: 批改教师:静批改时间: 一.实验目的 1.掌握用Matlab软件编程计算电力系统元件参数的方法. 2.通过对不同长度的电力线路的三种模型进行建模比较,学会选取根据电路要求选取模 型。 3.掌握多级电力网络的等值电路计算方法。 4.理解有名制和标幺制。 二.实验容 1.电力线路建模 有一回220kV架空电力线路,导线型号为LGJ-120,导线计算外径为15.2mm,三相导线水平排列,两相邻导线之间的距离为4m。试计算该电力线路的参数,假设该线路长度分别为60km,200km,500km,作出三种等值电路模型,并列表给出计算值。 2.多级电力网络的等值电路计算 部分多级电力网络结线图如图1-1所示,变压器均为主分接头,作出它的等值电路模型,并列表给出用有名制表示的各参数值和用标幺制表示的各参数值。 线路额定电压电阻 (欧/km) 电抗 (欧/km) 电纳 (S/km) 线路长度 (km) L1(架空线)220kv 0.08 0.406 2.81*10-6 200 L2(架空线)110kV 0.105 0.383 2.81*10-6 60 L3(架空线)10kV 0.17 0.38 忽略15 变压器额定容量P k(kw) U k% I o% P o(kW) T1 180MVA 893 13 0.5 175 T2 63MVA 280 10.5 0.61 60 三.实验设备 1.PC一台 2.Matlab软件 四.实验记录 1.电力线路建模 电阻电抗电纳电阻电抗电纳电阻电抗电纳

自动控制原理实验六线性系统的频域分析

实验六 线性系统的频域分析 一. 实验目的 (1)熟练掌握使用MA TLAB 命令绘制控制系统Nyquist 图的方法; (2)能够分析控制系统Nyquist 图的基本规律; (3)加深理解控制系统乃奎斯特稳定性判据的实际应用; (4)学会利用奈氏图设计控制系统; (5)熟练掌握运用MA TLAB 命令绘制控制系统伯德图的方法; (6)了解系统伯德图的一般规律及其频域指标的获取方法; (7)熟练掌握运用伯德图分析控制系统稳定性的方法; (8)设计超前校正环节并绘制Bode 图; (9)设计滞后校正环节并绘制Bode 图。 二. 实验原理及内容 1、频率特性函数)(ωj G 。 频率特性函数为: 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) 2、用MATLAB 作奈魁斯特图。 控制系统工具箱中提供了一个MA TLAB 函数nyquist( ),该函数可以用来直接求解Nyquist 阵列或绘制奈氏图。当命令中不包含左端返回变量时,nyquist ()函数仅在屏幕上产生奈氏图,命令调用格式为: nyquist(num,den) ; 作Nyquist 图, nyquist(num,den,w); 作开环系统的奈氏曲线, 3、奈奎斯特稳定性判据(又称奈氏判据) 反馈控制系统稳定的充分必要条件是当ω从-∞变到∞时,开环系统的奈氏曲线不穿过点(-1,j0)且逆时针包围临界点(-1,j0)点的圈数R 等于开环传递函数的正实部极点数。 4、用MATLAB 作伯德图 控制系统工具箱里提供的bode()函数可以直接求取、绘制给定线性系统的伯德图。 命令的调用格式为: [mag,phase,w]=bode(num,den) [mag,phase,w]=bode(num,den,w) 由于伯德图是半对数坐标图且幅频图和相频图要同时在一个绘图窗口中绘制,因此,要用到半对数坐标绘图函数和子图命令。 (1) 对数坐标绘图函数 利用工作空间中的向量x ,y 绘图,要调用plot 函数,若要绘制对数或半对数坐标图,只需要用相应函数名取代plot 即可,其余参数应用与plot 完全一致。 (2) 子图命令

相关文档
最新文档