济南大学matlab试题,计算机仿真技术期末考核.doc

计算机仿真技术期末考核

济 南 大 学 课 程 考 核( 封 面)

课程名称： 计算机仿真技术 学时 32 学分 2

考生学号： 20080302015 考生姓名 樊军河 班级 电气0803班 成绩

1、考核形式：考查

2、考核具体要求：

MATLAB仿真实验50分

课程感想或论文报告50分

任课教师：李实

2011 年 11 月 15 日

一. MATLAB上机仿真实验(每题10分，共50分)：

1. 已知传递函数为

在MATLAB中将传递函数转化为零极点和状态方程两种模型形式，要求写出MATLAB指令和MATLAB中得到的模型方程.

>> num=[25.04,5.008];

>> den=[1,5.03247,25.1026,5.008];

>> sys_tf=tf(num,den);

>> [A,B,C,D]=tf2ss(num,den); %将tf形式转换为ss形式

>> sys_ss=ss(A,B,C,D);

>> [z,p,k]=tf2zp(num,den); %将zpk形式转换为tf形式

>> sys_zpk=zpk(z,p,k)

Zero/pole/gain:

25.04 (s+0.2)

--------------------------------

(s+0.2078) (s^2 + 4.825s + 24.1)

>> num=[25.04,5.008];

>> den=[1,5.03247,25.1026,5.008];

>> sys_tf=tf(num,den);

>> [A,B,C,D]=tf2ss(num,den); %将tf形式转换为ss形式

>> sys_ss=ss(A,B,C,D)

a =

x1 x2 x3

x1 -5.032 -25.1 -5.008

x2 1 0 0

x3 0 1 0

b =

u1

x1 1

x2 0

x3 0

c =

x1 x2 x3

y1 0 25.04 5.008

d =

u1

y1 0

Continuous-time model.

>> [z,p,k]=tf2zp(num,den); %将zpk形式转换为tf形式

>> sys_zpk=zpk(z,p,k)

Zero/pole/gain:

25.04 (s+0.2)

--------------------------------

(s+0.2078) (s^2 + 4.825s + 24.1)

2. 已知状态空间模型为2u-2y形式：

在MATLAB中将状态方程转化为传递函数，要求写出MATLAB指令和得到的模型方程.

注意：2个输入和2个输出，得到4个传递函数：

G11=y1/u1, G12=y1/u2,G21=y2/u1,G22=y2/u2.

MATLAB指令为：

A=[0 1?; -25 -4]?; B=[1 1?; 0 1]?; C=[1 0?; 0 1]?; D=[0 0?; 0 0]?;

sys_ss=ss(A,B,C,D)?;

[num,den]=ss2tf(A,B,C,D,1);

sys_tf11=tf(num(1,:),den)

sys_tf21=tf(num(2,:),den)

[num,den]=ss2tf(A,B,C,D,2);

sys_tf12=tf(num(1,:),den)

sys_tf22=tf(num(2,:),den)

Transfer function y1/u1

s + 4

--------------

s^2 + 4 s + 25

Transfer function: y2/u1

-25

--------------

s^2 + 4 s + 25

Transfer function: y1/u2

s + 5

--------------

s^2 + 4 s + 25

Transfer function: y2/u2

s - 25

--------------

s^2 + 4 s + 25

3. 已知连续时间传递函数为

取采样时间为0.1s，求离散化的传递函数，要求写出MATLAB指令和MATLAB得到的离散化传递函数。

>> H=tf( [10 10], [1 5 25 6])

Transfer funct