Legge di potenza
# Legge di potenza

Una legge di potenza (power law) è una qualsiasi relazione del tipo:

    $\ f(x)=ax^{k}+o(x^{k}),\ x\to 0,$

dove a e k sono costanti e $o(x^{k})$ è una funzione asintoticamente piccola di $x^{k}$. k è di solito chiamato esponente di scala.

Leggi di potenza ricorrono nelle distribuzioni di probabilità di molti fenomeni fisici (ad esempio la magnitudo dei terremoti, il diametro dei crateri dei pianeti, la dimensione dei frammenti degli oggetti che si infrangono per urti, l'intensità delle esplosioni solari), sociali (il numero dei morti nelle guerre, la popolazione delle città, il numero di collegamenti ai siti web, il numero di citazioni) ed economici (la distribuzione della ricchezza, le vendite di libri e cd, ecc.); così come ricorrono in altri tipi di relazioni, come quella tra il tasso metabolico di una specie e la sua massa corporea (cosiddetta legge di Kleiber), o quella tra la forza di gravità e la distanza tra le masse.

Nel caso delle distribuzioni di probabilità, una distribuzione che obbedisce alla legge di potenza è denominata power law distribution, scale-free distribution (distribuzione a invarianza di scala), o anche distribuzione di Pareto - dal nome dell'economista Vilfredo Pareto, che per primo la individuò nella distribuzione del reddito - o infine legge di Zipf - dal linguista George Kingsley Zipf che la individuò studiando la frequenza d'uso delle parole nei testi.

La particolarità di questo tipo di distribuzioni sta proprio nell'assenza di una scala caratteristica dei fenomeni. Così, ad esempio, dire che la distribuzione del reddito rispecchia la legge di potenza, o che la distribuzione del reddito è una paretiana, significa dire che, se ogni quattro individui con reddito annuo pari a diecimila euro, ne esiste uno con reddito pari a ventimila, allora ci sarà una persona che guadagna 2 trilioni di euro per ogni quattro con reddito pari a 1 trilione.

## Distribuzioni a legge di potenza

Una distribuzione a legge di potenza, o distribuzione power-law, o distribuzione di Pareto, nella sua forma più generale ha la forma:

    $\ p(x)\propto L(x)x^{-\alpha }$

dove il simbolo $\propto$ indica la proporzionalità, ovvero l'eguaglianza a meno di un fattore moltiplicativo. α (> 1) è chiamato esponente della power-law e L(x) è una funzione "che varia lentamente", cioè qualsiasi funzione tale per cui

    $\ \lim _{x\rightarrow \infty }{\frac {L(t\,x)}{L(x)}}=1$

con t costante. Questa proprietà di L(x) segue direttamente dalla condizione che p(x) sia asintoticamente invariante alla scala.

Quando L(x) è una costante, la distribuzione diventa:

    $\ p(x)=Cx^{-\alpha }$

Inoltre, poiché per qualsiasi valore positivo di α la distribuzione diverge quando x tende a zero, è normale imporre un valore minimo $x_{\min }$. Dato l'esponente, quando x è una variabile continua si ha che la costante di normalizzazione C è data da:

    $\ C=1/\int _{x_{\min }}^{\infty }x^{-\alpha }\,\mathrm {d} x$

da cui:

    $p(x)={\frac {\alpha -1}{x_{\min }}}\left({\frac {x}{x_{\min }}}\right)^{-\alpha }$.

### Momenti della distribuzione

Il valore atteso della distribuzione a legge di potenza è:

    $\langle x\rangle =\int _{x_{\min }}^{\infty }xp(x)\mathrm {d} x={\frac {C}{2-\alpha }}\left[x^{2-\alpha }\right]_{x_{\min }}^{\infty }$

che è finito solo se α > 2.

Più in generale, il momento di ordine m è dato da:

    $\langle x^{m}\rangle =\int _{x_{\min }}^{\infty }x^{m}p(x)\,\mathrm {d} x={\frac {\alpha -1}{\alpha -1-m}}x_{\min }^{m}$

ed è finito solo se m < α -1.

### Rappresentare graficamente le distribuzioni a legge di potenza

Di solito le distribuzioni power-law sono rappresentate su un grafico log-log, cioè un grafico in cui entrambe le variabili sugli assi sono misurate in logaritmi. Trasformando in logaritmi si ha infatti:

    $\ \log p(x)=\log C-\alpha \log x$

e la relazione diventa lineare.

Inoltre, per evitare le fluttuazioni frazionarie nella coda della distribuzione dovute alla costanza dell'intervallo di classe nell'istogramma, si utilizza un intervallo di classe logaritmico, o, più frequentemente, si rappresenta la funzione di probabilità cumulata complementare, cioè:

    $P(x)=\Pr(X>x)=C\int _{x}^{\infty }p(X)\mathrm {d} X={\frac {\alpha -1}{x_{\min }^{-\alpha +1}}}\int _{x}^{\infty }X^{-\alpha }\mathrm {d} X=\left({\frac {x}{x_{\min }}}\right)^{(-\alpha +1)}$

che segue anch'essa una legge di potenza con esponente più piccolo.

### Percentili della distribuzione

La mediana della distribuzione è data da:

    $\ \int _{x_{1/2}}^{\infty }p(x)\mathrm {d} x={\frac {1}{2}}$

da cui:

    $\ x_{1/2}=2^{\frac {1}{\alpha -1}}x_{\min }$

Il k-esimo percentile è:

    $\ \int _{x_{k/100}}^{\infty }p(x)\mathrm {d} x={\frac {100-k}{100}}$

ovvero:

    $\ x_{k/100}=\left({\frac {100}{100-k}}\right)^{\frac {1}{\alpha -1}}x_{\min }$

### La regola 80/20

Se si considera la distribuzione della ricchezza, la frazione della popolazione la cui ricchezza eccede un determinato valore x è semplicemente:

    $\ P(x)=\left({\frac {x}{x_{\min }}}\right)^{1-\alpha }$

La percentuale della ricchezza totale in mano a queste persone sarà:

    $\ W(x)={\frac {\int _{x}^{\infty }x'p(x')\mathrm {d} x'}{\int _{x_{\min }}^{\infty }x'p(x')\mathrm {d} x'}}=\left({\frac {x}{x_{\min }}}\right)^{2-\alpha }$

da cui segue:

    $\ W=P^{\frac {\alpha -2}{\alpha -1}}$

Poiché l'esponente nel caso della distribuzione della ricchezza è circa 2,15, si ha:

    $\ W=0.2^{\frac {2,15-2}{2,15-1}}\approx 0,8$

Che è nota come legge 80/20 (o principio di Pareto), in base alla quale il 20% della popolazione detiene l'80% della ricchezza mondiale.

Va notato come, in base alla stessa distribuzione, segue come, ad esempio, il 50% della ricchezza mondiale sia detenuto solo da meno dell'1% della popolazione, e così via.

## Invarianza di scala e legge di potenza

Una distribuzione di probabilità è invariante alla scala se:

    $\ p(bx)=f(b)\ p(x)$

cioè se, cambiando la scala o l'unità di misura della variabile x di un fattore b, la distribuzione di probabilità rimane invariata a meno di una costante moltiplicativa f(b).

L'invarianza di scala è condizione necessaria e sufficiente affinché una distribuzione sia power-law.

Infatti, una distribuzione è invariante alla scala se è una distribuzione power-law. Dato

    $\ p(x)=Cx^{-\alpha }$

segue infatti:

    $\ p(bx)=C(bx)^{-\alpha }=b^{-\alpha }Cx^{-\alpha }=b^{-\alpha }p(x).$

Inoltre, una distribuzione è power-law se è invariante alla scala. Infatti, se una distribuzione è scale-free, per x = 1 si ha:

    $\ f(b)={\frac {p(b)}{p(1)}}$

e quindi:

    $\ p(bx)={\frac {p(b)}{p(1)}}p(x)$

Differenziando rispetto a b otteniamo:

    $\ xp'(bx)={\frac {p'(b)}{p(1)}}p(x)$

