Tilbake
6.5

6.5 Drill: Gram–Schmidt, projeksjon og minste kvadrater

Ortogonalitetsapparatet drillet som én kjede: identifiser indreproduktet → Gram–Schmidt → ortogonal projeksjon/avstand → minste kvadrater via normallikningene, ofte på samme matrise.

85 min
13 oppgaver
DrillGram–Schmidtprojeksjonminste kvadrater
Din fremgang i kapitlet
0 / 13 oppgaver

Løsningsoppskrift — ortogonalitetskjeden

Hver oppgave i sjanger H/I løses med samme faste kjede. Lær den utenat (kode E — ingen formelsamling på eksamen):

1. Identifiser indreproduktet. Er det standard prikkprodukt uv=uivi\mathbf u\cdot\mathbf v=\sum u_iv_i, et vektet wiuivi\sum w_iu_iv_i, eller et integral p,q=abp(t)q(t)dt\langle p,q\rangle=\int_a^b p(t)q(t)\,dt? Alle skalarprodukt, normer og vinkler skal regnes i dette produktet.

2. Gram–Schmidt til ortogonal basis:
v1=a1,vk=aki<kak,vivi,vivi.\mathbf v_1=\mathbf a_1,\qquad \mathbf v_k=\mathbf a_k-\sum_{i<k}\frac{\langle\mathbf a_k,\mathbf v_i\rangle}{\langle\mathbf v_i,\mathbf v_i\rangle}\mathbf v_i.
Skalér gjerne hver vk\mathbf v_k til heltall underveis (det endrer ikke retningen). Normaliser til slutt hvis oppgaven ber om ortonormal basis: ui=vi/vi\mathbf u_i=\mathbf v_i/\|\mathbf v_i\|.

3. Projeksjon og avstand med den ortogonale basisen:
projWy=iy,vivi,vivi,y=y^+z,avstand=z=yprojWy.\operatorname{proj}_W\mathbf y=\sum_i\frac{\langle\mathbf y,\mathbf v_i\rangle}{\langle\mathbf v_i,\mathbf v_i\rangle}\mathbf v_i,\qquad \mathbf y=\hat{\mathbf y}+\mathbf z,\qquad \text{avstand}=\|\mathbf z\|=\|\mathbf y-\operatorname{proj}_W\mathbf y\|.

4. Minste kvadrater via normallikningene:
ATAx^=ATb.A^{T}A\hat{\mathbf x}=A^{T}\mathbf b.
Entydig når kolonnene i AA er uavhengige; affin løsningsmengde ved rangdefekt. Prediksjonen er Ax^=projColAbA\hat{\mathbf x}=\operatorname{proj}_{\operatorname{Col}A}\mathbf b — samme projeksjon som i steg 3.

Under følger en gjennomregnet case der samme matrise går gjennom Gram–Schmidt, projeksjon og minste kvadrater, med sensor-margnotater, så et oppgavesett på eksamensnivå.

Forkunnskaper — drillen forutsetter hele Del 6:

- kap. 6.1 — indreprodukt, norm og ortogonalitet (standard og integral).
- kap. 6.2 — Gram–Schmidt og ortonormalisering.
- kap. 6.3 — ortogonal projeksjon, dekomposisjon, avstand og komplement.
- kap. 6.4 — minste kvadrater og normallikningene.
- kap. 3.3 — kolonnerom og lineær uavhengighet.

Sist du var her (nøkkelformler):
vk=aki<kak,vivi,vivi,projWy=iy,vivi,vivi,ATAx^=ATb.\mathbf v_k=\mathbf a_k-\sum_{i<k}\frac{\langle\mathbf a_k,\mathbf v_i\rangle}{\langle\mathbf v_i,\mathbf v_i\rangle}\mathbf v_i,\quad \operatorname{proj}_W\mathbf y=\sum_i\frac{\langle\mathbf y,\mathbf v_i\rangle}{\langle\mathbf v_i,\mathbf v_i\rangle}\mathbf v_i,\quad A^{T}A\hat{\mathbf x}=A^{T}\mathbf b.

