在matlab中，用ode23求解初值问题dy^2/dx^2+2dy/dx+y=cosx,y(0)=0,y'(0)=3/2 的数值解。

先化成多元一阶微分方程组,设 y=y1Dy=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,y