| 查看: 2818 | 回复: 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)); |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有138人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有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人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
zhangholmes: 金币+10, ★有帮助 2013-01-18 09:40:17
感谢参与，应助指数 +1
zhangholmes: 金币+10, ★有帮助 2013-01-18 09:40:17
| 内容已删除 |
» 本帖已获得的红花（最新10朵）
zhangholmes2楼2013-01-16 11:09:19
zhangholmes
至尊木虫 (著名写手)
逍遥门主
- 应助: 17 (小学生)
- 金币: 17040
- 散金: 50
- 红花: 12
- 帖子: 1242
- 在线: 226.3小时
- 虫号: 967701
- 注册: 2010-03-11
- 性别: GG
- 专业: 能源化工
3楼2013-01-16 13:30:11
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
4楼2013-01-16 14:02:18
zhangholmes
至尊木虫 (著名写手)
逍遥门主
- 应助: 17 (小学生)
- 金币: 17040
- 散金: 50
- 红花: 12
- 帖子: 1242
- 在线: 226.3小时
- 虫号: 967701
- 注册: 2010-03-11
- 性别: GG
- 专业: 能源化工
5楼2013-01-17 09:59:10
zhangholmes
至尊木虫 (著名写手)
逍遥门主
- 应助: 17 (小学生)
- 金币: 17040
- 散金: 50
- 红花: 12
- 帖子: 1242
- 在线: 226.3小时
- 虫号: 967701
- 注册: 2010-03-11
- 性别: GG
- 专业: 能源化工
6楼2013-01-17 10:00:43
投票: