1、内容简介
略
2、内容说明
3、仿真分析
clc
close all
clear
% 基本参数
C = 28;
S = 10;
B1 = 0.7;
B2= 0.3;
a1 = 10;
a2 = 2;
L = 10;
A1 = 8;
A2 = 5;
%% 组织规制型安全生产管理模式决策复制动态
% 第一种情况 组织采取任何安全生产管理模式均是稳定状态
z = (C-S)/(a1*(1-B1)*L-a2*(1-B2)*L)
t = 0:0.01:10;
Fy = zeros(1,length(t));
figure
plot(t,Fy)
xlabel t/时间
ylabel Fy
title 组织采取任何安全生产管理模式均是稳定状态
% 第二种情况，分2种情况
% 当z>....的时候，y=1是稳定点
z = 2;
y0 = 0.2; % 初始值
t0=0;
tend = 10; % 仿真时间
[t,y] = ode45(@hanshu,[t0,tend],y0); % 调用ode45 求解微分方程2-3
figure
plot(t,y)
xlabel t/时间
ylabel Fy
title y=1是演化稳定点
% 当z<....的时候，y=0是稳定点
z = 0.6;
y0 = 0.2; % 初始值
t0=0;
tend = 10; % 仿真时间
[t,y] = ode45(@hanshu1,[t0,tend],y0); % 调用ode45 求解微分方程2-3
figure
plot(t,y)
xlabel t/时间
ylabel Fy
title y=0是演化稳定点
4、参考论文
略
