| 查看: 869 | 回复: 6 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
空向量是因为我的自己的式子巧合出现的吗？ 已有1人参与
|
|||
|
程序老出错，据说出现空向量 怎么避免出现空向量？、 谢谢 【m.file】ODE45_fun： function dx=ODE45_fun(t,x) a1=1;a2=1;e1=9;e2=7;f1=56;f2=98;g1=76;g2=665; a=a1+a2; e=e1+e2; f=f1+f2; g=g1+g2; dx(1)=x(3)+x(2)-x(4)*x(2)+2*x(5)*x(2)+x(2)*x(2)+x(5)+e+f+g; dx(2)=x(1)+x(3)-x(4)*x(2)+x(2)*x(2)+x(5)*x(5)+x(5)*x(2)+e*f+g; dx(3)=x(2)+x(1)-x(5)*x(1)+x(4)*x(3)-x(3)-x(4)-e*g-f; x(6)=x(4)*x(5)+x(2)*x(3)+e*f+e*g; % 就是x(6) 而非dx(6) dx(4)=x(5)*a-x(3)*x(1)-x(3)-e*f*g*x(6) dx(5)=x(4)*a-x(2)*x(1)+e*f*x(1)+x(2)*x(2)*x(5)+x(2)*x(5)*x(5)+e*f+f*g dx=[dx(1);dx(2);dx(3);dx(4);dx(5);dx(6)]; 【m.file】DE45_main： tspan=[0,10]; x0=[0;0.1;2;0;0.1;2]; [t,x]=ode45('ODE45_fun',tspan,x0); data=[t,x]; save ODE45_data.txt data -ascii subplot(2,3,1),plot(t,x(1)) subplot(2,3,2),plot(t,x(2)) subplot(2,3,3),plot(t,x(3)) subplot(2,3,4),plot(t,x(4)) subplot(2,3,5),plot(t,x(5)) subplot(2,3,6),plot(t,x(6)) >> ODE45_main dx = 1.0e+010 * 0.0000 0.0000 -0.0000 -2.6146 dx = 1.0e+010 * 0.0000 0.0000 -0.0000 -2.6146 0.0000 ??? Attempted to access dx(6); index out of bounds because numel(dx)=5. Error in ==> ODE45_fun at 14 dx=[dx(1);dx(2);dx(3);dx(4);dx(5);dx(6)]; Error in ==> funfun\private\odearguments at 110 f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0. Error in ==> ode45 at 173 [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, ... Error in ==> ODE45_main at 3 [t,x]=ode45('ODE45_fun',tspan,x0); |
回复此楼» 猜你喜欢
40+中年人失业后，还能找到好工作吗？
已经有13人回复
-
用NBS溴化时溶剂除了CCl4外还有什么
已经有6人回复
-
申请科研项目纯靠实力能成功吗？
已经有11人回复
-
极限拉扯，如何抉择
已经有15人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求四面体体积公式和点的法向量求解。
已经有12人回复
- 求助一个向量范数的不等式问题!
已经有3人回复
- 使用mathtype里面的括号，长公式，括号内没法换行，一换行就变成维数加1的向量
已经有11人回复
- 欧式空间中定义的内积与向量2范数之间的关系
已经有3人回复
- 高手请进，求向量叉乘公式。
已经有5人回复
3楼2015-03-16 13:33:54
信彼南山
木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 4142.9
- 散金: 1221
- 红花: 16
- 帖子: 1178
- 在线: 233.5小时
- 虫号: 1133529
- 注册: 2010-10-27
- 专业: 导航、制导与传感技术
2楼2015-03-16 11:22:31
|
http://muchong.com/bbs/viewthread.php?tid=8647830&pid=10#pid10 已经看到了，谢谢耐心讲解 MATLAB程序中的式子已经按照图片修改过了 前期是想程序的可行性，还是想的简单了 运算后还是那些错误结果 dx=[dx(1);dx(2);dx(3);dx(4);dx(5);dx(6)] 确实是习惯写成dx(6)了，当时没注意 但确实有六个变量啊 写成：dx=[dx(1);dx(2);dx(3);dx(4);dx(5)] 那么x(6)怎么办呢？ 彻底不会写了  |
4楼2015-03-16 14:00:18
信彼南山
木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 4142.9
- 散金: 1221
- 红花: 16
- 帖子: 1178
- 在线: 233.5小时
- 虫号: 1133529
- 注册: 2010-10-27
- 专业: 导航、制导与传感技术
5楼2015-03-16 14:45:14
200