PDF-format
PDF-format
PDF-format
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
38<br />
løsninger til overbestemte ligningssystemer).<br />
Hvis man bare skriver<br />
> LeastSquares(data13,x);<br />
8.812137559 − 0.494268374915711906 x<br />
får man den linie, der bedst (i mindste kvadraters betydning) beskriver datapunkterne. Nu er en linie næppe nogen<br />
god interpolation af disse datapunkter, det kunne tænkes at et andengradspolynomium ville være bedre. Man kan til<br />
LeastSquares angive hvilken form den funktion, man ønsker at tilpasse med, skal have. Funktionen skal være lineær i<br />
sine parametre, men ikke nødvendigvis lineær i funktionsvariablen. Så vi kan forsøge med et andengradspolynomium:<br />
> kurve1:=LeastSquares(data13,x,curve=a*x^2+b*x+b);<br />
kurve1 := 2.833792084 − 0.313247450643479031 x 2 + 2.83379208407862526 x<br />
Et alternativ kunne være et andengradspolynomium i logaritmen til x i stedet:<br />
> kurve2:=LeastSquares(data13,x,curve=a*log(x)^2+b*log(x)+c);<br />
kurve2 := 3.775907628 + 10.5396548308231992 ln(x) − 4.93559635161722721 ln(x) 2<br />
Vi kan tegne de to kurver sammen med de oprindelige punkter for at sammenligne interpolationerne:<br />
> display(pointplot(data13,symbol=circle),plot([kurve1,kurve2],x=1..10)<br />
> );<br />
8<br />
6<br />
4<br />
2<br />
0<br />
2 4 6 8 10<br />
Som det sidste i dette afsnit skal kort nævnes hvordan man får sine data eksporteret som en tabulatorsepareret<br />
tekstfil. Det gør man med funktionen writedata, der er fuldstændig analog med readdata, inklusiv muligheder der er<br />
for at angive kolonne<strong>format</strong> (kommatal er standard):<br />
> writedata("nyfil.txt",data13);<br />
> restart;<br />
x