没有预定义间隔的分段函数

Question

我想制作一个 .m 文件，它代表一个分段函数并返回一个包含所有离散值计算的向量。

为了更清楚一点，我想要一个函数（我将其命名为 Iapp，并且与时间相关，因此 Iapp(t)）在前 100 秒内返回零，然后在 100-120 秒内返回 0.5，然后在 120 秒内返回零-220s 然后 0.5+0.2 持续 220-240s，然后继续。

我知道可以使用逻辑索引定义分段函数，但我的问题是我想要函数的时间间隔没有预定义。所以我不知道逻辑索引是如何工作的……如果时间间隔不是 120 的倍数，它就不起作用。

我尝试了以下方法：

function Vect_Iapp = Iapp_morceaux(tspan, h)
i = 1;
j = 1;
t = tspan(1):h:tspan(2);
while t(i) < tspan(2)         
    while(t(i)< (j*100 + (j-1)*20))
        Iapp(i) = 0;
        i = i + 1;
    end
    while (t(i)>j*100 && t(i) < j*100 + j*20)
        Iapp(i) = 0.5 + j*0.2;
        i = i + 1;
    end
    j = j + 1;
end    
Vect_Iapp = Iapp;
end

但算法并不总是像它应该的那样工作。关于如何定义此功能的任何想法？ 请注意，我还希望能够以某种方式为 tspan 提供一个标量值，并使函数只返回一个标量值。

Answer 1

恐怕我没有得到卷积循环的目的，但据我所知，您可以从函数参数中推断出第 k 个部分向量的长度 nk 及其值 xk 。那为什么不一开始就这样做，然后通过类似的方式创建每个部分向量

iappk = xk*ones(1,nk);

最后将它们连接在一起

Iapp = [iapp1 iapp2 iapp3 iapp4]

希望对你有帮助。

Answer 2

这个函数？

x=linspace(0,600,1000);
y=Iapp(x);
plot(x,y)

function y=Iapp(t)
    r=mod(t,120);
    c=floor(t/120);
    VAL1=0;
    VAL2=0.5 + 0.2*c;
    y=VAL1.*(r<=100) + VAL2.*(r>100);
end