24小时热门版块排行榜    

查看: 416  |  回复: 2

yangang0310

银虫 (初入文坛)

[求助] MATLAB 数据拟合的问题 已有1人参与

有个问题困扰为好多天了,希望虫友 能解答。
现有一模型 : y=a*f*x.^t+b*(f*x).^2+e*(f*x).^1.5
其中,f ,x是自变量,y是因变量,且都有数据。a,b,e,t是未知参数。
MATLAB里有个 fminsearch函数好像能解决这类问题。可是,我是个MATLAB菜鸟,希望广大虫友不吝赐教啊。
另 附下数据:

数据分为三组:分别为f=50,100,150.

f=50;x=[0.2:0.1:1.9]  y=[0.0148        0.0296        0.0501        0.0747        0.1031        0.1352        0.1709        0.2115        0.2575        0.3088        0.3644        0.4267        0.4952        0.5698   0.6572         0.7755        1.0057        1.3361]

f=100; x=[0.2:0.1:1.9]  y=[0.0368         0.0799        0.1371        0.2074        0.2909        0.3866        0.4948        0.6177        0.7582  0.916        1.09        1.285 1.4992 1.7222   1.9709  2.263        2.7734        3.5511]
f=150; x=[0.2:0.1:1.9]  y=[0.072  0.1568  0.2689        0.4411        0.5684        0.7533        0.9618        1.2058        1.4794        1.8289        2.1415        2.5103   2.9037        3.326   3.784        4.3331   5.18890  6.4422]

谢谢大家了!
回复此楼

» 猜你喜欢

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

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

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
1stOpt求解:
CODE:
VarConstant f=[50,100,150];
Variable x,y;
Function  y=a*f*x^t+b*(f*x)^2+e*(f*x)^1.5;
Data;
x=[0.2:0.1:1.9];
y=[0.0148,0.0296,0.0501,0.0747,0.1031,0.1352,0.1709,0.2115,0.2575,0.3088,0.3644,0.4267,0.4952,0.5698,0.6572,0.7755,1.0057,1.3361];
Data;
x=[0.2:0.1:1.9];
y=[0.0368,0.0799,0.1371,0.2074,0.2909,0.3866,0.4948,0.6177,0.7582,0.916,1.09,1.285,1.4992,1.7222,1.9709,2.263,2.7734,3.5511;
Data;
x=[0.2:0.1:1.9];
y=[0.072,0.1568,0.2689,0.4411,0.5684,0.7533,0.9618,1.2058,1.4794,1.8289,2.1415,2.5103,2.9037,3.326,3.784,4.3331,5.18890,6.4422];

均方差(RMSE):0.0449907733571378
残差平方和(SSE):0.10930516311276
相关系数(R): 0.9993394941954
相关系数之平方(R^2): 0.998679424658717
确定系数(DC): 0.997962342657576
F统计(F-Statistic): 3848.72254811153

参数                  最佳估算
--------------------        -------------
a        2.23611194227693E-5
t        9.77849410813785
b        2.80454370471318E-5
e        0.000491364528602323
2楼2014-08-18 11:37:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangang0310

银虫 (初入文坛)

引用回帖:
2楼: Originally posted by dingd at 2014-08-18 11:37:54
1stOpt求解:

VarConstant f=;
Variable x,y;
Function  y=a*f*x^t+b*(f*x)^2+e*(f*x)^1.5;
Data;
x=;
y=;
Data;
x=;
y=[0.0368,0.0799,0.1371,0.2074,0.2909,0.3866,0.4948,0.6177,0.7582,0.916,1.09,1 ...

非常之感谢,但是t有个参考值,一般在1~3之间,用1stOpt求解,是不是初值选取的不对啊?
另外,MATLAB里用fminsearch函数求解,能否给点意见,谢谢了
3楼2014-08-18 15:20:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yangang0310 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 我不理解! +19 Edward_pc 2026-08-26 27/1350 2026-09-02 12:08 by Ai19911123
[基金申请] 要骂人了,新模版改版就是要淡化问题凝练这种虚的东西,结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 16/800 2026-09-02 12:07 by Ai19911123
[基金申请] 科研人应该花精力去思考如何解决问题,而不是去凝练问题 +4 瞬息宇宙 2026-09-01 8/400 2026-09-02 10:35 by 瞬息宇宙
[基金申请] 学科评审组评审是指会评吗? +5 瞬息宇宙 2026-08-31 5/250 2026-09-02 10:13 by 雨冰共舞
[文学芳草园] 梦想 +5 myrtle 2026-08-26 8/400 2026-09-02 02:42 by Leogzhya
[基金申请] 面上合作单位盖章 +6 ssyjh 2026-08-27 9/450 2026-09-01 20:01 by huagongfeihu
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 17/850 2026-09-01 19:51 by 超级无敌华子
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
信息提示
请填处理意见