Octave nie porównuje zmiennych w odpowiedni sposób.

Question

Witam,

rozwiązując sobie jedno z zadań na Metody Numeryczne trafiłem na pewien problem, który nie do końca rozumiem z czego wynika. Jeden z ostatnich kroków zadania mówi, abym obliczył błąd rzeczywisty, a następnie sprawdził czy jest on na poziomie epsilonu maszynowego.

Błąd wychodzi mi taki: b_rz = 1.343834447586900e-05

A kiedy sprawdzam wartość eps pojawia się taka liczba: eps = 2.220446049250313e-16

Cały problem pojawia się w momencie kiedy wpisuję następującą linijkę:

b_rz <= eps

Wynikiem tego ciągle jest 0, a wydaje mi się, że powinno być 1.

Czy jest tu może osoba, która byłaby w stanie mi wytłumaczyć czemu tak się dzieje? Czy to błąd Octave, czy ja coś źle robię lub źle rozumiem działanie porównań w tym programie?

 

 

Answer 1

no dlatego, że EPS nie jest większy ? zauważ że b_rz jest e-05, jest większe niz e-16
to są wartości po przecinku więc potęgi -16 są znacznie mniejszymi liczbami niż potęgi liczb do -5 

0.0001 > niz 0.0000001 np. 
 

x <= y

True if x is less than or equal to y.

Jeśli analizujesz pełne wartości a nie ich odcięcia, to program działa prawidłowo.  

Comments

Dziękuję za odpowiedź, jest tak jak mówisz .. siedzę z tym od rana i myślenie się chyba już samo wyłącza.

Jednak skoro wynik który pokazuje Octave jest poprawny, oznacza to, że musiałem zrobić gdzieś błąd w obliczeniach, bo zgodnie z komentarzem który usłyszałem na zajęciach "błąd powinien być na poziomie epsilonu maszynowego". Czy byłbyś w stanie pomóc mi zlokalizować go gdybym wstawił polecenie ze wzorem i swoje obliczenia? Bo ja już siedzę i siedzę, ale błędu za nic w świecie znaleźć nie mogę.

---

postarać się mogę, nie obiecuje że pomogę wrzuć co masz- przeanalizujemy. Aktualnie wartości wyraźnie wskazują  

b_rz = 1.343834447586900e-05  ( e-05)

eps = 2.220446049250313e-16  (e-16)
w oby przypadkach  -> tutaj istotne jest ze epsilony sa ujemne 

 

---

OK, udało mi się jakoś ogarnąć. Dzięki mimo wszystko za poprzednią odpowiedź.