版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

wyjjf

新虫 (正式写手)

应助: 1 (幼儿园)
金币: 294.1
散金: 46
红花: 3
帖子: 742
在线: 160.5小时
虫号: 2588385
注册: 2013-08-09
专业: 海洋工程

[交流] MATLAB程序总是出错? 已有1人参与

编辑一MATLAB程序，总是出错，按提示信息改动了很多遍
还是不行
不明白了
麻烦帮帮看看怎么修改

编辑程序：
在Matlab下输入：edit，然后将下面两行百分号之间的内容，复制进去，保存?

CODE:

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%function dx=myfun(t,x)

syms a e f g a1 a2 e1 e2 f1 f2 g1 g2

dx = zeros(5,1);

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(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;

a=a1+a2;

e=e1+e2;

f=f1+f2;

g=g1+g2;

a1=1;

a2=1; 

e1=9; 

e2=7; 

f1=56; 

f2=98; 

g1=76; 

g2=665;

dx=dx(:);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

然后，在Matlab下面输入：

CODE:

clear all

clc

x0=[1,1,1,1,1,1];

t0=0:0.1:2;

[t,x]=ode45('myfun',[0,10],x0); %?%ode45会自动调整步长

plot(t,x)                         % ?

legend('x(1)','x(2)','x(3)','x(4)','x(5)','x(6)')

结果出错信息

第一：
dx(1)=x(3)-x(2)+x(4)*x(2)-2*x(5)*x(2)-x(2)*x(2)-x(5)-e-f-g;
有错吗？

二：
Error in ==> funfun\private\odearguments at 110
f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0.

即：
f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0.
[m,n] = size(f0);
if n > 1
error('MATLAB

dearguments:FoMustReturnCol',...
      [funstring(ode) ' must return a column vector.'])
elseif m ~= neq
  msg = sprintf('an initial condition vector of length %d.',m);
  error('MATLAB

dearguments:SizeIC', ['Solving ' funstring(ode) ' requires ' msg]);
End

三：
Error in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, ...
即：
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, ...
options, threshold, rtol, normcontrol, normy, hmax, htry, htspan, ...
dataType] = ...
odearguments(FcnHandlesUsed, solver_name, odeFcn, tspan, y0, ...
options, varargin);
nfevals = nfevals + 1;

出错信息1.jpg

2.jpg

3.jpg

[ Last edited by xiegangmai on 2015-3-13 at 16:40 ]

回复此楼

» 猜你喜欢

A期刊撤稿已经有3人回复
职称评审没过，求安慰已经有34人回复
垃圾破二本职称评审标准已经有17人回复
回收溶剂求助已经有6人回复
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有22人回复
申请26博士已经有5人回复
EST投稿状态问题已经有7人回复
毕业后当辅导员了，天天各种学生超烦已经有4人回复
求助文献已经有3人回复
投稿返修后收到这样的回复，还有希望吗已经有8人回复

» 本主题相关价值贴推荐，对您同样有帮助:

matlab优化报错已经有8人回复
matlab 运行ode45出错，不知道什么原因已经有4人回复
matlab编程约束条件问题已经有10人回复
matlab一运行就出现的问题已经有7人回复
XFree86-devel安装时出错已经有5人回复
关于激光器matlab程序已经有5人回复
请大侠指教matlab的这个错误提示到底是个神马意思？已经有12人回复
有个matlab程序，但是运行不下来，求高手已经有23人回复
matlab程序中调用m文件的问题已经有13人回复
目标函数出现错误已经有3人回复
请matlab高手看一下这个问题怎样编程才好，附程序代码，谢谢！已经有9人回复
matlab-出错提示Undefined method for input arguments type 'double'. 已经有5人回复
请教matlab在linux下的并行问题已经有18人回复
一个matlab 小问题~~ 已经有4人回复
大家爱帮我看看这段MATLAB程序为啥出错了已经有5人回复
Matlab中save问题已经有5人回复
matlab程序问题，总是出错，麻烦大家看看问题在那里？已经有6人回复
同样的MATLAB程序，为何在别人的电脑上运行不了已经有21人回复

1楼 2015-03-12 17:48:17

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

xiegangmai

版主 (职业作家)

我没头衔

仿真EPI: 3
应助: 157 (高中生)
贵宾: 5.735
金币: 40609
散金: 2394
红花: 63
沙发: 4
帖子: 4046
在线: 1330.2小时
虫号: 653607
注册: 2008-11-13
性别: GG
专业: 机械动力学
管辖: 仿真模拟

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

看了看，程序中有几个问题，供参考：
1. function dx=myfun(t,x)中，应先赋值：

CODE:

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;

% syms a e f g a1 a2 e1 e2 f1 f2 g1 g2

2. dx=myfun(t,x)中，x(1) ~ x(6)有6个变量，只有5个方程？无解！

赞一下

回复此楼

明德厚学、求是创新

2楼2015-03-13 16:48:49

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 wyjjf 的主题更新

返回列表