查看: 409|回复: 3|关注: 0

[已解决] 使用ode45出现odearguments 错误

[复制链接]

新手

16 麦片

财富积分


050


5

主题

15

帖子

0

最佳答案
  1. clear all
  2. close all
  3. tspan=[0:0.01:10];
  4. y0=[0;0;0;0;0;0];
  5. [t,y]=ode45('kuta1',tspan,y0);
  6. subplot(311);
  7. plot(t,y(:,1));
  8. xlabel('t');
  9. ylabel('x1');
  10. subplot(312);
  11. plot(t,y(:,3));
  12. xlabel('t');
  13. ylabel('x2');
  14. subplot(313);
  15. plot(t,y(:,5));
  16. xlabel('t');
  17. ylabel('x3');
复制代码

新手

16 麦片

财富积分


050


5

主题

15

帖子

0

最佳答案
 楼主| 发表于 2019-3-14 20:59:18 | 显示全部楼层

  1. function f=kuta1(t,y)
  2. f=zeros(6:1);
  3. k1=4e5;
  4. k2=1e4;
  5. k3=5e4;
  6. k4=8e5;
  7. c1=0;
  8. c2=0;
  9. c3=0;
  10. c4=1600;
  11. m1=1.5;
  12. m2=0.075;
  13. m3=45;
  14. kbc=8*10^7;
  15. ktc=8*10^7;
  16. a=0.00004;
  17. ht=0.00004;
  18. hb=0.00006;
  19. a1=(k1+k2+(kbc+ktc)*a)/m1;
  20. a2=k2/m1;
  21. a3=(c1+c2)/m1;
  22. a4=c2/m1;
  23. a5=k2/m2;
  24. a6=(k2+k3)/m2;
  25. a7=c2/m2;
  26. a8=(c2+c3)/m2;
  27. b1=kbc*a;
  28. b2=ktc*a;
  29. a9=b1/m1;
  30. a10=c4/m3;
  31. a11=(b2*ht-b1*hb)/m1;
  32. a12=(b1+k4)/m3;
  33. b3=(b1*hb)/m3;
  34. b5=b1/m3;

  35. f(1)=y(2);
  36. f(2)=a4*y(4)+a2*y(3)-a3*y(2)-a1*y(1)+a9*y(5)+a11;
  37. f(3)=y(4);
  38. f(4)=a7*y(2)+a5*y(1)-a8*y(4)-a6*y(3);
  39. f(5)=y(6);
  40. f(6)=b5*y(1)-a10*y(6)-a12*y(5)+b3;
  41. end
复制代码

新手

16 麦片

财富积分


050


5

主题

15

帖子

0

最佳答案
 楼主| 发表于 2019-3-14 21:00:50 | 显示全部楼层
报错
错误使用 odearguments (line 93)
KUTA1 必须返回列向量。

出错 ode45 (line 115)
  odearguments(FcnHandlesUsed, solver_name, ode, tspan, y0, options, varargin);

出错 oode45 (line 6)
[t,y]=ode45('kuta1',tspan,y0);

论坛优秀回答者

入门

154 麦片

财富积分


50500


0

主题

373

帖子

33

最佳答案
  • 关注者: 8
发表于 2019-3-18 22:44:17 | 显示全部楼层 |此回复为最佳答案
把futa1函数     f=zeros(6:1)换成f=zeros(6,1)就可以了
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /4 下一条

快速回复 返回顶部 返回列表