Próbuję użyć ode45
rozwiązać system ODE to:Matlab: Czy możliwe jest numeryczne rozwiązanie układu odów z mieszaniną warunków początkowych i końcowych?
[X,Y]= ode45(@sys,[0, T],y0);
gdzie
function dy = sys(t,y)
dy(1) = f_1(y)
dy(2) = f_2(y)
dy(3) = f_3(y)
end
Problemem jest to, że funkcja ode45
wymaga y0
być wartościami początkowymi [y_1(0), y_2(0), y_3(0)]
, natomiast w moim system, mam dostępne tylko wartości [y_2(0), y_3(0), y_3(T)]
.
Matematycznie, ten zestaw warunków początkowych/końcowych powinien wystarczyć do ustalenia systemu, ale czy jest jakiś sposób, w jaki mogę z nim pracować przez ode45
lub jakąkolwiek inną funkcję w MATLAB?
Dzięki!
Naprawdę jestem zainteresowany tym pytaniem, ale obawiam się, że nie mogę pomóc; Nigdy nie spotkałem się z tego rodzaju problemem przed ... Wiem, że 'ode45' może integrować się wstecz (po prostu użyj' tspan = [tend tstart] '), więc możesz wymyślić iteracyjny schemat, aby uzyskać' y_1' w taki sposób, że 'y_3 (0) 'i' y_3 (T) 'są spełnione. Nie trzeba dodawać, że będzie to * bardzo * powolne i dość niezgrabne, ale będzie to ** rozwiązanie **. Będę miał na to oko :) Czy mógłbyś zamieścić równania i warunki początkowe/końcowe? –
@RodyOldenhuis Myślę, że znalazłem sposób na rozwiązanie tego problemu. – Vokram