Reductie (kegelsnede)

De algemene vergelijking van een kegelsnede is van de vorm:

${\displaystyle ax^{2}+2bxy+cy^{2}+dx+ey+f=0}$

waarbij a, b en c niet alle drie nul zijn. Indien de coëfficiënt b niet nul is bevat de vergelijking een xy-term en zal een rotatie van het assenkruis nodig zijn om deze te verwijderen. Indien b nul is is geen rotatie vereist. Bovenstaande vergelijking kan in matrixvorm geschreven worden als:

${\displaystyle {\begin{bmatrix}x&y\end{bmatrix}}.{\begin{bmatrix}a&b\\b&c\end{bmatrix}}.{\begin{bmatrix}x\\y\end{bmatrix}}+{\begin{bmatrix}d&e\end{bmatrix}}.{\begin{bmatrix}x\\y\end{bmatrix}}+f=0.}$

of nog korter:

${\displaystyle X^{\tau }AX\,+\,LX\,+\,f\,=\,0}$

met:

${\displaystyle A\,=\,{\begin{bmatrix}a&b\\b&c\end{bmatrix}}\quad en\quad L\ =\,{\begin{bmatrix}d&e\end{bmatrix}}\quad en\quad X\ =\,{\begin{bmatrix}x\\y\end{bmatrix}}}$

De matrix A, die de informatie van de tweedegraadstermen bevat, toont nu waarom de xy-coëfficiënt in de algemene vergelijking als 2b werd geschreven. Op die manier kan die 2b gelijkmatig worden verdeeld over de twee niet-diagonaalelementen van de matrix, die daardoor symmetrisch wordt.

Een symmetrische matrix heeft een aantal belangrijke eigenschappen:

• Hij is steeds diagonaliseerbaar.
• Alle eigenwaarden zijn reëel.
• Eigenvectoren van verschillende eigenwaarden zijn orthogonaal (loodrecht).
• Bij meervoudige eigenwaarden kan nooit ontaarding optreden. Dit heeft hier tot gevolg dat, zelfs indien de twee eigenwaarden van A gelijk zijn, het steeds mogelijk is twee lineair onafhankelijke, en dus orthogonaliseerbare, eigenvectoren te vinden.

In alle gevallen kunnen er dus twee orthogonale eigenvectoren gevonden worden, die vervolgens na normering zelfs orthonormaal zijn.

Deze eigenschap wordt gebruikt om de rotatie concreet uit te voeren. Stel dat de matrix A twee eigenwaarden heeft, ${\displaystyle \lambda _{1}}$ en ${\displaystyle \lambda _{2}}$, met respectievelijke orthonormale eigenvectoren ${\displaystyle X_{1}}$ en ${\displaystyle X_{2}}$. De matrix P met als kolommen deze eigenvectoren is dan een orthogonale matrix die A diagonaliseert, met de matrix D als diagonaalvorm. D is een diagonaalmatrix met op de hoofddiagonaal de eigenwaarden in dezelfde volgorde als de eigenvectoren op de kolommen van P. Daarbij moet de zin van de orthonormale eigenvectoren zodanig gekozen worden dat de determinant van P gelijk aan 1 is, zodat de rotatie de rechtshandigheid van het assenkruis bewaart. Een eigenvector waarvan men de zin omdraait blijft steeds een eigenvector van dezelfde eigenwaarde zodat deze eventuele aanpassing steeds mogelijk is.

${\displaystyle P\,=\,[X_{1}\,X_{2}]\quad }$

${\displaystyle D\,=\,P^{-1}AP\,=\,{\begin{bmatrix}\lambda _{1}&0\\0&\lambda _{2}\end{bmatrix}}}$

Gezien P een orthogonale matrix is, bestaat zijn inverse en is die gelijk aan de getransponeerde van P. De eerste kolom bevat de richting van de x'-as uitgedrukt in het xy-assenkruis. De tweede kolom bevat de richting van de y'-as uitgedrukt in het xy-assenkruis. Door nu op drie plaatsen in de matrixvergelijking de eenheidsmatrix in de vorm ${\displaystyle I\,=\,P.P^{\tau }}$ in te voeren verkrijgt men:

