Re: [semi-OT] integrale di linea sinusoidale

From: Soviet_Mario <SovietMario_at_CCCP.MIR>
Date: Tue, 17 Jul 2018 00:33:45 +0200

Il 16/07/2018 20:37, Elio Fabri ha scritto:
> Giorgio Bibbiani ha scritto:
>> Io avevo capito tutt'altro, e la risposta che ho dato si
>> basava su
>> quanto avevo capito...
> In effetti non è facile capirlo.

Oh ... non ci sto più capendo molto.
Allora diciamo questo con un esempio non ambiguo
per me la sinusoide classica senza parametri (y=sen x, o
y=cos x) ha ampiezza DUE (da -1 a +1).
Se però dite (non l'ho capito ancora) che questo è il doppio
dell'ampiezza, che sarebbe quindi 1, allora va bene

> A me ha risposto che va tutto bene quello che ho detto e il
> calcolo
> che ho fatto, ma secondo me, da quello che dice, avrei
> dovuto prendere
> come ampiezza della sinusoide 27 mm.
> Tu dici che avevi capito tutt'altro, ma stai un po' a vedere
> :-)
>
> Avevo detto di aver fatto il calcolo con Maxima.

mmm, indagherò su questo SW.
Stavo provando a esplorare OCTAVE, ma trovandolo non
intuitivo, e sarebbe tra l'altro il classico "overkill" mi sa.
Se magari Maxima ha una sintassi più intuitiva, beh lo
assaggerò volentieri

> Copio qui la routine, che è chiara, salvo l'ultima riga.
>
> ======================================l: 200$
> h: 54$
> f : h*sin(2*%pi*x/l)$
> g: diff(f,x)$
> F: sqrt(1 + g^2)$
> quad_qag(F,x,0,l,1,'epsabs=0.01);
>
> ======================================I $ in fine riga sono
> terminatori che non danno risposta.
> Invece il ";" produce la risposta.
>
> "quad_qag" è la funzione d'integrazione numerica.
> Gli argomenti sono:
> - la funzione integranda
> - la var. d'integrazione
> - i due estremi
> - "1" indica il procedimento più semplice, largamente suff.
> nel nostro
> caso
> - epsabs è il max errore assoluto che si richiede.
>
> Questa mi risponde:
>
> [304.8898287907637,7.036287023254317*10^-6,105,0]
>
> Il primo numero è il risultato.
> Il secondo la stima dell'errore assoluto.
> Il terzo il n. d'intervalli usati.
> L'ultimo è il codice di errore: 0 significa OK.
>
> Questo era il calcolo con h=54.
> E ora lo ripeto con h=27.
> Ecco che cosa mi dice.
>
> [232.1788401312093,7.060734150132753*10^-10,105,0]
>
> Dunque proprio quello che avevi trovato tu!
>
> Osservazione finale.
> Non ho provato, ma si poteva usare un arco di parabola al
> posto della
> sinusoide.

ah non ci ho pensato. Avevo usato l'ellissi ma non era il
profilo migliore.

> L'integrale si esprime con f. elementari, e sono sicuro che il
> risultato cambierebbe pochissimo.
> Molto meno delle tolleranze di lavorazione del lattoniere :-)

specialmente dal momento che il lattoniere sono le mie mani,
forbicioni e martello :(

>


-- 
1) Resistere, resistere, resistere.
2) Se tutti pagano le tasse, le tasse le pagano tutti
Soviet_Mario - (aka Gatto_Vizzato)
Received on Tue Jul 17 2018 - 00:33:45 CEST

This archive was generated by hypermail 2.3.0 : Mon Jan 20 2025 - 04:22:55 CET