Bruce Lin - Conceptual Design and Modeling of a Fuel Cell Scooter for Urban Asia

= more smoothing.
%
% If I wanted to close the loop I would use this
% [numc,denc]=feedback(num,den,gclose,1);
s0=1.5; k=1; gclose=1;
num=k;
den=[1 s0];
v2=transpose(lsim(num,den,v,t))*s0;
v=v2; clear v2;
%% EXAMINE DRIVING CYCLE
%
% Evaluate driving cycle to calculate acceleration from finite

290

%
%
%
%
%
%
%
%
%
%

differences of velocity; also, determine what fraction of the
driving cycle is spent accelerating, decelerating, etc.
plustime
minustime
steadytime

if positive acceleration
deceleration
steady velocity, no acceleration

We also accumulate acceleration values to separate out
average of positive accelerations and average of
negative accelerations.

plustime=0; minustime=0; steadytime=0;
for i=2:1:cyclelength
a(i)=(v(i)-v(i-1))/timestep;
if a(i)>0,
aplus(i)=a(i);
aminus(i)=0;
plustime=plustime+1;
elseif a(i) --">