matlab机电耦合系统相位分岔图

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、参考论文

略