应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 1stOpt拟合求助
61/1返回列表
查看: 921  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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);
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
人就这一辈子，能够从事自己喜欢的事情，并且可以从中获取生活资本，这是非常幸运的事情。
1楼 2015-10-22 22:24:15
已阅   回复此楼   关注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 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 346求调剂[0856] +3 WayneLim327 2026-03-16 6/300 2026-03-19 11:21 by WayneLim327
[考研] 一志愿天大材料与化工（085600）总分338 +5 蔡大美女 2026-03-13 5/250 2026-03-19 10:44 by 是小刘呀～
[考研] 267一志愿南京工业大学0817化工求调剂 +10 SUICHILD 2026-03-12 10/500 2026-03-19 09:51 by Delta2012
[考研] 274求调剂 +6 S.H1 2026-03-18 6/300 2026-03-19 09:34 by 花店相见
[考研] 0703化学调剂，求各位老师收留 +10 秋有木北 2026-03-14 10/500 2026-03-19 05:52 by anny19840123
[考研] 295求调剂 +3 一志愿京区211 2026-03-18 5/250 2026-03-18 17:03 by zhaoqian0518
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +8 Ncdx123456 2026-03-13 9/450 2026-03-18 14:40 by haxia
[考研] 0703化学调剂 +4 pupcoco 2026-03-17 7/350 2026-03-18 12:14 by djl2006
[考研] 材料，纺织，生物（0856、0710），化学招生啦 +3 Eember. 2026-03-17 9/450 2026-03-18 10:28 by Eember.
[考研] 301求调剂 +9 yy要上岸呀 2026-03-17 9/450 2026-03-18 08:58 by 无际的草原
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 0703 物理化学调剂 +3 我可以上岸的对� 2026-03-13 5/250 2026-03-16 10:50 by 我可以上岸的对�
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 材料与化工 323 英一+数二+物化，一志愿：哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[基金申请] 现在如何回避去年的某一个专家，不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 307求调剂 +5 超级伊昂大王 2026-03-12 5/250 2026-03-13 15:56 by 棒棒球手
[考研] 工科材料085601 279求调剂 +8 困于星晨 2026-03-12 10/500 2026-03-13 15:42 by ms629
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[考研] 321求调剂（食品/专硕） +3 mxcz321 2026-03-12 6/300 2026-03-13 08:45 by xc321
信息提示
关闭
请填处理意见
关闭