Step response with initial condition

Question

Drishya Dinesh 2021-4-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/806347-step-response-with-initial-condition

评论： Paul 2021-4-20

How to obtain a step response starting from 10?

ie; the initial value of u_del is 10. Hence the response should start from 10.

k=0.2;
t2=400;
u_del=(0.0022*k*(s+0.06931)*(s^2+0.4852*s+0.1492))/((s+0.04833)*(s+0.004352)*(s^2+0.06012*s+0.01331));
figure
step(u_del,t2);
ylabel('Velocity,u (m/s)','fontsize',10);
title('Time Response');
grid 

2 个评论
显示无隐藏无

Paul 2021-4-19

编辑：Paul 2021-4-20

Do you mean the step response output starts at 10, i.e., y(t=0) = 10?

Or do you mean the step is applied at t = 0, but the output is delayed and doesn't start to respond until t = 10?

It sounds like you mean the former, but it's not really clear.

Drishya Dinesh 2021-4-20

The step response output should start at 10, i.e., y(t=0) = 10

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Star Strider 2021-4-19

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/806347-step-response-with-initial-condition#answer_678882

在 MATLAB Online 中打开

To set the step amplitude to 10, try this:

k=0.2;
t2=400;
s = tf('s');
u_del=(0.0022*k*(s+0.06931)*(s^2+0.4852*s+0.1492))/((s+0.04833)*(s+0.004352)*(s^2+0.06012*s+0.01331));
opts = stepDataOptions('StepAmplitude',10);                         % Set Amplitude
figure
step(u_del,t2,opts);                                                % Add ‘opts’ Argument
ylabel('Velocity,u (m/s)','fontsize',10);
title('Time Response');
grid 

.

4 个评论
显示 2更早的评论隐藏 2更早的评论

Star Strider 2021-4-19

在 MATLAB Online 中打开

My pleasure!

From the step documentation for time vector t:

step always applies the step input at t = 0, regardless of Ti.

So, it will be necessary to use lsim to do what you want —

k=0.2;
t2=400;
s = tf('s');
u_del=(0.0022*k*(s+0.06931)*(s^2+0.4852*s+0.1492))/((s+0.04833)*(s+0.004352)*(s^2+0.06012*s+0.01331));
Ts = 1;                                                             % Time Step (sec)
tv = linspace(0, 400, fix(400/Ts))*Ts;                              % Time Vector
u = tv>10;                                                          % Unit Step Input Starting AT t=10 sec
% [y,t] = lsim(u_del, u, tv);                                         % Outputs (If Desired)
figure
lsim(u_del, u, tv)                                                  % Simulate System & Plot
ylabel('Velocity,u (m/s)','fontsize',10);
title('Time Response');
grid k=0.2;
t2=400;
s = tf('s');
u_del=(0.0022*k*(s+0.06931)*(s^2+0.4852*s+0.1492))/((s+0.04833)*(s+0.004352)*(s^2+0.06012*s+0.01331));
Ts = 1;                                                             % Time Step (sec)
tv = linspace(0, 400, fix(400/Ts))*Ts;                              % Time Vector
u = tv>10;                                                          % Unit Step Input Starting AT t=10 sec
% [y,t] = lsim(u_del, u, tv);                                         % Outputs (If Desired)
figure
lsim(u_del, u, tv)                                                  % Simulate System & Plot
ylabel('Velocity,u (m/s)','fontsize',10);
title('Time Response');
grid 

That works correctly, at least as I understand what you want to do. If you want a different amplitude than 1 for the input, multiply ‘u’ by that value. If you want a different version of the step input, please be specific.

Drishya Dinesh 2021-4-19

在 MATLAB Online 中打开

Sir I'm getting an error

Index exceeds the number of array elements (20).
Error in Untitled18 (line 3)
s = tf('s');

Star Strider 2021-4-19

在 MATLAB Online 中打开

When I first tried to run your code, it threw:

Unrecognized function or variable 's'.

so I added the tf call, since that makes sense in the context of the Control System Toolbox, and the posted code.

I cannot guess as to what the exact problem is, only that ‘tf’ appears to have been assigned as a vector or other variable in code you did not post, and MATLAB is interpreting the 's' as the ASCII numeric equivalent, and throwing that error. The unposted code therefore overshadowed the very useful tf function with something else.

Ignore whatever exists before my posted code and run only my code as I posted it to get the desired result.

My posted code runs without error.

请先登录，再进行评论。

Answer 2

Paul 2021-4-20

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/806347-step-response-with-initial-condition#answer_679407

According to the question, the output should satisfy y(t=0) = 10. But the system, u_del, as specified will yield a step response that starts at y(t=0) = 0, in the absence of any initial conditions on the states of the system. So the next question is: what causes the step response to start at y(t=0) = 10? If it is initial conditions on the states of the system, then we have to use a state space representation for u_del, and even then there will be different initial conditions on states that will result in an initial output y(t=0) = 10 but different dynamics for y(t).

2 个评论
显示无隐藏无

Drishya Dinesh 2021-4-20

How to give initial conditions if the system u_del is represented in state space to obtain output with y(0)=10

Paul 2021-4-20

在 MATLAB Online 中打开

One of many, many possibilities would be:

>> u_delss=ss(u_del);
>> x0=[10/u_delss.c(1) 0 0 0];
>> t = 0:1:t2;
>> y = step(u_delss,t) + lsim(u_delss,0*t,t,x0);
>> plot(t,y,t(1:10:end),y(1:10:end),'o'),grid

Probably doesn't give the response you're looking for, but it is the response of the system u_delss to a unit step input with initial conditions x0 that result in y(t=0) = 10.

Did you just want something simple like:

[y,t] = step(u_del);
plot(t,y+10),grid

Note that this isn't the step response of u_del.

请先登录，再进行评论。

Step response with initial condition

2 个评论
显示无隐藏无

回答（2 个）

4 个评论
显示 2更早的评论隐藏 2更早的评论

2 个评论
显示无隐藏无

另请参阅

类别

标签

Community Treasure Hunt

Step response with initial condition

2 个评论 显示 无隐藏 无

回答（2 个）

4 个评论 显示 2更早的评论隐藏 2更早的评论

2 个评论 显示 无隐藏 无

另请参阅

类别

标签

Community Treasure Hunt

2 个评论
显示无隐藏无

4 个评论
显示 2更早的评论隐藏 2更早的评论

2 个评论
显示无隐藏无