Rekurencja pośrednia - czyli zabawy z rekurencją (część 4)

W pierwszych trzech częściach "Zabaw z rekurencją" skupialiśmy się na rekurencji bezpośredniej, tzn. na sytuacji, kiedy w ciele funkcji dochodzi do wywołania "siebie samej". Przebieg rekurencji bezpośredniej jest dość oczywisty, struktura wywołania, argumenty, jak też warunek stopu, są takie same dla wszystkich odwołań.

Rekurencja pośrednia

O rekurencji pośredniej mówimy w sytuacji"łańcucha wywołań". Przykładowo funkcja f(.) wywołuje funkcję g(.), następnie funkcja g(.) wywołuje f(.), zatem ponowne wywołanie funkcji f(.) realizowane jest bezpośrednio przez funkcję g(.), jednak pośrednio przez f(.), gdyż to f(.) wywołała g(.).

Typy rekurencji

Długość łańcucha nie musi być ograniczona, w rzeczywistości wywołania pośrednie mogą mieć nietrywialną strukturę, mogą "cofać się" do przednich elementów, "iść na skróty", "rozdzielać się", a w szczególności może dochodzić do wariantów mieszanych - tzn. wywołań bezpośrednich i pośrednich (różnego typu) w ramach jednej procedury. Dobrze to obrazuje poniższy schemat.

Rekurencja pośrednia

Aproksymacja funkcji sin(x) oraz cos(x) przy wykorzystaniu połączenia rekurencji bezpośredniej i rekurencji pośredniej

Przypomnijmy podstawowe tożsamości trygonometryczne dla wielokrotności kątów.

\sin(2x)=2\sin(x)\cos(x)

\cos(2x)=\cos^2(x)-\sin^2(x)

Równoważnie powyższe można zapisać jako

\sin(x)=2\sin\big(\frac{x}{2}\big)\cos\big(\frac{x}{2}\big)

\cos(x)=\cos^2\big(\frac{x}{2}\big)-\sin^2\big(\frac{x}{2}\big)

Zwróćmy uwagę, że znając rozwiązanie dla argumentu mniejszego \frac{x}{2} możemy podać rozwiązanie dla x – zatem tożsamości trygonometryczne są w istocie rekurencją z odwołaniami bezpośrednimi i pośrednimi! Funkcje \sin(x) w otoczeniu 0 można przybliżyć przez x, natomiast funkcję \cos(x) przez stałą wartość 1. Im mniejsze otoczenie 0 wybierzemy tym lepsza aproksymacja w zadanym przedziale, a w konsekwencji mniejszy błąd oszacowania w całości. Przyjęte wartości w otoczeniu 0 dają również pewny warunek stopu! Mamy więc wszystko co niezbędne do zastosowania strategii rekurencyjnej w aproksymacji.

Ustalmy stałą a>0 (reprezentującą otoczenie 0), następnie definiujemy dwie funkcje rekurencyjne

\text{s}(x)=\begin{cases}x&\text{dla}\quad |x|<a\\2\text{s}\big(\frac{x}{2}\big)\text{c}\big(\frac{x}{2}\big)&\text{dla}\quad |x|\geq a\end{cases}

\text{c}(x)=\begin{cases}1&\text{dla}\quad |x|<a\\\text{c}^2\big(\frac{x}{2}\big)-\text{s}^2\big(\frac{x}{2}\big)&\text{dla}\quad |x|\geq a\end{cases}

Podkreślmy ponownie, że funkcja \text{s}(x) wywołuje siebie bezpośrednio oraz wskazuje na funkcję \text{c}(x), która, oprócz bezpośredniego wywołania siebie samej, wskazuje ponownie na \text{s}(x). Jest to zatem ciekawa kombinacja rekurencji bezpośredniej z rekurencją pośrednią. Zapiszmy to w mXparser.

/* Definicja funkcji rekurencyjncyh */
Constant a = new Constant("a", 0.1);
Function s = new Function("s(x) =  if( abs(x) < a, x, 2*s(x/2)*c(x/2) )", a);
Function c = new Function("c(x) =  if( abs(x) < a, 1, c(x/2)^2-s(x/2)^2 )", a);

