rigid.m

来自「数字信号处理的MATLAB实现」· M 代码 · 共 5 行

function dy = rigid(t,y)
dy = zeros(3,1);    % 行向量
dy(1) = y(2) * y(3);
dy(2) = -y(1) * y(3);
dy(3) = -0.51 * y(1) * y(2);