Level-2 S-function: finding the correct input value and type

I used chatgpt to write a code for a s-function block for a system of 2 ODEs. It is suppose to have an input and two outputs. The two ODEs are as follows:

m, g, R, L are constants, T is the input, the outputs are supposed to be v and theta. The code, saved as myODESFunction.m, is as follows:

function sys = myODESFunction(t, x, u, flag)

m = 384986.21; % Mass
g = 9.81; % Gravity
R = 3.83; % drag constant
L = 59.367; % lift constant

switch flag
    case 0 % Initialization
        sizes = simsizes;
        sizes.NumContStates = 2; % Number of continuous states
        sizes.NumDiscStates = 0; % Number of discrete states
        sizes.NumOutputs = 2; % Number of outputs
        sizes.NumInputs = 1; % Number of inputs
        sizes.DirFeedthrough = 0; % No direct feedthrough
        sizes.NumSampleTimes = 1; % One continuous sample time
        
        sys = simsizes(sizes);

        x0 = [252.222; 0.1745]; % Initial state vector
        str = [];
        ts = [0 0]; % Continuous sample time
        sys = x0;
        sys(2) = str;
        sys(3) = ts;

    case 3 % Update
        sys = [];

    case 1 % Derivatives
        x1 = x(1); % v
        x2 = x(2); % theta
        u1 = u(1); % T

        % Define your differential equations here
        dx1_dt = (-m*g*sin(x2) - R*x1^2 + u1) / m;
        dx2_dt = (-m*g*cos(x2) + L*x1^2) / (m*x1);

        sys = [dx1_dt; dx2_dt];

    case 2 % Outputs
        sys = [x(2); x(1)]; % Output theta and velocity
end

But when I tried to use it in a Level-2 s-function block, by entering [252.222 0.1745] as the two intial conditions in my code into the "s-function parameters" field, it tells me I don't have enough input arguments. But what more could I add?

Answer 1

[Aiswarya]

I understand that you are trying to write a Level 2 S-function for a system of 2 ODEs. The error you are getting "Not enough input arguments" is because your function is a Level 1 S-function which you are trying to use with the Level-2 S-function block in Simulink. There are some issues with the chatgpt code for the Level 1 S-function as well.

If you want to write Level 2 S-function, you may refer to the following documentation: https://www.mathworks.com/help/simulink/sfg/writing-level-2-matlab-s-functions.html. We can also convert Level 1 S-function to Level 2 S-function by following the steps explained in detail in the following link : https://www.mathworks.com/help/simulink/sfg/maintaining-level-1-matlab-s-functions.html.

Based on your code, I have provided a Level 2 S-function code below:

function myODESFunction_Level2(block)
% Level-2 MATLAB S-Function for myODESFunction

setup(block);

%endfunction
function setup(block)

%% Register number of dialog parameters  (2 parameters, v and theta)
block.NumDialogPrms = 2;
%% Register number of input and output ports
block.NumInputPorts  = 1;
block.NumOutputPorts = 2;
%% Setup functional port properties to dynamically
%% inherited.
block.SetPreCompInpPortInfoToDynamic;
block.SetPreCompOutPortInfoToDynamic;

block.InputPort(1).Dimensions        = 1;
block.InputPort(1).DirectFeedthrough = false;

block.OutputPort(1).Dimensions       = 1;
block.OutputPort(2).Dimensions       = 1;
%% Set block sample time to continuous
block.SampleTimes = [0 0];

%% Setup Dwork
block.NumContStates = 2;
%% Set the block simStateCompliance to default (i.e., same as a built-in block)
block.SimStateCompliance = 'DefaultSimState';
%% Register methods
block.RegBlockMethod('InitializeConditions',    @InitConditions);  
block.RegBlockMethod('Outputs',                 @Output);  
block.RegBlockMethod('Derivatives',             @Derivative);  

%endfunction
function InitConditions(block)
%% Initialize Dwork
block.ContStates.Data(1) = block.DialogPrm(1).Data;
block.ContStates.Data(2) = block.DialogPrm(2).Data;
%endfunction
function Output(block)
block.OutputPort(1).Data = block.ContStates.Data(1);
block.OutputPort(2).Data = block.ContStates.Data(2);

%endfunction
function Derivative(block)
m = 384986.21; % Mass
g = 9.81; % Gravity
R = 3.83; % drag constant
L = 59.367; % lift constant
x1 = block.DialogPrm(1).Data;
x2 = block.DialogPrm(2).Data;
u1 =  block.InputPort(1).Data;
block.Derivatives.Data(1) = (-m*g*sin(x2) - R*x1^2 + u1) / m;
block.Derivatives.Data(2) = (-m*g*cos(x2) + L*x1^2) / (m*x1); 
%endfunction

Now you can use the Level 2 S-function block in Simulink and set the block parameters as shown in the figure below:

This is similar to a Level 2 S-function example provided in MATLAB (msfcn_limintm.m). This can be accessed by the command

edit('msfcn_limintm.m')

Hope this helps!