✏️Gjennomregnet eksamenscase (med sensor-margnotater)

La A=[10111213]A=\begin{bmatrix}1&0\\1&1\\1&2\\1&3\end{bmatrix} og b=[1335]\mathbf b=\begin{bmatrix}1\\3\\3\\5\end{bmatrix} (standard prikkprodukt i R4\mathbb R^4).

a) Bruk Gram–Schmidt på kolonnene i AA og finn en ortogonal basis for ColA\operatorname{Col}A.
b) Finn projColAb\operatorname{proj}_{\operatorname{Col}A}\mathbf b og avstanden fra b\mathbf b til ColA\operatorname{Col}A.
c) Finn minste kvadraters løsningen av Ax=bA\mathbf x=\mathbf b via normallikningene, og bekreft at Ax^A\hat{\mathbf x} er lik projeksjonen fra (b).

a) Gram–Schmidt. Kolonnene er a1=(1,1,1,1)\mathbf a_1=(1,1,1,1) og a2=(0,1,2,3)\mathbf a_2=(0,1,2,3).
v1=a1=(1,1,1,1).\mathbf v_1=\mathbf a_1=(1,1,1,1).
v2=a2a2v1v1v1v1=(0,1,2,3)64(1,1,1,1)=(32,12,12,32).\mathbf v_2=\mathbf a_2-\frac{\mathbf a_2\cdot\mathbf v_1}{\mathbf v_1\cdot\mathbf v_1}\mathbf v_1=(0,1,2,3)-\frac{6}{4}(1,1,1,1)=\left(-\tfrac32,-\tfrac12,\tfrac12,\tfrac32\right).
Skalér med 22 til heltall: v2=(3,1,1,3)\mathbf v_2=(-3,-1,1,3). Ortogonal basis: {(1,1,1,1),(3,1,1,3)}\{(1,1,1,1),(-3,-1,1,3)\}.
> Sensor: heltalls-oppskalering endrer ikke retningen og gir renere regning videre — men si eksplisitt at du skalerer.

b) Projeksjon og avstand. Med den ortogonale basisen:
projColAb=bv1v1v1v1+bv2v2v2v2.\operatorname{proj}_{\operatorname{Col}A}\mathbf b=\frac{\mathbf b\cdot\mathbf v_1}{\mathbf v_1\cdot\mathbf v_1}\mathbf v_1+\frac{\mathbf b\cdot\mathbf v_2}{\mathbf v_2\cdot\mathbf v_2}\mathbf v_2.
Her er bv1=1+3+3+5=12\mathbf b\cdot\mathbf v_1=1+3+3+5=12, v1v1=4\mathbf v_1\cdot\mathbf v_1=4; bv2=33+3+15=12\mathbf b\cdot\mathbf v_2=-3-3+3+15=12, v2v2=9+1+1+9=20\mathbf v_2\cdot\mathbf v_2=9+1+1+9=20. Da
projColAb=124(1,1,1,1)+1220(3,1,1,3)=3(1,1,1,1)+35(3,1,1,3)=(65,125,185,245).\operatorname{proj}_{\operatorname{Col}A}\mathbf b=\frac{12}{4}(1,1,1,1)+\frac{12}{20}(-3,-1,1,3)=3(1,1,1,1)+\tfrac35(-3,-1,1,3)=\left(\tfrac65,\tfrac{12}{5},\tfrac{18}{5},\tfrac{24}{5}\right).
> Sensor: bruk ortogonal basis FØR projeksjonsformelen. Bruker du a\mathbf a-vektorene direkte (de er ikke ortogonale), blir formelen feil.

