24小时热门版块排行榜    

查看: 2626  |  回复: 3
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

kittyzoe

木虫 (正式写手)

[求助] matlab 求解隐式常微分方程时用ode15i 函数出现的问题求助

在matlab中利用ode15i函数求解隐式常微分方程时,需要用decic函数先确定初值。求初值的语句如下:
%求初值
t0=0;
Y00=[a;0;b;0];
fix_Y00=ones(4,1);%题目中Y0全给出时,全为1
dY0=[0;0;0;0];
fix_dY0=zeros(4,1);%题目中dY0没给出,全为0
[Y02,dY02]=decic(@function,t0,Y00,fix_Y00,dY0,fix_dY0);

但是运行以后出现以下错误:
??? Error using ==> decic>sls at 131
Try freeing 4 fixed components.

Error in ==> decic at 80
    [dy,dyp] = sls(res,dfdy,dfdyp,neq,free_y,free_yp);

查了好久也不知道是出现什么原因了,不知高手们有没有遇到过这种问题呢?请指点一二啊。
回复此楼

» 猜你喜欢

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

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

WF1987

铜虫 (小有名气)

怎么不把所有程序都贴出来呢?
4楼2012-02-27 17:02:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

kittyzoe

木虫 (正式写手)

我刚尝试把第三句:fix_Y00=ones(4,1);%题目中Y0全给出时,全为1
改为:fix_Y00=zeros(4,1);
就不再出现错误,但是此时给定的Y0初值也没法固定了啊,不知道有什么办法解决么
2楼2011-11-09 09:15:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kittyzoe

木虫 (正式写手)

当出现如下警告的时候是不是说明方程没办法求解呢?
>> [t,Y]=ode15i(@function,Tspan,Y0,Yp0,options);
Warning: Failure at t=8.015625e-003.  Unable to meet integration tolerances without reducing
the step size below the smallest value allowed (2.847722e-017) at time t.
> In ode15i at 416
>>
3楼2011-11-09 09:24:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 生物学学硕求调剂 +6 小羊睡着了? 2026-03-23 9/450 2026-03-24 20:21 by 小羊睡着了?
[考研] 307求调剂 +5 超级伊昂大王 2026-03-24 5/250 2026-03-24 15:46 by 星空星月
[考研] 277分求调剂,跨调材料 +3 考研调剂lxh 2026-03-24 3/150 2026-03-24 13:52 by JourneyLucky
[考研] 材料与化工328分调剂 +4 。,。,。,。i 2026-03-23 4/200 2026-03-24 11:03 by 544594351
[考研] 一志愿哈工大,085400,320,求调剂 +3 gdlf9999 2026-03-24 3/150 2026-03-24 10:08 by 搏击518
[考研] 350求调剂 +6 weudhdk 2026-03-19 6/300 2026-03-23 15:47 by tangyuan0840221
[考研] 求老师收我 +3 zzh16938784 2026-03-23 3/150 2026-03-23 12:56 by ztnimte
[考研] 280分求调剂 一志愿085802 +4 PUMPT 2026-03-22 7/350 2026-03-22 22:13 by 星空星月
[考研] 319求调剂 +4 小力气珂珂 2026-03-20 4/200 2026-03-22 15:53 by ColorlessPI
[考研] 260求调剂 +3 朱芷琳 2026-03-20 4/200 2026-03-22 15:12 by 朱芷琳
[基金申请] 山东省面上项目限额评审 +4 石瑞0426 2026-03-19 4/200 2026-03-22 08:50 by Wei_ren
[考研] 286求调剂 +10 Faune 2026-03-21 10/500 2026-03-21 23:34 by 314126402
[考研] 求调剂 +3 13341 2026-03-20 3/150 2026-03-21 18:28 by 学员8dgXkO
[考研] 一志愿深大,0703化学,总分302,求调剂 +4 七月-七七 2026-03-21 4/200 2026-03-21 18:20 by 学员8dgXkO
[考研] 279求调剂 +5 红衣隐官 2026-03-21 5/250 2026-03-21 14:59 by lature00
[考研] 一志愿重庆大学085700资源与环境专硕,总分308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-21 00:39 by JourneyLucky
[考研] 330求调剂 +4 小材化本科 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[考研] 材料学硕297已过四六级求调剂推荐 +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[考研] 一志愿西南交通 专硕 材料355 本科双非 求调剂 +5 西南交通专材355 2026-03-19 5/250 2026-03-20 21:10 by JourneyLucky
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
信息提示
请填处理意见