| 查看: 685 | 回复: 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); |
回复此楼» 猜你喜欢
自荐读博
已经有5人回复
-
求个博导看看
已经有16人回复
-
上海工程技术大学张培磊教授团队招收博士生
已经有4人回复
-
上海工程技术大学【激光智能制造】课题组招收硕士
已经有5人回复
-
求助院士们，这个如何合成呀
已经有4人回复
-
临港实验室与上科大联培博士招生1名
已经有9人回复
-
需要合成515-64-0，50g，能接单的留言
已经有4人回复
-
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
-
带资进组求博导收留
已经有10人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求四面体体积公式和点的法向量求解。
已经有12人回复
- 求助一个向量范数的不等式问题!
已经有3人回复
- 使用mathtype里面的括号，长公式，括号内没法换行，一换行就变成维数加1的向量
已经有11人回复
- 欧式空间中定义的内积与向量2范数之间的关系
已经有3人回复
- 高手请进，求向量叉乘公式。
已经有5人回复
信彼南山
木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 4142.9
- 散金: 1221
- 红花: 16
- 帖子: 1178
- 在线: 233.5小时
- 虫号: 1133529
- 注册: 2010-10-27
- 专业: 导航、制导与传感技术
2楼2015-03-16 11:22:31
3楼2015-03-16 13:33:54
|
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
6楼2015-03-16 16:07:39
|
编辑MATLAB程序： x y z b c 分别对应x(1——5)， Y暂时没求 ODE45_fun【m文件】： function dx=ODE45_fun(t,x) dx = zeros(5,1); % a column vector 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; dx(4)=x(5)*a+x(3)*x(1)+x(3)+e*f*g*(x(4)*x(5)+x(2)*x(3)+e*f+e*g) 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=[dx(1);dx(2);dx(3);dx(4);dx(5);dx(6)]; ODE45_main【m文件】： tspan=[0,1e-5]; x0=[0;0.1;2;0;0.1]; [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)) 以上的程序没问题 【问】 （1）由于Y=x(4)*x(5)+x(2)*x(3)+e*f+e*g Y=x(4)*x(5)+x(2)*x(3)+e*f+e*g,该怎么添加到上面的程序里？ （2）tspan=[0,1e-5]; 表示在MATLAB里算了多少次？或需要多长时间呢？ 谢谢  六元方程组怎么用MATLAB编辑程序？.jpg |
7楼2015-03-17 09:26:04