发布网友 发布时间:2024-05-03 21:20
共1个回答
热心网友 时间:2024-08-11 06:25
先化成多元一阶微分方程组,设
y=y1
Dy=y2
则有
Dy1=y2
Dy2=cos(x)-2y2-y1
建立m文件vdpol.m
function yp = vdpol(t,y);
yp=[y(2);cos(t)-2*y(2)-y(1)];
在命令窗口执行
[t,y]=ode23('vdpol',[0,20],[0,3/2]);
y1=y(:,1); % 原方程的解
y2=y(:,2);
plot(t,y1,t,y2,'--') % y1(t),y2(t) 曲线图
pause,
plot(y1,y2),grid,
% 相轨迹图,即y2(y1)曲线