Is it possible to create a stairstep function?

Question

Tyler Bikaun 2017-2-13

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/324847-is-it-possible-to-create-a-stairstep-function

编辑： Walter Roberson 2017-2-14

Hi all,

I am trying to create a function that emulates what is seen below, I don't think it may be possible with stairstep but I am unsure. Does anyone know how this could be done? As I want an equation that I can have an input x value and then output a y value from the function.

% Stairstep function 
x = [0:15];
y = [19,19,19,19,19,25,25,25,25,25,25,19,19,19,19,19];
figure
stairs(x,y)
ylim([0,30]);

Much appreciated

3 个评论
显示 1更早的评论隐藏 1更早的评论

John Chilleri 2017-2-14

在 MATLAB Online 中打开

Not sure I understand, unless you're speaking of exactly the x and y you have above. In that case,

function y = foo(x)
    vec = [ones(1,5)+18 ones(1,6)+24 ones(1,5)+18];
    y = vec(x+1);
end

which accepts any x in [0,15].

Tyler Bikaun 2017-2-14

在 MATLAB Online 中打开

I am using ODE45 to evaluate some equations, x being one of the variables within ODE45 will be used to evaluate a constant that is multiplied onto one of the state space equations within ODE45.

For example something along these lines

function xdot = ode(t,x)
 xdot(1) = x(2)
 xdot(2) = y_function*x(1);

The y_function is of the shape of a stairstep as shown above where the actual range I'll be using for x is approximately 0:0.001:100 (I don't need the stairstep explicitly, just need a way to call data from a function that is of the same nature).

I've done this using Simulinks 1D lookup table, but I don't want to use Simulink for this, however, I cannot find a way to do it using ODE45 and MATLAB for some reason.

请先登录，再进行评论。

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

Answer 1

Walter Roberson 2017-2-14

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/324847-is-it-possible-to-create-a-stairstep-function#answer_254668

编辑：Walter Roberson 2017-2-14

Although it would be possible to program that, you should not use it in the way you propose. The ode*() routines require that the values have continuous derivatives, which you would not satisfy with a stair-steps function.

If your discontinuities are according to fixed time-points, then you need to break up your ode45 call at the timepoints, start at one end, ode45, use the result for the initial conditions of the next ode45 call and so on.

If your discontinuities are according to x rather than t, then you need to use event functions to tell it to terminate the integration at that point, and loop back and do the next ode45 call using the results from each as the boundary conditions for the next.

ode45 typically notices the discontinuity of a stairs function and considers it a singularity at which it is unable to meet integration tolerances so it gives up.