vtb9_3ex

来自「vtoolbox3.5是振动数值模拟工具箱」· 代码 · 共 11 行

function [zd]=rfunc(z,u,t)
%    function for    2
%                  dx      dx
%                m --  + c -- +k x = f(t)
%                    2     dt
%                  dt                     dx
%    where m=1,k=1,c=.1, and z(1)=x, z(2)=--
%                                         dt
zd=[z(2);-.1*z(2)-z(1)+u(2)];
end;