Примеры S-функций на основе М-файлов

S-функции многотактовых дискретных и дискретно-непрерывных систем

S-функции могут описывать дискретно-непрерывные системы, а также многотактовые дискретные системы, которые имеют различные периоды квантования. В общем случае, рекомендуется разрабатывать S-функции, ориентированные на описание моделей конкретного типа. Это упрощает обслуживание и формирование S-функций, а также делает их более приспособленными к многократному использованию в других проектах.

Приведем примеры, показывающие, как создать S-функции, описывающие непрерывную, дискретную системы, а также модель, задаваемую алгебраическим уравнением. Другие примеры, входящие в комплект поставки системы Simulink, можно найти в директории MATLAB/toolbox/simulink/simdemos.

Пример1: S-функция на основе М-файла для модели, заданной алгебраическим уравнением .

S-функция имеет два входа, один выход и не является моделью в переменных состояния. Поскольку между входом и выходом отсутствуют динамические звенья, то данная S-функция определяет прямую передачу между входом и выходом. Структурная схема системы Simulink, включающая в себя S-функцию, представлена на рисунке 10.

Рисунок 10 – Структурная схема системы Simulink

В диалоговом окне блока Block Parameters: S-Function в текстовом поле S-function name содержится имя функции s_xmpl. Код S-функции представлен ниже

function[sys,xO,str,ts]=s_xmpl(t, x, u, flag)

switch flag,

case 0,

[sys,x0,str,ts]=mdlInitializeSizes;

case 1,

sys=mdlDerivatives(t,x,u);

case 2,

sys=mdlUpdate(t,x,u);

case 3,

sys=mdl0utputs(t,x,u);

case 4,

sys=mdlGetTimeOfNextVarHit(t,x, u);

case 9,

sys=mdlTerminate(t,x,u);

otherwise, %0шибочный вход.

error(['опущен flag= ',num2str(flag) ]);

end

function[sys,xO,str,ts]=mdlInitializeSizes ()

size=simsizes;

sizes.NumContStates= 0;

sizes.NumDiscStates= 0;

sizes.NumOutputs= 1;

sizes.Numlnputs= 2;

sizes.DirFeedthrough=l;

sizes.NumSampleTimes=l;

sys = simsizes (sizes);

х0 = []; % Устанавливает начальное значение для всех

% переменных состояния.

str = []; % str всегда пустая матрица

ts = [0 0]; % Инициализирует матрицу периодов квантования.

%***********************************************************

function sys = mdlDerivatives(t, x, u)

sys = []; % Пустой, т.к. S-функция не содержит переменных

% состояния непрерывной части системы.

%***********************************************************

function sys = mdlUpdate(t, x, u)

sys =[]; % Пустой, т.к. модель не содержит переменных

% состояния дискретной части системы.

%***********************************************************

function sys = mdlOutputs(t,x,u)

sys = [u(l)+u(2).^2];

%***********************************************************

% mdlGetTimeOfNextVarHit

function sys=mdlGetTimeOfNextVarHit(t,x,u)

% Возвращает следующее значение периода квантования для

% данного блока. Результат – абсолютное время.

% Эта функция используется, только когда

% период квантования определен как переменная величина.

sampleTime=[];

%****************************************************************

% mdlTerminate

function sys = mdlTerminate(t,x,u)

sys=[];

Дата добавления: 2014-01-11