控制系统计算机仿真作业

  • 格式:doc
  • 大小:122.61 KB
  • 文档页数:29

下载文档原格式

  / 29
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

兰州理工大学

《控制系统计算机仿真》

上机报告Ⅰ

院系:电气工程与信息工程学院

班级:14级自动化3班

姓名:孙悦

学号: 1405220323

时间: 2017 年 6 月 15 日

电气工程与信息工程学院

《控制系统计算机仿真》上机实验任务书Ⅰ(2017)

一、上机实验内容及要求

1.matlab软件

要求利用课余时间熟悉掌握matlab软件的基本数值运算、基本符号运算、基本程序设计方法及常用的图形命令操作。

2.各章节仿真实验内容及要求

具体实验内容及要求请详见上机实验报告。

二、上机实验时间安排及相关事宜

1.依据课程教学大纲要求,上机实验学时共16学时,学生须在每次上机之前

做好相应的准备工作,以确保在有限的机时内完成仿真实验要求的内容;

2.实验完成后按规定完成相关的仿真实验报告;

3.仿真实验报告请按有关样本制作并A4打印,侧面装订,作为成绩评定的一

部分。

自动化系《控制系统计算机仿真》课程组

2017年3月

一、Matlab 基础操作

1-1用MATLAB 语言求下列系统的状态方程、传递函数、零极点增益和部分分式形式的模型参数,并分别写出其相应的数学模型表达式:

(1)2450351024

247)(23423+++++++=s s s s s s s s G

程序如下:

num=[7,24,24]

den=[10,35,50,24]

[A,B,C,D]=tf2ss(num,den)

系统的状态方程:

A =

-3.5000 -5.0000 -2.4000

1.0000 0 0

0 1.0000 0

B =

1

C =

0.7000 2.4000 2.4000

D =

零极点增益形式:

[Z,P,K]=tf2zp(num,den)

Z =

-1.7143 + 0.6999i

-1.7143 - 0.6999i

P =

-1.2973 + 0.9838i

-1.2973 - 0.9838i

-0.9053

K =

0.7000

部分分式:

[R,P,H]=residue(num,den)

R =

-0.0071 - 0.2939i

-0.0071 + 0.2939i

0.7141

P =

-1.2973 + 0.9838i

-1.2973 - 0.9838i

-0.9053

H =

[]

(2)u X X ⋅⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡+⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡------------=⋅022475.025.075.125.1125.15.025.025.025.125.425.25.025

.1525.2 []X y 2020=

A=[2.25 -5 -1.25 -0.5;2.25 -4.25 -1.25 -0.25;0.25 -0.5 -1.25 -1;1.25 -1.75 -0.25 -0.75]

A =

2.2500 -5.0000 -1.2500 -0.5000

2.2500 -4.2500 -1.2500 -0.2500

0.2500 -0.5000 -1.2500 -1.0000

1.2500 -1.7500 -0.2500 -0.7500 >> B=[4;2;2;0]

B =

4

2

2

>> C=[0 2 0 2]

C =

0 2 0 2

>> D=0

D =

零极点增益形式:

>> [Z,P,K]=ss2zp(A,B,C,D)

Z =

-1.0000 + 1.2247i

-1.0000 - 1.2247i

-1.5000

P =

-0.5000 + 0.8660i

-0.5000 - 0.8660i

-1.5000 + 0.0000i

-1.5000 - 0.0000i

K =

4.0000

传递函数形式:

>> num=[0 4 14 22 15]

num =

0 4 14 22 15

>> den=[1 4 6.25 5.25 2.25]

den =

1.0000 4.0000 6.2500 5.2500

2.2500

部分分式:

>> [R ,P ,H]=residue(num,den)

R =

4.0000

-0.0000

0.0000 - 2.3094i

0.0000 + 2.3094i

P =

-1.5000

-1.5000

-0.5000 + 0.8660i

-0.5000 - 0.8660i

H =

[]

1-2 用殴拉法matlab 编程实现下列系统的输出响应)(t y 在10≤≤t 上,1.0=h 时的数值解。

y y -=',1)0(=y

要求保留4位小数,并将结果以图形的方式与真解t e t y -=)(比较。

t=0:0.1:1

t =