24小时热门版块排行榜    

Znn3bq.jpeg
查看: 962  |  回复: 5

dingx2005

金虫 (正式写手)

[求助] 1stOpt拟合求助 已有1人参与

为了求出参数a,b,c,d。编写了下边的程序,运行时老提示方程不正确。新手不知如何改进,还望大侠指点
1stOpt拟合求助
Title "step 2";
Parameter a,b,c,d;
ConstStr  1=2*b*m*n^2+sqrt(b*d)*m*n+n;
DataSet;
       X, Y =
         0.6        0.51
         0.7        0.477
         0.8        0.461
         0.9        0.449
         1        0.443
EndDataSet;
MinFunction sum(i=1:5)(-y+a*b*m*n^2+c*sqrt(b*d)*m*n)+
                    sum(i=1:5)(-x+m+b*m*(n)^2+sqrt(b*d)*m*n);
回复此楼
人就这一辈子,能够从事自己喜欢的事情,并且可以从中获取生活资本,这是非常幸运的事情。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
这样写:
CODE:
Title "step 2";
Parameter a,b,c,d;
DataSet;
       X, Y =
         0.6        0.51
         0.7        0.477
         0.8        0.461
         0.9        0.449
         1        0.443
EndDataSet;
MinFunction sum(i=1:5,y)(-y+a*b*m*n^2+c*sqrt(b*d)*m*n)+sum(i=1:5,x)(-x+m+b*m*(n)^2+sqrt(b*d)*m*n);
            1=2*b*m*n^2+sqrt(b*d)*m*n+n;

只是结果不对,无穷小,参数有范围限制吧?
2楼2015-10-23 19:04:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingx2005

金虫 (正式写手)

引用回帖:
2楼: Originally posted by dingd at 2015-10-23 19:04:44
这样写:

Title "step 2";
Parameter a,b,c,d;
DataSet;
       X, Y =
         0.6        0.51
         0.7        0.477
         0.8        0.461
         0.9        0.449
         ...

谢谢您的回复
参数都是大于零
人就这一辈子,能够从事自己喜欢的事情,并且可以从中获取生活资本,这是非常幸运的事情。
3楼2015-10-23 22:29:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

不唯一吧:

1:
目标函数值(最小): -6.34
a: 0
b: 1.43801264441738
c: 6.93475569356464
d: 11.8339281140646
m: 0
n: 0.999999999999951

约束函数:
1: 1-(2*b*m*n^2+sqrt(b*d)*m*n+n) = 4.862776848E-014

2:
目标函数值(最小): -6.34
a: 0.690915185725316
b: 0.0244748821230291
c: 0.0396681396914308
d: 0.0143748852975545
m: 1.58458487214934E-16
n: 0.999999999999912

约束函数:
1: 1-(2*b*m*n^2+sqrt(b*d)*m*n+n) = 8.792966355E-014
4楼2015-10-23 22:45:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingx2005

金虫 (正式写手)

谢谢,我也不知道有几个解。但是m,n不是一个值,每组x,y对应一组m,n。
我想再请教一下我的函数处理的算法对不对:
我想用实验数据拟合以下函数,从而求得未知数a,b,c,d.
y=a*b*m*n^2+c*sqrt(b*d)*m*n
x=m+b*m*(n)^2+sqrt(b*d)*m*n
2*b*m*n^2+sqrt(b*d)*m*n+n=1
在编程时将函数处理为MinFunction sum(i=1:5,y)(-y+a*b*m*n^2+c*sqrt(b*d)*m*n)+sum(i=1:5,x)(-x+m+b*m*(n)^2+sqrt(b*d)*m*n);1=2*b*m*n^2+sqrt(b*d)*m*n+n;
不知道这种处理合适不合适?
人就这一辈子,能够从事自己喜欢的事情,并且可以从中获取生活资本,这是非常幸运的事情。
5楼2015-10-24 10:40:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingx2005

金虫 (正式写手)

引用回帖:
4楼: Originally posted by dingd at 2015-10-23 22:45:22
不唯一吧:

