MATLAB基础与控制系统CAD仿真练习1(19)

控制系统计算机辅助设计,matlab仿真的一些小练习

其源程序见zuoye31.m中%（2）对应的部分。 （3）A的行列式为：46； A的迹为：14； A的秩为：3

A的特征值为：2.0000 2.3944 9.6056

A的特征多项式为：1.0000*x^3-14.0000*x^2+47.0000*x-46.0000

此外，在Command Window中输入命令行：

A=[7 3 -2;3 4 -1;-2 -1 3];

q=sym(A); dxs=poly(q)

也可以得到A的特征多项式：

dxs =

x^3-14*x^2+47*x-46

2．考虑如下Lorenz模型，其微分方程为： 1(t) x1(t) x2(t)x3(t)x 2(t) x2(t) x3(t)x

3(t) x2(t)x1(t) x2(t) x3(t)x

10, 28, 8/3其中参数为：，试写出该方程的S函数表示，并分

别通过Matlab 编程和Simulink对其进行仿真分析并绘制图形。

解答：

（1）用该方程的S函数绘制图形： 方法一、

根据题目所给Lorenz模型的相关信息修改S-Function的模板，得到M-file并命名为Lorenz.m：

function [sys,x0] = Lorenz(t,x,u,flag) switch flag, case 0,

sys=[3;0;3;0;0;0]; %3个连续状态;0个离散状态;3个输出变量;0个...

... 输入变量; 没有直接馈通，即输出y中没有输入参数u;采用时