Question

用matlab解：y=dsolve('D2y+exp(-t)*Dy-sin(sqrt(t))*y=exp(-t)*sin

Answer 1

此微分方程没有显式解,建议用数值解法

function dyy=xielei(t,y) % 保存为 xielei.m 文件

dyy=[y(2);-exp(-t)*y(2)+sin(sqrt(t))*y(1)+exp(-t)*sin(sqrt(t))];

--------------------------------------------

运行以下：

clear;clc

y0=[1;-2];

[t,y]=ode45('xielei',[0 5],y0);

plot(t,y(:,1),t,y(:,2))