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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 拜求！matlab拟合动力学参数
152/2返回列表上一页12
查看: 2327  |  回复: 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
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