24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 367  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 277跪求调剂 +3 1915668 2026-03-27 4/200 2026-03-27 11:39 by zhshch
[考研] 0703化学/290求调剂/本科经历丰富/工科也可 +7 丹青奶盖 2026-03-26 8/400 2026-03-27 11:24 by 丹青奶盖
[考研] 求调剂 +3 张zz111 2026-03-27 4/200 2026-03-27 11:07 by 帕尔马拉特
[考研] 324求调剂 +5 hanamiko 2026-03-26 5/250 2026-03-27 10:33 by wangjy2002
[考研] 一志愿哈尔滨工业大学材料与化工方向336分 +3 辰沐5211314 2026-03-26 3/150 2026-03-27 10:24 by 尽舜尧1
[考研] 材料调剂 +8 匹克i 2026-03-23 8/400 2026-03-27 08:11 by hypershenger
[考研] 求调剂,一志愿 南京航空航天大学大学 ,080500材料科学与工程学硕 +4 @taotao 2026-03-26 5/250 2026-03-27 08:10 by hypershenger
[考研] 一志愿太原理工安全工程300分,求调剂 +3 0857求调剂. 2026-03-24 3/150 2026-03-27 01:18 by 求知若渴lz
[考研] 一志愿北京化工大学材料与化工(085600)296求调剂 +9 稻妻小编 2026-03-26 9/450 2026-03-26 16:16 by 不吃魚的貓
[考研] 0856求调剂 +8 zhn03 2026-03-25 9/450 2026-03-26 13:42 by zzll406
[考研] 调剂310 +3 温柔的晚安 2026-03-25 4/200 2026-03-25 23:16 by peike
[考研] 一志愿上海交大生物与医药专硕324分,求调剂 +6 jiajunX 2026-03-22 6/300 2026-03-25 23:05 by licg0208
[考研] 考研调剂 +6 来好运来来来 2026-03-21 7/350 2026-03-25 22:43 by 418490947
[考研] 299求调剂 +7 某某某某位 2026-03-21 8/400 2026-03-25 20:34 by 热情沙漠
[考研] 【2026考研调剂】制药工程 284分 求相关专业调剂名额 +4 袁奂奂 2026-03-25 8/400 2026-03-25 14:32 by lbsjt
[考研] 293求调剂 +7 加一一九 2026-03-24 7/350 2026-03-25 12:02 by userper
[考研] 一志愿吉林大学材料与化工303分求调剂 +4 为学666 2026-03-24 4/200 2026-03-25 11:27 by BruceLiu320
[考研] 318求调剂 +5 plum李子 2026-03-21 8/400 2026-03-25 09:26 by aa331100
[考研] 080500求调剂 +3 zzzzfan 2026-03-24 3/150 2026-03-24 16:38 by barlinike
[考研] A区线材料学调剂 +5 周周无极 2026-03-20 5/250 2026-03-20 21:33 by laoshidan
信息提示
请填处理意见