/* Wskazanie, ze 's' korzysta z 'c', a 'c' korzysta z 's' */
s.addDefinitions(c);
c.addDefinitions(s);

Oczekujemy, że im mniejszy parametr a>0 tym lepsza aproksymacja funkcji \sin(x) oraz \cos(x) przez odpowiednio \text{s}(x) oraz \text{c}(x). Poniżej wykresy dla a=0.5 oraz a=0.01.


/* Dane do wykresu s(x) vs sin(x) */
for (double x = -MathConstants.PI; x <= MathConstants.PI; x=x+0.02)
mXparser.consolePrintln("[ " + x +", " + MathFunctions.sin(x) + ", " + s.calculate(x) + " ],");

/* Dane do wykresu c(x) vs cos(x) */
for (double x = -MathConstants.PI; x <= MathConstants.PI; x=x+0.02)
mXparser.consolePrintln("[ " + x +", " + MathFunctions.cos(x) + ", " + c.calculate(x) + " ],");

Wniosek - proste zapisy rekurencyjne dają złożone wyniki! 🙂

Pozdrowienia,

Mariusz Gromada

Pamiętaj o:

import org.mariuszgromada.math.mxparser.*;
import org.mariuszgromada.math.mxparser.mathcollection.*;

Kod:

Zobacz również:

  1. Polowanie na czarownice - czyli zabawy z rekurencją (część 1)
  2. Prędkość ucieczki do nieskończoności - czyli zabawy z rekurencją (część 2)
  3. Naiwny test pierwszości - czyli zabawy z rekurencją (część 3)

Naiwny test pierwszości - czyli zabawy z rekurencją (część 3)

Liczby pierwsze

Jednym z najprostszych testów pierwszości jest weryfikacja czy dana liczba n posiada dzielnik z przedziału (2, \sqrt{n}) - takie podejście nazywane jest metodą naiwną - i niestety charakteryzuje się dużą złożonością obliczeniową. Nawet przy wykorzystaniu Sita Eratostenesa złożoność obliczeniowa sięga \frac{\sqrt{n}}{\log{n}}. Jednak w cyklu "Zabawy z rekurencją" nie bardzo zwracamy uwagę na złożoność 🙂 , bardziej chodzi o zobrazowanie jak całe algorytmy mogą być łatwo zapisane w postaci krótkich matematycznych funkcji rekurencyjnych - zatem do dzieła 🙂

Rekurencyjne poszukiwanie dzielników

Naszym zadaniem będzie zdefiniowanie funkcji zwracającej 1 jeśli podana liczba n jest liczbą pierwszą oraz 0 w przeciwnym wypadku. Zacznijmy jednak od podania funkcji weryfikującej czy liczba posiada dzielniki.

{\small\text{CzyDzielnik}(n, a, b)=}

{\small=\begin{cases}0&\text{dla}\quad a>b\\1&\text{dla}\quad n \mod a=0\\ \text{CzyDzielnik}(n, a+1, b)&\text{w inn. przyp.}\end{cases}}

Powyższa funkcja zwraca 1 jeśli liczba n posiada dzielnik z przedziału (a,b), oraz 0 w przeciwnym wypadku. Następnie definiujemy wyrażenie reprezentujące naiwny test pierwszości.

{\small\text{CzyPierwsza}(n)=}

{\small=\begin{cases}0&\text{dla}\quad n<2\\ \neg\text{CzyDzielnik}(n,2,\sqrt{n})&\text{dla}\quad n>=2\end{cases}}

Rolą funkcji "CzyPierwsza" jest jedynie "wprawienie algorytmu w ruch" oraz zwrócenie negacji wyniki funkcji "CzyDzielnik". Proste prawda? 🙂 Sprawdźmy więc w mXparser czy to faktycznie działa.

/* Definicja funkcji rekurencyjnych */
Function CzyDzielnik = new Function("CzyDzielnik(n, a, b) = if( a>b, 0, if( n%a = 0, 1, CzyDzielnik(n, a+1, b) ) )");
Function CzyPierwsza = new Function("CzyPierwsza(n) = if( n<2, 0, ~CzyDzielnik(n, 2, sqrt(n)) )", CzyDzielnik);

