《Matlab/Simulink与控制系统仿真》程序指令总结

Matlab_Simulink_BookExample11. 非线性系统11.1 Matlab 函数例题 Eg11_2

书中详细实例代码可见：Github

Matlab_Simulink_BookExample

图书：《Matlab/Simulink与控制系统仿真》

11. 非线性系统

11.1 Matlab 函数

ode45()求解微分方程

odefuntspant0tfinal

例题 Eg11_2

% Page307：已知二阶线性系统的微分方程，绘制系统的相平面图和零输入响应曲线clear;clc;% Eg11_2_test 是系统微分方程的描述函数% 初始化状态变量为 [0,1]，计算时间为 [0,20][t, x] = ode45('Eg11_2_test', [0,20], [0,1]);% 初始化状态变量为 [1,1]，计算时间为 [0,20][t1, x1] = ode45('Eg11_2_test', [0,20], [1,1]);% 绘制相轨迹plot(x(:,1),x(:,2),'-', x1(:,1),x1(:,2),'.')grid% 添加图标题title('相轨迹图')% 绘制时域响应曲线plot(t,x(:,2),'-', t1,x1(:,2),'.')% 添加图标题title('时间响应曲线')

% Eg11_2 程序的子函数代码function xdot = Eg11_2_test(t, x)% 导数关系式xdot = [-2*x(2); x(1)];