1:
目标函数值(最小): -6.34
a: 0
b: 1.43801264441738
c: 6.93475569356464
d: 11.8339281140646
m: 0
n: 0.999999999999951

约束函数:
1: 1-(2*b*m*n^2+sqrt(b*d)*m*n+n) = 4.86277684 ...

谢谢,我也不知道有几个解。但是m,n不是一个值,每组x,y对应一组m,n。
我想再请教一下我的函数处理的算法对不对:
我想用实验数据拟合以下函数,从而求得未知数a,b,c,d.
y=a*b*m*n^2+c*sqrt(b*d)*m*n
x=m+b*m*(n)^2+sqrt(b*d)*m*n
2*b*m*n^2+sqrt(b*d)*m*n+n=1
在编程时将函数处理为MinFunction sum(i=1:5,y)(-y+a*b*m*n^2+c*sqrt(b*d)*m*n)+sum(i=1:5,x)(-x+m+b*m*(n)^2+sqrt(b*d)*m*n);1=2*b*m*n^2+sqrt(b*d)*m*n+n;
不知道这种处理合适不合适?
人就这一辈子,能够从事自己喜欢的事情,并且可以从中获取生活资本,这是非常幸运的事情。
6楼2015-10-24 10:42:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dingx2005 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 297求调剂 +27 GENJIOW 2026-04-07 30/1500 2026-04-09 23:20 by wolf97
[考研] 085400 328分 求调剂 +8 喂你一个大橙子 2026-04-09 12/600 2026-04-09 19:20 by zl8213662
[考研] 288求调剂 +11 代fish 2026-04-09 11/550 2026-04-09 17:59 by qluyyj
[考研] 312求调剂 +3 李鸿飞飞 2026-04-06 3/150 2026-04-09 17:32 by wp06
[考研] 085501机械专硕 302分 不挑专业求调剂 +5 汪某. 2026-04-09 5/250 2026-04-09 15:38 by 蒋皓禹
[考研] 材料工程322 +18 哈哈哈吼吼吼哈 2026-04-07 19/950 2026-04-09 10:44 by cymywx
[考研] 求调剂 +8 吃口冰激凌 2026-04-07 8/400 2026-04-09 08:03 by 5268321
[考研] 调剂 +22 不逢春 2026-04-07 23/1150 2026-04-09 08:01 by Sammy2
[考研] 285求调剂 +12 AZMK 2026-04-05 18/900 2026-04-08 20:43 by 逆水乘风
[考研] 336求调剂,一志愿中科大 +9 墨彧 yuyu 2026-04-06 9/450 2026-04-08 11:24 by 想读书的菌菌
[考研] 材料工程专业日语生求调剂 +9 111623 2026-04-07 9/450 2026-04-07 23:31 by 一只好果子?
[考研] 生物学调剂 可调剂到生物与医药 +3 李政莹 2026-04-06 3/150 2026-04-06 19:02 by macy2011
[考研] 工科370求调剂 +3 溏心煎鸡蛋 2026-04-05 3/150 2026-04-06 10:55 by 这是一个无聊的
[考研] 0703求调剂383分 +9 W55j 2026-04-03 9/450 2026-04-06 06:50 by houyaoxu
[考研] 求调剂 +7 张.1 2026-04-05 7/350 2026-04-05 20:40 by 啵啵啵0119
[考研] 材料调剂 +9 革微桂 2026-04-04 9/450 2026-04-05 08:27 by 544594351
[考研] 考研调剂 +11 小sun要好运 2026-04-04 11/550 2026-04-05 08:02 by qlm5820
[考研] 400分求调剂 +3 尴尬且挠头 2026-04-04 3/150 2026-04-04 08:41 by jp9609
[考研] 一志愿中国石油大学化学工程323分求调剂 +4 化工专硕323分 2026-04-03 6/300 2026-04-03 22:12 by dongzh2009
[考研] 求调剂 +8 akdhjs 2026-04-03 8/400 2026-04-03 18:17 by 戴维ING
信息提示
请填处理意见