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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 拜求！matlab拟合动力学参数
152/2返回列表上一页12
查看: 3100  |  回复: 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
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 厉娜 +3 张丹枫1991 2026-11-09 5/250 2026-11-09 16:53 by 张丹枫1991
[考博] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 12/600 2026-11-09 13:26 by Equinoxhua
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 15/750 2026-11-09 13:19 by Equinoxhua
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +5 5sJFDKp2kwuI 2026-11-06 16/800 2026-11-09 11:47 by GwrQVVR6Oenx
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 5sJFDKp2kwuI 2026-11-06 15/750 2026-11-09 11:41 by GwrQVVR6Oenx
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 10/500 2026-11-09 08:00 by IhnEQauGmT3K
[考研] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 12/600 2026-11-09 06:32 by l48hFd8Nwae2
[公派出国] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-09 00:14 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 23:01 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
[教师之家] 延安大学回应“副教授留绝笔信”：确已离世，正核实举报内容 +3 yexuqing 2026-11-04 5/250 2026-11-07 13:22 by 苏东坡二世
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
信息提示
关闭
请填处理意见
关闭