3楼: Originally posted by 月只蓝 at 2014-11-27 20:29:38
你说的这种情况比较少见。
麻烦给出完整代码。

4楼: Originally posted by brucexmc at 2014-11-27 21:20:06
非常感谢您的回复，我的调用函数代码为：
clear all;
clc;
options=odeset('reltol',1e-5);
=ode45(@fun,,,options);
plot(t,y(:,1),t,y(:,2));

函数fun的定义为
function dx=fun(t,x)

if t<1
    ...

5楼: Originally posted by 月只蓝 at 2014-11-27 22:07:05
MATLAB求解结果无误。
代码中
fdata=fopen('data.txt','a');
fprintf(fdata,'%f\n',t);
其中data.txt数据是什么。...

6楼: Originally posted by brucexmc at 2014-11-27 23:55:13
我是把时间t写入到一个名为“data”的txt文件中了
...

7楼: Originally posted by 月只蓝 at 2014-11-28 08:44:03
可以确定，当前代码没问题，计算结果可靠。t y的计算均不存在问题。...

9楼: Originally posted by shikang999 at 2014-11-29 13:51:10
1、ode45计算流程数值算法书上已经有，这里不详细叙述
2、你函数里记录的t和你外面输出的那个t不存在线性相关性（特别是一些算法做了技巧性处理，比如加入自适应或者复化分段自控误差等等），所以不影响你外面t按大 ...