/* Obliczenie i wyświetlenie wartości */
mXparser.consolePrintln( "CzyPierwsza(1) = " + CzyPierwsza.calculate(1) + ", czas oblicz. = " + CzyPierwsza.getComputingTime() + " s." );
mXparser.consolePrintln( "CzyPierwsza(2) = " + CzyPierwsza.calculate(2) + ", czas oblicz. = " + CzyPierwsza.getComputingTime() + " s." );
mXparser.consolePrintln( "CzyPierwsza(3) = " + CzyPierwsza.calculate(3) + ", czas oblicz. = " + CzyPierwsza.getComputingTime() + " s." );
mXparser.consolePrintln( "CzyPierwsza(4) = " + CzyPierwsza.calculate(4) + ", czas oblicz. = " + CzyPierwsza.getComputingTime() + " s." );
mXparser.consolePrintln( "CzyPierwsza(5) = " + CzyPierwsza.calculate(5) + ", czas oblicz. = " + CzyPierwsza.getComputingTime() + " s." );
mXparser.consolePrintln( "CzyPierwsza(6) = " + CzyPierwsza.calculate(6) + ", czas oblicz. = " + CzyPierwsza.getComputingTime() + " s." );
mXparser.consolePrintln( "CzyPierwsza(7) = " + CzyPierwsza.calculate(7) + ", czas oblicz. = " + CzyPierwsza.getComputingTime() + " s." );
mXparser.consolePrintln( "CzyPierwsza(8) = " + CzyPierwsza.calculate(8) + ", czas oblicz. = " + CzyPierwsza.getComputingTime() + " s." );
mXparser.consolePrintln( "CzyPierwsza(9) = " + CzyPierwsza.calculate(9) + ", czas oblicz. = " + CzyPierwsza.getComputingTime() + " s." );
mXparser.consolePrintln( "CzyPierwsza(10) = " + CzyPierwsza.calculate(10) + ", czas oblicz. = " + CzyPierwsza.getComputingTime() + " s." );

+ wynik

CzyPierwsza(1) = 0.0, czas oblicz. = 0.08 s.
CzyPierwsza(2) = 1.0, czas oblicz. = 0.03 s.
CzyPierwsza(3) = 1.0, czas oblicz. = 0.026 s.
CzyPierwsza(4) = 0.0, czas oblicz. = 0.022 s.
CzyPierwsza(5) = 1.0, czas oblicz. = 0.038 s.
CzyPierwsza(6) = 0.0, czas oblicz. = 0.015 s.
CzyPierwsza(7) = 1.0, czas oblicz. = 0.028 s.
CzyPierwsza(8) = 0.0, czas oblicz. = 0.015 s.
CzyPierwsza(9) = 0.0, czas oblicz. = 0.053 s.
CzyPierwsza(10) = 0.0, czas oblicz. = 0.011 s.

Wygląda na to, że obliczenia są poprawne! Teraz możemy zweryfikować ile jest liczb pierwszych w podanym przedziale, definiując

\pi(n)=\sum_{i=1}^n \text{CzyPierwsza}(i)

/* Definicja wyrażenia sumującego wynik funkcji CzyPierwsza */
Expression pi100 = new Expression("sum(i, 1, 100, CzyPierwsza(i) )");
pi100.addFunctions(CzyPierwsza);

/* Obliczenie i wyświetlenie wyniku */
mXparser.consolePrintln( "Liczba liczb pierwszych w przedziale (1,100) = " + pi100.calculate());

+ wynik

Liczba liczb pierwszych w przedziale (1,100) = 25.0

Pozdrowienia,

Mariusz Gromada

Pamiętajcie, że uruchamiając kody mXparsera należy dodać w nagłówku:

import org.mariuszgromada.math.mxparser.*;

Kod:

 

Zobacz również:

  1. Polowanie na czarownice - czyli zabawy z rekurencją (część 1)
  2. Prędkość ucieczki do nieskończoności - czyli zabawy z rekurencją (część 2)
  3. Rekurencja pośrednia - czyli zabawy z rekurencją (część 4)

