Dear Everyone
I’m working on my Phd thsesis on measure synchronization under nonlinear dynamics.am currently stuck with the average bare energy and interaction energy code for my work.i have gotten the Time series,phase plot,poincare section,frquency and power spectrum of my work.

am trying to reproduce the avareage bare energy code for ”Exploring the route to measure synchronization in non-linearly coupled Hamiltonian systems” journal. I have done the phase plot,time series ,frequency and power spectrum of this jounal.
this is the code i try to use for the average bare energy of this journal and its not working.

function dxdt = exploring(t,x,y)
dxdt = zeros(4,1);
dxdt(1) = x(3);
dxdt(2) = x(4);
dxdt(3) = -x(1) + y*x(1)*x(2^2);
dxdt(4) = -x(2) + y*x(1)^2*x(2);

%% %% Average bear energy
clear all; close all; clc; %#ok<CLALL>
y=[0 0.07 0.15 0.52];
m=length(y);
tspan=[500:0.001:1400];
T=0:0.01:1400;
L=10;
E=0.2;
p1=0.2;
p2=sqrt(2*E – (p1)^2);
E1=((p1)^2 + y*(p1)^2)/2;
E2=((p2)^2 + y*(p2)^2)/2;
R =(1./ T) * integral((E1*E2),0,T);
for i=1:m
[t,R] = ode45(@(t,R) exploring(t,R,y(i)),tspan,[0 0 p1 p2]);
end
figure(i)
plot(t,R(:,1),’-b’);
xlabel(‘t’);
ylabel(‘X’);
hold on
plot(t,R(:,2),’-r’)
xlabel(‘t’);
ylabel(‘x1,x2’);%ylim([-1 1]);
I need your assistance, support and suggestion.
i would welcome any code with related to the average bare energy and intereaction energy for any oscillator.
ThanksDear Everyone
I’m working on my Phd thsesis on measure synchronization under nonlinear dynamics.am currently stuck with the average bare energy and interaction energy code for my work.i have gotten the Time series,phase plot,poincare section,frquency and power spectrum of my work.

am trying to reproduce the avareage bare energy code for ”Exploring the route to measure synchronization in non-linearly coupled Hamiltonian systems” journal. I have done the phase plot,time series ,frequency and power spectrum of this jounal.
this is the code i try to use for the average bare energy of this journal and its not working.

function dxdt = exploring(t,x,y)
dxdt = zeros(4,1);
dxdt(1) = x(3);
dxdt(2) = x(4);
dxdt(3) = -x(1) + y*x(1)*x(2^2);
dxdt(4) = -x(2) + y*x(1)^2*x(2);

%% %% Average bear energy
clear all; close all; clc; %#ok<CLALL>
y=[0 0.07 0.15 0.52];
m=length(y);
tspan=[500:0.001:1400];
T=0:0.01:1400;
L=10;
E=0.2;
p1=0.2;
p2=sqrt(2*E – (p1)^2);
E1=((p1)^2 + y*(p1)^2)/2;
E2=((p2)^2 + y*(p2)^2)/2;
R =(1./ T) * integral((E1*E2),0,T);
for i=1:m
[t,R] = ode45(@(t,R) exploring(t,R,y(i)),tspan,[0 0 p1 p2]);
end
figure(i)
plot(t,R(:,1),’-b’);
xlabel(‘t’);
ylabel(‘X’);
hold on
plot(t,R(:,2),’-r’)
xlabel(‘t’);
ylabel(‘x1,x2’);%ylim([-1 1]);
I need your assistance, support and suggestion.
i would welcome any code with related to the average bare energy and intereaction energy for any oscillator.
Thanks Dear Everyone
I’m working on my Phd thsesis on measure synchronization under nonlinear dynamics.am currently stuck with the average bare energy and interaction energy code for my work.i have gotten the Time series,phase plot,poincare section,frquency and power spectrum of my work.

am trying to reproduce the avareage bare energy code for ”Exploring the route to measure synchronization in non-linearly coupled Hamiltonian systems” journal. I have done the phase plot,time series ,frequency and power spectrum of this jounal.
this is the code i try to use for the average bare energy of this journal and its not working.

function dxdt = exploring(t,x,y)
dxdt = zeros(4,1);
dxdt(1) = x(3);
dxdt(2) = x(4);
dxdt(3) = -x(1) + y*x(1)*x(2^2);
dxdt(4) = -x(2) + y*x(1)^2*x(2);

%% %% Average bear energy
clear all; close all; clc; %#ok<CLALL>
y=[0 0.07 0.15 0.52];
m=length(y);
tspan=[500:0.001:1400];
T=0:0.01:1400;
L=10;
E=0.2;
p1=0.2;
p2=sqrt(2*E – (p1)^2);
E1=((p1)^2 + y*(p1)^2)/2;
E2=((p2)^2 + y*(p2)^2)/2;
R =(1./ T) * integral((E1*E2),0,T);
for i=1:m
[t,R] = ode45(@(t,R) exploring(t,R,y(i)),tspan,[0 0 p1 p2]);
end
figure(i)
plot(t,R(:,1),’-b’);
xlabel(‘t’);
ylabel(‘X’);
hold on
plot(t,R(:,2),’-r’)
xlabel(‘t’);
ylabel(‘x1,x2’);%ylim([-1 1]);
I need your assistance, support and suggestion.
i would welcome any code with related to the average bare energy and intereaction energy for any oscillator.
Thanks signal processing MATLAB Answers — New Questions

​