Dekomposisjon b=b^+z\mathbf b=\hat{\mathbf b}+\mathbf z med z=bb^=(15,35,35,15)\mathbf z=\mathbf b-\hat{\mathbf b}=\left(-\tfrac15,\tfrac35,-\tfrac35,\tfrac15\right). Avstand:
z=1+9+9+125=2025=45=25.\|\mathbf z\|=\sqrt{\tfrac{1+9+9+1}{25}}=\sqrt{\tfrac{20}{25}}=\sqrt{\tfrac45}=\tfrac{2}{\sqrt5}.

c) Minste kvadrater. Normallikningene ATAx^=ATbA^{T}A\hat{\mathbf x}=A^{T}\mathbf b:
ATA=[46614],ATb=[1224].A^{T}A=\begin{bmatrix}4&6\\6&14\end{bmatrix},\qquad A^{T}\mathbf b=\begin{bmatrix}12\\24\end{bmatrix}.
Løsning: det=5636=20\det=56-36=20, x^=120[14664][1224]=120[2424]=(65,65).\displaystyle \hat{\mathbf x}=\frac{1}{20}\begin{bmatrix}14&-6\\-6&4\end{bmatrix}\begin{bmatrix}12\\24\end{bmatrix}=\frac{1}{20}\begin{bmatrix}24\\24\end{bmatrix}=\left(\tfrac65,\tfrac65\right).
> Sensor: navngi normallikningene. Her tilsvarer x^=(c0,c1)\hat{\mathbf x}=(c_0,c_1) regresjonslinja y=65+65xy=\tfrac65+\tfrac65x for punktene (0,1),(1,3),(2,3),(3,5)(0,1),(1,3),(2,3),(3,5).

Bekreftelse: Ax^=65(1,1,1,1)+65(0,1,2,3)=(65,125,185,245)=projColAbA\hat{\mathbf x}=\tfrac65(1,1,1,1)+\tfrac65(0,1,2,3)=\left(\tfrac65,\tfrac{12}{5},\tfrac{18}{5},\tfrac{24}{5}\right)=\operatorname{proj}_{\operatorname{Col}A}\mathbf b ✓. Kjeden lukker seg: minste kvadrater er projeksjon. \blacksquare

Oppgavesett (eksamensnivå, stigende)

Blandet drill av hele kjeden. Gjør oppgavene uten fasit først; alle svar er eksakte.

📝Oppgave 1
Sjanger H

I R3\mathbb R^3 med standard prikkprodukt: la u=(2,1,2)\mathbf u=(2,1,2) og v=(1,2,0)\mathbf v=(1,-2,0). a) Regn u\|\mathbf u\|. b) Er u\mathbf u og v\mathbf v ortogonale? c) Finn vinkelen mellom dem hvis de ikke er det.

📝Oppgave 2
Sjanger H

Finn en ortogonal basis for Span{(1,1,1),(2,0,1)}\operatorname{Span}\{(1,1,1),(2,0,1)\} ved Gram–Schmidt.

📝Oppgave 3
Sjanger H

Bruk Gram–Schmidt på a1=(1,0,1)\mathbf a_1=(1,0,1), a2=(1,1,0)\mathbf a_2=(1,1,0), a3=(0,1,1)\mathbf a_3=(0,1,1), og oppgi en ortogonal basis (skalert til heltall).

📝Oppgave 4
Sjanger H

Den ortogonale basisen {(1,1,1),(1,1,0)}\{(1,1,1),(1,-1,0)\} skal gjøres ortonormal. Finn den ortonormale basisen.

📝Oppgave 5
Sjanger H

La W=Span{(1,2,1),(1,0,1)}W=\operatorname{Span}\{(1,2,1),(1,0,-1)\} og y=(1,2,4)\mathbf y=(1,2,4). Finn projWy\operatorname{proj}_W\mathbf y, dekomposisjonen y=y^+z\mathbf y=\hat{\mathbf y}+\mathbf z, og avstanden fra y\mathbf y til WW.

📝Oppgave 6
Sjanger H

