odedc

Аргументы

y0

вещественный вектор-столбец (начальные условия), y0=[y0c;y0d], где y0d имеет nd элементов.

nd

целое число, размер y0d

stdel

вещественный вектор с одним или двумя элементами stdel=[h, delta] (с delta=0 в качестве значения по умолчанию).

t0

вещественный скаляр (начальное время).

t

вещественный вектор-(строка), моменты времени, в которые вычислено yt.

f

внешняя функция, то есть функция или символьная строка или список с последовательностью вызова: yp = f(t,yc,yd,flag).

a list

Эта форма внешней функции используется для передачи параметров в функцию. Она должна быть в следующем виде:

list(f, p1, p2,...)

где синтаксис функции f теперь

yp = f(t, yc, yd, flag, p1, p2,...)

f по-прежнему возвращает значение функции в виде функции (t, yc, yd, flag, p1, p2,...), а p1, p2,... являются параметрами функции.

символьная строка

она должна ссылаться на имя модуля C или fortran, в предположении, что <f_name> является указанным именем.

• Последовательность вызова Fortran должна быть

  <f_name>(iflag, nc, nd, t, y, ydp)

  double precision t, y(*), ydp(*)

  integer iflag, nc, nd

• C-синтаксис должен быть

  void <f_name> (int *iflag, int *nc, int *nd, double *t, double *y, double *ydp)

Как в случае Fortran, так и в случае C входными аргументами являются:

• iflag = 0 или 1

• nc = число непрерывных состояний yc

• nd = число дискретных состояний yd

• t = время

• y = [yc; yd; param]. param может быть использован для получения дополнительных аргументов, которые были указаны при вызове odedc (y = odedc([y0c; y0d], nd, stdel, t0, t, list('fexcd', param)))

• В качестве выходного аргумента ydp, процедура должна вычислять ydp[0:nc-1] = d/dt ( yc(t) ) для iflag=0 и ydp[0:nd-1] = yd(t+) для iflag=1.

Описание

y = odedc([y0c;y0d], nd, [h,delta], t0, t, f) вычисляет решение смешанной дискретно-непрерывной системы. Состояние дискретной системы ydk вложено в кусочно-постоянную функцию времени yd(t) в виде: yd(t) = ydk для t в [tk = h*(delta+k) , tk+1 = h*(delta+k+1) ).

Теперь моделируемые уравнения имеют вид, для t в [ t_k , t_k+1 ):

dyc/dt = f(t, yc(t), yd(t), 0)
yc(t0) = y0c

и в моменты времени tk дискретная переменная yd обновляется как

yd(tk+) = f(yc(tk-), yd(tk-), 1).

Заметьте, что, используя определение yd(t), последнее уравнение даёт

ydk = f(tk, yc(tk-), yd(tk-1), 1)

(yc непрерывно во времени : yc(tk-) = yc(tk)).

Параметры вызова f фиксированы: ycd = f(t,yc,yd,flag); эта функция должно возвращать либо производную вектора yc, если flag=0, либо обновлять yd, если flag=1.

ycd = dot(yc) должна быть вектором того же размера, что и yc, если flag=0 и ycd = update(yd) должна быть вектором того же размера, что и yd, если flag=1.

t -- это вектор моментов времени, в которых вычисляется решение y.

y -- это вектор y = [y(t(1)), y(t(2)),...].

Эта функция может быть вызвана с теми же самыми необязательными параметрами, что и функция ode (указанные nd и stdel даны в последовательности вызова как второй и третий параметры). В частности, могут быть установлены флаги интегрирования, допуски. Необязательные параметры могут быть установлены с помощью функции odeoptions.

Внешние процедуры могут быть динамически скомпонованы (см. link).

Примеры

//Линейная система с переключающимися входами
deff('xdu=phis(t,x,u,flag)','if flag==0 then xdu=A*x+B*u; else xdu=1-u;end');
x0=[1;1];A=[-1,2;-2,-1];B=[1;2];u=0;nu=1;stdel=[1,0];u0=0;t=0:0.05:10;
xu=odedc([x0;u0],nu,stdel,0,t,phis);x=xu(1:2,:);u=xu(3,:);
nx=2;
plot2d(t, x', [1:nx], '161')
plot2d2('onn',t',u',[nx+1:nx+nu],'000');
//Внешняя процедура Fortran (см. fydot2.f):
norm(xu-odedc([x0;u0],nu,stdel,0,t,'phis'),1)

//Дискретизированная обратная связь
//
//      |    xcdot=fc(t,xc,u)
//  (система)      |
//      |      y=hc(t,xc)
//
//
//      |     xd+=fd(xd,y)
//  (обр. связь)   |
//      |      u=hd(t,xd)
//
deff('xcd=f(t,xc,xd,iflag)',...
  ['if iflag==0 then '
   '  xcd=fc(t,xc,e(t)-hd(t,xd));'
   'else '
   '  xcd=fd(xd,hc(t,xc));'
   'end']);
A=[-10,2,3;4,-10,6;7,8,-10];B=[1;1;1];C=[1,1,1];
Ad=[1/2,1;0,1/20];Bd=[1;1];Cd=[1,1];
deff('st=e(t)','st=sin(3*t)')
deff('xdot=fc(t,x,u)','xdot=A*x+B*u')
deff('y=hc(t,x)','y=C*x')
deff('xp=fd(x,y)','xp=Ad*x + Bd*y')
deff('u=hd(t,x)','u=Cd*x')
h=0.1;t0=0;t=0:0.1:2;
x0c=[0;0;0];x0d=[0;0];nd=2;
xcd=odedc([x0c;x0d],nd,h,t0,t,f);
norm(xcd-odedc([x0c;x0d],nd,h,t0,t,'fcd1')) // Быстрое вычисление (см. fydot2.f)
plot2d([t',t',t'],xcd(1:3,:)');
scf(2);
plot2d2("gnn",[t',t'],xcd(4:5,:)');
scf(0);

Смотрите также

• ode — программа решения обыкновенных дифференциальных уравнений
• odeoptions — установка опций для программ решения ОДУ
• link — dynamic linker
• csim — simulation (time response) of linear system
• external — объект Scilab'а, внешняя функция или подпрограмма