請問您是研究的自激振動問題嗎，或者說是非線性動力學、非線性振動的有關問題？
vdpol 這個方程和MATLAB程序我都看過。
vdpol 就是英文:Van de pol（範德波爾）的縮寫，而ydot是個輸出變量組。是個數值數組，意思是把非線性微分方程的近似數值解y(t)和y'(t)的值都賦予到y這個數組中。
源代碼如下：
function ydot=vdpfun(t,y,flag,u,x0,w0,v,w)
ydot=[y(2);u*(x0^2-y(1)^2)*y(2)-y(1)*w0^2-v*cos(w*t)];
end
您還可以利用ode23這個命令求解這個微分方程從而畫出速度y'關於控制參數u（類似於阻尼的東西）變化的分岔圖：