math – Programmare con Python

Fattoriale

04/03/202304/03/2023 di admin

Puoi usare la funzione math.factorial() Definizione iterativa Per n=5: 5! = 5*4*3*2*1 = 120 Il fattoriale come funzione di n Con for con while utilizzo

def fatt(n):

prod=1

for i in range(2,n+1):

prod=i*prod

return prod

def fatt(n):

prod=1

i=2

while(i <= n):

prod=i*prod

i =i+1

return prod

n=5

f5=fatt(n)

print(n, f5) # 5 120

Definizione ricorsiva oppure Per N=5

5!=5*4!

4!=4*3!

3!=3*2!

2!=2*1!

1!=1

2!=2*1=2

3!=3*2=6

4!=4*6=24

5!=5*24=120

Codifica

def fatt(n):

if(n == 0) or (n == 1):

return 1

else:

return n*fatt(n-1)

Approssimazione All’aumentare del valore di n il valore del fattoriale diventa astronomico! Per esempio, . Se è sufficiente un … Leggi tutto

Angoli e funzioni

30/08/202228/08/2022 di admin

Visualizza la tabella della funzione seno e coseno con angolo in [0°, 360°) e passo=15°

import math

import prettytable as pt

t=pt.PrettyTable()

t.field_names = ["Gradi", "Radianti", "Seno", "Coseno"] # Etichette

t.align = "r" # Tutti i campi allineati a destra

t.float_format = ".10" # 10 cifre decimali

START = 0

STOP = 360

STEP = 15

for gradi in range(START, STOP, STEP):

radianti = math.radians(gradi)

seno = math.sin(radianti)

coseno = math.cos(radianti)

t.add_row([gradi, radianti, seno, coseno])

print(t)

Con il doppio bordo? Aggiungi la riga di codice

1	t.set_style(pt.DOUBLE_BORDER)

Vedi la discussione teorica Utilizza i moduli per fractions, mantenere formali i risultati delle operazioni con le frazioni matplotlib, rappresentare graficamente la distribuzione delle probabilità math, radice quadrata

from fractions import Fraction as F # Fraction()

import matplotlib.pyplot as plt # plot()

import math # sqrt()

N = 13*[0] # Conteggio degli esiti del lancio di due dadi

for d1 in range(1, 7):

for d2 in range(1,7):

dadi=d1+d2

N[dadi] += 1

X = list(range(2, 13))

P = [] # Probabilità: frequenze assolute

for x in X:

p=F(N[x], 36)

P.append(p)

plt.bar(X, P) # Grafico a barre verticali

plt.grid()

plt.show()

print(" x | p(x)") # Tabella

print("----|------")

for x, p in zip(X, P):

print(" %2i | %s" %(x, p))

Risultato

x | p(x)

----|------

2 | 1/36

3 | 1/18

4 | 1/12

5 | 1/9

6 | 5/36

7 | 1/6

8 | 5/36

9 | 1/9

10 | 1/12

11 | 1/18

12 | 1/36

Il resto del codice esegue i calcoli

media = 0

media_quadrati = 0

for p, x in zip(P, X):

media += p*x

media_quadrati += p*x**2

varianza_1 = media_quadrati-media**2

print("Media = %5s %10.6f" %(media , media ))

print("Media dei quadrati = %5s %10.6f" %(media_quadrati, media_quadrati))

print("Varianza (1) = %5s %10.6f" %(varianza_1 , varianza_1 ))

scarto_medio_assoluto = 0

devianza = 0

varianza_2 = 0

for p, x in zip(P, X):

scarto_medio_assoluto += p*abs(media-x)

devianza += (media-x)**2

varianza_2 += p*(media-x)**2

print("Scarto medio assoluto = %5s %10.6f" %(scarto_medio_assoluto, scarto_medio_assoluto))

print("Devianza = %5s %10.6f" %(devianza , devianza ))

print("Varianza (2) = %5s %10.6f" %(varianza_2 , varianza_2 ))

deviazione_standard = math.sqrt(varianza_2)

deviazione_standard_relativa = deviazione_standard/abs(media)

print("Deviazione standard =", deviazione_standard )

print("Deviazione standard rel. =", deviazione_standard_relativa)

Risultati

Media = 7 7.000000

Media dei quadrati = 329/6 54.833333

Varianza (1) = 35/6 5.833333

Scarto medio assoluto = 35/18 1.944444

