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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 拜求！matlab拟合动力学参数
152/2返回列表上一页12
查看: 2825  |  回复: 14
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

小小莓邮

铁虫 (小有名气)
送红花一朵
引用回帖:
8楼: Originally posted by pdl9527 at 2014-02-19 15:16:14
随便编了一下，拟合结果有点问题，估计是k的初值k0取得有问题或者数据太少或者实验数据有问题。。。
程序如下，其中k0是我随便设的初值，k1=k1,k2=k1',k3=k8,k4=k4*,k5=k5,k6=f



function piadatfit2
clear ...

可能是我数据太少了吧！这个不行就后面接着补，虽然没整出来，但觉着还是很厉害，学习了！
赞一下 回复此楼
11楼2014-02-19 17:42:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

天花板.

新虫 (初入文坛)

【答案】应助回帖

感谢参与，应助指数 +1
引用回帖:
7楼: Originally posted by 钟静 at 2014-02-19 13:10:49
我还是觉得数据太少了，哪怕可以拟合出来，结果的意义也不大了。其中还有一个导数，拟合的时候是要把这个算出来吧，如果间隔是一个小时的话求出的速率有代表性不？...

数据预处理一下，包括插值以及利用插值数据进行求导得到r_RH.我把模型简化了一下，其中有些参数是乘积然后放在一起的，这显然可以作为一个参数来看待。-y'=c1/(c1+c2*c_RH)*c2*c_RH*sqrt(c3*c_ROOH)+c4*c_ROOH.
1.新建parameter_fix函数
function F=parameter_fix(c,xdata)
F=c(1)./(c(1)+c(2)*xdata(1,)*c(2).*xdata(1,.*sqrt(c(3).*xdata(2,)+c(4)*xdata(2,;
F=-F;
2.数据预处理及参数优化
clear
clc
close all;
t=[0 1 2 3 4 5 6 8 10];
c_RH=[8.32016 7.37673 6.74433 5.92613 5.68806 5.01584 6.09249 5.34827 4.9961];
c_ROOH=[0 0.21115 0.7379 1.10562 1.29586 1.58327 1.36404 1.84065 2.00986];
dt=0.1;
t_x=0:dt:10;
y_RH=interp1([t,10.1],[c_RH,4.9961-(5.34827-4.9961)],[t_x,t_x(size(t_x,2))+dt]);
dy_RH=y_RH(2:size(t_x,2)+1)-y_RH(1:size(t_x,2));
ydata=interp1(t_x,dy_RH,t);
xdata=[c_RH;c_ROOH];
result=lsqcurvefit(@(c,xdata) parameter_fix(c,xdata),[1;1;1;1],xdata,ydata,[-inf;-inf;0;-inf],[inf;inf;inf;inf]);
赞一下 回复此楼
12楼2014-02-19 17:52:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

天花板.

新虫 (初入文坛)
引用回帖:
12楼: Originally posted by 天花板. at 2014-02-19 17:52:03
数据预处理一下，包括插值以及利用插值数据进行求导得到r_RH.我把模型简化了一下，其中有些参数是乘积然后放在一起的，这显然可以作为一个参数来看待。-y'=c1/(c1+c2*c_RH)*c2*c_RH*sqrt(c3*c_ROOH)+c4*c_ROOH.
1 ...

那个笑脸是个代码。。。尴尬居然会这样，是“: )"，引号里面的部分，识别成笑脸了。
赞一下 回复此楼
13楼2014-02-19 17:55:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

13488143026

金虫 (著名写手)
引用回帖:
8楼: Originally posted by pdl9527 at 2014-02-19 17:16:14
随便编了一下，拟合结果有点问题，估计是k的初值k0取得有问题或者数据太少或者实验数据有问题。。。
程序如下，其中k0是我随便设的初值，k1=k1,k2=k1',k3=k8,k4=k4*,k5=k5,k6=f



function piadatfit2
clear ...

k0是什么。。。
回复此楼
14楼2015-04-07 20:10:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

飞行系

新虫 (正式写手)
引用回帖:
8楼: Originally posted by pdl9527 at 2014-02-19 15:16:14
随便编了一下，拟合结果有点问题，估计是k的初值k0取得有问题或者数据太少或者实验数据有问题。。。
程序如下，其中k0是我随便设的初值，k1=k1,k2=k1',k3=k8,k4=k4*,k5=k5,k6=f



function piadatfit2
clear ...

朋友 能不能留个QQ向你请教点问题
赞一下 回复此楼
15楼2015-05-11 21:51:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小小莓邮 的主题更新
152/2返回列表上一页12
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 jzg3xk1po8 2026-07-31 6/300 2026-08-01 15:46 by y4XuShGyllvT
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 jzg3xk1po8 2026-07-31 11/550 2026-08-01 15:43 by y4XuShGyllvT
[基金申请] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 jzg3xk1po8 2026-07-31 11/550 2026-08-01 15:43 by y4XuShGyllvT
[基金申请] 小木虫看见有人已经知道结果了 +17 1234567wang 2026-07-25 22/1100 2026-08-01 15:43 by iaeyuan
[教师之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 jzg3xk1po8 2026-07-31 8/400 2026-08-01 15:42 by y4XuShGyllvT
[基金申请] 面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！ +16 低垂的野花 2026-07-31 21/1050 2026-08-01 15:35 by jklily
[基金申请] 2026年国自然面上资助率 +15 布布和一二 2026-07-30 19/950 2026-08-01 15:18 by 40862026
[基金申请] 微信指数没变化，科研之友没阅读 +14 wangze12014 2026-07-28 18/900 2026-07-31 18:35 by lzchaihuo
[教师之家] 基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 8/400 2026-07-31 18:13 by beefly
[高分子] HXDI做水性聚氨酯乳液，是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了，估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[考博] 申请2027年材料/化学类博士 10+3 考研老狗? 2026-07-25 5/250 2026-07-30 19:33 by 白色木鱼
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
[基金申请] 今年的WR进展到哪一步了？ +5 wsgjhwz 2026-07-25 9/450 2026-07-30 12:59 by 可淡不可忘
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 时间戳他又来了 +14 晓晓爱翠翠 2026-07-26 17/850 2026-07-29 16:30 by 呵呵二哥
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
[基金申请] 同事接到电话了，我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
[基金申请] 准备明年的基金了 +5 Tide man 2026-07-27 5/250 2026-07-28 09:25 by newfuzzy1
[分析] 水蒸气蒸馏法乳化现象 +3 月光娇娇 2026-07-26 3/150 2026-07-26 23:25 by caiyun
信息提示
关闭
请填处理意见
关闭