| 查看: 1118 | 回复: 2 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
求大神---MATLAB-ode23tb-微分方程组--运行失败
|
||
|
有一个包含117个方程的常微分方程组。 用MATLAB ode23tb求解。 问题是这样的：第一次我把tspan设为[0 12*60] （为方便描述，单位为分钟，12为小时数）。得到运算结果[T0 Y0] (设s为计算步数，size(T0)=[s 1] size(Y0)=[s 117]) 现在，我要把这12个小时分成几个时间段，在每个时间段上要改变一个初始值（想象一下在不同时间段给小鼠注射不同药物，然后观察体内某些物质浓度变化）。 我是这样做的，如果是分成两段，第一段运行，初始值y0，tspan=[0 6*60]； 第二段运行时，取第一段Y0中最后一个行向量YL，并改动其中某一项，然后作为初始值y1，在运行下6个小时。这时候运行就失败了(不收敛)！并且只取YL不作改动，运行也失败！ 请问这是什么情况？ |
回复此楼» 猜你喜欢
346求调剂[0856]
已经有7人回复
-
一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂！
已经有3人回复
-
310求调剂
已经有3人回复
-
机械专硕299求调剂至材料
已经有4人回复
-
070300化学319求调剂
已经有7人回复
-
08工科 320总分 求调剂
已经有6人回复
-
一志愿天津大学化学工艺专业（081702）315分求调剂
已经有12人回复
-
307求调剂
已经有3人回复
-
265求调剂
已经有3人回复
-
301求调剂
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 微分方程组求参数问题，求高人指教，感谢
已经有11人回复
- MATLAB微分方程参数拟合问题，求大神
已经有7人回复
- 一阶常微分方程组，matlab进行数值求解
已经有3人回复
- 请教mathematic问题？关于解微分方程组的和画图的一些问题。
已经有19人回复
- matlab解一阶常微分方程组
已经有4人回复
- matlab有限差分法求解复杂的偏微分方程组
已经有19人回复
- matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较
已经有4人回复
- Matlab求解偏微分方程组
已经有22人回复
- 请问如何求解二元一阶微分方程组
已经有24人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 高金求助matlab解微分方程组
已经有12人回复
- matlab解微分方程组
已经有15人回复
- LABVIEW四阶龙格库塔法求解一阶微分方程组
已经有12人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】请教非齐次常微分方程组的解析解法
已经有4人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 【求助】有限元法求解二阶偏微分方程组
已经有6人回复
3楼2013-07-11 10:50:48
|
%------ first dose-------- [T0 Y0]=ode23tb(@model,[0 6*60],y0); y01=Y0(length(Y0),  ;y01(1)=y01(1)+0.3;% 0.3-> value of the 2nd dose tt=T0(1:length(T0)-1); yy=Y0(1:length(Y0)-1); % ------second dose -------- [T1 Y1]=ode23tb(@model,[0 6*60],y01); T1=T1+6*60; T=[tt;T1]; Y=[yy;Y1]; figure(1) plot(T/60,Y(:,50));% 50-> system output |
2楼2013-07-11 10:47:33
