Создать Simulink модель функции f(t), используя подсистемы.

sum(u)

2*pi

function [sys,x0,str,ts] = S_DUE(t,x,M,flag,J,UgSk0) 
% S-функция S_DUE Динамических Уравнений Ейлера 
% Реализует динамику вращательного движения твердого тела, 
% отыскивая вектор абсолютной угловой скорости тела 
% по заданному вектору моментов внешних сил, 
% действующих на тело 
% ВХОД блока: 
% M - вектор проекций момента внешних сил на оси 
% X, Y i Z связанной с телом системы координат 
% ВЫХОД блока: 
% y - вектор из шести элементов: первые три - проекции 
% абсолютной угловой скорости om тела на указанные оси, 
% последние три - проекции на те же оси 
% углового ускорения тела 
% Входные ПАРАМЕТРЫ S-функции: 
% J - матрица моментов инерции тела в указанных осях; 
% UgSk0 - вектор начальных значений проекций 
% угловой скорости тела 
% Лазарев Ю.Ф., Украина, 18-12-2001 
IJ=inv(J); % вычисление обратной матрицы моментов инерции 
switch flag, 
case 0 
[sys,x0,str,ts] = mdlInitializeSizes(UgSk0); 
case 1, 
sys = mdlDerivatives(t,x,M,J,IJ); 
case 3, 
sys = mdlOutputs(x,M,J,IJ); 
case 9 
sys = [0]; 
end 
% Конец процедуры 
%======================================================= 
function [sys,x0,str,ts] = mdlInitializeSizes(UgSk0,Ug0) 
sizes = simsizes; 
sizes.NumContStates = 3; 
sizes.NumDiscStates = 0; 
sizes.NumOutputs = 6; 
sizes.NumInputs = 3; 
sizes.DirFeedthrough = 1; 
sizes.NumSampleTimes = 1; 
sys = simsizes(sizes); 
x0 = UgSk0; 
str = []; 
ts = [0 0]; 
% Конец процедуры mdlInitializeSizes 
%========================================== 
function z = mdlDerivatives(t,x,M,J,IJ) 
% ВХОДНОЙ вектор M является вектором проекций моментов внешних сил, 
% действующих на космический аппарат соответственно по осям X Y Z 
% х(1)=om(1); х(2)=om(2); х(3)=om(3); 
% z(1)=d(om(1))/dt; z(2)=d(om(2))/dt; z(3)=d(om(3))/dt; 
% |Jx Jxy Jxz| 
% J = |Jxy Jy Jyz| - матрица моментов инерции КА 
% |Jxz Jyz Jz | 
om=x(1:3); 
z=IJ*(M-cross(om,J*om)); % ДИНАМИЧЕСКИЕ уравнения Ейлера 
% Конец процедуры mdlDerivatives 
%================================= 
function y = mdlOutputs(x,M,J,IJ) 
y(1:3)=x; 
om=x(1:3); 
zom=IJ*(M-cross(om,J*om)) ; % Определения УСКОРЕНИЙ 
y(4:6)=zom; 
% Конец процедуры mdlOutputs