1
初值问题
范德波德莫是界说范德波尔方程的函数。
type vanderpoldemo
2
按“Enter”键。
如图1所示。
3
方程被写当作一个二阶一阶赋格。对于参数Mu的分歧值,对它们进行计较。为了加速积分速度,我们按照参数μ的值来选择合适的求解器。
对于μ=1,肆意一个MATLAB的ODE解算器都能有用地求解范德波尔方程。下面利用的ODE45求解器就是这样一个例子。方程在[0,20]域内求解。
tspan = [0, 20];
y0 = [2; 0];
Mu = 1;
ode = @(t,y) vanderpoldemo(t,y,Mu);
[t,y] = ode45(ode, tspan, y0);
% Plot of the solution
plot(t,y(:,1))
xlabel('t')
ylabel('solution y')
title('van der Pol Equation, \mu = 1')
4
按“Enter”键。
得图2所示。