Prędkość ucieczki do nieskończoności - czyli zabawy z rekurencją (część 2)

Dziś ciekawostka w nawiązaniu do wpisu z dnia 20 października 2015 roku "Liczba PI ukryta w zbiorze Mandelbrota" - przeanalizujmy równanie rekurencyjne dla liczb rzeczywistych:

x_n=\begin{cases}x_{n-1}^2+\frac{1}{4}+\epsilon,&\text{dla}\quad n>0\\0,&\text{dla}\quad n=0\end{cases}

Zbliżanie się do "ostrza" zbioru Mandelbrota

Powyższe wyrażenie powstaje na bazie równania zbioru Mandelbrota z_n=z_{n-1}^2+c (w liczbach zespolonych), jeśli ograniczymy się do prostej rzeczywistej (dlatego użyłem zapisu x_n) i będzie nas interesowało zbliżanie się elementu x_1=\frac{1}{4}+\epsilon do "ostrza" (ang. "cusp") zbioru o współrzędnych (\frac{1}{4},0).

Mandelbrot - Ostrze

Szybkość ucieczki do nieskończoności

Ustalając odpowiednio małe \epsilon>0 decydujemy jak bardzo chcemy się zbliżyć do "ostrza", teraz zadanie polega na znalezieniu pierwszego n, dla którego x_n>=2. Takie minimalne n jest dobrą miarą prędkości ucieczki x_n do nieskończoności w zależności od wybranego \epsilon. Na marginesie dodam, że zbiór Juli dla równania Mandelbrota jest na powyższym obrazku oznaczony kolorem czarnym, który reprezentuje punkty "nieuciekające" do nieskończoności w trakcie nieskończonej iteracji (ta tematyka jest sama w sobie bardzo ciekawa i zapewne kiedyś coś napiszę o atraktorach).

Rekurencja na rekurencji

Zapiszmy nasze zadanie wykorzystując rekurencję w celu poszukiwania rozwiązania:

f(n,\epsilon)=\begin{cases}f(n+1,\epsilon),&\text{dla}\quad x_n<2\\n,&\text{dla}\quad x_n>=2\end{cases}

Niestety zdefiniowaliśmy rekurencję na rekurencji - to zły znak dla wydajności - cóż - sprawdźmy to używając mXparsera.

/* Definicja funkcji rekurencyjnej */
Function x = new Function("x(n,c) = if( n>0, x(n-1,c)^2+0.25+c, 0 )");
Function f = new Function("f(n,c) = if( x(n,c) >= 2, n, f(n+1, c) )", x);

/* Obliczenia i wyświetlenie wyniku */
mXparser.consolePrintln( "c = 0.01" + ", f(0,c) = " + f.calculate(0, 0.01) + ", czas = " + f.getComputingTime() + " s" );
mXparser.consolePrintln( "c = 0.0001" + ", f(0,c) = " + f.calculate(0, 0.0001) + ", czas = " + f.getComputingTime() + " s" );
mXparser.consolePrintln( "c = 0.000001" + ", f(0,c) = " + f.calculate(0,0.000001) + ", czas = " + f.getComputingTime() + " s" );
mXparser.consolePrintln( "c = 0.00000001" + ", f(0,c) = " + f.calculate(0,0.00000001) + ", czas = " + f.getComputingTime() + " s" );

+ wyczekiwany wynik


c = 0.01, f(0,c) = 30.0, czas = 0.224 s
c = 0.0001, f(0,c) = 312.0, czas = 1.532 s
c = 0.000001, f(0,c) = 3140.0, czas = 37.343 s
c = 0.00000001, f(0,c) = 31414.0, czas = 4068.338 s


Wzorzec prędkości ucieczki

Wow - jaki przedziwny wzorzec liczby iteracji wymaganych, aby przekroczyć 2!! Dostajemy coś, co przypomina \pi, jednak wymaga postawienia przecinka w odpowiednim miejscu! Można również zauważyć, że 100-krotne zmniejszenie \epsilon zwiększa niezbędną liczbę iteracji około 10-krotnie. Zmniejszając \epsilon otrzymamy liczbą coraz bardziej przypominającą \pi 🙂

