用matlab解数值题,利用matlab解决数值解的问题

利用matlab解决数值解的问题

利用matlab解决数值解的问题 实验一 操作一：如下操作file-->new-->function-m-files，创建一个my.m文件，并保存在适当的路径中(如：bin) function y=my(t,x) a=1; b=0.3; y=[a*x(1)*x(2)-b*x(1),-a*x(1)*x(2)] ; 操作二：在matlab窗口的操作 ts=0:50; x0=[0.02,0.98]; [t,x]=ode45( my ,ts,x0) plot(t,x(:,1),t,x(:,2)),grid,pause 运行结果： t = 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 x = 0.0200 0.9800 0.0390 0.9525 0.0732 0.9019 0.1285 0.8169 0.2033 0.6927 0.2795 0.5438 0.3312 0.3995 0.3444 0.2839 0.3247 0.2027 0.2863 0.1493 0.2418 0.1145 0.1986 0.0917 0.1599 0.0767 0.1272 0.0665 0.1004 0.0593 0.0787 0.0543 0.0614 0.0507 0.0478 0.0480 0.0371 0.0460 0.0287 0.0445 0.0223 0.0434 0.0172 0.0426 0.0133 0.0419 0.0103 0.0415 0.0079 0.0411 0.0061 0.0408 0.0047 0.0406 0.0036 0.0404 0.0028 0.0403 0.0022 0.0402 0.0017 0.0401 0.0013 0.0400 0.0010 0.0400 0.0008 0.0400 0.0006 0.0399 0.0005 0.0399 0.0004 0.0399 0.0003 0.0399 0.0002 0.0399 0.0002 0.0399 0.0001 0.0399 0.0001 0.0399 0.0001 0.0399 0.0001 0.0399 0.0000 0.0398 0.0000 0.0398 0.0000 0.0398 0.0000 0.0398 0.0000 0.0398 0.0000 0.0398 0.0000 0.0398 plot(x(:,2),x(:,1)),grid 实验二：(一次函数) 操作一：如下操作file-->new-->function-m-files，创建一个my.m文件，并保存在适当的路径中(如：bin) function y=wyp(t,x) a=5; y=a*x(1)-2; 操作二：在matlab窗口的操作 ts=0:50; x0=0.02; [t,x]=ode45( wyp ,ts,x0) plot(t,x),grid,pause 运行结果： t = 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 x = 1.0e+042 * 0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.