05-23.txt

来自「MATLABR2007教程以及i习题和答案」· 文本 代码 · 共 12 行

function dydt = vdp1(t,y)
dydt = [y(2); (1-y(1)^2)*y(2)-y(1)];


>> [t,y] = ode45(@vdp1,[0 20],[2; 0]);

>> plot(t,y(:,1),'-',t,y(:,2),'--')
>> title('Solution of van der Pol Equation, \mu = 1');
>> xlabel('time t');
>> ylabel('solution y');
>> legend('y_1','y_2')