数值法求解微分方程.doc

数值法求解微分方程P168 第27题建立一个M的函数文件function xdot=exf(t,x)u=1-2*t;xdot=0 1;1 -t*x+0 1*u;再新建一个M文件t0=0;tf=pi;x0t=0.1;0.2;t,x=ode23(exf,t0,tf,x0t)y1=x(:,1)y2=x(:,2)plot(t,y1,y,t,y2,m:)legend(y,d(y)grid on运行结果：t = 0 0.0145 0.0873 0.2014 0.3259 0.4621 0.6121 0.7778 0.9621 1.1482 1.2767 1.4053 1.5188 1.6706 1.8601 2.0891 2.3569 2.6546 2.9687 3.1416x = 0.1000 0.2000 0.1030 0.2158 0.1214 0.2883 0.1598 0.3798 0.2116 0.4479 0.2756 0.4847 0.3485 0.4813 0.4245 0.4282 0.4939 0.3168 0.5388 0.1599 0.5513 0.0320 0.5466 -0.1069 0.5272 -0.2358 0.4780 -0.4127 0.3788 -0.6347 0.2037 -0.8952 -0.0742 -1.1797 -0.4677 -1.4646 -0.9691 -1.7290 -1.2793 -1.8586y1 = 0.1000 0.1030 0.1214 0.1598 0.2116 0.2756 0.3485 0.4245 0.4939 0.5388 0.5513 0.5466 0.5272 0.4780 0.3788 0.2037 -0.0742 -0.4677 -0.9691 -1.2793y2 = 0.2000 0.2158 0.2883 0.3798 0.4479 0.4847 0.4813 0.4282 0.3168 0.1599 0.0320 -0.1069 -0.2358 -0.4127 -0.6347 -0.8952 -1.1797 -1.4646 -1.7290 -1.8586t = 0 0.0145 0.0873 0.2014 0.3259 0.4621 0.6121 0.7778 0.9621 1.1482 1.2767 1.4053 1.5188 1.6706 1.8601 2.0891 2.3569 2.6546 2.9687 3.1416x = 0.1000 0.2000 0.1030 0.2158 0.1214 0.2883 0.1598 0.3798 0.2116 0.4479 0.2756 0.4847 0.3485 0.4813 0.4245 0.4282 0.4939 0.3168 0.5388 0.1599 0.5513 0.0320 0.5466 -0.1069 0.5272 -0.2358 0.4780 -0.4127 0.3788 -0.6347 0.2037 -0.8952 -0.0742 -1.1797 -0.4677 -1.4646 -0.9691 -1.7290 -1.2793 -1.8586y1 = 0.1000 0.1030 0.1214 0.1598 0.2116 0.2756 0.3485 0.4245 0.4939 0.5388 0.5513 0.5466 0.5272 0.4780 0.3788 0.2037 -0.0742 -0.4677 -0.9691 -1.2793y2 = 0.2000 0.2158 0.2883 0.3798 0.4479 0.4847 0.4813 0.4282 0.3168 0.