1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
| function lancer2(obj,event)
data2=guidata(gcbf);
P = str2num(get(data2.pulsationspropres1,'string'));
W = str2num(get(data2.condinitiale1,'string'));
K = str2num(get(data2.matricedinter,'string'));
wp=str2num(get(data2.pulsationP,'string'));
Kp=str2num(get(data2.KP,'string'));
X=size(P);
n=X(1,2);
w0 = [n,Kp,W,P];
for l=1:n-1
for s=l+1:n
w0=[w0,K(l,s)];
end
end
w0(1,2*n+3+n*(n-1)/2)=wp;
tspan=0:0.001:5;
options = odeset('RelTol',1e-12);
[t,w]=ode45(@kuramuto2,tspan,w0',options);
A=size(t);
A(1,1);
M=zeros(size(t));
for m=3:n+2
M=M+exp(i*w(:,m));
end
r=abs(M/n);
ksi=angle(M/n);
if (data2.Vrp==1)
figure('units','pixels',...
'position',[250 250 500 500],...
'color',[0.925 0.913 0.687],...
'numbertitle','off',...
'name','graphe de r',...
'menubar','figure',...
'tag','graphe');
sli=uicontrol('Style','slider',...
'Min' ,0,'Max',Kp+50, ...
'Position',[10,220,130,20], ...
'Value', Kp,...
'SliderStep',[0.001 0.1], ...
'CallBack',@slider(w0,n),...
'tag','slider');
data3=guihandles(gcf);
guidata(gcf,data3);
plot(t,r);
end |
