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
| function [Time w]=vitesse (e,T,n)
w=zeros(1,n+1);
t=0 ; % temps initial
Time=zeros(1,n+1);
for i=1:n+1
Time(i+1)=t+i*T/n;
w(i+1)=w(i)*(i+1)*T/n+e;
t=t+T/n;
end
end
function [Time Cm]=couple (T,n,e,J,a)
Cm=zeros(1,n+1);
t=0 ; % temps initial
Time=zeros(1,n+1);
e=24000;
T=0.3*3600;
n=10000000;
[Time w]=vitesse (e,T,n);
for i=1:n+1
Time(i) = t+(i-1)*T/n;
t = t+i*T/n;
Cm(i+1)=J/T/n*((i+1).*vitesse(i+1)+T/n*a*(vitesse(i+1)).^2+e);
end
end
clear all
close all
clc
e=24000;
T=0.3*3600;
n=10000;
J=120;
a=6;
[Time w]=vitesse(e,T,n);
[Time Cm]=couple(T,n,e,J,a);
plot(Time,w)
hold on
plot(Time,Cm) |
Partager