| 查看: 2879 | 回复: 5 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
zhangholmes至尊木虫 (著名写手)
逍遥门主
|
[求助]
怎样用一组参数同时拟合两个曲线--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)); |
回复此楼» 猜你喜欢
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有123人回复
-
物理学I论文润色/翻译怎么收费?
已经有250人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有5人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 使用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人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
4楼2013-01-16 14:02:18
zhangholmes
至尊木虫 (著名写手)
逍遥门主
- 应助: 17 (小学生)
- 金币: 17045
- 散金: 50
- 红花: 12
- 帖子: 1242
- 在线: 226.3小时
- 虫号: 967701
- 注册: 2010-03-11
- 性别: GG
- 专业: 能源化工
送鲜花一朵 3楼2013-01-16 13:30:11
zhangholmes
至尊木虫 (著名写手)
逍遥门主
- 应助: 17 (小学生)
- 金币: 17045
- 散金: 50
- 红花: 12
- 帖子: 1242
- 在线: 226.3小时
- 虫号: 967701
- 注册: 2010-03-11
- 性别: GG
- 专业: 能源化工
5楼2013-01-17 09:59:10
zhangholmes
至尊木虫 (著名写手)
逍遥门主
- 应助: 17 (小学生)
- 金币: 17045
- 散金: 50
- 红花: 12
- 帖子: 1242
- 在线: 226.3小时
- 虫号: 967701
- 注册: 2010-03-11
- 性别: GG
- 专业: 能源化工
6楼2013-01-17 10:00:43
