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
| figNumber=figure( ...
'Name','Pendule simple', ...
'NumberTitle','off', ...
'Visible','off', ...
'DoubleBuffer','on', ...
'Color', [1 1 1], ...
'BackingStore','off');
axes( ...
'Units','normalized', ...
'Position',[0.05 0.05 0.75 0.90], ...
'YGrid','on','XGrid','on',...
'NextPlot','add');
figure(figNumber);
L=10; % longueur du fil
A=pi/4; % amplitude maximale
x=[0 0];
y=[0 0];
y0=10;
h=plot(x,y);
set(h,'color',[1 0 0]);
set(h,'lineWidth',1);
h1=plot(-L*sin(-A),-L*cos(-A)+y0,'.');
set(h1,'MarkerSize',60);
set(h1,'color',[1 1 0]);
set(h1,'lineWidth',38);
axis([-10 10 -10 10])
g=9.81;
T=2*pi*sqrt(L/g);amort=0.1;
pas=0.01;
%t=0;
%=======================================
for t =0:pas:100
teta=-A*cos(2*pi/T*t)*exp(-amort*t);
x=[0 -L*sin(teta)];
y=[y0 -L*cos(teta)+y0];
set(h,'Xdata',x,'YData',y)
set(h1,'Xdata',x(2),'YData',y(2));
drawnow
end |
