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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 拜求！matlab拟合动力学参数
152/2返回列表上一页12
查看: 2439  |  回复: 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
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 304求调剂 +6 司空. 2026-03-18 6/300 2026-03-18 23:03 by 星空星月
[考研] 085601专硕，总分342求调剂，地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 311求调剂 +6 26研0 2026-03-15 6/300 2026-03-18 14:43 by haxia
[考研] 297求调剂 +8 戏精丹丹丹 2026-03-17 8/400 2026-03-18 14:30 by laoshidan
[考研] 070300化学319求调剂 +6 锦鲤0909 2026-03-17 6/300 2026-03-18 13:22 by Iveryant
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[考研] 考研化学学硕调剂，一志愿985 +4 张vvvv 2026-03-15 6/300 2026-03-17 17:15 by ruiyingmiao
[考研] 材料与化工专硕调剂 +5 heming3743 2026-03-16 5/250 2026-03-17 14:03 by 勇敢太监王公公
[论文投稿] 有没有大佬发小论文能带我个二作 +3 增锐漏人 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
[考研] 326求调剂 +4 诺贝尔化学奖觊� 2026-03-15 7/350 2026-03-16 17:11 by 诺贝尔化学奖觊�
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 0856专硕279求调剂 +5 加油加油！? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 255求调剂 +3 李嘉慧， 2026-03-12 4/200 2026-03-14 16:58 by 有只狸奴
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
信息提示
关闭
请填处理意见
关闭