应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » MATLAB 数据拟合的问题
31/1返回列表
查看: 431  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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]

谢谢大家了!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-08-18 10:47:47
已阅   回复此楼   关注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
赞一下(1人) 回复此楼
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 的主题更新
31/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-01 10:43 by BCdvBAEpt3zm
[论文投稿] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 10:35 by BCdvBAEpt3zm
[基金申请] 论文的挂名怎么评判贡献？基金申请中，排名第一的挂名文章 +3 biily 2026-10-25 4/200 2026-11-01 10:32 by biily
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 10:27 by BCdvBAEpt3zm
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 10:27 by BCdvBAEpt3zm
[论文投稿] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 X5bHxOULxFqb 2026-10-30 26/1300 2026-11-01 10:22 by BCdvBAEpt3zm
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-01 07:01 by YXM1FDCSE8Bn
[基金申请] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 X5bHxOULxFqb 2026-10-30 21/1050 2026-11-01 06:46 by YXM1FDCSE8Bn
[基金申请] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-01 05:48 by YXM1FDCSE8Bn
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-01 05:48 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-01 05:43 by YXM1FDCSE8Bn
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 05:41 by YXM1FDCSE8Bn
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 04:04 by YXM1FDCSE8Bn
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 10/500 2026-11-01 03:51 by YXM1FDCSE8Bn
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 4/200 2026-11-01 00:41 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 X5bHxOULxFqb 2026-10-30 12/600 2026-11-01 00:15 by YXM1FDCSE8Bn
[考博] 2027博士申请-临床检验诊断学、病原生物学、微生物方向 50+3 yyh@ 2026-10-28 4/200 2026-10-31 12:47 by 可淡不可忘
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +8 zuesli 2026-10-27 14/700 2026-10-30 22:44 by 好运连连公主
[文学芳草园] 闲聊 +3 myrtle 2026-10-25 6/300 2026-10-29 15:28 by myrtle
[考博] 申博求助 20+3 Katrinae 2026-10-26 4/200 2026-10-27 09:56 by 北京莱茵编辑
信息提示
关闭
请填处理意见
关闭