Rysowanie gęstości
Sposób na rysowanie wykresu funkcji dwuwymiarowej.
Jak dodać paczkę do OCTAVE
Istnieją na to dwa sposoby:
Instalacja z Octave Forge. W repozytorium Octave Forge wyszukaj paczkę, którą chcesz zainstalować np. "control". Aby dokonać instalacji uruchom Octave i wpisz: "pkg install -forge control". Paczka powinna automatycznie ściągnąć się z repozytorium i zainstalować.
Instalacji można dokonać również ręcznie, ściągając plik: "control-***.tar.gz" (gwiazdki oznaczają wersję) i wpisując "pkg install control-***.tar.gz". Uwaga! Octave musi być otwarte w katalogu w którym znajduje się ściągnięta paczka.
Zainstalowane paczki można wyświetlić poleceniem: "pkg list". Paczki pozbyć się można poleceniem "pkg uninstall control"
Jak pomnożyć wielomiany w Octave
W środowisku Octave wielomian jest reprezentowany przez współczynniki funkcji wielomianowej poczynając od najbardziej znaczącego. Na przykład wektor: p=[3 -2 4 0 1]; reprezentuje funkcję wielomianową: p(x)=3*x^4-2*x^3+4*x^2+1 Możemy teraz wyświetlić sobie funkcje wielomianową następującym poleceniem: polyout(p,'x') 3*x^4 - 2*x^3 + 4*x^2 + 0*x^1 + 1 Mnożenie dwóch wielomianów p,q możemy wykonać poleceniem: r = conv(p,q);. Dzielenie natomiast poleceniem: [b,r] = deconv(y,a);. Gdzie b-wynik dzielenia wielomianów, r-reszta z dzielenia.
Jak numerycznie rozwiązać równanie różniczkowe?
Octave pozwala na rozwiązanie wektorowych równań różniczkowych w postaci dx/dt = f(t,x).
- Pierwszym krokiem jest definicja funkcji f(x,t). Załóżmy, że f(x,t)=sin(t).
octave:1> function x_pirm = f(x,t) > x_prim = sin(t); > endfunction
- Definiujemy czas rozwiązania: t=[0:0.01:6*pi]; lub t=linspace(0,6*pi,2000);
- Definiujemy warunki początkowe x_0=1;
- Rozwiązujemy równanie poleceniem y=lsode("f",x_0,t);
- Rysujemy rozwiązanie: plot(t,y)
Jak znaleźć zero funkcji?
Do rozwiązywanie nie liniowych równań postaci f(x)=0 służy funkcja fsolve.
- Najpierw musimy zdefiniować funkcję lub ich układ.
octave:1> function retval = wielo(x) > retval = x.^3-x.^2 +4; > endfunction
- Kolejnym krokiem jest wybór wartości początkowych x_0. Niech x_0 = -1.
- Możemy spróbować rozwiązać równanie poleceniem: fsolve(@wielo,x_0). Powinniśmy otrzymać
Warto zwrócić uwagę, że nie zawsze polecenie fsolve znajduje poprawne rozwiązanie. Jeśli położymy x_0 = 5, rozwiązanie okaże się błędne.
Jak zmienić grubość linii na wykresie?
Służy do tego parametr 'LineWidth' w funkcji plot. Na przykład: plot(x,y,'LineWidth',2)
Brak komentarzy:
Prześlij komentarz