La W=Span{(1,1,1),(1,1,1)}W=\operatorname{Span}\{(1,1,1),(1,-1,1)\} og y=(4,1,3)\mathbf y=(4,1,3). Finn projWy\operatorname{proj}_W\mathbf y og avstanden fra y\mathbf y til WW.

📝Oppgave 7
Sjanger I

Finn minste kvadraters løsningen av Ax=bA\mathbf x=\mathbf b med A=[101112]A=\begin{bmatrix}1&0\\1&1\\1&2\end{bmatrix}, b=[312]\mathbf b=\begin{bmatrix}3\\1\\2\end{bmatrix}, og minste-kvadraters-feilen.

📝Oppgave 8
Sjanger I

Finn den rette linja y=c0+c1xy=c_0+c_1x som best tilpasser (0,1),(1,1),(2,2),(3,2)(0,1),(1,1),(2,2),(3,2).

📝Oppgave 9
Sjanger I

Finn parabelen y=c0+c1x+c2x2y=c_0+c_1x+c_2x^2 som best tilpasser (2,0),(1,1),(0,0),(1,4),(2,6)(-2,0),(-1,1),(0,0),(1,4),(2,6).

📝Oppgave 10
Sjanger H

P2\operatorname{P}_2 med indreproduktet p,q=11p(t)q(t)dt\langle p,q\rangle=\int_{-1}^{1}p(t)q(t)\,dt: bruk Gram–Schmidt på {1,t,t2}\{1,t,t^2\} og finn en ortogonal basis.

📝Oppgave 11
Sjanger H

Med samme indreprodukt p,q=11pqdt\langle p,q\rangle=\int_{-1}^{1}pq\,dt: finn projeksjonen av f(t)=t2f(t)=t^2 ned på W=Span{1,t}W=\operatorname{Span}\{1,t\}, og avstanden fprojWf\|f-\operatorname{proj}_W f\|.

📝Oppgave 12

(Kjedet — sjanger H + I.) La A=[11121314]A=\begin{bmatrix}1&1\\1&2\\1&3\\1&4\end{bmatrix} og b=[1345]\mathbf b=\begin{bmatrix}1\\3\\4\\5\end{bmatrix}. a) Gram–Schmidt på kolonnene i AA. b) projColAb\operatorname{proj}_{\operatorname{Col}A}\mathbf b. c) Minste kvadraters løsningen av Ax=bA\mathbf x=\mathbf b, og bekreft Ax^=projColAbA\hat{\mathbf x}=\operatorname{proj}_{\operatorname{Col}A}\mathbf b.

📝Oppgave 13

(Kald bank — rangdefekt.) La A=[120121122123]A=\begin{bmatrix}1&2&0\\1&2&1\\1&2&2\\1&2&3\end{bmatrix} (kolonne 2 =2×=2\,\times kolonne 1) og b=(1,2,2,5)\mathbf b=(1,2,2,5). Finn alle minste kvadraters løsninger av Ax=bA\mathbf x=\mathbf b.

Symbol- og formelliste

Formelkort for kjeden — repetisjon

Flashcard-/repetisjonsstoff — hopp trygt over ved førstegangslesing; tidsanslaget over gjelder kjernestoffet. Kortene under komprimerer kjeden.

Løsningsoppskriften for ortogonalitetskjeden

Den faste rekkefølgen for sjanger H/I: (1) identifiser indreproduktet; (2) Gram–Schmidt til ortogonal basis; (3) projeksjon/dekomposisjon/avstand med den ortogonale basisen; (4) minste kvadrater via normallikningene. Samme matrise går ofte gjennom hele kjeden.

Steg 1 — identifiser indreproduktet

