Включение функции, зависящей от времени в численном difeq

Question

В попытке решить дифференциальное уравнение показано ниже

function dr = ball(t,b) 

T2 = 1; 
T1 = T2/2; 
d = 0; 
R = 0.2; 

dr = zeros(3,1); 
dr(1) = (-1/T2)*b(1)-d*b(2); 
dr(2) = (-1/T2)*b(2) + d*b(1) + R*b(3); 
dr(3) = (-1/T1)*b(3) - R*b(2) ; 

end 

Я сделал это с помощью ode45. Теперь я хочу сделать это, когда мое R не является функцией t. Как бы включить это в настройку дифференциального уравнения.

Answer 1

Вы можете использовать вход t, который дает вам текущий временной интервал и определяет R, что бы вы ни хотели.

function dr = ball(t,b) 

T2 = 1; 
T1 = T2/2; 
d = 0; 
R = sin(t); %// for example 

dr = zeros(3,1); 
dr(1) = (-1/T2)*b(1)-d*b(2); 
dr(2) = (-1/T2)*b(2) + d*b(1) + R*b(3); 
dr(3) = (-1/T1)*b(3) - R*b(2) ; 

end