Using feedback() with time delay

Dear Matlab community,

I would expect the same step response in the plot in the following code, but the step response of G_y_r is unstable, while T_CL is stable.

The only difference is that T_CL is built using the feedback() function, while G_y_r is built manually.

Does someone have a hint why G_y_r has two internal delays (“Internal delays (seconds): 2.5e-05 2.5e-05”), while T_CL has only one?

It seems that feedback() is doing something special here?

Also, G_y_r has more internal state variables than T_CL.

It would be great if someone could clarify this!

Thanks and best regards,

%
C = tf([1.5, 47.12], [1.0, 0.0]);
G_delay = tf([1], [1], 'OutputDelay', 1.0 / 40.0e3);
P = tf([6773, 5.131e7], [1.0, 1.123e4, 5.193e7]);
H = tf([5.988e19], [1.0, 348900, 4.591e10, 2.7e15, 5.988e19]);
% build closed-loop using feedback()
T_CL = minreal(feedback(P * C * G_delay, H));
% build closed-loop manually
G_y_r = minreal((P * C * G_delay) / (1.0 + P * C * G_delay * H));
%
T_CL, G_y_r
%
figure;
hold on;
step(T_CL, '-b');
step(G_y_r, '*r');
grid on;
legend('T_{CL}', 'G_{y, r}', ...
  'Location', 'best');

NOTE:-

Matlabsolutions.com provide latest MatLab Homework Help,MatLab Assignment Help for students, engineers and researchers in Multiple Branches like ECE, EEE, CSE, Mechanical, Civil with 100% output.Matlab Code for B.E, B.Tech,M.E,M.Tech, Ph.D. Scholars with 100% privacy guaranteed. Get MATLAB projects with source code for your learning and research.

Your code is just fine except for one tiny but a curicial point that is internally built-in solver step size related issue of the control toolbox. In fact, there is only one time delay. Here is the correct code that addresses solver step size in computing the systems’ step responses. Now, your simulation results of the both systems converge perfectly well.

clc; clear variables
%
C = tf([1.5, 47.12], [1.0, 0.0]);
G_delay = tf(1, 1, 'OutputDelay', 1.0 / 40.0e3);
P = tf([6773, 5.131e7], [1.0, 1.123e4, 5.193e7]);
H = tf(5.988e19, [1.0, 348900, 4.591e10, 2.7e15, 5.988e19]);
% build closed-loop using feedback()
T_CL = minreal(feedback(P * C * G_delay, H));
% build closed-loop manually
G_y_r = minreal((P * C * G_delay) / (1.0 + P * C * G_delay * H));
%
disp('System T_CL:  '); T_CL  %#ok
disp('System G_y_r: '); G_y_r %#ok
%
figure;

SEE COMPLETE ANSWER CLICK THE LINK

Using feedback() with time delay

Using feedback() with time delay