如何在 Matlab 中使用变化的输入/时间相关的输入函数求解 DAE？

我正在使用 ode15i 求解器求解 DAE 问题。我有 8 个变量和 8 个方程，系统很复杂，到目前为止唯一可用的求解器是 ode15i。

我使用过该指南：http://se.mathworks.com/help/symbolic/set-up-your-dae-problem.html http://se.mathworks.com/help/symbolic/set-up-your-dae-problem.html。但是，本指南无法帮助我解决不同输入的问题。

我的系统有一个与时间相关的输入，一个函数。函数本身很简单，但问题是DAE系统中的时间t是符号形式的，而我的输入函数无法求解相应的输入值，因为它无法计算符号时间值t。

这是我的 DAE 系统的输入函数：

    function [delta] = delta(t)
    t=vpa(t)
       if t <= 1.01
       delta = 0;
       elseif t <= 3.61
       delta = 0.1222;
       elseif t <= 4.33
       delta = 0;
       elseif t <= 7.21
       delta = -0.1222;
       else
       delta = 0;
    end

戴伊问题：

syms v1(t) r1(t) r2(t) r3(t) r4(t) fii(t) sig(t) psi(t);
tspan = 0:0.01:10;
....
a11 = delta(t)-(v1(t)+r1(t)*l_1_1)/u;
....
eqs = [...]
vars [...]
f = daeFunction(eqs,vars);
y0est = zeros(8,1);
yp0est = zeros(8,1);
opt = odeset('RelTol',10.0^(-7),'Abstol',10.0^(-7));
[y0,yp0] = decic(f,0,y0est, [], yp0est, [], opt);
[t,y] = ode15i(f,tspan, y0, yp0, opt);

我收到的错误如下：

Conversion to logical from sym is not possible.
 Error in delta (line 3)
    if t <= 1.01
 Error in Nonlin_painonsiirto_DUO2 (line 181)
 a11 = delta(t)-(v1(t)+r1(t)*l_1_1)/u;

如果 delta 输入是常数或 (t) 的三角函数，则系统可以工作，例如：

delta = 0.0175;
delta = sin(t)*0.0175;

我也尝试过不使用 vpa(t) 命令和使用命令 double(t) 的 delta 函数，但它什么也没做。我还尝试使用系统的时间向量 tspan=0:0.01:10 作为输入，它作为 ode15i 的时间跨度给出：

delta(tspan)

但是，然后它尝试计算整个向量 tspan，这会导致错误，因为矩阵维度不一致。

希望这里的问题是可以理解的，谢谢。

-Jere