24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zhangholmes

至尊木虫 (著名写手)

逍遥门主

应助: 17 (小学生)
金币: 17040
散金: 50
红花: 12
帖子: 1242
在线: 226.3小时
虫号: 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));

回复此楼

» 猜你喜欢

津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名已经有0人回复
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有266人回复
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）已经有6人回复
【修正版】单环用磁约束低速电子实验方案（简化版）已经有0人回复
桂林理工大学物理学专业招收调剂，还有三个名额！！！已经有22人回复
考博自荐已经有4人回复
山东大学第二批博士研究生招生已经有0人回复
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学已经有0人回复
《电磁学》教材推荐已经有1人回复
【急招】合肥工大核聚变材料计算方向2026级工程博士生已经有4人回复

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

使用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

【答案】应助回帖

如果用的是非正版软件，就无法保证了。

赞一下

回复此楼

4楼2013-01-16 14:02:18

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

查看全部 6 个回答

zhangholmes

至尊木虫 (著名写手)

逍遥门主

应助: 17 (小学生)
金币: 17040
散金: 50
红花: 12
帖子: 1242
在线: 226.3小时
虫号: 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的回帖