chinamiracle 发表于 2008-5-13 03:01

请分析下doc ode45 中示例程序产生错误的问题

function ha
clear all,clc
options = odeset('RelTol',1e-4,'AbsTol',);
= ode45(@rigid,,,options);
plot(t,y(:,1),'-',t,y(:,2),'-.',t,y(:,3),'.')

function dy = rigid(t,y)
dy = zeros(3,1);    % a column vector
dy(1) = y(2) * y(3);
dy(2) = -y(1) * y(3);
dy(3) = -0.51 * y(1) * y(2);
我发现好怪了,一下想不通。这是ode45帮助文档里的例子程序,我把它分开存放运行可以做出图形,而同时存放 在一个m文件里面却出错这是为什么?我觉得它应该可以的啊

ch_j1985 发表于 2008-5-13 07:12

原帖由 chinamiracle 于 2008-5-13 03:01 发表 http://www.chinavib.com/forum/images/common/back.gif
function ha
clear all,clc
options = odeset('RelTol',1e-4,'AbsTol',);
= ode45(@rigid,,,options);
plot(t,y(:,1),'-',t,y(:,2),'-.',t,y(:,3),'.')

function dy = ri ...

在我电脑上运行没有报错,能说明一下报的什么错吗?

eight 发表于 2008-5-13 11:30

原帖由 chinamiracle 于 2008-5-13 03:01 发表 http://www.chinavib.com/forum/images/common/back.gif
function ha
clear all,clc
options = odeset('RelTol',1e-4,'AbsTol',);
= ode45(@rigid,,,options);
plot(t,y(:,1),'-',t,y(:,2),'-.',t,y(:,3),'.')

function dy = ri ... 发帖前请多看看置顶帖,清晰地描述你遇到的问题

chinamiracle 发表于 2008-5-13 12:03

哦呵呵,是我的问题,昨天晚上傻了哈一下没注意 ,函数名没一致
页: [1]
查看完整版本: 请分析下doc ode45 中示例程序产生错误的问题