【调剂】北京石油化工学院2024年16个专业接受调剂
查看: 811  |  回复: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 博士论文被抄袭 +15 和尚敲小木鱼 2024-04-22 24/1200 2024-04-23 13:33 by 和尚敲小木鱼
[考博] 申博成果界定是根据Jcr分区还是中科院分区 +4 我属驴核动力驴 2024-04-22 4/200 2024-04-23 12:18 by lyfbangong
[教师之家] 国内科研环境差导致国内科研水平差的根源 +8 zju2000 2024-04-18 8/400 2024-04-23 10:40 by Quakerbird
[考研] 0854-0855调剂 +4 shangannum1 2024-04-21 7/350 2024-04-23 09:46 by alkb
[教师之家] 刚博士毕业,入职普通一本高校 +16 保持肌肤 2024-04-20 16/800 2024-04-23 08:31 by zxf984
[论文投稿] 一直找不到审稿人 +4 lizhengke06 2024-04-21 5/250 2024-04-22 23:37 by topedit
[论文投稿] 国家重点研发计划论文可以发MDPI里的期刊吗 +6 lin_2018 2024-04-17 6/300 2024-04-22 17:32 by 20094958
[考博] 25申博记录贴 +6 我属驴核动力驴 2024-04-18 12/600 2024-04-22 15:50 by 诩雾
[教师之家] 期末给学生划重点都是什么话术啊 +13 luokereng 2024-04-20 14/700 2024-04-22 14:58 by luokereng
[访问学者] CSC的访问学者申请,没有个评审意见,也不知道怎么改,还有必要申请吗 +3 flyingship 2024-04-20 3/150 2024-04-22 13:52 by louistom_521
[论文投稿] 研究光催化的,好中的三四区 20+3 sl.0117 2024-04-20 3/150 2024-04-22 09:53 by bobvan
[找工作] 工作问题求助 +3 向前?@ 2024-04-18 4/200 2024-04-22 09:40 by zhaojiang427
[论文投稿] 编辑是选国外的好还是国内的好。 +8 lizhengke06 2024-04-20 8/400 2024-04-22 08:58 by cuiyunjian
[考博] 山东大学化学与化工学院刘磊课题组博士研究生招生启事 +6 sdorganic 2024-04-17 8/400 2024-04-21 21:49 by lmr1212
[公派出国] CSC在职申请求问 +8 Amber0919 2024-04-16 8/400 2024-04-20 18:33 by 371122gyd
[考研] 本科211 277分求调剂 +3 0854通信工程 2024-04-18 5/250 2024-04-20 17:55 by 0854通信工程
[考博] 24/25读博求博导 +4 宝23 2024-04-16 4/200 2024-04-19 17:18 by wangzhe_bs
[论文投稿] CCS Chemistry投稿求助 30+3 wfqtriumph 2024-04-17 4/200 2024-04-19 08:56 by Bletilla
[基金申请] 国自然青年基金只有一篇能第一标注,其他都是第二标注能结题吗? 10+4 lice_1987 2024-04-18 6/300 2024-04-19 08:29 by kudofaye
[考研] 工学322,英语六级求调剂 +3 domax 2024-04-17 10/500 2024-04-17 21:25 by 东方花园小万
信息提示
请填处理意见