Czasy wykonania

Tak jak oczekiwaliśmy - czasy wykonania są koszmarne - głównie z dwóch powodów:

  1. wykorzystanie generycznych funkcji rekurencyjnych (możliwe dowolne definicje, jednak kosztem wydajności - konieczność klonowania obiektów, brak możliwości zapamiętywania wartości z poprzednich przebiegów)
  2. definicji rekurencji na rekurencji,

mXparser dostarcza również nieco wydajniejszą wersję rekurencji - tzw. argumenty rekurencyjne z indexem (klasa RecursiveArgument). Wykorzystanie  RecursiveArgument jest ograniczone do jednego indexu, co umożliwia przechowywanie rezultatów z poprzednich iteracji. Zyskujemy na wydajności, tracąc na elastyczności i przejrzystości kodu. Działanie RecursiveArgument zaprezentuję w kolejnej części.

Pozdrowienia,

Mariusz Gromada

BTW: uruchamiając kody mXparsera nie zapomnijcie o dodaniu

import org.mariuszgromada.math.mxparser.*;

Zobacz również:

  1. Polowanie na czarownice - czyli zabawy z rekurencją (część 1)
  2. Naiwny test pierwszości - czyli zabawy z rekurencją (część 3)
  3. Rekurencja pośrednia - czyli zabawy z rekurencją (część 4)
  4. Liczba PI ukryta w zbiorze Mandelbrota

Polowanie na czarownice - czyli zabawy z rekurencją (część 1)

Rozpocznijmy branżowym humorem

Okres średniowiecza, kobieta winna uprawiania magii, kara straszna - spalenie na stosie! Nadszedł dzień, tłum gawiedzi, czarownica na stosie, płomienie, wiedźma krzyczy - więcej drewna! Więcej drewna! Tłum zdziwiony, mimo wszystko spełnia ostatnie życzenie opętanej. Wiedźma nie przerywa - jeszcze więcej drewna! Więcej drewna! Z oddali dobiega nagły i stanowczy sprzeciw - STOP! Czarownica chce przepełnić stos!

🙂

Czym jest rekurencja?

Zazwyczaj o rekurencji myślimy jako o procesie podziału zadania na mniejsze, następnie podziału na jeszcze mniejsze, i jeszcze mniejsze ... dochodząc do zadań, dla których rozwiązanie jest znane. Od tego momentu zaczyna się składanie "mniejszych" rozwiązań w "większe", następnie tych większych w jeszcze większe, ... i w jeszcze większe ... kończąc na rozwiązaniu zadania początkowego. Dla przykładu zapiszmy funkcję n! w postaci rekurencyjnej.

n!=\begin{cases}n\cdot(n-1)!&\text{dla}\quad n>0\\1,&\text{dla}\quad n=0\end{cases}

W celu zobrazowania reprezentacja powyższego podana w mXparser:

/* Definicja funkcji rekurencyjnej */
Function silnia = new Function("s(n) = if( n>0, n*s(n-1), 1 )");

/* Obliczenia i wyświetlenie wyniku */
System.out.println( "n = 0, s(n) = " + silnia.calculate(0) );
System.out.println( "n = 1, s(n) = " + silnia.calculate(1) );
System.out.println( "n = 2, s(n) = " + silnia.calculate(2) );
System.out.println( "n = 3, s(n) = " + silnia.calculate(3) );
System.out.println( "n = 4, s(n) = " + silnia.calculate(4) );
System.out.println( "n = 5, s(n) = " + silnia.calculate(5) );

+ wynik:

n = 0, s(n) = 1.0
n = 1, s(n) = 1.0
n = 2, s(n) = 2.0
n = 3, s(n) = 6.0
n = 4, s(n) = 24.0
n = 5, s(n) = 120.0

Wynik jest zgodny z oczekiwanym. Innym przykład rekurencji to iterowany operator sumowania, niech

A_n=a_1+a_2+\ldots+a_n=\sum_{i=1}^n a_i

