用MATLAB写y”-0.01(y')^2+2y=sint ,y(0）=1,y'(0)=1,0≦t≦3.作y^(t)的图像

首先，将二阶微分方程转化为一组一阶微分方程：

令y1=y，y2=y'，则有：

y1'=y2

y2'=-0.01*y2^2*2*y1+sin(t)

将初始条件代入，即可用ode45函数求解：

function dydt = odeFcn(t,y)
dydt = [y(2);-0.01*y(2)^2*2*y(1)+sin(t)];
end

[t,y] = ode45(@odeFcn,[0 3],[1 1]);

plot(t,y(:,2))
xlabel('t')
ylabel('y''(t)')

最后得到y'(t)的图像如下：