4 175
muokkausta
Ero sivun ”Gnuplot” versioiden välillä
Siirry navigaatioon
Siirry hakuun
→Käyrän sovitus pistejoukkoon: virherajojen huomioiminen
Ape (keskustelu | muokkaukset) pEi muokkausyhteenvetoa |
(→Käyrän sovitus pistejoukkoon: virherajojen huomioiminen) |
||
| Rivi 129: | Rivi 129: | ||
Määritellään ensin funktio f(x): | Määritellään ensin funktio f(x): | ||
f(x)=A*x+B | f(x)=A*x+B | ||
y=A*x+B siis on suoran yhtälö, joten olemme sovittamassa suoraa. Gnuplot selvittää vakiot A ja B. | <tt>y=A*x+B</tt> siis on suoran yhtälö, joten olemme sovittamassa suoraa. Gnuplot selvittää vakiot A ja B. | ||
Varsinainen sovitus tehdään komennolla | Varsinainen sovitus tehdään komennolla | ||
| Rivi 150: | Rivi 150: | ||
B -0.946 1.000 | B -0.946 1.000 | ||
</pre> | </pre> | ||
Eli pistejoukkoon sovitetun suoran yhtälö on 3,66307*x+55,6119. Kulmakertoimen virhe on +/- 0,3708 ja vakiotermin +/- 2,34. | Eli pistejoukkoon sovitetun suoran yhtälö on <tt>3,66307*x+55,6119</tt>. Kulmakertoimen virhe on <tt>+/- 0,3708</tt> ja vakiotermin <tt>+/- 2,34</tt>. | ||
Nyt kun sovitus on tehty, voidaan sovitettu kuvaaja piirtää pistejoukkoon komennolla | Nyt kun sovitus on tehty, voidaan sovitettu kuvaaja piirtää pistejoukkoon komennolla | ||
| Rivi 158: | Rivi 158: | ||
plot "mittaus.dat" with xyerrorbars,f(x) | plot "mittaus.dat" with xyerrorbars,f(x) | ||
Tässä siis oikeastaan piirretään kaksi kuvaajaa, jotka erotetaan toisistaan pilkulla. | Tässä siis oikeastaan piirretään kaksi kuvaajaa, jotka erotetaan toisistaan pilkulla. | ||
=== Virherajojen huomioiminen sovituksessa === | |||
Oletuksena Gnuplot ei huomioi pisteille asetettuja virhetermejä käyrää sovittaessaan. Sovitukseen voidaan kuitenkin asettaa jokaista pisteparia vastaamaan tietty painokerroin joka kuvaa pisteen virherajoja. Kaksiulotteisessa sovituksessa annettava virhetermi on y-koordinaatin keskihajonta, ts. Gnuplot asettaa pisteparin painokertoimeksi <tt>1/s^2</tt>, missä s on annettu keskihajonta. | |||
Käytännössä tämä tapahtuu siten, että sovituskomennossa kerrotaan Gnuplotille missä tiedoston sarakkeessa virhetermi on. Jos edellä olevassa esimerkissä päättäisimme, että voimme jättää x-koordinaatin virheen huomioon ottamatta, kertoisimme virheen olevan tiedoston 4. sarakkeessa jolloin sovitus tapahtuisi komennolla | |||
fit f(x) "mittaus.dat" using 1:2:4 via A,B | |||
Tässä siis parametri <tt>using 1:2:4</tt> kuvaa tiedoston rakenteen: ensimmäisessä sarakkeessa on x-koordinaatti, toisessa y-koordinaatti ja neljännessä y-koordinaatin keskihajonta. | |||
Lisätietoja sovitustoiminnosta löytyy Gnuplotin ohjeista komennolla | |||
help fit | |||
==Aiheesta muualla== | ==Aiheesta muualla== | ||