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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 1stOpt拟合求助
61/1返回列表
查看: 922  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 307求调剂 +5 冷笙123 2026-03-17 5/250 2026-03-19 14:12 by kakakapanpan
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 10/500 2026-03-19 13:23 by houyaoxu
[考研] 材料080500调剂求收留 +4 一颗meteor 2026-03-13 4/200 2026-03-19 10:32 by 30660438
[考研] 一志愿吉林大学材料学硕321求调剂 +6 Ymlll 2026-03-18 9/450 2026-03-19 10:28 by 星空星月
[考研] 0703化学调剂 +8 妮妮ninicgb 2026-03-15 12/600 2026-03-19 05:55 by anny19840123
[考研] 280求调剂 +6 咕噜晓晓 2026-03-18 7/350 2026-03-18 11:25 by 无际的草原
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] 一志愿苏州大学材料工程（085601）专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 材料与化工专硕调剂 +5 heming3743 2026-03-16 5/250 2026-03-17 14:03 by 勇敢太监王公公
[考研] 278求调剂 +3 Yy7400 2026-03-13 3/150 2026-03-17 08:24 by laoshidan
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[考研] 311求调剂 +3 冬十三 2026-03-13 3/150 2026-03-13 20:41 by JourneyLucky
[硕博家园] 085600 260分求调剂 +3 天空还下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空还下雨么
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-13 3/150 2026-03-13 10:43 by houyaoxu
信息提示
关闭
请填处理意见
关闭