Zaloguj się
Blog na Matlab.pl
Forum polskich użytkowników
 
UŻYTKOWNICY GRUPY PROFIL Zaloguj się, by sprawdzić wiadomości FAQ
 



Napisz nowy temat     Odpowiedz do tematu Zobacz poprzedni temat :: Zobacz następny temat

Matlab portret fazowy
Forum MATLAB Strona Główna-> MATLAB
Post Wysłany: 23 Maj 2012, Sro 7:18 pm Temat postu: Matlab portret fazowy Odpowiedz z cytatem
 
AUTOR:
dexter90
Początki pisania


Dołączył: 14 Maj 2012
Posty: 17


Ogląda profil użytkownika Wyślij prywatną wiadomość
Witam

W jaki sposób mogę zrobić portret fazowy?? Szukam w internecie, ale mało znajduje konkretów. Mam na myśli portret fazowy coś ala modelu Lotki Voltery


 

Post Wysłany: 23 Maj 2012, Sro 7:45 pm Temat postu: Odpowiedz z cytatem
 
AUTOR:
mc2
Może pisać książki


Dołączył: 19 Maj 2009
Posty: 2056


Ogląda profil użytkownika Wyślij prywatną wiadomość Wyślij email Odwiedź stronę autora
Po prostu użyj komendy plot. Reszta to kwestia przygotowania danych, np:
Kod:
[t,y]=ode45(@(t,y)diag([1 - .01*y(2), -1 + .02*y(1)])*y,[0 15],[20 20])
plot(y(:,1),y(:,2))
 

Post Wysłany: 23 Maj 2012, Sro 8:21 pm Temat postu: Odpowiedz z cytatem
 
AUTOR:
dexter90
Początki pisania


Dołączył: 14 Maj 2012
Posty: 17


Ogląda profil użytkownika Wyślij prywatną wiadomość
Kod:
[t,y]=ode45(@(t,y)diag([-0.02*y(1)+y(1)*y(1)*y(2)+0.01, 0.02*y(1)-y(1)*y(1)*y(2)])*y,[0 15],[0 0.2])
plot(y(:,1),y(:,2))


Nieznacznie przerobiłem Twój kod, chciałem wklepać tam swoje równanie, ale mi nie wyszło. Nie ukrywam, że nie mam rzetelnej wiedzy, dopiero wypożyczyłem książke o Matlabie i mam zamiar się w nią wgłębić...


 

Post Wysłany: 23 Maj 2012, Sro 9:48 pm Temat postu: Odpowiedz z cytatem
 
AUTOR:
mc2
Może pisać książki


Dołączył: 19 Maj 2009
Posty: 2056


Ogląda profil użytkownika Wyślij prywatną wiadomość Wyślij email Odwiedź stronę autora
Zrób sobie przykłady z tego poradnika, wtedy Ci się rozjaśni:
www.igf.fuw.edu.pl/ZOI/Matlab/Lessons/MatlabII/tm2.pdf

Ja w przykładzie użyłem może trochę skomplikowanego zapisu z wykorzystaniem funkcji anonimowej, ale pytałeś o diagram fazowy, a nie o metodę rozwiązania układu równań różniczkowych Smile
 

Post Wysłany: 23 Maj 2012, Sro 10:00 pm Temat postu: Odpowiedz z cytatem
 
AUTOR:
dexter90
Początki pisania


Dołączył: 14 Maj 2012
Posty: 17


Ogląda profil użytkownika Wyślij prywatną wiadomość
Hmmm nic tam nie mogę znaleźć na temat portretu fazowego.


 

Post Wysłany: 24 Maj 2012, Czw 6:55 am Temat postu: Odpowiedz z cytatem
 
AUTOR:
mc2
Może pisać książki


Dołączył: 19 Maj 2009
Posty: 2056


Ogląda profil użytkownika Wyślij prywatną wiadomość Wyślij email Odwiedź stronę autora
Chyba się nie rozumiemy Smile
W poradniku masz opisany sposób rozwiązywania równań różniczkowych.
Portret fazowy, to zwykły wykres w matlabie, który wykonasz za pomocą polecenia plot.
Ty musisz najpierw rozwiązać jakieś równanie różniczkowe, a następnie narysować odpowiedni wykres, czyli nie rysujesz np. zmian populacji w czasie, tylko zmianę jednej populacji względem liczebności drugiej populacji.
 

Post Wysłany: 24 Maj 2012, Czw 8:03 am Temat postu: Odpowiedz z cytatem
 
AUTOR:
dexter90
Początki pisania


Dołączył: 14 Maj 2012
Posty: 17


Ogląda profil użytkownika Wyślij prywatną wiadomość
Równanie różniczkowe normlanie już mam rozwiązane, jak mówie nie wiem jak to ugryźć - brak doświadczenia :/

Oto mój układ równań:

Kod:
x=y(1);
y=y(2);
dxdy(1,1)=-0.02*x+x*x*y+0.01;
dxdy(2,1)=0.02*x-x*x*y;


Bardzo proszę o pomoc. Jeżeli mam go rozwiązać najpierw analitycznie proszę pisać.

EDIT:

Znalazłem skrypt, być może innym forumowiczom pomoże, natomiast jakoś nie satysfakcjonuje mnie ten wykres..

http://math.uchicago.edu/~emstrawb/Phase_how_to.pdf


 

Forum MATLAB Strona Główna-> MATLAB
Wyświetl posty z ostatnich:   

Napisz nowy temat     Odpowiedz do tematu Zobacz poprzedni temat :: Zobacz następny temat

Wszystkie czasy w strefie CET (Europa)

Skocz do:  

Statystyki forum:



Od dnia 08.06.2006 forum odwiedzano 20376293
Najwięcej użytkowników 90 było obecnych 29 Maj 2010, Sob 7:37 pm

Aktualnie online:




Najnowsze posty na forum:
Macierz zerojedynkowa o różnych wierszach  (26 Października 2014, Nie 12:19 am)
rownanie różniczkowe  (25 Października 2014, Sob 11:26 pm)
Równanie różniczkowe warunki początkowe - schemat  (25 Października 2014, Sob 8:15 pm)
Mini edytor tekstu  (25 Października 2014, Sob 6:24 pm)
Rekurencja w Matlabie  (24 Października 2014, Pią 1:15 pm)
Wygenerowanie sygnału PRBS, a następnie obróbka.  (24 Października 2014, Pią 11:25 am)
Wygenerowanie sygnału PRBS, a następnie obróbka.  (24 Października 2014, Pią 11:24 am)
error using => plus  (23 Października 2014, Czw 11:01 pm)
Podziałka a etykieta - problem  (23 Października 2014, Czw 1:42 pm)
Interpolacja funkcji  (21 Października 2014, Wto 8:09 pm)
Twoje prawa:
Nie możesz pisać nowych tematów
Nie możesz odpowiadać w tematach
Nie możesz zmieniać swoich postów
Nie możesz usuwać swoich postów
Nie możesz głosować w ankietach
Nie możesz załączać plików na tym forum
Nie możesz ściągać plików na tym forum