24小时热门版块排行榜    

查看: 1434  |  回复: 7
【悬赏金币】回答本帖问题,作者xingyaoliu将赠送您 5 个金币

xingyaoliu

铁虫 (小有名气)

[求助] 1stopt公式参数拟合 已有3人参与

各位虫友,我尝试用1stopt对下面的公式进行拟合,想求解出里面的6个参数,但是运行一直出错,本人是新手,望不吝赐教。谢谢。
Parameter E;A;N;B;C;D;
Variable x,y;
y=(1+E*x)exp[(A+N*E)*x+B*x^2+C*x^3+D*x^4]
data
0 0
5 0.35
10 0.7
20 1.85
50 6
100 12
200 54
回复此楼
橡胶硫化
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lambert2014

捐助贵宾 (著名写手)

2楼2019-07-13 21:30:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xingyaoliu

铁虫 (小有名气)

引用回帖:
2楼: Originally posted by lambert2014 at 2019-07-13 01:30:22
我知道matlab拟合法

可以赐教下吗
橡胶硫化
3楼2019-07-13 22:25:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lambert2014

捐助贵宾 (著名写手)

%
x1=temperature(:,1);
x2=temperature(:,4);
x3=temperature(:,5);  
y=temperature(:,3);   % 目标 获得 f=g(x1,x2,x3)的函数

   X = [ones(size(x1))  x1  x2  x1.*x3   x2.*x3 ]; % 选取的函数具有形式: y=b0+b1*x1+b2*x2+b3*x1*x3+b4*x2*x3
       [b,bint,r,rint,stats] = regress(y,X);
b'    % 上述展开系数 b0-b4
stats(1)    % 相关系数  反映拟合函数的精度  接近1为优

x1fit = min(x1)max(x1)-min(x1))/20:max(x1);
x2fit = min(x2)max(x2)-min(x2))/20:max(x2);
x3fit = min(x3)max(x3)-min(x3))/20:max(x3);

YFIT =  b(1)+b(2)*X1FIT+b(3)*(X1FIT.*X3FIT)+ b(4)*(X2FIT.*X3FIT); %目标函数形式
4楼2019-07-15 13:43:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lambert2014

捐助贵宾 (著名写手)

【答案】应助回帖

具体百度regress函数。 刚才发的是我的算例函数
5楼2019-07-15 13:44:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

小毛病不少:
CODE:
Parameter E,A,N,B,C,D;
Variable x,y;
Function y=(1+E*x)*exp((A+N*E)*x+B*x^2+C*x^3+D*x^4);
data;
0 0
5 0.35
10 0.7
20 1.85
50 6
100 12
200 54

模型公式有问题:“(A+N*E)”,A和N两个参数完全可以用一个表示即可,否则过拟合(A和N不唯一)。
6楼2019-07-18 10:47:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雨夜孤星

金虫 (正式写手)

【答案】应助回帖

函数前要加Function ,要不然它识别不了啊,函数中需要加乘号的地方也不能吝啬,比如说exp前面。。。
Parameter E;A;N;B;C;D;
Variable x,y;
Function y=(1+E*x)*exp[(A+N*E)*x+B*x^2+C*x^3+D*x^4];
data;
0 0
5 0.35
10 0.7
20 1.85
50 6
100 12
200 54
稳住,我们能赢
7楼2019-07-20 12:09:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xingyaoliu

铁虫 (小有名气)

引用回帖:
2楼: Originally posted by lambert2014 at 2019-07-13 01:30:22
我知道matlab拟合法

还请赐教
橡胶硫化
8楼2019-07-25 22:33:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xingyaoliu 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 274求调剂 +6 cgyzqwn 2026-03-01 11/550 2026-03-05 23:39 by ,!?
[考研] 一志愿东北大学080500,英一数二343分求调剂 +5 Seele_v5 2026-03-04 5/250 2026-03-05 13:26 by Iveryant
[基金申请] 本子写完了,给DS兄弟看了,得了92分 +6 Doma 2026-03-01 10/500 2026-03-05 10:40 by ms629
[考研] 一志愿清华深研院材料专硕294分,专业课111分,本科中南大学材料,有六级,有工作经验 +3 H14528 2026-03-04 3/150 2026-03-05 10:35 by ms629
[考研] 化学 0703求调剂 总分293 一志愿211 +9 土土小虫 2026-03-03 11/550 2026-03-05 09:26 by vfdbghn
[考研] 求调剂院校 +6 云朵452 2026-03-02 12/600 2026-03-04 23:17 by 云朵452
[考研] 一志愿武汉理工大学-085602-总分296分-求调剂 +7 紫川葡柚 2026-03-04 7/350 2026-03-04 21:04 by kakakapanpan
[考研] 0703化学 一志愿东北大学 279分求调剂 +4 聂乐平 2026-03-03 4/200 2026-03-04 16:31 by zhukairuo
[考博] 26申博 求博导 +3 爱读书的小帅 2026-02-28 5/250 2026-03-03 22:17 by 远方13579
[基金申请] 请问大家,研究风险与应对措施那里, 大家都怎么写呢 ? +3 cauasen 2026-03-02 3/150 2026-03-03 19:43 by Tide man
[考研] 292求调剂 +3 sgbl 2026-03-03 3/150 2026-03-03 17:30 by barlinike
[考研] 一志愿华中科技大学,化学专业344分,求调剂 +6 邢xing1 2026-03-02 6/300 2026-03-03 15:46 by stewie_jz
[考研] 化工专硕348,一志愿985求调剂 +8 弗格个 2026-02-28 11/550 2026-03-03 14:55 by relum
[考研] 299求调剂 +5 kkcoco25 2026-03-02 9/450 2026-03-03 12:55 by 公瑾逍遥
[论文投稿] 通讯作者写谁,问题是你意想不到的问题 15+3 阿尔法啊 2026-03-01 3/150 2026-03-03 09:13 by 北京莱茵润色
[考研] 338求调剂 +5 18162027187 2026-03-02 6/300 2026-03-03 07:52 by njzyff
[考研] 298求调剂 +10 人间唯你是清欢 2026-02-28 14/700 2026-03-02 22:49 by 人间唯你是清欢
[考研] 290分材料工程085601求调剂 数二英一 +8 llx0610 2026-03-02 9/450 2026-03-02 22:09 by 无际的草原
[考研] 材料085601调剂 +5 多多子. 2026-03-02 5/250 2026-03-02 19:15 by zhukairuo
[基金申请] 刚录用,没有期刊号,但是在线可看的论文可以放为代表作吗 10+3 arang1 2026-03-01 3/150 2026-03-01 16:43 by babero
信息提示
请填处理意见