Devianza = 110 110.000000

Varianza (2) = 35/6 5.833333

Deviazione standard = 2.41522945769824

Deviazione standard rel. = 0.3450327796711771

Riepilogo collections SymPy I moduli … permettono di … collections, conteggiare gli esiti di … Leggi tutto

Analisi di tre dadi

10/10/202209/07/2022 di admin

Vedi la discussione teorica la codifica per due dadi Modifica le 2 righe di codice iniziali

...

DADI = [d1+d2+d3 for d1 in DADO for d2 in DADO for d3 in DADO] # Lancio di 3 dadi

...

P = [sp.Rational(count, 216) for count in COUNT.values()] # Probabilità = Frequenze assolute

...

Il resto del codice rimane lo stesso

import collections # Counter()

import matplotlib.pyplot as plt

import sympy as sp # Rational(), sqrt()

DADO = list(range(1, 7))

DADI = [d1+d2+d3 for d1 in DADO for d2 in DADO for d3 in DADO]

COUNT = collections.Counter(DADI);

X = COUNT.keys()

P = [sp.Rational(count, 216) for count in COUNT.values()]

plt.bar(X, P)

plt.grid()

plt.show()

print(" x | p(x)")

print("----|------")

for x, p in zip(X, P):

print(" %2i | %s" %(x, p))

#-------------------------------------------------------------------------------

media = 0

media_quadrati = 0

for p, x in zip(P, X):

media += p*x

media_quadrati += p*x**2

varianza_1 = media_quadrati-media**2

print("Media = %5s %10.6f" %(media , media ))

print("Media dei quadrati = %5s %10.6f" %(media_quadrati, media_quadrati))

print("Varianza (1) = %5s %10.6f" %(varianza_1 , varianza_1 ))

#-------------------------------------------------------------------------------

scarto_medio_assoluto = 0

devianza = 0

varianza_2 = 0

for p, x in zip(P, X):

scarto_medio_assoluto += p*abs(media-x)

devianza += (media-x)**2

varianza_2 += p*(media-x)**2

print("Scarto medio assoluto = %5s %10.6f" %(scarto_medio_assoluto, scarto_medio_assoluto))

print("Devianza = %5s %10.6f" %(devianza , devianza ))

print("Varianza (2) = %5s %10.6f" %(varianza_2 , varianza_2 ))

#-------------------------------------------------------------------------------

deviazione_standard = sp.sqrt(varianza_2)

deviazione_standard_relativa = deviazione_standard/abs(media)

print("Deviazione standard = %11s %10.6f" %(deviazione_standard , deviazione_standard ))

print("Deviazione standard rel. = %11s %10.6f" %(deviazione_standard_relativa, deviazione_standard_relativa))

Per ottenere

x | p(x)

----|------

3 | 1/216

4 | 1/72

5 | 1/36

6 | 5/108

7 | 5/72

8 | 7/72

9 | 25/216

10 | 1/8

11 | 1/8

12 | 25/216

13 | 7/72

14 | 5/72

15 | 5/108

16 | 1/36

17 | 1/72

18 | 1/216

Media = 21/2 10.500000

Media dei quadrati = 119 119.000000

Varianza (1) = 35/4 8.750000

Scarto medio assoluto = 29/12 2.416667

Devianza = 340 340.000000

Varianza (2) = 35/4 8.750000

Deviazione standard = sqrt(35)/2 2.958040

Deviazione standard rel. = sqrt(35)/21 0.281718

Riepilogo

Analisi di un dado

12/07/202208/07/2022 di admin

Vedi la discussione Le sei facce di un dado sono equiprobabili Sia X la variabile casuale “punti realizzati lanciando un dado”, allora Il modulo fractions permette di operare con le frazioni in modo simbolico, cioè sono trattate formalmente fino alla fine dei calcoli. fractions

import fractions # Fractions

import math # sqrt()

X = list(range(1, 7)) # Esiti : 1, 2, 3, ...

P = 6*[fractions.Fraction(1,6)] # probabilità: 1/6, 1/6, 1/6, ...

print(" x | p(x)")

print("---|------")

for x, p in zip(X, P):

print(" %i | %s" %(x, p))

Risultato

x | p(x)

---|------

1 | 1/6

2 | 1/6

3 | 1/6

4 | 1/6

5 | 1/6

6 | 1/6

Con i calcoli successivi sarà possibile leggere i … Leggi tutto