matlab求解微分方程结果不对
问题是这样的（见附图）然后我的程序如下function Untitled
clear all;clc;
f=@(t)(2*sin(t)*(t<(4*pi)) + 0);
g=@(t)(0+cos(t).*(t>=((7*pi)/2)));
function dy = rigid(t,y)
dy = zeros(2,1);
dy(1) = y(2)-f(t);
dy(2) = y(1)*g(t)-y(2);
end
options = odeset('RelTol',1e-4,'AbsTol',);
% = ode45(@rigid, , ,options)
sol = ode45(@rigid, , ,options);
x=linspace(0,20,200000);
y=deval(sol,x);
res=y(1,:)+y(2,:);
idx=find(abs(res-0)<1e-4) %相加，当和小于误差运行范围的时候可以认为它就为0
xx=x(idx) %算出在x中的下标
F=@(t)(f(t)+g(t));
r=[]; %得出的解反代入方程求值，得到的值保存在r中。如果解正确，r中的值应该非常接近0
for i=1:size(xx,2)
r = ; %将解的值依次带入。当前问题在于r中的值都和0差的很远
end
r
end问题就是最后算出来的解再带入方程进行检验得出的结果和0差的好远，都到了1.73几。
我也不知道哪里用错了，但猜测可能是由于微分方程里出现了f(t)，g(t)的缘故
求教大家我的程序是哪里出现了问题，该如何改呢？
谢谢
算了，就知道这种论坛一般也没人
页:
[1]