Radbyte när vi kör LU-faktorisering!
Hallå!
Har stött på ett litet mindre problem, när man ska göra radbyte, vilket jag hade för mig att man gör efter varje "iteration" (dvs. efter varje gauss-elimination, IFALL ena elementet är större än andra), men det blir bara fel för min del. När jag jämför mina beräkningar och lösning med facit verkar det som jag på något konstigt sätt också fått rätt men kan inte se hur:
Det jag stör mig på är att jag tycker att |-0,5|>0 och därför borde vi göra ett till radbyte, men så verkar inte vara fallet. Varför? Vart går det snett?
uppsalairaniern skrev :Hallå!
Har stött på ett litet mindre problem, när man ska göra radbyte, vilket jag hade för mig att man gör efter varje "iteration" (dvs. efter varje gauss-elimination, IFALL ena elementet är större än andra), men det blir bara fel för min del. När jag jämför mina beräkningar och lösning med facit verkar det som jag på något konstigt sätt också fått rätt men kan inte se hur:
Det jag stör mig på är att jag tycker att |-0,5|>0 och därför borde vi göra ett till radbyte, men så verkar inte vara fallet. Varför? Vart går det snett?
Man gör ett radbyte om du har t.ex. 1 som första element och du vill ts bort 4 under. Det betyder att du måste multiplicera raden med 4 och sedan subtrahera från den undre raden. Det betyder att felen blir större. Därför vill du alltid multiplicera med ett tal 0<x<1. Så du byter rad på 1 och 4 så får du multiplicera 4 med 1/4 och sedan subtrahera. Då har du en stabil algoritm.