che, per b = 1 diventa:

    $\ xp'(x)={\frac {p'(1)}{p(1)}}p(x)$

che è un'equazione differenziale la cui soluzione è:

    $\ p(x)=p(1)\ x^{p(1)/p'(1)}$

che è una distribuzione power-law.

## Invarianza di scala e cutoff

Nonostante siano molti i fenomeni che presentano distribuzioni ad invarianza di scala per alcuni intervalli, sono rari i casi in cui questo vale lungo tutto il supporto.

Così, ad esempio, nonostante la legge Gutenberg-Richter venga di solito citata come un esempio di distribuzione a legge di potenza, la distribuzione reale della magnitudo dei terremoti, visto il limite costituito dall'energia totale racchiusa nella crosta terrestre, smette di essere scalabile quando ci si avvicina ad esso.

Per tenere conto di questo, si usa spesso introdurre un cutoff esponenziale nella distribuzione originale, adottando la forma seguente:

    $p(x)\propto x^{-\alpha }e^{-\lambda x}.$

In questa distribuzione, il saggio di decadimento esponenziale $\ e^{-\lambda x}$ alla fine tenderà a prevalere sull'invarianza di scala per valori via via più grandi di x.

## Distribuzioni lognormali e a legge di potenza

Data una variabile casuale logonormale, la cui funzione di probabilità è:

    $p(x)={\frac {1}{x}}\exp \left(-{\frac {\left(\ln x-\mu \right)^{2}}{2\sigma ^{2}}}\right)$

prendendo il logaritmo si ottiene:

    $\ln p(x)=-{\frac {1}{2\sigma ^{2}}}(\ln x)^{2}+\left({\frac {\mu }{\sigma ^{2}}}-1\right)\ln x-{\frac {\mu ^{2}}{2\sigma ^{2}}}$

Questa è una funzione quadratica di ln x, che quindi presenta una curvatura verso il basso. Tuttavia, può approssimarsi ad una linea per piccoli tratti e l'approssimazione è tanto migliore quanto più grande è la deviazione standard.

Data questa forma, di fatto risulta molto difficile in pratica distinguere una power-law con un cut-off da una lognormale.

