Top Banner
V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA, July 25 th , 2016 Joint work with Junxing Yang, Radu Grosu, and Scott A. Smolka
31

V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

Aug 17, 2020

Download

Documents

dariahiddleston
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Page 1: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

V-FormationasOptimalControl

AshishTiwari

SRIInternational,MenloPark,CA,USA

BDA,July25th,2016

JointworkwithJunxingYang,RaduGrosu,andScottA.Smolka

Page 2: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

Outline

• Introduction

• TheV-FormationProblem

• ModelPredictiveControlforV-Formation

• ExperimentalResults

• Conclusions&FutureWork

Page 3: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

Outline

• Introduction

• TheV-FormationProblem

• ModelPredictiveControlforV-Formation

• ExperimentalResults

• Conclusions&FutureWork

Page 4: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

V-Formation

• FlocksofbirdsorganizethemselvesintoV-formations

Eurasian Cranes migrating in a V-formation (HamidHajihusseini, Wikipedia)

Reason:Savesenergyasbirdsbenefit fromupwashregion;providesclearvisualfieldwithvisibilityoflateralneighbors

Page 5: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

ReachingaV-Formation

• Rule-basedApproach:ØCombinationsofdynamicalflightrules asdrivingforcesØNotcompletelysatisfying

• ViewasaDistributedControlProblem:ØFlockwantstogetintoanoptimalconfigurationthatprovidesbestview,energybenefit,andstability

• OurApproach:ØUsesModel-Predictive Control(MPC)ØWhichusesParticle-SwarmOptimization(PSO)

Page 6: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

Reynolds’Rules

Reynolds (1987)presentedthreerulesforgeneratingV-formations:

Alignment Cohesion SeparationAlignment:steertowardstheaverageheadingoflocalflockmatesCohesion:steertomovetowardtheaveragepositionoflocalflockmatesSeparation:steertoavoidcrowdinglocalflockmates

Page 7: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

ExtendedReynoldsModel

Reynolds’modelwasextendedbyadditionalrules:• Arulethatforcesabirdtomovelaterallyawayfromanybirdthatblocksitsview(Flake(1998)).

• Dragreductionrule:computingthe induceddraggradientandsteeringalongthisgradient(Dimock&Selig(2003)).

Nathan&Barbosa’smodel(2008):• Coalescing:seekproximityofnearestbird• Gap-seeking:seeknearestpositionaffordingclearview• Stationingrule:movetoupwashofaleadingbird

Page 8: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

ARule-basedAttempt

DesignedrulesthatgenerateaV-formation• Drivebirdstowardstheoptimalupwashpositionw.r.t.thenearestbirdinfront;unsatisfactory solution

Page 9: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

Outline

• Introduction

• TheV-FormationProblem

• ModelPredictiveControlforV-Formation

• ExperimentalResults

• Conclusions&FutureWork

Page 10: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

TheV-FormationProblem

Assumeageneric2-ddynamicmodelofaflockof

birds

xi(t+1)=xi(t)+vi(t+1)

vi(t+1)=vi(t)+ai(t)

Goal:findbestaccelerationsai(t) ateachtimestep

thatwillfinallyleadtoaV-formation.

Thisisadistributedcontrolproblem

Page 11: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

WhatisaV-Formation?

Wewantaformationthatachievestheoptimumvaluesforthefollowingthreefitnessmetrics:1. VelocityMatching2. ClearView3. UpwashBenefit

Page 12: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

VelocityMatching(VM)

s=stateofthen-birds=npositions,nvelocitiesVM(s)=normalizedsumofpairwisevelocitydifferenceVM(s)=0ifallbirdshavethesamevelocityVM(s)increasesasthevelocitiesgetmoremismatched

VMisminimizedwhenallbirdshaveequalvelocity.

Velocitynotmatched Velocitymatched

Page 13: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

ClearView(CV)

• Accumulatethepercentageofthebird’sviewthatisblocked• CV(s)=0ifeverybirdhasa100%clearview• CV(s)increasesasmoreoftheviewofanybirdisblocked

(b) i’s view is completelyblocked by j and k.Clear view: 1

Page 14: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

UpwashBenefit(UB)

• AGaussian-likemodelofupwash anddownwash• UB(s)=sumofupwash benefiteachbirdgetsfromeveryother• UB(s)=1ifn-1birdsgetsmaxpossibleUBbenefit• UB(s)increasesasbirdsgetlesserupwash benefit

Page 15: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

FitnessFunction

Fitnessofastateisasum-of-squares

combinationofVM,CVandUB

F(s)=(VM(s)-VM(s*))2+(CV(s)-CV(s*))2+(UB(s)-UB(s*))2

• stateachieving optimal fitnessvalue (i.e.,aV-

formation)

Page 16: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

TheV-FormationProblem

