尝试用matlab代码求解Bouc-Wen模型

好的，下面是如何用 MATLAB 代码求解 Bouc-Wen 模型的示例代码：

function dydt = BoucWen(t,y,k,n,A,omega)
% Bouc-Wen Model
dydt = zeros(2,1);
dydt(1) = y(2);
dydt(2) = -k * (abs(y(2)))^n * sign(y(2)) - omega^2 * y(1) + A * cos(omega*t);
end