mrdivide tf, ss

modules/control_system/functions/@ss/mrdivide.m

@@ -0,0 +1,15 @@
+%=============================================================================
+% Copyright (c) 2023-present Allan CORNET (Nelson)
+%=============================================================================
+% This file is part of the Nelson.
+%=============================================================================
+% LICENCE_BLOCK_BEGIN
+% SPDX-License-Identifier: LGPL-3.0-or-later
+% LICENCE_BLOCK_END
+%=============================================================================
+function varargout = mrdivide(varargin)
+  sys1 = ss(varargin{1});
+  sys2 = ss(varargin{2});
+  varargout{1} = sys1 * inv(sys2);
+end
+%=============================================================================

modules/control_system/functions/@tf/mrdivide.m

@@ -0,0 +1,72 @@
+%=============================================================================
+% Copyright (c) 2023-present Allan CORNET (Nelson)
+%=============================================================================
+% This file is part of the Nelson.
+%=============================================================================
+% LICENCE_BLOCK_BEGIN
+% SPDX-License-Identifier: LGPL-3.0-or-later
+% LICENCE_BLOCK_END
+%=============================================================================
+function varargout = mrdivide(varargin)
+  narginchk(2, 2);
+  nargoutchk(0, 1);
+  sysA = tf(varargin{1});
+  sysB = tf(varargin{2});
+  if ~issiso(sysA) && ~issiso(sysB)
+    error(_('SISO LTI model expected.'));
+  end
+  Ts = mrdivide_timesample(sysA.Ts, sysB.Ts);
+
+  numeratorA = sysA.Numerator{1}{1};
+  denominatorA = sysA.Denominator{1}{1};
+  numeratorB = sysB.Numerator{1}{1};
+  denominatorB = sysB.Denominator{1}{1};
+
+  if isempty(numeratorA) || isempty(numeratorB) || isempty(denominatorA) || isempty(denominatorB)
+    error(_('Matrix dimensions must agree.'));
+  end
+
+  if isstatic(sysA) && isstatic(sysA)
+    numerator = numeratorA / numeratorB;
+    sys = tf(numerator, 1, Ts);
+  else
+    numerator = conv(numeratorA, denominatorB);
+    denominator = conv(denominatorA, numeratorB);
+    sys = tf(numerator, denominator, Ts);
+  end 
+
+  if isa(sysA, 'tf') && strcmp(sysA.Variable,'z^-1')
+    sys.Variable = 'z^-1';
+  elseif isa(sysB, 'tf') && strcmp(sysB.Variable,'z^-1')
+    sys.Variable = 'z^-1';
+  end
+
+  UserData = [];
+  if ~isempty(sysA.UserData) && ~isempty(sysB.UserData)
+    UserData = sysB.UserData;
+  else
+    if ~isempty(sysA.UserData)
+      UserData = sysA.UserData;
+    end
+    if ~isempty(sysA.UserData)
+      UserData = sysA.UserData;
+    end
+  end
+  if ~isempty(UserData)
+    sys.UserData = UserData;
+  end
+  varargout{1} = sys;
+end
+%=============================================================================
+function Ts = mrdivide_timesample(TsA, TsB)
+  Ts = TsA;
+  if Ts ~= TsB
+    if ((Ts > 0) && (TsB > 0))
+      error(_('Sampling times must agree.'));
+    end
+    if (Ts < 0)
+      Ts = TsB;
+    end
+  end
+end    
+%=============================================================================

modules/control_system/tests/test_ss_mrdivide.m

@@ -0,0 +1,16 @@
+%=============================================================================
+% Copyright (c) 2023-present Allan CORNET (Nelson)
+%=============================================================================
+% This file is part of the Nelson.
+%=============================================================================
+% LICENCE_BLOCK_BEGIN
+% SPDX-License-Identifier: LGPL-3.0-or-later
+% LICENCE_BLOCK_END
+%=============================================================================
+A = [1 2; 3 4];
+B = [1 0; 0 1];
+C = [1 1; 1 1];
+D = [1 2; 3 4];
+sys1 = ss(A, B, C, D);
+sys2 = ss(A, B, C, D);
+R = sy1 / sys2;

modules/control_system/tests/test_tf_mrdivide.m

@@ -0,0 +1,35 @@
+%=============================================================================
+% Copyright (c) 2023-present Allan CORNET (Nelson)
+%=============================================================================
+% This file is part of the Nelson.
+%=============================================================================
+% LICENCE_BLOCK_BEGIN
+% SPDX-License-Identifier: LGPL-3.0-or-later
+% LICENCE_BLOCK_END
+%=============================================================================
+sysA = tf(3);
+sysB = tf(4);
+R = sysA / sysB;
+REF = tf(0.75);
+assert_isequal(R, REF);
+%=============================================================================
+sysA = tf(3);
+R = sysA / 3;
+REF = tf(1);
+assert_isequal(R, REF);
+%=============================================================================
+sysA = tf(3);
+R = 3 / sysA;
+REF = tf(1);
+assert_isequal(R, REF);
+%=============================================================================
+num = [3 4];
+den = [3 1 5];
+Ts = 0.2;
+sys = tf(num, den, Ts);
+R = sys / sys;
+num_REF = [9  15  19  20];
+den_REF = [9  15  19  20];
+REF = tf(num_REF, den_REF, Ts);
+assert_isequal(R, REF);
+%=============================================================================