diaaa

Untitled

Oct 30th, 2020
602
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. for nig = 1:NoOfShifts
  2.    
  3.     if nig == NoOfShifts
  4.         DragTimeDistance(1,nig) = vpaintegral( fMig(GBGeometricRatios(nig)) /sigFx(GBGeometricRatios(nig)),...
  5.             symVx , fNe2V(Ne(1) , GBGeometricRatios(nig))/3.6 , 0.99*MaxSpeed/3.6 , 'MaxFunctionCalls',Inf);
  6.        
  7.         DragTimeDistance(2,nig) = vpaintegral( symVx * fMig(GBGeometricRatios(nig)) /sigFx(GBGeometricRatios(nig)),...
  8.             symVx , fNe2V(Ne(1) , GBGeometricRatios(nig))/3.6 , 0.99*MaxSpeed/3.6 , 'MaxFunctionCalls',Inf);
  9.         break
  10.     end
  11.    
  12.     solvx = vpasolve(fFx(GBGeometricRatios(nig))==fFx(GBGeometricRatios(nig+1)),symVx);
  13.     UpShiftSpeeds(nig) = double(solvx(end))*3.6;
  14.    
  15.     solvx = vpasolve(fFx(GBProgressiveRatios(nig))==fFx(GBProgressiveRatios(nig+1)),symVx);
  16.     UpShiftSpeedsProg(nig) = double(solvx(end))*3.6;
  17.    
  18.     DragTimeDistance(1,nig) = vpaintegral( fMig(GBGeometricRatios(nig)) /sigFx(GBGeometricRatios(nig)),...
  19.         symVx , fNe2V(Ne(1) , GBGeometricRatios(nig))/3.6 , UpShiftSpeeds(nig)/3.6 , 'MaxFunctionCalls',Inf);
  20.    
  21.     DragTimeDistance(2,nig) = vpaintegral( symVx * fMig(GBGeometricRatios(nig)) /sigFx(GBGeometricRatios(nig)),...
  22.         symVx , fNe2V(Ne(1) , GBGeometricRatios(nig))/3.6 , UpShiftSpeeds(nig)/3.6 , 'MaxFunctionCalls',Inf);
  23.    
  24. end
RAW Paste Data