|
Post by opperbanana on Jan 31, 2006 20:29:54 GMT 1
as ge de ortogonale projectie wil zoeken, kan je dit dan het best doen met de singuliere waarde ontbinding die je net hebt berekent dus : x(hoetje) = pseudoinverse*b met dan de projectie op kolom van A: U(gereduceerd)*U(gereduseerd+getransponeerd)*b
ik hoop dat jullier er aan uit kunnen
|
|
|
Post by Niels Bosmans on Jan 31, 2006 22:33:24 GMT 1
hoe vind ge vraag 2d) eigenlijk? Int boek staat wel uitgelegd als xTx=1 maar niet als xTx=4...
|
|
quetzl
~Machines~
~Machine 2~
Posts: 168
|
Post by quetzl on Jan 31, 2006 22:39:28 GMT 1
das gewoon (0,0,2,0,0) gebruiken voor x ipv (0,0,1,0,0) 4a is dan uw oplossing
|
|
|
Post by Niels Bosmans on Jan 31, 2006 22:46:45 GMT 1
mmm ja ok maar dr is zeker niet zo een kant en klare regel voor da was k aant zoeken . Merci nog!
|
|
|
Post by opperban on Feb 1, 2006 10:20:53 GMT 1
ma oe zoekte gelle die orthogonale projectie?
|
|
|
Post by Niels Bosmans on Feb 1, 2006 10:44:04 GMT 1
ja daar heb ik ook nog een probleem mee met die orthogonale projectie:
Ik heb dus eerst de kleinste kwadraten oplossing van Ax=B gezocht.
ATA: [4 0 0 0 4] [0 4 0 4 0] [0 0 a^2 0 0] [0 4 0 4 0] [4 0 0 0 4]
ATB: [-2 ] [ 2 ] [3a] [ 2 ] [-2 ]
Als ik dan ATAx=ATB oplos kom ik als x(hoedje) een oplossing die afhankelijk is van 2 vrije onbekenden.
Nu neem ik de orthogonale projectie van B op ColA. Als ik dan het stelsel Ax=B(hoedje) oplos zou ik dezelfde kleinste kwadraat oplossing moeten uitkomen... dat is niet zo.
voor die projectie heb ik UUTB berekend met U de 2 eerste genormaliseerde eigenvectoren uit de eigenwaardenontbinding horende bij eigenwaarden a en 4.
dit zijn: [0] [ 1/2 ] [0] [-1/2] [1] [ 0 ] [0] [-1/2] [0] [ 1/2 ]
je kan gemakkelijk nagaan dat dit een basis is voor Col A btw. De projectie van B op ColA is dan UUTB en dit is: [-1/2] [ 1/2] [ 3 ] [ 1/2] [-1/2 ]
Als ik dan Ax=B(hoedje) oplos dan kom ik een x(hoedje) afhankelijk van 3 vrije onbekenden.
Ook bij Ax=C en de projectie van C op ColA heb ik hetzelfde probleem...
Kan iemand helpen aub?
|
|
|
Post by Zeus on Feb 1, 2006 13:48:55 GMT 1
Wat zijn om te beginnen de eigenwaarden bij deze vraag, ik heb a, 2, 1, 1 en 0. Vond die twee op het eerste zicht verdacht en het zou zot zijn om met rekenfouten verder te rekenen... alvast bedankt!
|
|
|
Post by Zeus on Feb 1, 2006 14:19:00 GMT 1
hm mijn eigenwaarden kunnen zeker niet kloppen superdomm fout gemaakt...
|
|
|
Post by Niels Bosmans on Feb 1, 2006 14:25:03 GMT 1
de eigenwaarden zijn: a,4,0,0,0 klopt 100% zeker, vraag maar aan maple
|
|
quetzl
~Machines~
~Machine 2~
Posts: 168
|
Post by quetzl on Feb 1, 2006 14:28:38 GMT 1
Ge kunt da allemaal op zicht zien , gewoon een beetje trainen op symmetrische matrix'n
|
|
|
Post by yentlswolfs on Feb 1, 2006 14:34:09 GMT 1
Uw Ata klopt niet niels. Waar jij nullen hebt staan heb ik -4 staan.(behalve op de rij en de kolom met a op, daar staan gewoon nullen en a².
Hoe vinde die eigenwaarden, met de spilmethode bekom a,1,2,3+i en 3-i . Wat dus zeker niet kan, want een symmetrische matrix heeft reele eigenwaarden. Wsl mag je de spilmethode dus niet toepassen, maar hoe doe je het dan? Met rij-operaties gaat toch niet omdat je met termen 1-lambda zit?
|
|
|
Post by Alex Prinsier on Feb 1, 2006 14:36:14 GMT 1
das gewoon (0,0,2,0,0) gebruiken voor x ipv (0,0,1,0,0) 4a is dan uw oplossing die (0,0,2,0,0) heb ik ook, maar vanwaar komt die 4a? De maximale waarde blijft bij mij a, niet 4a.
|
|
|
Post by Niels Bosmans on Feb 1, 2006 14:52:37 GMT 1
Uw Ata klopt niet niels. Waar jij nullen hebt staan heb ik -4 staan.(behalve op de rij en de kolom met a op, daar staan gewoon nullen en a². Hoe vinde die eigenwaarden, met de spilmethode bekom a,1,2,3+i en 3-i . Wat dus zeker niet kan, want een symmetrische matrix heeft reele eigenwaarden. Wsl mag je de spilmethode dus niet toepassen, maar hoe doe je het dan? Met rij-operaties gaat toch niet omdat je met termen 1-lambda zit? zucht ik ga mijzelf ff slaan man wat een rekenfouten heb ik hier weer! over uw spilprobleem. h1.ripway.com/nigel87/Thematrix.doc
|
|
quetzl
~Machines~
~Machine 2~
Posts: 168
|
Post by quetzl on Feb 1, 2006 15:05:09 GMT 1
doe eens (0,0,2,0,0)*A*(0,0,2,0,0)t, ge zult het wel zien
|
|
|
Post by yentlswolfs on Feb 1, 2006 15:20:21 GMT 1
K merci niels.
|
|