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: 1957


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: 1957


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: 1957


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 18393638
Najwięcej użytkowników 90 było obecnych 29 Maj 2010, Sob 7:37 pm

Aktualnie online:




Najnowsze posty na forum:
funkcja z dwoma argumętami  (15 Kwietnia 2014, Wto 6:08 pm)
Przykładowe kolokwium z MatLaba nr 2  (14 Kwietnia 2014, Pon 11:16 pm)
Kolokwium z MatLaba  (14 Kwietnia 2014, Pon 10:04 pm)
Problem z przekazywaniem danych do funkcji.  (14 Kwietnia 2014, Pon 9:34 pm)
Odczytytwanie wartości z wykresu  (14 Kwietnia 2014, Pon 6:05 pm)
odległości między punktami  (14 Kwietnia 2014, Pon 5:16 pm)
Scilab Kombinacje Kolumn Macierzy  (14 Kwietnia 2014, Pon 5:00 pm)
Czy pochodna jest liczbą  (13 Kwietnia 2014, Nie 5:39 pm)
Kinematyka odwrotna.  (13 Kwietnia 2014, Nie 5:36 pm)
Usunięcie tekstu 'ans'  (13 Kwietnia 2014, Nie 4:21 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