# Step response with initial condition

77 ビュー (過去 30 日間)
Drishya Dinesh 2021 年 4 月 19 日 8:46
コメント済み: Paul 2021 年 4 月 20 日 4:08
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 件のコメント表示非表示 1 件の古いコメント
Drishya Dinesh 2021 年 4 月 20 日 2:06
The step response output should start at 10, i.e., y(t=0) = 10

サインインしてコメントする。

### 回答 (3 件)

VBBV 2021 年 4 月 19 日 9:14
k = 1:10;
t2 = linspace(1,400,length(k))
u_del=@(k,s) (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(1:10,1:10),t2);
ylabel('Velocity,u (m/s)','fontsize',10);
title('Time Response');
grid
try this
##### 1 件のコメント表示非表示 なし
Drishya Dinesh 2021 年 4 月 19 日 9:20
But it does'nt start from 10 at t=0.

サインインしてコメントする。

Star Strider 2021 年 4 月 19 日 12:30
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
ylabel('Velocity,u (m/s)','fontsize',10);
title('Time Response');
grid
.
##### 4 件のコメント表示非表示 3 件の古いコメント
Star Strider 2021 年 4 月 19 日 15:31
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.

サインインしてコメントする。

Paul 2021 年 4 月 20 日 2:43
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 件のコメント表示非表示 1 件の古いコメント
Paul 2021 年 4 月 20 日 4:08
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.

サインインしてコメントする。

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by