Łatwo zauważyć, że

A_n=\begin{cases}a_n+A_{n-1},&\text{dla}\quad n>1\\a_1,&\text{dla}\quad n=1\end{cases}

Jak więc widzicie, rekurencja jest powszechna, często będąc nieco innym sposobem patrzenia na iteracje.

Formalna definicja rekurencji

O rekurencji mówimy jeśli metoda (funkcja, zachowanie, obiekt) może być opisana przez:

  1. elementy bazowe / rozwiązania bazowe;
  2. zestaw reguł, które redukuję (sprowadzają) każdy inny przypadek do (w kierunku) elementów bazowych.

Rekurencja
Powyższe określenie jest szerokie, ale takie być musi, bo i typów rekurencji jest wiele.

Rekurencja jako złożenie funkcji

Jednym (ale nie jedynym) sposobem zapisu ogólnych równań rekurencyjnych jest złożenie funkcji:

f_n=\begin{cases}F\big(f_{n-1},f_{n-2},\ldots,f_{n-k}\big)&\text{dla}\quad n>k\\f_1,f_2,\ldots,f_k&\text{el. baz. dla}\quad n<=k\end{cases}

Dobrą ilustracją powyższego jest ciąg Fibonacciego:

f_n=\begin{cases}0&\text{dla}\quad n=0\\1&\text{dla}\quad n=1\\f_{n-1}+f_{n-2}&\text{dla}\quad n>1\end{cases}

Zapiszmy ciąg Fibonacciego w mXparser:

/* Definicja funkcji rekurencyjnej */
Function fib = new Function("fib(n) = if( n>1, fib(n-1)+fib(n-2), if(n=1,1,0) )");

/* Obliczenia i wyświetlenie wyniku */
System.out.println( "fib(0) = " + fib.calculate(0) );
System.out.println( "fib(1) = " + fib.calculate(1) );
System.out.println( "fib(2) = " + fib.calculate(2) );
System.out.println( "fib(3) = " + fib.calculate(3) );
System.out.println( "fib(4) = " + fib.calculate(4) );
System.out.println( "fib(5) = " + fib.calculate(5) );

+ rezultat:

fib(0) = 0.0
fib(1) = 1.0
fib(2) = 1.0
fib(3) = 2.0
fib(4) = 3.0
fib(5) = 5.0

Rekurencja w roli pętli "For"

Podane wyżej przykłady zapisów rekurencyjnych (n!, suma n-pierwszych wyrazów ciągu, ciąg Fibonacciego) są tak naprawdę rekurencyjną realizacją pętli "for" - znamy przecież dokładnie liczbę niezbędnych operacji do wykonania, a i same operacje są raczej łatwe oraz czytelne - zatem zagnieżdżenie ich w pętli "for" nie powinno spowodować utraty przejrzystości kodu.

Poszukiwanie rozwiązania - czyli rekurencja w roli pętli "While/Until"

W metodach numerycznych często stosuje się strategie rekurencyjne - w takiej sytuacji, będąc w kroku n, weryfikujemy czy propozycja rozwiązania n spełnia kryterium stopu (np. jakość oszacowania), jeśli tak - kończymy z rozwiązaniem n, jeśli nie - przechodzimy do badania propozycji rozwiązania n+1. Procedurę rozpoczynamy od kroku 0 (zerowego).

Przykład: znając definicję silni chcemy znaleźć pierwsze n, dla którego n! >= 100 - takie zadanie formalnie możemy zapisać jako:

S_n=\begin{cases}S_{n+1},&\text{dla}\quad n!<100\\n,&\text{dla}\quad n!>=100\end{cases}

n_{100} = S(0)

Reprezentacja w mXparser:

/* Definicja funkcji rekurencyjnej */
Function S = new Function("S(n) = if( n! < 100, S(n+1), n )"); /* Obliczenia i wyświetlenie wyniku */ System.out.println( "Pierwsze n, że n! >= 100 to n = " + S.calculate(0) );

+ wynik:

Pierwsze n, że n! >= 100 to n = 5.0

Rekurencja bezpośrednia