Assumeageneric2-ddynamicmodelofaflockof

birds

xi(t+1)=xi(t)+vi(t+1)

vi(t+1)=vi(t)+ai(t)

Goal:findbestaccelerationsai(t) ateachtimestep

thatwillfinallyleadtoastatewithminimumF(s)

Thisisadistributedcontrolproblem

Page 17: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

Outline

• Introduction

• TheV-FormationProblem

• ModelPredictiveControlforV-Formation

• ExperimentalResults

• Conclusions&FutureWork

Page 18: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

ModelPredictiveControl(1)

Ateachtimet,considerhowthemodelwillbehaveinthe

nextTstepsunderdifferentchoicesforthecontrolinputs

• Useamodel thatrepresentsthebehavioroftheplant

Useanoptimizationsolvertofindthebestcontrolinputs

overthisfinitepredictionhorizon

Onlyapplythefirstoptimalcontrolaction

Repeatatt+1

Page 19: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

ModelPredictiveControl(2)

• Attimet+1,updatemodelstatewithnewmeasurementsoftheplant.• Repeattheoptimizationwithnewstates.

AdiscreteMPCscheme(Wikipedia):horizon=p,currenttime=k

Page 20: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

ModelPredictiveControlforV-Formation(1)

Birdi attimetsolvesthefollowingoptimizationproblem:

a*i(t),…,a*i(t+T-1)=argmin ai(t),…,ai(t+T-1)F(sNi(t+T-1))

• sNi(t) :stateattimet consistingofpositionsandvelocitiesofbird’sneighbors

• Centralized controlifNi includesallbirds• F:fitness function.• T:predictionhorizon.

Page 21: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

ModelPredictiveControlforV-Formation(2)

• Subjecttoconstraints:• Modeldynamics:Stateupdatesofeachbirdare

governedbythemodeldynamics

• Boundedvelocitiesandaccelerations:Thevelocitiesareupper-bounded byaconstant,andtheaccelerationsare

upper-bounded byafactorofthevelocities

• Finally,birdi usestheoptimalaccelerationforbird

itfoundfortime.

Page 22: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

ParticleSwarmOptimization(1)

TheoptimizationproblemissolvedusingPSO• Inspiredbysocialbehaviorofbirdflockingorfishschooling.

• Initializeapopulation(swarm)ofcandidatesolutions(particles)thatmovearoundinthesearch-space.

• Eachparticlekeepstrackofthebestsolutionithasachievedsofar(pbest)andthebestsolutionobtainedsofarbyanyparticleintheneighborsoftheparticle(gbest).

Page 23: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

ParticleSwarmOptimization(2)

• Repeatedlyupdatetheparticle’svelocityandpositionby:vi(t+1)=wvi(t)+c1 r1 (pbesti– xi(t))+c2 r2 (gbesti – xi(t))

xi(t+1)=xi(t)+vi(t+1)

where

w:inertiaweightr1,r2 :randomnumbersin(0,1)sampledeveryiterationc1,c2 :constantlearningfactors

• Terminatewhenmaximumiterationsordesiredfitnesscriteriaisattained.

Page 24: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

DistributedMPCProcedure

Ateverytimestep:• EachbirdlooksatitsneighborsØPlaysseveralscenariosinitsheadtofindthebestconfiguration thattheneighborhoodcanreachin3 steps

ØThebirdthenappliesthefirstmove ofthatsolutiontoupdateitsposition

Inthenexttimestep,eachbirdupdatesitsknowledgeoftheneighbors (positionsandvelocities),whichmaynotbethesameofwhatthatbirdpredictedforitsneighbors

Page 25: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

Outline

• Introduction

• TheV-FormationProblem

• ModelPredictiveControlforV-Formation

• ExperimentalResults

• Conclusions&FutureWork

Page 26: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

ExperimentalResults(1)

Page 27: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

ExperimentalResults(2)

Page 28: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

Outline

• Introduction

• TheV-FormationProblem

• ModelPredictiveControlforV-Formation

• ExperimentalResults

• Conclusions&FutureWork

Page 29: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

Conclusions

• Usedistributedcontrolinsteadofbehavioralrulesto

achieveV-formation.

• IntegrateMPCwithPSOtosolvetheoptimization

problem.

Page 30: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

OngoingandFutureWork

• Deploytheapproachtoactualplants(drones).

• Collisionavoidance.

• ImprovesuccessrateofconvergingtoV-formation.

• UseSMCtoquantifytheprobabilityofsuccess.

• Energyconsumptionandleaderselection.

Page 31: V-Formation as Optimal Controlnavlakha/BDA2016/Material/JunxingYang/... · 2016. 9. 11. · V-Formation as Optimal Control Ashish Tiwari SRI International, Menlo Park, CA, USA BDA,

Thankyou!