OldClaudio” post=79168Hai quattro file corrispondenti a 4 passi di tempo in progressione geometrica che contengono la posizione nello spazio e le tre componenti della velocità nello spazio.
Vuoi tracciare l’errore fra i moduli delle posizioni e delle velocità all’ultimo campione di ciascun file, visto che tutti e quattro i file si riferiscono ad un intervallo complessivo di 50s.
Dunque vorresti fare un diagramma con quattro ascisse in progressione geometrica che in ordinata abbiano l’errore finale (sicuro la differenza fra i moduli e non il modulo delle differenze?). Secondo te questo diagramma dovrebbe rappresentare un diagramma lineare in cui l’errore diminuisce linearmente con il passo scelto in progressione geometrica.Mah; non metto lingua nell’approccio concettuale che evidentemente non ho capito, perché per quel che ho capito (ed evidentemente non l’ho capito) esso è errato.
Ma veniamo al punto che chiedi; invece di spremere le meningi a trovare un accrocchio che ti permetta di scrivere del codice complicatissimo epe fare una cosa che a mano si fa in pochi secondi, io ti suggerirei di fare a mano.
Prepari un quinto file dove nella prima colonna metti il passo di tempo e nelle colonne 2-4 le tre coordinate della posizione simulata e nelle colonne 5-7 le tre coordinate della posizione calcolata analiticamente; nelle colonne 8-10 e 11-13 fai lo stesso per le tre componenti di velocità simulate e calcolate; Il file contiene solo quattro righe e i valori da mettere nelle colonne sono trasferiti per copia e incolla dai file che hai già. Poi con questo quinto file disegni quello che credo di aver capito che tu voglia disegnare.
Che l’errore diminuisca al diminuire del passo è un fatto noto; che diminuisca linearmente con il passo non credo che si possa vedere da un diagramma lineare errore vs passo; forse con entrambe le scale logaritmiche, ma non certo con le scale lineari specialmente se il passo è in progressione geometrica. Stai attento all’errore algoritmico che potrebbe consegnarti funzioni di errore per passi piccoli costituite da rumore “deterministico”.
Naturalmente se non o capito il tuo problema, è evidente che ho fatto commenti e ho dato suggerimenti che non stanno in piedi. 😐 😯 🙄
Hai ragione che sono forse andato a finire all’ufficio “complicazione affari semplici” e che fare a mano è più veloce. Farò così.
Per quanto riguarda l’approccio concettuale mi sembra giusto. Se ho che per la soluzione analitica il raggio a 50s è r(50s), per una soluzione numerica con un passo di tempo di 1s il raggio a 50s è r_1(50s) e se per una soluzione numerica con un passo di tempo di 0.1s è r_0.1(50s) un grafico che riporta i seguenti valori
PassoDiTempo Errore
0s r(50s)-r(50s)
0.1s r_0.1(50s)-r(50s)
1s r_1(50s)-r(50s)
è una retta (per via dell’errore di troncatura). Penso di essermi spiegato male precedentemente.