szukanie zaawansowane
 [ Posty: 18 ]  Przejdź na stronę Poprzednia strona  1, 2
Autor Wiadomość
Mężczyzna Offline
PostNapisane: 5 gru 2017, o 20:07 
Moderator

Posty: 3938
Lokalizacja: Kraków PL
  1. SlotaWoj napisał(a):
    Wartość drugiego równania zależy od pochodnej wartości czwartego równania. Nie wszystkie implementacje algorytmu Rungego-Kuty dopuszczają taką możliwość w sensie definicyjnym i proceduralnym.
    Może mieć to związek z tym co pisałem uprzednio.
  2. Dlaczego ostatnie współczynniki modelu są inne niż te, które wyliczyłem? Co się zmieniło?
    Dla jakich wartości współczynników otrzymałeś w WorkingModel pokazane przebiegi?
  3. Jeśli będę miał czas, to może dzisiaj zamodeluję Twój układ w Pascalu i zrobię wykresy w Excelu.


-- czwartek, 7 grudnia 2017, 02:21 --

Tak jak pisałem wcześniej, zamodelowałem w Pascalu Twój układ w postaci:

    \begin{cases}\ s_1'=s_2 \\
\ s_2'=-6,153846s_1+0,2{\red{s_4'}}-0,346154s_4 \\
\ s_3'=s_4 \\
\ s_4'=2,197802s_1+0,123626s_4-6,185714\end{cases}

z warunkami początkowymi:

    \begin{cases}
\ s_1(0)=0 \\
\ s_2(0)=0 \\
\ s_3(0)=\pi/3 \\
\ s_4(0)=0\end{cases}

zaimportowałem otrzymane wyniki do Excela i zrobiłem wykresy. Są one tak rażąco niepodobne do tego, co wyświetlił WorkingModel, więc nawet ich nie będę tu publikował.

  1. Podaj mi proszę dokładną postać równań modelu, które wprowadziłeś do WorkingModel.

    Obrazek

  2. Dlaczego w MathCADzie stała \pi/3 jest jako druga współrzędna wektora y , a nie jako trzecia?
  3. Dlaczego wektor pochodnych jest oznaczona przez D(x,\varphi) , a nie przez D(x,x',\varphi,\varphi') ? Skąd MathCAD wie, że x i \varphi są wektorami?
  4. Skąd w ogóle się wziął ten wektor? Powyżej widzę ostatnie równanie modelu, zakodowane „po bożemu”, tzn. \ddfrac{s_4}{t}=f(t,\mathbf{s}) (t nie jest używane).
  5. Co to jest to drugie równanie w wektorze D ? Skąd się ono wzięło?
  6. Dlaczego w ww. równaniu występuje pierwsza pochodna z \varphi , a nie druga?
  7. Dlaczego wywołałeś w MathCADzie metodę rkfixed z parametrami t_\text{start}=-180 i t_\text{stop}=0,5 , a warunki początkowe masz dla t=0 ?
  8. Nie będę zmieniał tego, co napisałem w pkt. 3 i 4, ale teraz zauważyłem, że masz błąd w definicji wektora pochodnych. Ma ona mieć postać D(t,\textbf{s}) , gdzie t to zmienna niezależna, czyli czas, a \textbf{s} to wektor stanu, czyli \left[x,x',\varphi,\varphi'\right]^\text{T} .
Góra
Mężczyzna Offline
PostNapisane: 7 gru 2017, o 18:02 
Użytkownik

Posty: 7
Lokalizacja: Katowice
@slotawoj
Możesz mi podesłać na priv ten wykres zerknę chociaż jak on mniej więcej wygląda.

Ad. 1 Program "Working model" działa na modelach fizycznych. To oznacza, że żadnych równań początkowych nie wprowadzam. Układ w nim projektuje się za pomocą brył i więzów oraz ich właściwości. Równanie, które wypisałem na samym początku było sprawdzone z literaturą, więc pochodne i zmienne powinny być okej (wykres mniej więcej podobny)

Ad. 2 Jest to mój błąd wprowadzę tą poprawkę.
Ad. 3 MathCad nie przyjmuje pochodnych w opisie. Spróbuję pochodną oznaczyć inną zmienną i wrzucić do D.
Ad. 4 Wektor jest opisany twoimi równaniami jednak zastąpiłem wszystkie zmienne S wcześniejszymi parametrami, aby mógł ruszyć dalej.
Ad. 5 Jak wyżej. Podstawione na bezczela s4.
Ad. 6 Mój błąd.
Ad. 7 Funkcji rkfixed nie znam (jak i MathCad'a). Skorzystałem z wcześniejszego opisu Janusza błędnie interpretując iż chodzi o kąt. Poprawię to.
Ad. 8 Spróbuję tak zrobić i wrócić do wcześniejszego wektora D.
Góra
Mężczyzna Offline
PostNapisane: 7 gru 2017, o 20:24 
Moderator

Posty: 3938
Lokalizacja: Kraków PL
Cavanus napisał(a):
Ad.5 Jak wyżej. Podstawione na bezczela s4.
Sprawdziłem. To nie jest równanie na s_2' , ani nawet proporcjonalne do niego (przeskalowane).
Góra
Utwórz nowy temat Odpowiedz w temacie  [ Posty: 18 ]  Przejdź na stronę Poprzednia strona  1, 2


 Zobacz podobne tematy
 Tytuł tematu   Autor   Odpowiedzi 
 Do użytkowników i znawców MathCAD'a!  Anonymous  2
 Wykres powierzchni f(x,y)=x*y  michael_94  4
 Mathematica 5.2 wykres krzywej w 3D  Alfredq  0
 różne wyniki z takiej samej formuły (MathCad)  kefirek  2
 program do równań (istnieje taki wogóle)?  matekleliczek  1
 
Atom [Regulamin Forum] [Instrukcja LaTeX-a] [Poradnik] [F.A.Q.] [Reklama] [Kontakt]
Copyright (C) Karpatka.pl