1、内容简介
略
638-可以交流、咨询、答疑
2、内容说明
略
3、仿真分析
clc
close all
clear
%% parameter
global C
C = 48;
tspan = [0 5]; % 仿真时间
x0 = [0.1 0.1 0.1 0.1]; % 变量初始值
[T,X] = ode45(@model_diff, tspan, x0); % 调用求解器
X3_dot = [0;diff(X(:,3))./diff(T)]; % 求解x3的微分数值
clear
close all
clc
%% parameter
% 绘制分叉图
global C
tspan = [0 5]; % 仿真时间
x0 = [0.1 0.1 0.1 0.1]; % 变量初始值
C_span=10:0.05:60;
figure
for i = 1:length(C_span)
i
disp(['C:',num2str(C_span(i))])
C = C_span(i);
[T,X] = ode45(@model_diff, tspan, x0); % 调用求解器
x3 = X(:,3);
x3max=getmax(x3); % 通过极大值求解分叉图
plot(C,x3max,'b.','markersize',4)
hold on
end
title('映射分岔图')
xlabel('C'),ylabel('|x3|')
hold off
4、参考论文
略
