05-23.txt

《MATLAB R2007基础教程》 清华大学出版社,包含电子课件,实例原文件,习题答案,都是基于matlab实现

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')