版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

zhangholmes

至尊木虫 (著名写手)

逍遥门主

应助: 17 (小学生)
金币: 16720
散金: 50
红花: 12
帖子: 1242
在线: 226.1小时
虫号: 967701
注册: 2010-03-11
性别: GG
专业: 能源化工

[求助] 怎样用一组参数同时拟合两个曲线--matlab

各位大侠，小弟由于研究需要，要进行参数估计，遇到了一个难题：怎样把涉及到2条曲线的同一组参数估计出来，请各位大侠鼎力相助。方法是先用fmincon函数拟合，得到的结果作为lsqnonlin函数的初值进行运算，最后得到参数值。请大侠能给我编程的m文件最好。
a0=[1 1 1 1 1 1 1 1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1]; %参数初值
lb=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0];
ub=[+inf +inf +inf +inf +inf +inf +inf +inf +inf +inf +inf +inf +inf +inf +inf +inf]; %参数上下限
x=[0.51333,0.48393,0.45268,0.42557,0.40512,0.38257,0.36301,0.34194,0.32613,0.31155,0.29826,0.28502,0.27131,0.25321,0.23865,0.21536,0.19563,0.17771];
y1=[0.50536,0.52537,0.54375,0.56052,0.57522,0.59563,0.61684,0.63566,0.64712,0.65980,0.67163,0.68357,0.69210,0.70533,0.71154,0.71536,0.71907,0.72182];
y2=[2.53630,2.19802,1.89159,1.59469,1.35026,1.14236,0.93376,0.78531,0.66720,0.58129,0.49402,0.41640,0.34402,0.27812,0.23651,0.18564,0.15012,0.12590];
%y1,y2对应两条曲线r1,r2的实验值
k1=a (1).*exp ((-a(8).*exp (-a(16)./4098.802).*((0.133+0.173*x)./(0.483-0.346*x))-a (9))./4098.802);
k2=a (2).*exp ((-a(8).*exp (-a(16)./4098.802).*((0.133+0.173*x)./(0.483-0.346*x))-a (10))./4098.802);
k3=a (3).*exp ((-a(8).*exp (-a(16)./4098.802).*((0.133+0.173*x)./(0.483-0.346*x))-a (11))./4098.802);
k4=a (4).*exp ((-a(8)*exp (-a(16)./4098.802).*((0.133+0.173*x)./(0.483-0.346*x))-a (12))./4098.802);
k5=a (5).*exp ((-a(8)*exp (-a(16)./4098.802).*((0.133+0.173*x)./(0.483-0.346*x))-a (13))./4098.802);
k6=a (6).*exp ((-a(8)*exp (-a(16)./4098.802).*((0.133+0.173*x)./(0.483-0.346*x))-a (14))./4098.802);
k7=a (7).*exp ((-a(8)*exp (-a(16)./4098.802).*((0.133+0.173*x)./(0.483-0.346*x))-a (15))./4098.802); %要求的参数是a(1)-a(16)这16个参数
r1=-(x.*k1.*k2.*k3.*(-k5-k7))./(-(x.*k1.*k2-k3.*(-x.*k1-k2-x.*k4-k6)).*(-k5-k7)-x.*k3.*k4.*(-x.*k1+k7));
r2=-(x^2).*k1.*k3.*k4.*k5./((x.*k1.*k2-k3.*(-x.*k1-k2-x.*k4-k6)).*(-k5-k7)-x.*k3.*k4.*(-x.*k1+k7));

回复此楼

» 猜你喜欢

投稿chemical physical letters不送审？已经有2人回复
南开大学物理学院张书辉副教授招收凝聚态物理理论方向博士生、硕士生已经有22人回复
物理学I论文润色/翻译怎么收费? 已经有253人回复
杂志要求提供official language-editing certificate 已经有4人回复
深圳大学2026年秋博士招生-物理学-活性胶体方向-高永祥课题组已经有16人回复
论物质与能量的统一模型及物理现象解释已经有0人回复
基于基元I统一理论的数学相关应用推导已经有0人回复
基元I统一理论：宇宙本质、层级演化与修炼文明的本源规律已经有1人回复
基元I理论下三大核心空间现象精准推导与细节解析已经有0人回复
基于基元 I 统一理论的反重力理论推导已经有0人回复
基于基元I统一理论的量子力学本源推导已经有1人回复

» 本主题相关价值贴推荐，对您同样有帮助:

使用DSLF拟合吸附等溫線時，怎麼樣才能讓得到的參數為identifiable(參數估計值唯一) 已经有17人回复
matlab数据拟合问题求助已经有7人回复
origin拟合出的不同曲线如何在一个图中导出，谢了已经有9人回复
1stopt或matlab如何求解以下的非线性方程并拟合出相应曲线？已经有13人回复
MATLAB能否将一条曲线拆分成两条曲线已经有3人回复
matlab 拟合反应动力学参数结果很差。大家帮忙看一下已经有14人回复
请教matlab反应动力学参数估计遇到的问题，谢谢已经有15人回复
matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较已经有4人回复
matlab拟合方程参数时初值的选择已经有15人回复
如何用matlab画出两组实验测定数据概率密度曲线，求两条曲线重叠部分面积已经有8人回复
如何做出拟合曲线求值已经有24人回复
高手帮忙拟合一组曲线，解决另加分！已经有19人回复
origin多组数据如何用一条曲线拟合已经有11人回复
有关origin曲线拟合的问题已经有7人回复
origin多组数据用一条曲线拟合已经有20人回复
【求助】如何同时拟合两条曲线或者两个函数！【已解决】已经有7人回复
【求助】如何使得matlab拟合的曲线强制经过一个点已经有8人回复
【求助】使用MATLAB怎么实现拟合力场参数的程序？【已解决】已经有15人回复
请教关于曲线拟合的问题已经有1人回复
【求助】在MATLAB中如何将离散点拟合成曲面图已经有18人回复
一条曲线如何进行分段拟合已经有13人回复

探索真理比占有真理更可贵

1楼 2013-01-16 10:04:16

已阅关注TA 给TA发消息送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

计算强帖: 4
应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.7小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
zhangholmes: 金币+10, ★有帮助 2013-01-18 09:40:17

内容已删除

赞一下

回复此楼

» 本帖已获得的红花（最新10朵）

zhangholmes

2楼2013-01-16 11:09:19

已阅关注TA 给TA发消息送TA红花 TA的回帖

zhangholmes

至尊木虫 (著名写手)

逍遥门主

应助: 17 (小学生)
金币: 16720
散金: 50
红花: 12
帖子: 1242
在线: 226.1小时
虫号: 967701
注册: 2010-03-11
性别: GG
专业: 能源化工

送鲜花一朵

引用回帖:

2楼: Originally posted by dingd at 2013-01-16 11:09:19
用1stOpt倒是很方便，关键是不依赖初值，只是你的模型似乎可以化简，否则有过拟合现象，参数组值难以唯一：

ConstStr For(i=1:7)(k=a*exp((-a*exp (-a/4098.802)*((0.133+0.173*x)/(0.483-0.346*x))-a)/4098.802 ...

您好，我按您给的代码试了一下，提示函数表达式定义有错，是什么原因呢，我用的是1.5版本的

赞一下

回复此楼

探索真理比占有真理更可贵

3楼2013-01-16 13:30:11

已阅关注TA 给TA发消息送TA红花 TA的回帖