Wszystkie omawiane wyżej typy rekurencji polegają na wywołaniu z ciała funkcji "siebie samej", dlatego należą do bardziej ogólnej klasy nazywanej rekurencją bezpośrednią.

Wydajność

Implementacje na bazie rekurencji są bardzo czytelne, o często minimalnym rozmiarze kodu. Jednak coś za coś - tracimy bardzo dużo na złożoności obliczeniowej (powtarzane operacje, dzielenie zadań, operacje na stosie) oraz wymogach pamięci (głównie struktura stosu) - to właśnie dlatego czarownica błagała o drewno - licząc na przerwanie procesu z tytułu przepełnienia stosu 🙂

Cdn 🙂

Pozdrowienia,

Mariusz Gromada

BTW: jeśli planujecie uruchamiać kody mXparsera nie zapomnijcie o dodaniu

import org.mariuszgromada.math.mxparser.*;

Zobacz również:

  1. Prędkość ucieczki do nieskończoności - czyli zabawy z rekurencją (część 2)
  2. Naiwny test pierwszości - czyli zabawy z rekurencją (część 3)
  3. Rekurencja pośrednia - czyli zabawy z rekurencją (część 4)

 

mXparser - wersja 1.0.2 dostępna do pobrania

mXparser - wysoce elastyczny parser (interpreter) wyrażeń matematycznych dla JAVA oraz C# .NET

mXparser - parser matematyczny

mXparser - wersja 1.0.2 dostępna do pobrania

Zmiana w stosunku do 1.0.1 to poprawa kontroli składni w funkcji definiowanej przez użytkownika. Do wersji 1.0.1 użycie argumentu rekurencyjnego w funkcji definiowanej przez użytkownika powodowało zwiększenie liczby oczekiwanych parametrów wywołania funkcji. Wartość argumentu rekurencyjnego jest wyznaczana automatycznie, stąd w kontroli składni (weryfikując liczbę podanych parametrów) taki argument należy pominąć. Od wersji 1.0.2 liczba parametrów funkcji definiowanej przez użytkownika opiera się na bazie liczby argumentów pomniejszonej o argumenty rekurencyjne.

Przykład

Poniżej przykład demonstrujący działanie wprowadzonej zmiany - proszę przesunąć tekst w celu uwidocznienia komentarzy.


import org.mariuszgromada.math.mxparser.*;

public class mXparserTests {
  public static void main(String[] args) {
    RecursiveArgument z = new RecursiveArgument("z","z(n-1)^2+c", "n"); /* Definicja argumentu rekurencyjnego */
    z.addBaseCase(0, 0);                                                /* Definicja elementu początkowego */
    Constant c = new Constant("c", 0.3);                                /* Deklaracja stałej występującej w równaniu rekurencyjnym */
    z.addConstants(c);                                                  /* Wskazanie argumentowi oczekiwanej stałej 'c' */
    Function diff1 = new Function("diff1", "z(n) - z(n-1)", "n");       /* Definicja funkcji o parametrze n opartej również na argumencie rekurencyjnym */
    diff1.addArguments(z);                                              /* Wskazanie funkcji oczekiwanego argumentu 'z' */
    System.out.println(diff1.calculate(10));                            /* Obliczenie wartości funkcji + wyświetlenie wyniku */
  }
}

Rezultat


0.17863353779409574

 

Pozdrowienia,

Mariusz Gromada

mXparser - wysoce elastyczny parser (interpreter) wyrażeń matematycznych dla JAVA oraz C# .NET

mXparser - parser matematyczny

Serdecznie zapraszam do zapoznania się z wysoce elastycznym interpreterem wyrażeń matematycznych. Oprogramowanie jest mojego autorstwa, powstało w 2010 roku i wtedy zostało opublikowane w serwisie SourceForge.net. Z racji, że teraz posiadam stronę o odpowiedniej tematyce, zdecydowałem się przygotować dedykowany opis, który znajdziecie pod tym linkiem. Dostępne są również tutorial oraz specyfikacja API.

Pozdrowienia,

Mariusz Gromada

Pobierz mXparser - parser matematyczny

Download mXparser

Download mXparser