我有以下代码要执行:
= [0;0;0;0];= 0;
for k = 1 : n
if (sw(k) == 0 && X(2,k)> 0.7 ) || (sw(k) == 1 && X(1,k) >0)
sw(k+1) = 1;
X(:,k+1) = X(:,k) +h*(A1*X(:,k) + B1*u);
else
sw(k+1) =0;
X(:,k+1) = X(:,k) +h*(A*X(:,k) + B*u);
end
end
这可以。我想使用内置的 ode 套件解决相同的 ODE。我如何检查和在每次迭代后查看要求解两个方程中的哪一个。