${\displaystyle X^{\tau }P.P^{\tau }AP.P^{\tau }X\,+\,LP.P^{\tau }X\,+\,f\,=\,0}$

Indien men vervolgens stelt:

${\displaystyle X'\,=\,{\begin{bmatrix}x'\\y'\end{bmatrix}}\,=\,P^{\tau }X}$

bekomt men:

${\displaystyle X'^{\tau }DX'\,+L.P.X'\,+\,f\,=\,0}$

of voluit geschreven:

${\displaystyle \lambda _{1}x'^{2}\,+\,\lambda _{2}y'^{2}\,+\,d'x'\,+\,e'y'\,+\,f\,=\,0}$

Hieruit blijkt dat de kruiselingse term in x.y verdwenen is door de coördinaatstransformatie. Meetkundig betekent de rotatie van het assenkruis dat de assen zodanig gedraaid worden dat een van beide evenwijdig is aan de hoofdas van de kegelsnede, en de andere as er loodrecht op staat. De figuur horend bij het voorbeeld (zie verder) toont dit effect. Door de rotatie zullen in het algemeen wel de eventueel aanwezige lineaire coëfficiënten gewijzigd worden.

De aarde van de kegelsnede kan aan de hand van de determinant van de matrix A bepaald worden:

${\displaystyle |A|\,=\,ac-b^{2}\,=\lambda _{1}\lambda _{2}}$

Voor een ellips is deze positief, voor een hyperbool negatief en voor een parabool nul. Bij een ellips hebben beide eigenwaarden dus hetzelfde teken, bij een hyperbool een tegengesteld teken, en bij een parabool is één eigenwaarde nul.

Een verschuiving is nodig indien:

• een onbekende x' of y' na rotatie nog zowel kwadratisch als lineair voorkomt
• een onbekende voorkomt samen met een constante term f (dit kan optreden bij een parabool)

Deze twee voorwaarden dienen in deze volgorde toegepast te worden.

• In het eerste geval kan de lineaire term worden opgeheven op volgende manier, uitgelegd op een getalvoorbeeld:

${\displaystyle 2x'^{2}+12x'\,=\,2(x'^{2}+6x)\,=\,2(x'^{2}+6x'+9)-18\,=\,2(x'+3)^{2}-18}$

Eerst wordt de kwadratische coëfficiënt (2) voorop gezet. Vervolgens wordt binnen de haakjes het kwadraat van de helft van de daar aanwezige lineaire term toegevoegd (9), en wordt deze toevoeging buiten de haakjes weer tenietgedaan (−18). Tussen de haakjes staat nu een merkwaardig product dat als één kwadraat kan worden geschreven. De verschuiving is dan:

${\displaystyle x''\,=\,x'+3}$

De lineaire term is dus verdwenen door de verschuiving. Op dezelfde manier kan een term in y' weggeschoven worden door een analoge verschuiving in de y-richting.

• In het tweede geval wordt de constante term gewoon opgenomen in de nog resterende lineaire term. Bijvoorbeeld:

${\displaystyle 2x'+10\,=\,2(x'+5)\,=\,2x''}$

Na deze bewerkingen staat de vergelijking van de kegelsnede in een standaardvorm.

• Voor een ellips met halve lange as a:

${\displaystyle \quad {\frac {x''^{2}}{a^{2}}}+{\frac {y''^{2}}{b^{2}}}=1\quad }$

of

${\displaystyle \quad {\frac {y''^{2}}{a^{2}}}+{\frac {x''^{2}}{b^{2}}}=1}$

• Voor een hyperbool:

${\displaystyle \quad {\frac {x''^{2}}{a^{2}}}-{\frac {y''^{2}}{b^{2}}}=1\quad }$

of

${\displaystyle \quad {\frac {y''^{2}}{b^{2}}}-{\frac {x''^{2}}{a^{2}}}=1}$

• Voor een parabool:

${\displaystyle \quad y''^{2}=2px''}$

of

${\displaystyle \quad x''^{2}=2py''}$

Gezien tijdens de berekening een aantal keuzes moeten gemaakt worden, zoals bijvoorbeeld de nummering van de eigenwaarden, is de gereduceerde vergelijking niet uniek. Echter, bij elke keuze bekomt men ook telkens een ander xy-assenkruis, en dit compenseert het verschil in formule.

De hyperbool uit het voorbeeld getekend in de originele assen. De assen na rotatie zijn in het blauw weergegeven. De vervolgens verschoven assen in het zwart. In het verschoven zwarte assenstelsel heeft de hyperbool een standaardvergelijking

De kegelsnede K:

${\displaystyle 7x^{2}\,+\,50xy\,+\,7y^{2}\,-\,46{\sqrt {2}}x-\,82{\sqrt {2}}y\,+\,160\,=\,0}$

heeft als matrixvorm:

${\displaystyle {\begin{bmatrix}x&y\end{bmatrix}}.{\begin{bmatrix}7&25\\25&7\end{bmatrix}}.{\begin{bmatrix}x\\y\end{bmatrix}}+{\begin{bmatrix}-46{\sqrt {2}}&-82{\sqrt {2}}\end{bmatrix}}.{\begin{bmatrix}x\\y\end{bmatrix}}+160=0.}$

Aan de determinant van A, die negatief is, is te zien dat dit een hyperbool voorstelt. De twee eigenwaarden zullen dus verschillend van nul zijn, met tegengesteld teken. De eigenwaarden met hun respectievelijke genormeerde eigenvectoren zijn:

${\displaystyle \lambda _{1}\,=\,32\quad X_{1}\,=\,{\begin{bmatrix}1/{\sqrt {2}}\\1/{\sqrt {2}}\end{bmatrix}}}$

${\displaystyle \lambda _{2}\,=\,-18\quad X_{2}\,=\,{\begin{bmatrix}-1/{\sqrt {2}}\\1/{\sqrt {2}}\end{bmatrix}}}$

De rotatiematrix P en de diagonaalvorm D zijn dus:

${\displaystyle P\,=\,{\begin{bmatrix}1/{\sqrt {2}}&-1/{\sqrt {2}}\\1/{\sqrt {2}}&1/{\sqrt {2}}\end{bmatrix}}\quad en\quad D\ =\,{\begin{bmatrix}32&0\\0&-18\end{bmatrix}}}$

Na rotatie wordt de kegelsnede:

${\displaystyle 32x'^{2}\,-\,18y'^{2}\,-\,128x'\,-36y'\,+\,160\,=\,0}$

Voor elke kwadratische term is nog steeds ook een lineaire term aanwezig, en dus dient er verschoven te worden zowel in x' als in y'. Na vorming van de merkwaardige producten wordt dit:

${\displaystyle 32(x'-2)^{2}\,-\,18(y'+1)^{2}\,=\,-50}$

en ten slotte, na verschuiving:

${\displaystyle 18y''^{2}\,-\,32x''^{2}\,=\,50}$

Een mogelijke gereduceerde vorm is dus:

${\displaystyle {\frac {y''^{2}}{(3/5)^{2}}}\,-\,{\frac {x''^{2}}{(4/5)^{2}}}\,=\,1}$

Dit is een hyperbool die in het x"y"-assenkruis op de y"-as ligt. De toppen liggen in dit assenkruis in (0,±3/5) en de brandpunten in (0,±1).

Bij een hyperbool kan ook een asymptoot als y-as gekozen worden, en de oorsprong weer op het snijpunt van de asymptoten. De gereduceerde vorm is dan

${\displaystyle y=gx+h/x}$