Aloha :)
Ein Polynom 3-ten Grades hat die Form: y(x)=ax3+bx2+cx+d. Da das Ausgleichspolynom den Punkt (0∣1) sicher enthalten soll, legen wir, wegen y(0)=d, die Konstante d=1 fest:y(x)=ax3+bx2+cx+1Du kannst nun 2 Vektoren bestimmen. Der eine Vektor enthält die gemessenen y-Werte, nennen wir sie mal ym, der andere Vektor enthält die gemäß des Ausgleichspolynoms y(xm) berechneten Werten für die Messorte xm. Das sieht dann etwa so aus:
ym=⎝⎜⎜⎜⎜⎜⎛−5−11,528⎠⎟⎟⎟⎟⎟⎞;y=⎝⎜⎜⎜⎜⎜⎛y(−2)y(−1)y(0)y(2)y(3)⎠⎟⎟⎟⎟⎟⎞=⎝⎜⎜⎜⎜⎜⎛−8a+4b−2c+1−a+b−c+118a+4b+2c+127a+9b+3c+1⎠⎟⎟⎟⎟⎟⎞Die Idee ist nun, die Parameter a,b,c so zu bestimmen, dass der Vektor y dem Vektor ym "möglichst nahe" kommt. Das heißt, der Betrag der Differenz der beiden Vektoren soll minimal sein. Wir definieren dafür eine Funktion f(a,b,c), die es zu minimieren gilt:f(a,b,c)=∥y−ym∥=∥∥∥∥∥∥∥∥∥∥∥⎝⎜⎜⎜⎜⎜⎛−8a+4b−2c+6−a+b−c+2−0,58a+4b+2c−127a+9b+3c−7⎠⎟⎟⎟⎟⎟⎞∥∥∥∥∥∥∥∥∥∥∥→!MinimumAnstatt den Betrag des Vektors zu minimieren, können wir mit demselben Ergebnis auch das Quadrat minimieren. Das spart uns beim Berechnen die Wurzel. Sei daher F(a,b,c)=f2(a,b,c):
F(a,b,c)=(−8a+4b−2c+6)2+(−a+b−c+2)2+(−0,5)2+F(a,b,c)=(8a+4b+2c−1)2+(27a+9b+3c−7)2F(a,b,c)=858a2+484ab+228ac−494a+114b2+52bc−82b+18c2−74c+90,25Die partiellen Ableitungen müssen verschwinden:0=∂a∂F=2(858a+242b+114c−247)0=∂b∂F=2(242a+114b+26c−41)0=∂c∂F=2(114a+26b+18c−37)Wir haben also ein lineares Gleichungssysem gefunden:
⎝⎜⎜⎜⎛a858242114b24211426c1142618=2474137⎠⎟⎟⎟⎞⇒⎝⎛abc⎠⎞≈⎝⎛0,437198−0,6050720,160628⎠⎞
Damit sind wir am Ziel:y(x)=0,437198x3−0,605072x2+0,160628x+1
Plotlux öffnen f1(x) = 0,437198x3-0,605072x2+0,160628x+1P(0|1)P(-2|-5)P(-1|-1)P(0|1,5)P(2|2)P(3|8)Zoom: x(-3…4) y(-6…9)