matlab求解微分方程结果不对

yzh07137

问题是这样的（见附图）
- y& D( W0 u0 M/ d" z6 W, p2 D
! h+ p" z. t8 d- p) I6 `' _然后我的程序如下

1. function Untitled
   3 n. f! Q& e; N: \
2.         clear all;clc;; i7 M, x% ^\" `' q# _& x
   
3.         f=@(t)(2*sin(t)*(t<(4*pi)) + 0);
   , d3 d! I* _& F% t: F  Y. T; ?! Q
4.         g=@(t)(0+cos(t).*(t>=((7*pi)/2)));
   , K( q: U- d# }
5.             function dy = rigid(t,y)
   7 F! p: a+ X  \2 u
6.                         dy = zeros(2,1);4 D% Y; h, k: w5 X
   
7.                         dy(1) = y(2)-f(t);
   9 O8 K% ~0 F3 ~# v, j\" y
8.                         dy(2) = y(1)*g(t)-y(2);
   & n4 Y4 G& f# Z
9.             end5 q. P0 J* ]* v9 Z
   
10.             options = odeset('RelTol',1e-4,'AbsTol',[1e-5 1e-5]);( B! S, ^- X9 [: \
    
11.             %[T, y] = ode45(@rigid, [0 20], [1 2],options)6 {( l0 _8 ?) z
    
12.             sol = ode45(@rigid, [0 20], [1 2],options);: J, t\" p! O  D7 P$ B
    
13.             x=linspace(0,20,200000);
    ( F+ \) p9 u5 m+ g
14.             y=deval(sol,x);
    7 L9 Z# O2 X- c\" X\" U5 m. F: b
15.             res=y(1,:)+y(2,:);/ h' S, C3 ]/ Z. I; |
    
16.             idx=find(abs(res-0)<1e-4) %相加，当和小于误差运行范围的时候可以认为它就为0; r6 u/ @# @2 V( H/ `1 _% z7 X  @
    
17.             xx=x(idx) %算出在x中的下标) u. D# w% H5 _( w
    
18.             F=@(t)(f(t)+g(t));
    * Y2 o: B3 n1 p5 u& T
19.             r=[]; %得出的解反代入方程求值，得到的值保存在r中。如果解正确，r中的值应该非常接近05 G0 B4 s0 p! \0 ?4 r: s
    
20.             for i=1:size(xx,2); Y# ?) \, r* u
    
21.                     r = [r F(xx(i))]; %将解的值依次带入。当前问题在于r中的值都和0差的很远
    * F2 `& k; N/ j+ o  g+ N. V, P
22.             end; ]- O$ E, p! B9 @
    
23.             r
    ; J9 ]8 ?\" g: z+ d4 x1 h
24. end

复制代码

问题就是最后算出来的解再带入方程进行检验得出的结果和0差的好远，都到了1.73几。
; Y" V! F1 h- z9 Z* ~( d) t我也不知道哪里用错了，但猜测可能是由于微分方程里出现了f(t)，g(t)的缘故
8 h' y# t4 M7 A5 D3 b1 F: j
求教大家我的程序是哪里出现了问题，该如何改呢？
谢谢
" N% G9 z5 T: p- F

yzh07137

算了，就知道这种论坛一般也没人