StartAlgorytmyProcedury numeryczneCiąg Sturma
Baza Wiedzy
Wersja offline serwisu przeznaczona na urządzenia z systemem Android.

Darowizny
darowiznaWspomóż rozwój serwisu


Nagłówki RSS
Kanał artykuły
Kanał implementacje
Kanał komentarze
Kanał forum


Bookmarki









Sonda
Implementacji w jakim języku programowania poszukujesz?
 
Ciąg Sturma
Ocena użytkowników:+++++ / 3
SłabyŚwietny 
Wpisany przez Tomasz Lubiński
poniedziałek, 08 sierpnia 2005 21:07
Ciąg p0(x), p1(x),..., pn(x) nazywamy ciągiem Sturma jeśli:
- p0(x) ma tylko pojedyncze miejsca zerowe,
- sign p1(a)=-sign p0'(a) dla wszystkich rzeczywistych zer a wielomianu p0(x),
- dla i=1,2,...,n-1 pi+1(a)*pi-1(a)<0, jeśli a oznacza zera rzeczywiste wielomianu pi(a).
Liczba miejsc zerowych wielomianu p(x) w przedziale <a,b) jest równa |W(b)-W(a)|, gdzie W(x) oznacza liczbę zmian znaku w punkcie x w ciągu Sturma p0(x), p1(x),..., pn(x) zbudowanym w oparciu o p(x).
Sposób konstrukcji ciągu jest następujący:
p0(x)=p(x)
p1(x)=-p'(x)

Pozostałe to reszta z dzielenia przedostatniego przez ostatni pomnożona przez stałą mniejszą od zera.
Przykład:
Skonstruujemy ciąg Sturma wielomianu p(x)=x3-2x2-5x+5
p0(x)=p(x)=x3-2x2-5x+5
p1(x)=-p'(x)=-3x2+4x+5
Image
Resztę należy pomnożyć przez liczbę mniejszą od zera - dla ułatwienia rachunków niech będzie to -9
p2(x)=38x-35
Image
Resztę należy pomnożyć przez liczbę mniejszą od zera - dla ułatwienia rachunków niech będzie to -1444
p3(x)=-8865


x -2 0 1
p0(x) - + -
p1(x) - + +
p2(x) - - +
p3(x) - - -
Liczba zmian znaków W(x) 0 1 2


Z tabelki możemy odczytać m.in, że p(x) ma w przedziale <-2;0) jedno miejsce zerowe bo |W(0)-W(-2)| natomiast przedziale <-2;1) ma dwa miejsca zerowe bo |W(1)-W(-2)|



Autor Język programowania Komentarz Otwórz Pobierz Ocena
Tomasz Lubiński Delphi/Pascal Borland Delphi 5
Implementacja w Delphi/Pascal
Implementacja w Delphi/Pascal
++++- / 2
 
Dodaj własną implementację tego algorytmu
  • Zaloguj się na stronie
Plik:
Język programowania:
Komentarz:
  By móc dodać implementacje zaloguj się na stronie



Poprawiony: wtorek, 21 czerwca 2011 17:36

Dodaj komentarz

Kod antysapmowy
Odśwież