版块导航
正在加载中...
客户端APP下载
论文辅导
调剂小程序
登录
注册
帖子
帖子
用户
本版
应《网络安全法》要求,自2017年10月1日起,未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用,请尽快对帐号进行手机号验证,感谢您的理解与支持!
24小时热门版块排行榜
>
论坛更新日志
(4214)
>
导师招生
(666)
>
考研
(541)
>
虫友互识
(281)
>
文献求助
(218)
>
休闲灌水
(145)
>
考博
(80)
>
招聘信息布告栏
(70)
>
硕博家园
(55)
>
博后之家
(54)
>
论文投稿
(49)
>
基金申请
(37)
>
找工作
(29)
>
教师之家
(25)
>
绿色求助(高悬赏)
(23)
>
公派出国
(14)
小木虫论坛-学术科研互动平台
»
计算模拟区
»
计算模拟
»
lsqnonlin函数拟合微分方程组参数拟合问题
5
1/1
返回列表
查看: 2820 | 回复: 10
只看楼主
@他人
存档
新回复提醒
(忽略)
收藏
在APP中查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖
yanyl09
禁虫
(初入文坛)
本帖内容被屏蔽
» 猜你喜欢
2026年循环经济功能材料国际会议(ICFMCE 2026)
已经有0人回复
2026年第五届电气、电子与信息工程国际会议(ISEEIE 2026)
已经有0人回复
物理学I论文润色/翻译怎么收费?
已经有76人回复
哈尔滨理工大学物理系招收物理学考研调剂
已经有5人回复
0702一志愿吉大B区求调剂
已经有5人回复
求调剂
已经有0人回复
0702一志愿吉大B区求调剂有论文
已经有0人回复
请问还有没有用Latex写文章的小伙伴们?
已经有0人回复
光学工程学硕调剂信息
已经有23人回复
欢迎加入课题组
已经有0人回复
» 本主题相关价值贴推荐,对您同样有帮助:
origin 8.0 对数函数拟合
已经有13人回复
请问一个有理函数的曲线拟合过程如何做,给出数据,请求拟合,
已经有14人回复
origin8 双曲函数拟合一问
已经有3人回复
【求助】催化反应动力学
已经有5人回复
【求助】matlab用lsqnonlin非线性拟合的问题【已解决】
已经有8人回复
【求助】如何同时拟合两条曲线或者两个函数!【已解决】
已经有7人回复
【求助】二元函数拟合
已经有12人回复
【求助】紧急求助非线性拟合问题
已经有3人回复
【讨论】怎么对数据进行特定函数的拟合
已经有4人回复
origin7.5建立函数拟合为什么会这样?
已经有11人回复
请教关于origin中自定义函数非线性回归拟合问题!
已经有10人回复
1楼
2011-12-05 14:54:16
已阅
关注TA
给TA发消息
送TA红花
TA的回帖
yanyl09
禁虫
(初入文坛)
本帖内容被屏蔽
高级回复
3楼
2011-12-06 09:02:08
已阅
关注TA
给TA发消息
送TA红花
TA的回帖
查看全部 11 个回答
dingd
铁杆木虫
(职业作家)
计算强帖: 4
应助: 1641
(讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.7小时
虫号: 291104
注册: 2006-10-28
【答案】应助回帖
微分方程拟合推荐使用1stOpt,更强大方便些。
赞
一下
回复此楼
2楼
2011-12-05 17:01:18
已阅
关注TA
给TA发消息
送TA红花
TA的回帖
yanyl09
禁虫
(初入文坛)
本帖内容被屏蔽
4楼
2011-12-06 15:58:49
已阅
关注TA
给TA发消息
送TA红花
TA的回帖
dbb627
荣誉版主
(著名写手)
专家经验: +4
计算强帖: 12
应助: 289
(大学生)
贵宾: 0.589
金币: 24640.4
散金: 551
红花: 61
沙发: 1
帖子: 1246
在线: 1794.8小时
虫号: 149791
注册: 2005-12-29
性别: GG
专业: 污染控制化学
管辖:
计算模拟
【答案】应助回帖
yanyl09(金币+10): 非常感谢,公式的修改对我启发很大,我再试试,谢谢!!! 2011-12-06 22:37:09
我觉得你的模型不太对
我该了下,你试试,其实第一个生物量还是有问题,你没有考虑死亡率,估计改后会好些
CODE:
function AA
k0=[0.1 2 2];
x0=[0.01 54.026 0.2435];
t1=0:6:48;
tspan=0:6:48;
yexp=load('sample_file.txt')';
lb=[-10 -10 -10];ub=[10 10 10];
[k,resnorm,resid,exitflag,output,lambda,jacobian]=lsqnonlin(@ObjFunc4LNL8,k0,[0 0 0],[],[],x0,yexp);
ci=nlparci(k,resid,jacobian);
y1=[yexp(:,1)]';
y2=[yexp(:,2)]';
y3=[yexp(:,3)]';
[t4plot x4plot]=ode45(@kineticseqs8,[tspan(1) tspan(end)],x0,[],k);
plot(t1,y1,'bo',t1,y2,'g+',t1,y3,'r*',t4plot,x4plot,'k-'),
legend('Exp-biomass','Exp-glucose','Exp-butanol','Model'),
xlabel('time(h)'),ylabel('biomass dry weight(g/l),glucose(g/l),butanol(g/l)')
fprintf('Estimated Parameter\n'),
fprintf('\tk=%.4f %.4f\n',k(1),ci(1,2)-k(1))
fprintf('\tk=%.4f %.4f\n',k(2),ci(2,2)-k(2))
fprintf('\tk=%.4f %.4f\n',k(3),ci(3,2)-k(3))
[t,x]=ode45(@kineticseqs8,tspan,x0,[],k);
figure,plot(x,resid,'*')
xlabel('biomass dry wight(g/l)、glucose(g/l)、butanol(u/ml)'),
ylabel('Residuals(g/l,u/ml)')
function dxdt=kineticseqs8(t,x,k)
% dxdt=[k(1)*(1-x(2)./12)*x(1)
% -k(2)*x(1).^2
% k(3)*k(1)*(1-x(2)./12)*x(1)];
dxdt=[k(1)*x(1)*(1-x(1)/3)
-k(2)*x(1)*x(2)
k(3)*x(2)*x(1)];
function f=ObjFunc4LNL8(k,x0,yexp)
x0=[0.01 54.026 0.2435];
tspan=[0:6:48]';
[t1 x]=ode45(@kineticseqs8,tspan,x0,[],k);
y(:,1)=x(:,1);y(:,2)=x(:,2);y(:,3)=x(:,3);
f1=y(:,1)-yexp(:,1);f2=y(:,2)-yexp(:,2);f3=y(:,3)-yexp(:,3);
f=[f1 f2 f3];
赞
一下
回复此楼
» 本帖已获得的红花(最新10朵)
yanyl09
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
5楼
2011-12-06 17:55:00
已阅
关注TA
给TA发消息
送TA红花
TA的回帖
查看全部 11 个回答
最具人气热帖推荐
[查看全部]
作者
回/看
最后发表
[
考研
]
285化工学硕求调剂(081700)
+9
柴郡猫_
2026-03-12
9/450
2026-03-17 10:18
by
Sammy2
[
考研
]
考研调剂
+3
淇ya_~
2026-03-17
5/250
2026-03-17 09:25
by
Winj1e
[
考研
]
08工科 320总分 求调剂
+3
梨花珞晚风
2026-03-17
3/150
2026-03-17 09:04
by
晗晗调剂上岸
[
考研
]
278求调剂
+3
Yy7400
2026-03-13
3/150
2026-03-17 08:24
by
laoshidan
[
考研
]
一志愿985,本科211,0817化学工程与技术319求调剂
+5
Liwangman
2026-03-15
5/250
2026-03-16 17:10
by
我的船我的海
[
考研
]
304求调剂
+4
ahbd
2026-03-14
4/200
2026-03-16 16:48
by
我的船我的海
[
考研
]
294求调剂
+3
Zys010410@
2026-03-13
4/200
2026-03-15 10:59
by
zhq0425
[
考研
]
中科大材料与化工319求调剂
+3
孟鑫材料
2026-03-14
3/150
2026-03-14 20:10
by
ms629
[
考研
]
328求调剂
+3
5201314Lsy!
2026-03-13
6/300
2026-03-14 15:31
by
hyswxzs
[
考研
]
266求调剂
+4
学员97LZgn
2026-03-13
4/200
2026-03-14 08:37
by
zhukairuo
[
考研
]
0805,333求调剂
+3
112253525
2026-03-10
3/150
2026-03-13 23:42
by
JourneyLucky
[
考研
]
341求调剂
+3
番茄头---
2026-03-10
3/150
2026-03-13 23:07
by
JourneyLucky
[
考研
]
(081700)化学工程与技术-298分求调剂
+12
11啦啦啦
2026-03-11
35/1750
2026-03-13 21:25
by
JourneyLucky
[
考研
]
281求调剂
+9
Koxui
2026-03-12
11/550
2026-03-13 20:50
by
Koxui
[
考研
]
289求调剂
+3
李政莹
2026-03-12
3/150
2026-03-13 11:02
by
求调剂zz
[
考研
]
070303一志愿西北大学学硕310找调剂
+3
d如愿上岸
2026-03-13
3/150
2026-03-13 10:43
by
houyaoxu
[
考研
]
270求调剂 085600材料与化工专硕
+3
YXCT
2026-03-11
3/150
2026-03-13 10:13
by
houyaoxu
[
考博
]
2026年博士申请
+3
QwQwQW10
2026-03-11
3/150
2026-03-12 17:58
by
gxch43
[
考研
]
420求调剂
+4
莫向外求11
2026-03-10
6/300
2026-03-12 14:41
by
ruiyingmiao
[
考研
]
279求调剂
+3
莫xiao
2026-03-10
4/200
2026-03-11 08:06
by
斩魂滴兔子!
信息提示
关闭
请填处理意见
关闭
确定