《Matlab/Simulink与控制系统仿真》程序指令总结
Matlab_Simulink_BookExample11. 非线性系统11.1 Matlab 函数例题 Eg11_2书中详细实例代码可见:Github
Matlab_Simulink_BookExample
图书:《Matlab/Simulink与控制系统仿真》
11. 非线性系统
11.1 Matlab 函数
ode45()
求解微分方程
odefun
tspan
t0
tfinal
例题 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)];