Alt (skalarprodukt, norm, vinkel, projeksjon) regnes i det oppgitte indreproduktet: standard uv=uivi\mathbf u\cdot\mathbf v=\sum u_iv_i, vektet wiuivi\sum w_iu_iv_i, eller integral p,q=abp(t)q(t)dt\langle p,q\rangle=\int_a^b p(t)q(t)\,dt. Å bruke standard prikkprodukt når et annet er oppgitt, er den klassiske fellen.

Gram–Schmidt-formelen
Fra en basis {a1,,an}\{\mathbf a_1,\dots,\mathbf a_n\} bygges en ortogonal basis:

v1=a1,vk=aki<kak,vivi,vivi.\mathbf v_1=\mathbf a_1,\qquad \mathbf v_k=\mathbf a_k-\sum_{i<k}\frac{\langle\mathbf a_k,\mathbf v_i\rangle}{\langle\mathbf v_i,\mathbf v_i\rangle}\mathbf v_i.

Trekk fra projeksjonen på det allerede ortogonaliserte (vi\mathbf v_i, ikke ai\mathbf a_i).

Heltalls-oppskalering

Underveis i Gram–Schmidt kan hver vk\mathbf v_k ganges med et positivt tall uten at retningen (eller ortogonaliteten) endres. Å skalere bort brøker gir renere regning i de neste stegene. Normaliser til slutt hvis ortonormal basis kreves: ui=vi/vi\mathbf u_i=\mathbf v_i/\|\mathbf v_i\|.

Ortogonal projeksjon med ortogonal basis
Med en ortogonal basis {v1,,vp}\{\mathbf v_1,\dots,\mathbf v_p\} for WW:

projWy=i=1py,vivi,vivi.\operatorname{proj}_W\mathbf y=\sum_{i=1}^p\frac{\langle\mathbf y,\mathbf v_i\rangle}{\langle\mathbf v_i,\mathbf v_i\rangle}\mathbf v_i.

Formelen krever ortogonal basis — er spennvektorene ikke ortogonale, må du Gram–Schmidte først.

Ortogonal dekomposisjon og avstand
Enhver y\mathbf y splittes entydig som y=y^+z\mathbf y=\hat{\mathbf y}+\mathbf z med y^=projWyW\hat{\mathbf y}=\operatorname{proj}_W\mathbf y\in W og z=yy^W\mathbf z=\mathbf y-\hat{\mathbf y}\in W^{\perp}. Avstanden fra y\mathbf y til WW er lengden av z\mathbf z:

avstand=z=yprojWy.\text{avstand}=\|\mathbf z\|=\|\mathbf y-\operatorname{proj}_W\mathbf y\|.

Normallikningene
Minste kvadraters løsningen x^\hat{\mathbf x} av et overbestemt Ax=bA\mathbf x=\mathbf b finnes fra

ATAx^=ATb,A^{T}A\hat{\mathbf x}=A^{T}\mathbf b,

og prediksjonen er Ax^=projColAbA\hat{\mathbf x}=\operatorname{proj}_{\operatorname{Col}A}\mathbf b. Navngi dem alltid i besvarelsen.

Affin løsning ved rangdefekt

Er kolonnene i AA lineært avhengige, er ATAA^{T}A singulær og normallikningene har uendelig mange løsninger — en affin mengde x^0+NulA\hat{\mathbf x}_0+\operatorname{Nul}A. Alle gir samme prediksjon Ax^A\hat{\mathbf x}. Anta aldri entydig løsning uten å ha sjekket at kolonnene er uavhengige.

Dette kapitlet er skrevet av Anthropics toppmodeller (Claude Opus og Claude Fable) og er foreløpig ikke manuelt gjennomgått — kvalitetskontrollen gjøres av uavhengige KI-agenter, og innmeldte feil rettes fortløpende. Funnet en feil? Meld fra, så retter vi den. Les mer om hvordan innholdet lages.

Skolesaga er en uavhengig læringsressurs og er ikke tilknyttet eller godkjent av Norges teknisk-naturvitenskapelige universitet. Dette er ikke offisielt studiemateriell. Les mer.