| 查看: 1708 | 回复: 5 | |||
| 【悬赏金币】回答本帖问题，作者木虫小生将赠送您 50 个金币 | |||
[求助]
拟合曲线，参数不恒定已有3人参与
|
|||
|
2*((1-x)-(1-x)^(-2))*(-0.6494-0.6082*((1-x)^2+2*(1-x)^(-1)-3)+0.8259*((1-x)^2+2*(1-x)^(-1)-3)^2)+3000*((1-x)^2+0.5*(1-x)^(-4))*(k0+k1*b1*(1-exp(-x/b1/3000))+k2*b2*(1-exp(-x/b2/3000))+k3*b3*(1-exp(-x/b3/3000)))*(a1+a2*3000*((1-x)-(1-x)^(-2))) 2*((1-x)-(1-x)^(-2))*(-0.6494-0.6082*((1-x)^2+2*(1-x)^(-1)-3)+0.8259*((1-x)^2+2*(1-x)^(-1)-3)^2)+2260*((1-x)^2+0.5*(1-x)^(-4))*(k0+k1*b1*(1-exp(-x/b1/2260))+k2*b2*(1-exp(-x/b2/2260))+k3*b3*(1-exp(-x/b3/2260)))*(a1+a2*2260*((1-x)-(1-x)^(-2))) 2*((1-x)-(1-x)^(-2))*(-0.6494-0.6082*((1-x)^2+2*(1-x)^(-1)-3)+0.8259*((1-x)^2+2*(1-x)^(-1)-3)^2)+1510*((1-x)^2+0.5*(1-x)^(-4))*(k0+k1*b1*(1-exp(-x/b1/1510))+k2*b2*(1-exp(-x/b2/1510))+k3*b3*(1-exp(-x/b3/1510)))*(a1+a2*1510*((1-x)-(1-x)^(-2))) 上面是所用模型，区别是1510，2260，3000三个值不同，想得到k0,k1,k2,k3,b1,b2,b3,a1,a2的恒定值，或者说一个平均值，使用的是matlab工具箱拟合的，参数变化太大，取平均不行，所以着急想问问可有解决方法，这是一个过拟合，解决这样的过拟合应该怎么办，请大家帮帮忙，非常感谢！@月只蓝@beefly |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：libolin3@tal.com - 附件 1 : 1.xlsx
- 附件 2 : 2.xlsx
- 附件 3 : 3.xlsx
2019-01-01 21:14:43, 30.56 K
2019-01-01 21:14:43, 30.88 K
2019-01-01 21:14:43, 31.83 K
» 猜你喜欢
法国博士招聘
已经有0人回复
-
法国博士招聘
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有60人回复
-
法国博士招聘
已经有0人回复
-
法国博士招聘
已经有0人回复
-
超导三元凸包图
已经有0人回复
-
最近遇到这样一个问题
已经有10人回复
-
论文接收，散金币了！
已经有175人回复
-
vasp算居里温度
已经有0人回复
-
郑州大学光电器件方向招收一名博士
已经有6人回复
lipenggg
专家顾问 (正式写手)
-
专家经验: +65 - 应助: 122 (高中生)
- 金币: 4841.3
- 红花: 43
- 帖子: 722
- 在线: 435.7小时
- 虫号: 3929201
- 注册: 2015-06-17
- 性别: GG
- 专业: 流体力学
- 管辖: 计算模拟
2楼2019-01-02 10:40:47
独孤神宇
版主 (知名作家)
- 应助: 488 (硕士)
- 贵宾: 0.008
- 金币: 29512.8
- 散金: 802
- 红花: 120
- 沙发: 1
- 帖子: 5333
- 在线: 850.2小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
月只蓝: 金币+10, 感谢帮助！ 2019-01-03 21:57:49
感谢参与，应助指数 +1
月只蓝: 金币+10, 感谢帮助！ 2019-01-03 21:57:49
|
QuickReg = 20; VarConstant W=[3000,2260,1510]; Variable x,y; Parameter k0,k1,k2,k3,b1,b2,b3,a1,a2; Function y=2*((1-x)-(1-x)^(-2))*(-0.6494-0.6082*((1-x)^2+2*(1-x)^(-1)-3)+0.8259*((1-x)^2+2*(1-x)^(-1)-3)^2) +W*((1-x)^2+0.5*(1-x)^(-4))*(k0+k1*b1*(1-exp(-x/b1/W))+k2*b2*(1-exp(-x/b2/W)) +k3*b3*(1-exp(-x/b3/W)))*(a1+a2*W*((1-x)-(1-x)^(-2))); Data; 0 -0.102524 0.000016 0.027285 ............... Data; ............... ..................... Data; .................... ................ ************************ 均方差(RMSE): 1.24392092366486 残差平方和(SSR): 5573.51603012112 相关系数(R): 0.998413934604383 相关系数之平方(R^2): 0.996830384812205 修正R平方(Adj. R^2): 0.996825312550693 确定系数(DC): 0.964529222376253 F统计(F-Statistic): 31306.0930058816 参数 最佳估算 -------------------- ------------- k0 12.2461642539881 k1 -23258503.8196115 k2 25440146.8111859 k3 20353209.282357 b1 3.83420353685459E-5 b2 5.34485605257666E-6 b3 5690563.09328754 a1 6.80432042395383E-6 a2 6.25937243062727E-10 |
3楼2019-01-02 11:07:56
hzlhm
至尊木虫 (著名写手)
- 应助: 372 (硕士)
- 金币: 16128.8
- 红花: 51
- 帖子: 2560
- 在线: 571.9小时
- 虫号: 1517335
- 注册: 2011-11-30
- 性别: GG
- 专业: 常微分方程与动力系统
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
独孤神宇: 金币+5, 鼓励交流 2019-01-03 14:29:23
感谢参与，应助指数 +1
独孤神宇: 金币+5, 鼓励交流 2019-01-03 14:29:23
|
用数据1拟合，得到 a1=6.7524，a2=-0.00070073，b1=-13017.6577，b2=2.4265e-06，b3=-36478.0183，k0=5.4223e-06，k1=2.5056，k2=18.2815，k3=3.2361 解决思路：自定义拟合函数y=func(a,x)，其内容：y1=。。。，y2=。。。，y3=。。。，y=(y1+y2+y3)/3；用nlinfit（）函数拟合。 A=xlsread('1.xlsx'); %数据 x=A(:,1); y=A(:,2); x0=[5.4892 -0.0010972 2.512 3.7974e-05 1.2398 2.9089e-05 1.7848 5.2093 2.5155] [a,r,J] = nlinfit(x,y,@func,x0);  123.jpg |
4楼2019-01-03 13:32:18
james阿布
新虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 1148.7
- 散金: 70
- 红花: 7
- 沙发: 3
- 帖子: 1409
- 在线: 40.1小时
- 虫号: 9642089
- 注册: 2018-07-27
- 性别: GG
- 专业: X线与CT、电子与离子束、放
5楼2019-01-04 18:55:07
南木126
木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 2289.9
- 散金: 1198
- 红花: 7
- 帖子: 1510
- 在线: 165.1小时
- 虫号: 2023859
- 注册: 2012-09-23
- 性别: GG
- 专业: 影像医学与生物医学工程其
6楼2019-06-24 12:09:03
100