9 svar
102 visningar
Liddas är nöjd med hjälpen!
Liddas 305
Postad: 21 jan 2020 Redigerad: 21 jan 2020

Newton Raphson metod!

Hej, jag har gjort denna för hand, men när jag gör det i matlab stämmer det inte riktigt, nu står det f(y), det är vad jag approximerar dvs var linjen skär y-axeln , men approximationen går till på samma sätt som om man gjorde det med f(x).

 

Så här ser koden ut i matlab

sprintf('Beräkna roten för ekvationen e^y+y^3')
y=-0.5 %Första iterations gissning
i=0
delta= 100
while (delta > 1e-6)
    y=exp(y)+y^3                        %fx
    dy=exp(y)+3*(y^2)              %f'x
    NewY=y-(y/dy)                     %beräknar nya värdet på x
    delta = abs(y-NewY)
    y=NewY                                 %uppdatera värdet på y
    i = i+1
end

Laguna 10333
Postad: 22 jan 2020

Glömde du termen -0,7731 i sista uträkningen?

Liddas 305
Postad: 22 jan 2020 Redigerad: 22 jan 2020
Ja det har jag gjort, den sista blir -0,7728. Men koden gör fel på första dy också. Men när jag sätter in värdet g(y) där y = -0,7728 blir funktionsvärdet 0,0000336, vilket är ungerfär noll, det står dock inget i uppgiften hur många iterations-steg man ska göra.
Laguna 10333
Postad: 22 jan 2020

Koden gör y - y/dy, det är inte rätt. 

Liddas 305
Postad: 22 jan 2020 Redigerad: 22 jan 2020

Okej, men när jag gör det för hand då är

det y-(y/dy).  Lite osäker på hur matlab tolkar det där.

Liddas 305
Postad: 22 jan 2020

Undra  om man ens kan skriva y=exp(y)+y^3 . Kanske måste döpa om den till g(y)

Laguna 10333
Postad: 22 jan 2020
Liddas skrev:

Okej, men när jag gör det för hand då är

det y-(y/dy).  Lite osäker på hur matlab tolkar det där.

Nej, när du gör det för hand tar du y - f(y)/f'(y).

Liddas 305
Postad: 22 jan 2020 Redigerad: 22 jan 2020

Ja men om du kollar i koden är ju y=f(y) och dy=f’(y)  det är iaf menat så, det som kanske blir fel det är att jag har y=f(y) vilken kanske blir en konflikt / fel

Laguna 10333
Postad: 22 jan 2020

Gå igenom ett steg i koden mycket noggrant med penna och papper, så ser du vad som blir annorlunda.

Det bästa är kanske att formulera om uppgiften så att y är en funktion av x, som vanligt.

Liddas 305
Postad: 22 jan 2020 Redigerad: 22 jan 2020

    x=exp(y)+y^3 %fx
    dx=exp(y)+3*(y^2) %f'x
    NewY=y-(x/dx) %beräknar nya värdet på x

Ändrade om dom där, vilket jag (vi) också va inne på, blir lite knepigt med x och y, men  då blev det rätt .. Men tacksam för din hjälp på  vägen laguna!

Jag märker också att jag slarvar väldigt mycket, ska jobba på det !

Svara Avbryt
Close