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 |
Partager