24小时热门版块排行榜    

查看: 1105  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-29 21:44 by 4qO8zhfjr2AX
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-29 21:33 by 4qO8zhfjr2AX
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-29 21:22 by 4qO8zhfjr2AX
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 ASdOkHsho7FD 2026-08-28 11/550 2026-08-29 20:18 by 63GmZBFV29mG
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 17:30 by 4FFAWE8HcgUD
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 14:13 by jCd0dEvKHShX
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +5 gdd2018 2026-08-28 10/500 2026-08-29 14:10 by Jacob678
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 ASdOkHsho7FD 2026-08-28 9/450 2026-08-29 14:03 by jCd0dEvKHShX
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 14:02 by jCd0dEvKHShX
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-29 11:18 by jCd0dEvKHShX
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见