解微分方程与代数方程的混合方程组出错求助

rabbie · 发表于 2007-11-21 15:44

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

我要解一个33阶的方程组，前11个是微分方程，如
y1'=...
y2'=...
后面22个则都是代数方程。
我参考help文件里面的做法，设
M=[eye(11),zeros(11,22);zeros(22,33)]
然后把M放进option里，然后
[t,y] = ode23t(@f,tspan,y0,options);
function out = f(t,y)
out=[...]
这样来求解，可是解不出来，程序报错：
??? Error using ==> minus
Matrix dimensions must agree.
请问这是为什么呀？？解微分方程与代数方程混合的方程组有什么更好的办法吗？？

[ 本帖最后由 eight 于 2007-11-22 10:01 编辑 ]

xjzuo · 发表于 2007-11-21 23:47

请将原问题贴出。

eight · 发表于 2007-11-22 10:01

原帖由 rabbie 于 2007-11-21 15:44 发表

 登录/注册后可看大图

我要解一个33阶的方程组，前11个是微分方程，如
y1'=...
y2'=...
后面22个则都是代数方程。
我参考help文件里面的做法，设
M=[eye(11),zeros(11,22);zeros(22,33)]
然后把M放进option里，然后
[t,y] = o ...

请看看本版置顶帖关于出错问题整理的帖子

账号		自动登录	找回密码
密码			我要加入

[编程技巧] 解微分方程与代数方程的混合方程组出错求助

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。