MARIO DI BERNARDO · A NSA EL MO DO Numero verd e(va lid oso p r 'It a) 800.422.433 I nf orm azi...

23
M ARIO DI BERNARDO PROFESSOR OF SYSTEMS & CONTROL ENGINEERING UNIVERSITY OF NAPLES FEDERICO II & SSM SCHOOL FOR ADVANCED STUDIES

Transcript of MARIO DI BERNARDO · A NSA EL MO DO Numero verd e(va lid oso p r 'It a) 800.422.433 I nf orm azi...

Page 1: MARIO DI BERNARDO · A NSA EL MO DO Numero verd e(va lid oso p r 'It a) 800.422.433 I nf orm azi eWe b Mo b i le Pro ge t Ed i t o ri al i A rc h iv PR ODO T T I A N SA HOME Ult ima

MARIO DI BERNARDO

PROFESSOR OF SYSTEMS & CONTROL ENGINEERING

UNIVERSITY OF NAPLES FEDERICO II

& SSM – SCHOOL FOR ADVANCED STUDIES

Page 2: MARIO DI BERNARDO · A NSA EL MO DO Numero verd e(va lid oso p r 'It a) 800.422.433 I nf orm azi eWe b Mo b i le Pro ge t Ed i t o ri al i A rc h iv PR ODO T T I A N SA HOME Ult ima

Who am I?

• I am Professor of Systems and

Control Engineering

• My group has 15+ years

expertise in modeling and

controlling large-scale

complex systems

• From drones and robots to

bacteria and cells…

2

Page 3: MARIO DI BERNARDO · A NSA EL MO DO Numero verd e(va lid oso p r 'It a) 800.422.433 I nf orm azi eWe b Mo b i le Pro ge t Ed i t o ri al i A rc h iv PR ODO T T I A N SA HOME Ult ima

A bit of history

• Following the COVID-19 outbreak

in Lombardy, a strict national

lockdown was enforced

in March 2020

• All regions were shut down

concurrently independently of their

epidemic status

• A hot debate ensued

3

Page 4: MARIO DI BERNARDO · A NSA EL MO DO Numero verd e(va lid oso p r 'It a) 800.422.433 I nf orm azi eWe b Mo b i le Pro ge t Ed i t o ri al i A rc h iv PR ODO T T I A N SA HOME Ult ima

What intervention strategies?

4

19th April 2020

Social: Rss FacebookTwitterLinkedinYoutubeInstagramEDIZIONI Mediterraneo Europa-Ue NuovaEuropa America Latina Brasil English ANSAcheck

PRIMOPIANO Ambiente • ANSA2030 • ANSA ViaggiArt • Eccellenze • Industry 4.0 • Legalità • Lifestyle • Mare • Motori • Salute • Scienza • Sisma • Terra&Gusto

•••••

Mappa

Disclaimer

Privacy

Copyright

Modifica

consenso Cookie

Fai di ANSA.it

la tua

homepage

Certificazione ISO 9001. I "processi di

Produzione, distribuzione e pubblicazione in

formato multimediale di notizie giornalistiche"

ANSA sono certificati in conformità alla normativa

internazionale UNI EN ISO 9001:2015.

Politica per la Qualità

informazione pubblicitaria

informazione pubblicitaria

Podcast

Cronaca Politica Economia Mondo Cultura Tecnologia SportRegioni FOTO VIDEO Tutte le sezioni

AGENZIA ANSA - periodicità quotidiana - Iscrizione al Registro della Stampa presso il Tribunale di Roma n. 212/1948

P.I. IT00876481003 - © Copyright ANSA - Tutti i diritti riservati

CANALI ANSA

Mobile RSS

Meteo Cinema

Finanza Codici

Sconto

SERVIZIAZIENDA

ANSA

ANSA NEL MONDO

Numero verde (valido solo per l'Italia)

800.422.433

Informazione Web eMobile

ProgettiEditoriali

Archivi

PRODOTTI ANSA

HOME

• Ultima Ora

• Cronaca

• Politica

• Economia

• Mondo

• Cultura

• Cinema

• Tecnologia

• Sport

• Calcio

• FOTO

• VIDEO

• PODCAST

• Magazine

• Speciali

• Meteo

ECONOMIA

• Borsa

• Industry 4.0

• Professioni

• Real Estate

• PMI

• Ambiente &

Energia

• Motori

• Mare

• Aziende ed

Emergenza

Covid19

REGIONI

• Abruzzo

• Basilicata

• Calabria

• Campania

• Emilia Romagna

• Friuli Venezia Giulia

• Lazio

• Liguria

• Lombardia

• Marche

• Molise

• Piemonte

• Puglia

• Sardegna

• Sicilia

• Toscana

• Trentino-Alto

Adige/Suedtirol

• Umbria

• Valle d'Aosta

• Veneto

MONDO

• Europa

• Nord America

• America

Latina

• Africa

• Medio Oriente

• Asia

• Oceania

• Dalla Cina

• Europa-Ue

CULTURA

• Cinema

• Moda

• Teatro

• TV

• Musica

• Libri

• Arte

• Un Libro al

giorno

• Un Film al

giorno

TECNOLOGIA

• Hi-Tech

• Internet &

Social

• TLC

• Software&App

SPORT

• Calcio

• Formula 1

• Moto

• Golf

• Basket

• Tennis

• Nuoto

• Vela

• Sport Vari

2030

ANSA

2030

AE.

AMBIENTE

&

ENERGIA

IA.

OSSERVATORIO

INTELLIGENZA

ARTIFICIALE

MA.

MARE

ST.

SCIENZA

&

TECNICA

SB.

SALUTE &

BENESSERE

AV.

ANSA

VIAGGIART

M.

MOTORI

TG.

TERRA

&

GUSTO

LS.

LIFESTYLE

L.

LEGALITÀ

&

SCUOLA

4.0

INDUSTRY

4.0

ECCELLENZE

ITALIANE

Fai laricerca

Il mondo inImmagini

Vai allaBorsa

Vai alMeteo

CorporateProdottiPolitica

ANSA.it Politica Coronavirus, De Luca: 'Pronto a chiudere i confini della Campania'

Redazione ANSA

MILANO

18 aprile 2020

14:59

NEWS

Stampa

Scrivi alla redazione

SuggerisciFacebookTwitterLinkedin Mail

Si riunirà in giornata la cabina di regia tra il governo, le Regioni

e i Comuni, convocata dal ministro per gli Affari regionali Francesco

Boccia per coordinare le decisioni e valutare le scelte da compiere in

vista della "fase 2" della ripartenza. Alla riunione dovrebbe prendere

parte il ministro della Salute Roberto Speranza.

"E' assolutamente impossibile oggi dire che il 4 maggio è una data

sicura" per la riapertura di tutte le attività oggi sospese per il

coronavirus, ha detto Vito Crimi, capo politico di M5s a In viva voce,

su Radio Rai. "Credo che dobbiamo aspettare per vedere come

evolve la situazione, dobbiamo sentire il parere del comitato tecnico-

scientifico, il parere delle parti sociali ed economiche, e la politica a

quel punto decide avendo come obiettivo prioritario la tutela della

salute; una riapertura affrettata avrebbe conseguenze

eccessivamente negative".

CAMPANIA - C'è "chi preme per affrettare la ripresa di tutto ma

dobbiamo avere grande senso di responsabilità". Lo ha detto il

presidente della Giunta regionale della Campania, Vincenzo De

Luca, nel corso di una conferenza per fare il punto

sull'emergenza Coronavirus. "Se dovessimo avere corse in avanti

in regioni dove c'è il contagio così forte, la Campania chiuderà i suoi

confini. Faremo una ordinanza per vietare l'ingresso dei cittadini

provenienti da quelle regioni", ha proseguito De Luca.

"Se manteniamo il rigore che abbiamo avuto finora, a metà maggio

potremo dire di aver sconfitto in Campania il Coronavirus", ha detto

De Luca. Ha analizzato i dati del contagio dei primi giorni di aprile

sottolineando che a questi ritmi si "sarà fuori a metà maggio". Nei

prossimi giorni, ha proseguito sempre De Luca, saranno intensificati

i controlli nelle strutture che ospitano gli anziani.

"Stiamo già distribuendo 3 milioni di mascherine perché dal 3/4

maggio l'uso delle mascherine sarà obbligatorio", ha detto ancora.

Le mascherine (in confezione da due), arriveranno, ha spiegato De

Luca, nelle case dei campani, grazie ad una intesa con le Poste

Italiane. Al momento è stata avviata la distribuzione alle categorie a

rischio.

LOMBARDIA - Sulle Rsa "credo proprio che non abbiamo

assolutamente sbagliato niente". Così il governatore Attilio

Fontana al Tgr Lombardia sul ricovero nelle strutture di pazienti

Covid. "Abbiamo fatto una scelta - ha detto Fontana -, l'Ats doveva

controllare le condizioni delle delibera, ossia l'isolamento in singoli

reparti e dipendenti dedicati esclusivamente a quei pazienti, e, sulla

base delle risultanze tecniche, abbiamo portato avanti il

provvedimento". Lo rifareste? "Certamente, in quel periodo

drammatico che stavamo vivendo. Abbiamo liberato posti in

ospedale".

"Nell'ipotesi in cui l'evoluzione del virus dovesse andare in senso

positivo e ci fossero le condizioni, noi il 4 maggio dovremo essere

pronti per la riapertura, purché non prescinda mai dalla sicurezza dei

nostri cittadini e lavoratori", aveva detto il governatore in mattinata in

collegamento con Mattino Cinque. "La condizione ineludibile per

parlare di riapertura è che ci sia il via libera della scienza. Se la

scienza ci dirà bisogna stare chiusi staremo chiusi, però allo stesso

tempo non possiamo farci trovare impreparati", ha aggiunto

Fontana.

Per quanto riguarda il capitolo delle inchieste sulle Rsa, le Case di

riposo e degenza per anziani, dove il coronavirus, anche a detta

dell'Oms ha compiuto "un massacro", Fontana, a Mattino5, risponde:

"Io aspetto con estrema serenità l'esito". E alla domanda 'se pensa

che dalle indagini della magistratura sui contagi e i decessi per

Coronavirus nella case di riposo possano emergere 'errori gravi' da

parte della giunta lombarda', ha risposto: "Noi abbiamo fatto una

delibera che è stata proposta dai nostri tecnici. Sono stati i nostri

esperti che ci hanno detto che a determinate condizioni - e cioè che

esistessero dei reparti assolutamente isolati dal resto della struttura

e addetti dedicati esclusivamente ai malati Covid - la cosa si poteva

fare".

"Noi abbiamo fatto questa proposta - ha proseguito - e le case di

riposo che avevano queste condizioni hanno aderito. Non bisogna

dimenticare che la scelta è stata fatta perché non avevamo più posti

negli ospedali per ricoverare la gente che non poteva più essere

curata a casa. Però sono stati i nostri tecnici che ci hanno fatto la

proposta, che hanno valutato le condizioni delle singole case di

riposo e noi ci siamo adeguati".

La Lombardia "in questo momento ha una condizione

epidemiologica di particolare gravità - ha detto Walter Ricciardi,

rappresentante italiano nell'executive board dell'Oms e consulente

del ministro della Salute al giornale ilcaffeonline.it. -. Certamente in

miglioramento ma di particolare gravità. Tra i paesi europei è la

regione che in questo momento ha maggiori problemi. Non l'Italia,

ma la Lombardia. Mi sembra non saggio immaginare delle aperture

a prescindere da valutazioni oggettive delle condizioni

epidemiologiche".

VENETO - "Il trend positivo sui positivi ci fa preparare bene alla

ripartenza: il lockdown non esiste più", ha detto il governatore del

Veneto Luca Zaia. "Stiamo lavorando alacremente - spiega - per

presentare un progetto di messa in sicurezza per la ripartenza". "Se

dipendesse da me riaprirei tutto il 4 maggio con gradualità e senso

di responsabilità".

EMILIA ROMAGNA - Per quanto riguarda la 'Fase 2' il presidente

della regione Emilia Romagna Stefano Bonaccini rinvia invece al

governo l'ultima parola. "È il governo che ha l'ultima parola sulla

ripartenza delle fabbriche. Le regioni possono solo fare proposte

sulle filiere strategiche nel loro territorio", afferma in un'interbista a

'La Repubblica'. Sulla ripresa "bisogna continuare a gestire

un'emergenza sanitaria che resta difficile. Non ho intenzione di fare

polemiche, voglio invece dare una mano per arrivare presto a un

Piano per il Paese. Di cui c'è bisogno". Su una possibile ripartenza

in base alle aree geografiche, aggiunge: "Penso non sia il momento

delle divisioni. Serve unità. Né mi permetto io di indicare quale sia la

strada giusta per le altre regioni. Il governo ha giustamente l'ultima

parola. Per parte nostra vogliamo dare una mano a definire un Piano

Paese: costruire accordi territoriali per una ripartenza sicura. In

Emilia Romagna abbiamo deciso con le parti sociali che metteremo

a punto un progetto per far ripartire gradualmente le filiere a valenza

internazionale e i cantieri delle opere pubbliche".

TRENTINO ALTO ADIGE - Il governatore altoatesino Arno

Kompatscher ha chiesto, durante una videoconferenza delle

Province e Regioni autonome con il viceministro all'Economia Laura

Castelli, la sospensione del contributo della Provincia autonoma al

risanamento del disavanzo statale. Si tratta di un contributo annuo di

476 milioni di euro. Kompatscher ha giustificato questa richiesta con

il mancato gettito fiscale causa emergenza coronavirus. "Per

garantire servizi essenziali, prevedere bond regionali", ha aggiunto.

Salute, Sicurezza sul lavoro Attilio Fontana

RIPRODUZIONE RISERVATA © Copyright ANSA

CONDIVIDI

Notizie Correlate

Sindaci M5s al governo, ora serve il 'Cura Comuni'

In una residenza psichiatrica ligure 3 morti e 38 positivi su

40

Coronavirus, Pressing di quattro Regioni. 'Riapertura il 4

maggio'

Video

Coronavirus, De Luca: "Pronto a chiudere i

confini della Campania"

Video

Coronavirus, Zaia: "Per noi il lockdown

non esiste piu'"

Coronavirus, De Luca: 'Pronto a chiuderei confini della Campania''Se c'è corsa ad aprire tutto in regioni ancora esposte'. Fontana: 'Sulle Rsa non

abbiamo sbagliato niente, liberati posti in ospedale'

informazione pubblicitaria

informazione pubblicitaria

informazione pubblicitaria

informazione pubblicitaria

+CLICCA PER

INGRANDIREVincenzo De Luca © ANSA/

Suggerisci

Facebook

Twitter

Altri

TI POTREBBERO INTERESSARE ANCHE:

Azioni Amazon: con soli 200€ puoi

ottenere una rendita mensile fissa.

Scoprilo!Markets Guide

AD

Nuova CLIO ZEN da 169€/mese.

Blocca l’offerta in VIDEO LIVE

CHAT.Renault

AD

Il gioco City Builder "da

giocare". Nessuna…Forge Of Empires

AD

Noa si è lasciata morire.

Papa, una sconfitta per tutti…

Queste sono le case di riposo

a Naples dove i tuoi

genitori…Casa di riposo - Sponsored…

AD

Artista fa causa a Federica

Panicucci

I prezzi delle auto a noleggio

potrebbero stupirticonfronto autonoleggio | sponsored listing

AD

Offerte energia elettrica - quale

conviene? Confrontale oraOfferte energia elettrica | Link…

AD

Burioni: "Se il coronavirus ha 3% di

mortalità ed è così diffusa è una…

Heather Parisi, per 7 anni

picchiata con calci e pugni …

Assicurazione auto senior: il

costo potrebbe sorprendertiCar Insurance For Seniors |…

AD

Suora accusa dolori, in

ospedale scopre di essere…

[Rapporto Forrester] Scopri i

vantaggi dell'utilizzo delle…IBM

AD

Sterminò famiglia, muore incidente

auto

VIDEO ANSA

21 APRILE, 18:02

MACRON TELEFONA AL PAPA: "LA

ASPETTIAMO IN FRANCIA"

21 aprile, 17:57

Concluso il varo della

penultima campata del

ponte di Genova

21 aprile, 17:24

Harry e Meghan rompono

i ponti coi tabloid: "Basta

bugie"

› tutti i video

17:50

17:38

15:58

15:56

15:48

15:42

15:32

15:29

15:29

15:19

Raddoppio pensione invalidità 2020

Venezia: Barbera, così preparo la Mostra

Conte, per nuovo dl non meno di 50 mld

Conte, no a compromessi al ribasso

Conte, su Mes ultima parola a Parlamento

Conte, aperti a opposizione responsabile

Per l'Earth Day 7 ore live streaming con

National Geographic

Conte, modifiche anche su distanziamento

Conte, riaperture omogenee in Italia

Conte, mascherine finché manca vaccino

› Tutte le news

ULTIMA ORA HOME

SkyTG24, 17th April 2020

Il Fatto Quotidiano, 20th April 2020

Corriere, 21st April 2020

Page 5: MARIO DI BERNARDO · A NSA EL MO DO Numero verd e(va lid oso p r 'It a) 800.422.433 I nf orm azi eWe b Mo b i le Pro ge t Ed i t o ri al i A rc h iv PR ODO T T I A N SA HOME Ult ima

5

Our research questions

9th October 2020

Page 6: MARIO DI BERNARDO · A NSA EL MO DO Numero verd e(va lid oso p r 'It a) 800.422.433 I nf orm azi eWe b Mo b i le Pro ge t Ed i t o ri al i A rc h iv PR ODO T T I A N SA HOME Ult ima

Our research questions

• Do we need to shut down all regions

concurrently independently from their

epidemic status?

• Can we envisage regional lockdown

strategies to mitigate the COVID19

epidemic…

• … where each region automatically

enforces or releases social distancing

and other measures according to some

indicators?

• Are such regional measures effective and

do they cause lesser economic costs?

6

Page 7: MARIO DI BERNARDO · A NSA EL MO DO Numero verd e(va lid oso p r 'It a) 800.422.433 I nf orm azi eWe b Mo b i le Pro ge t Ed i t o ri al i A rc h iv PR ODO T T I A N SA HOME Ult ima

Modeling Italy as a complex system

Figure 1: Network structure of the epidemic model representing the flux f i j between neighboring regions, each

represented by a vertex of the graph.

As done in [?], we focus on equations (1a-1b) and start by recasting them in vector notation. Namely

we define S( t) = [S1 S2 . . . SM]T and I ( t) = [ I1 I2 . . . IM]T . Then, to find the basic reproduction numbers

(regional and national), we study local stability of the equilibrium point (S∗ , I ∗ ) with S∗ = [N1 N2 . . . NM]T

and I∗ = [0 0 . . . 0]T corresponding to an healthy population in all nodes. Specifically, letting dS∶= S− S∗ and

dI ∶= I − I ∗ be small perturbations away from the equilibrium of interest, linearizing we obtain:

dSi =∂ j i

S

∂Si

dSi +∂ j i

S

∂ I1dI1 +

∂ j iS

∂ I2dI2 + . . .+

∂ j iS

∂ IMdIM, (8)

d Ii =∂ j i

I

∂Si

dSi +∂ j i

I

∂ I1dI1 +

∂ j iI

∂ I2dI2 + . . .+

∂ j iI

∂ IMdIM, (9)

where

j iS = −

M

j=1

M

k=1

r j( t)bf i j ( t)Si

f kj ( t) Ik

Npj

j iI =

M

j=1

M

k=1

r j ( t)bf i j ( t)Si

f kj ( t) Ik

Npj

− (a + g) Ii

and, therefore,∂ j i

S

∂Si

=M

j=1

M

k=1

r jbf i j f kj

I ∗kN

pj

= 0 (10)

while, for q ≠ i,

∂ j iI

∂ Iq=

M

j=1

M

k=1

r jbf i jS∗i

f qj

Npj

=M

j=1

r jbf i j f qjNi

Npj

(11)

and, for q = i,

∂ j iI

∂ Ii=

M

j=1

M

k=1

r jbf 2i j

S∗iN

pj

− (a + g) =M

j=1

r jbf 2i j

Ni

Npj

− (a +g) (12)

4

7

Page 8: MARIO DI BERNARDO · A NSA EL MO DO Numero verd e(va lid oso p r 'It a) 800.422.433 I nf orm azi eWe b Mo b i le Pro ge t Ed i t o ri al i A rc h iv PR ODO T T I A N SA HOME Ult ima

A data-driven model of each region

• We derived a mathematical model of the

epidemic spread in each region…

• …and used data to calibrate and estimate all

model parameters and unknowns

8

Page 9: MARIO DI BERNARDO · A NSA EL MO DO Numero verd e(va lid oso p r 'It a) 800.422.433 I nf orm azi eWe b Mo b i le Pro ge t Ed i t o ri al i A rc h iv PR ODO T T I A N SA HOME Ult ima

Constructing a network model

Figure 1: Network structure of the epidemic model representing the flux f i j between neighboring regions, each

represented by a vertex of the graph.

As done in [?], we focus on equations (1a-1b) and start by recasting them in vector notation. Namely

we define S( t) = [S1 S2 . . . SM]T and I ( t) = [ I1 I2 . . . IM]T . Then, to find the basic reproduction numbers

(regional and national), we study local stability of the equilibrium point (S∗ , I ∗ ) with S∗ = [N1 N2 . . . NM]T

and I∗ = [0 0 . . . 0]T corresponding to an healthy population in all nodes. Specifically, letting dS∶= S− S∗ and

dI ∶= I − I ∗ be small perturbations away from the equilibrium of interest, linearizing we obtain:

dSi =∂ j i

S

∂Si

dSi +∂ j i

S

∂ I1dI1 +

∂ j iS

∂ I2dI2 + . . .+

∂ j iS

∂ IMdIM, (8)

d Ii =∂ j i

I

∂Si

dSi +∂ j i

I

∂ I1dI1 +

∂ j iI

∂ I2dI2 + . . .+

∂ j iI

∂ IMdIM, (9)

where

j iS = −

M

j=1

M

k=1

r j( t)bf i j ( t)Si

f kj ( t) Ik

Npj

j iI =

M

j=1

M

k=1

r j ( t)bf i j ( t)Si

f kj ( t) Ik

Npj

− (a + g) Ii

and, therefore,∂ j i

S

∂Si

=M

j=1

M

k=1

r jbf i j f kj

I ∗kN

pj

= 0 (10)

while, for q ≠ i,

∂ j iI

∂ Iq=

M

j=1

M

k=1

r jbf i jS∗i

f qj

Npj

=M

j=1

r jbf i j f qjNi

Npj

(11)

and, for q = i,

∂ j iI

∂ Ii=

M

j=1

M

k=1

r jbf 2i j

S∗iN

pj

− (a + g) =M

j=1

r jbf 2i j

Ni

Npj

− (a +g) (12)

4

• Once each region was

modelled, we introduced

fluxes among regions and

assessed their effects

9

Page 10: MARIO DI BERNARDO · A NSA EL MO DO Numero verd e(va lid oso p r 'It a) 800.422.433 I nf orm azi eWe b Mo b i le Pro ge t Ed i t o ri al i A rc h iv PR ODO T T I A N SA HOME Ult ima

Estimating the fluxes

10

We estimated fluxes by using

publicly available data on trains, air

traffic and capacity, key ferries routes

Page 11: MARIO DI BERNARDO · A NSA EL MO DO Numero verd e(va lid oso p r 'It a) 800.422.433 I nf orm azi eWe b Mo b i le Pro ge t Ed i t o ri al i A rc h iv PR ODO T T I A N SA HOME Ult ima

Identification method

11

Constrained Nonlinear

Optimization

fmincon (Optimization Toolbox)

Page 12: MARIO DI BERNARDO · A NSA EL MO DO Numero verd e(va lid oso p r 'It a) 800.422.433 I nf orm azi eWe b Mo b i le Pro ge t Ed i t o ri al i A rc h iv PR ODO T T I A N SA HOME Ult ima

Breakpoint identification

• The procedure self-determines the breakpoints

where notable parameter changes are observed

12

parameters

change

Breakpoints identified using

Chow statistical test

fitted

Parameters Change

(Econometrics Toolbox)

Page 13: MARIO DI BERNARDO · A NSA EL MO DO Numero verd e(va lid oso p r 'It a) 800.422.433 I nf orm azi eWe b Mo b i le Pro ge t Ed i t o ri al i A rc h iv PR ODO T T I A N SA HOME Ult ima

Identification of the remaining parameters

13

+

Constrained Least-

squares identification

lsqlin (Optimization Toolbox)

Page 14: MARIO DI BERNARDO · A NSA EL MO DO Numero verd e(va lid oso p r 'It a) 800.422.433 I nf orm azi eWe b Mo b i le Pro ge t Ed i t o ri al i A rc h iv PR ODO T T I A N SA HOME Ult ima

Example: Lombardy

𝑯%

0,15

1,23

1,21

0,63

𝝆𝜷 𝜶 𝝍 𝜻 𝑻

0,499 0,020 0,050 0,026 27 Feb – 2 Mar

0,268 0,038 0,032 0,034 2 Mar – 21 Mar

0,081 0,043 0,027 0,030 21 Mar – 8 Apr

0,047 0,052 0,018 0,019 13 Apr - 21 Apr

27/02/2020

RED ZONE

Page 15: MARIO DI BERNARDO · A NSA EL MO DO Numero verd e(va lid oso p r 'It a) 800.422.433 I nf orm azi eWe b Mo b i le Pro ge t Ed i t o ri al i A rc h iv PR ODO T T I A N SA HOME Ult ima

Example: Lombardy

𝝆𝜷 𝜶 𝝍 𝜻 𝑻

0,499 0,020 0,050 0,026 27 Feb – 2 Mar

0,268 0,038 0,032 0,034 2 Mar – 21 Mar

0,081 0,043 0,027 0,030 21 Mar – 8 Apr

0,047 0,052 0,018 0,019 13 Apr - 22 Apr

𝑯%

0,15

1,23

1,21

0,63

8/03/2020

Page 16: MARIO DI BERNARDO · A NSA EL MO DO Numero verd e(va lid oso p r 'It a) 800.422.433 I nf orm azi eWe b Mo b i le Pro ge t Ed i t o ri al i A rc h iv PR ODO T T I A N SA HOME Ult ima

Regional mitigation strategies

• Goal: introduce automatic feedback

measures to avoid saturation of the

regional health services while mitigating

the epidemic spread

• When its regional health system becomes

under pressure, a region enforces a

lockdown and prevents travel to/from

other regions

• e.g. when the number of patients

requiring ICU treatment crosses a certain

threshold

16

Page 17: MARIO DI BERNARDO · A NSA EL MO DO Numero verd e(va lid oso p r 'It a) 800.422.433 I nf orm azi eWe b Mo b i le Pro ge t Ed i t o ri al i A rc h iv PR ODO T T I A N SA HOME Ult ima

How we implemented our solutions

• We tested the robustness of our

controller via a Monte Carlo

simulation campaign

• All parameters' sets are generated

using a Latin hypercube sampling

17

lhsdesign (Statistics and Machine LearningToolbox)

Widely scattered

parameters

Page 18: MARIO DI BERNARDO · A NSA EL MO DO Numero verd e(va lid oso p r 'It a) 800.422.433 I nf orm azi eWe b Mo b i le Pro ge t Ed i t o ri al i A rc h iv PR ODO T T I A N SA HOME Ult ima

MATLAB was fundamental

18

Analyze the results

Run Monte Carlo

Simulations

Generatethe set of

Parameters

lhsdesign function

from Statistics and

Machine Learning

Toolbox.

parfor loop from

Parallel Computing

Toolbox

computing

control metrics

using MATLAB

Page 19: MARIO DI BERNARDO · A NSA EL MO DO Numero verd e(va lid oso p r 'It a) 800.422.433 I nf orm azi eWe b Mo b i le Pro ge t Ed i t o ri al i A rc h iv PR ODO T T I A N SA HOME Ult ima

A regional feedback strategy

Page 20: MARIO DI BERNARDO · A NSA EL MO DO Numero verd e(va lid oso p r 'It a) 800.422.433 I nf orm azi eWe b Mo b i le Pro ge t Ed i t o ri al i A rc h iv PR ODO T T I A N SA HOME Ult ima
Page 21: MARIO DI BERNARDO · A NSA EL MO DO Numero verd e(va lid oso p r 'It a) 800.422.433 I nf orm azi eWe b Mo b i le Pro ge t Ed i t o ri al i A rc h iv PR ODO T T I A N SA HOME Ult ima

Comparison with nationwide lockdown

Page 22: MARIO DI BERNARDO · A NSA EL MO DO Numero verd e(va lid oso p r 'It a) 800.422.433 I nf orm azi eWe b Mo b i le Pro ge t Ed i t o ri al i A rc h iv PR ODO T T I A N SA HOME Ult ima

Take home message

• Our work convincingly showed the benefits of regional rather than national NPIs to deal with the pandemic

• This is precisely the idea behind the strategy adopted in Italy since 3rd November 2020

• The use of the optimization algorithms provided by MATLAB Optimization Toolbox allowed us to quickly test our hypothesis with minimum effort

• The tools from MATLAB Parallel Computing Toolbox allowed us to easily write optimized code to be launched in parallel reducing noticeably the execution time

Page 23: MARIO DI BERNARDO · A NSA EL MO DO Numero verd e(va lid oso p r 'It a) 800.422.433 I nf orm azi eWe b Mo b i le Pro ge t Ed i t o ri al i A rc h iv PR ODO T T I A N SA HOME Ult ima

https://sites.google.com/site/dibernardogroup/ @mdiberna