| 查看: 1468 | 回复: 9 | |||
[交流]
MATLAB求助 已有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);  六元方程组怎么用MATLAB编辑程序？.jpg  六元方程组怎么用MATLAB编辑程序？.jpg |
回复此楼» 猜你喜欢
极限拉扯，如何抉择
已经有15人回复
-
申请科研项目纯靠实力能成功吗？
已经有8人回复
-
中科院海洋所招收海洋材料相关领域博士研究生及研究助理
已经有4人回复
-
40+中年人失业后，还能找到好工作吗？
已经有3人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
|
一个老外给的回答 What is x0=[0:0.1:2:0:0.1:2]; 原来写的是x0=[0:0.1:2:0:0.1:2]; supposed to mean? On my version of matlab (2015a) it does not give me any warning but it creates an empty vector. I suspect it's the cause of the warning and ultimately the cause of the error since you're passing empty for the initial condition. 希望论坛的大神多多帮忙！ |
2楼2015-03-14 17:24:26
3楼2015-03-14 17:25:43
信彼南山
木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 4142.9
- 散金: 1221
- 红花: 16
- 帖子: 1178
- 在线: 233.5小时
- 虫号: 1133529
- 注册: 2010-10-27
- 专业: 导航、制导与传感技术
4楼2015-03-15 01:16:14
信彼南山
木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 4142.9
- 散金: 1221
- 红花: 16
- 帖子: 1178
- 在线: 233.5小时
- 虫号: 1133529
- 注册: 2010-10-27
- 专业: 导航、制导与传感技术
5楼2015-03-15 02:11:05
6楼2015-03-15 07:20:43
7楼2015-03-15 07:53:24
信彼南山
木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 4142.9
- 散金: 1221
- 红花: 16
- 帖子: 1178
- 在线: 233.5小时
- 虫号: 1133529
- 注册: 2010-10-27
- 专业: 导航、制导与传感技术
8楼2015-03-16 00:20:57
9楼2015-03-16 06:54:50
信彼南山
木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 4142.9
- 散金: 1221
- 红花: 16
- 帖子: 1178
- 在线: 233.5小时
- 虫号: 1133529
- 注册: 2010-10-27
- 专业: 导航、制导与传感技术
10楼2015-03-16 11:18:05
200