Question

MATLAB微分方程问题 dsolve('D2y-t*Dy=1-t^2','y(0)=0.1,Dy(0)=0.2','t')

dsolve('D2y-t*Dy=1-t^2','y(0)=0.1,Dy(0)=0.2','t')

function xdot=exf(t,x)
u=(1-t^2);
xdot=[0 1;t 0]*x+[0 1]'*u

>> clf
>> t0=0;tf=2;x0t=[0.1 0.2];
>> [t,x]=ode23('exf',[t0,tf],x0t)

y=x(:,1),y2=（t.^2）/2 - (2^(1/2)*pi^(1/2)*erf((2^(1/2)*t*i)/2)*i)/10 + 1/10;就是这一步，一直提醒我错误使用erf，是因为 dsolve('D2y-t*Dy=1-t^2','y(0)=0.1,Dy(0)=0.2','t')
求解出来的有虚数吗，那这道题如何用数值积分方法求解呢

Answer 1

问题出在你定义exf(t,x)有错误。

要改成下列形式才对。

function xdot=exf(t,x)

xdot=[x(2);1-t^2+t*x(2)];

end

运行结果