0 Daumen
653 Aufrufe

Aufgabe:

Gegeben ist die folgende Messreihe:

x-2-1023
y-5-11,528


Berechnen Sie mit Hilfe der Differentialrechnung für Funktionen mehrerer Veränderlicher durch Minimierung der Summe der quadratischen Abweichungen das beste Ausgleichspolynom dritten Grades für die Betrachtung y = y(x) , das den Punkt (0, 1) enthält.


Ansatz:

Wäre für einen kompletten Lösungsweg dankbar.

Verstehe nicht wie ich das betrachten soll, eine Herleitung der Funktion wäre auch schonmal ein guter Anfang und eventuell wie ich das dann später mit den gegebenen Punkt mache.

Avatar von

1 Antwort

0 Daumen
 
Beste Antwort

Aloha :)

Ein Polynom 3-ten Grades hat die Form: y(x)=ax3+bx2+cx+dy(x)=ax^3+bx^2+cx+d. Da das Ausgleichspolynom den Punkt (01)(0|1) sicher enthalten soll, legen wir, wegen y(0)=dy(0)=d, die Konstante d=1d=1 fest:y(x)=ax3+bx2+cx+1y(x)=ax^3+bx^2+cx+1Du kannst nun 2 Vektoren bestimmen. Der eine Vektor enthält die gemessenen yy-Werte, nennen wir sie mal ymy_m, der andere Vektor enthält die gemäß des Ausgleichspolynoms y(xm)y(x_m) berechneten Werten für die Messorte xmx_m. Das sieht dann etwa so aus:

ym=(511,528);y=(y(2)y(1)y(0)y(2)y(3))=(8a+4b2c+1a+bc+118a+4b+2c+127a+9b+3c+1)\vec y_m=\begin{pmatrix}-5\\-1\\1,5\\2\\8\end{pmatrix}\quad;\quad \vec y=\begin{pmatrix}y(-2)\\y(-1)\\y(0)\\y(2)\\y(3)\end{pmatrix}=\begin{pmatrix}-8a+4b-2c+1\\-a+b-c+1\\1\\8a+4b+2c+1\\27a+9b+3c+1\end{pmatrix}Die Idee ist nun, die Parameter a,b,ca,b,c so zu bestimmen, dass der Vektor y\vec y dem Vektor ym\vec y_m "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)f(a,b,c), die es zu minimieren gilt:f(a,b,c)=yym=(8a+4b2c+6a+bc+20,58a+4b+2c127a+9b+3c7)!Minimumf(a,b,c)=\left\|\vec y-\vec y_m\right\|=\left\|\begin{pmatrix}-8a+4b-2c+6\\-a+b-c+2\\-0,5\\8a+4b+2c-1\\27a+9b+3c-7\end{pmatrix}\right\|\stackrel{!}{\to}\text{Minimum}Anstatt 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)=f^2(a,b,c):

F(a,b,c)=(8a+4b2c+6)2+(a+bc+2)2+(0,5)2+F(a,b,c)=(-8a+4b-2c+6)^2+(-a+b-c+2)^2+(-0,5)^2+F(a,b,c)=(8a+4b+2c1)2+(27a+9b+3c7)2\phantom{F(a,b,c)}=(8a+4b+2c-1)^2+(27a+9b+3c-7)^2F(a,b,c)=858a2+484ab+228ac494a+114b2+52bc82b+18c274c+90,25\phantom{F(a,b,c)}=858 a^2 + 484 a b + 228 a c - 494 a + 114 b^2 + 52 b c - 82 b + 18 c^2 - 74 c + 90,25Die partiellen Ableitungen müssen verschwinden:0=Fa=2(858a+242b+114c247)0=\frac{\partial F}{\partial a}=2 (858 a + 242 b + 114 c - 247)0=Fb=2(242a+114b+26c41)0=\frac{\partial F}{\partial b}=2 (242 a + 114 b + 26 c - 41)0=Fc=2(114a+26b+18c37)0=\frac{\partial F}{\partial c}=2 (114 a + 26 b + 18 c - 37)Wir haben also ein lineares Gleichungssysem gefunden:

(abc=8582421142472421142641114261837)(abc)(0,4371980,6050720,160628)\left(\begin{array}{r}a & b & c & =\\\hline 858 & 242 & 114 & 247 \\ 242 & 114 & 26 & 41 \\ 114 & 26 & 18 & 37\end{array}\right)\quad\Rightarrow\quad\begin{pmatrix}a\\b\\c\end{pmatrix}\approx\begin{pmatrix}0,437198\\-0,605072\\0,160628\end{pmatrix}

Damit sind wir am Ziel:y(x)=0,437198x30,605072x2+0,160628x+1y(x)=0,437198x^3-0,605072x^2+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)


Avatar von 153 k 🚀

Ein anderes Problem?

Stell deine Frage