Я новичок в Modelica, и мне интересно, можно ли написать своего рода уравнение динамического программирования. Предположим, что время дискретизировано целым числом i, и в моем конкретном приложении x является булевым значением, а f является булевой функцией от x.
x(t_i) = f(x(t_{i+d}))
Где d может быть положительным или отрицательным целым числом. Конечно, я бы инициализировал x соответственно, либо true, либо false.
Любая помощь или ссылки будут с благодарностью!
Возможно. В Modelica дискретизация во времени обычно выполняется компилятором, вам нужно позаботиться об уравнениях (непрерывная динамика). В противном случае, если вы хотите генерировать события в дискретные моменты времени, вы можете сделать это с помощью операторов when. Я предлагаю вам взглянуть на Introduction to Object-Oriented Modeling and Simulation with OpenModelica ( Формат PDF, 6,6 МБ) – более свежий учебник (2012 г.) Питера Фрицсона. Существует раздел, посвященный дискретным событиям и гибридным системам, который должен пояснить, как реализовать ваши уравнения в Modelica. Ниже вы можете найти пример из этого руководства о модели прыгающего мяча, поскольку вы можете видеть, что дискретизация во времени не учитывается при написании ваших динамических уравнений. Таким образом, непрерывная модель мяча v=der(s), a=der(v) и дискретная часть внутри предложения when, которая обрабатывает контакт с землей:
model BouncingBall "the bouncing ball model"
parameter Real g=9.81; //gravitational acc.
parameter Real c=0.90; //elasticity constant
Real height(start=10),velocity(start=0);
equation
der(height) = velocity;
der(velocity)=-g;
when height<0 then
reinit(velocity, -c*velocity);
end when;
end BouncingBall;
Надеюсь, это поможет, Марко
Насколько я понимаю ваш вопрос, вы хотите использовать последние n оценки x для определения следующего значения x. Если это так, этот код показывает, как это сделать:
model BooleanHistory
parameter Integer n=10 "How many points to keep";
parameter Modelica.SIunits.Time dt=1e-3;
protected
Boolean x[n];
function f
input Integer n;
input Boolean past[n-1];
output Boolean next;
algorithm
next :=not past[1]; // Example
end f;
initial equation
x = {false for i in 1:n};
equation
when sample(0,dt) then
x[2:n] = pre(x[1:(n-1)]);
x[1] = f(n, x[2:n]);
end when;
end BooleanHistory;
