24小时热门版块排行榜    

查看: 1564  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

danfreda

铜虫 (小有名气)

[求助] 1stopt非线性拟合求助 已有2人参与

需拟合的方程如下(亦附件图片所示):
y = A*2669/1.5043^i * [1 + B*(A/1.5043^i) + C*(A/1.5043^i)^2 + D*(A/1.5043^i)^3 + E*(A/1.5043^i)^4]
需拟合的变量:A,B,C,D,E
i 为控制变量
i 从1到5变化,对应 5个y 值(拟合五个参数用5个yi值可以求解吧?如果不行可舍弃参数E及式子中含E的部分),y值如data部分所示。

编写1stopt程序如下,问题是无法运行出结果,求大神指导纠错:
Title "test";
Parameters a,b,c,d,e;
Variable y[Output];
StartProgram;
var i:interger;
begin
     for i:= 0 to DataLength -1 do begin
     y:=a*2669/(1.5043^i)*[1+b*(a/(1.5043)^i)+c*(a/(1.5043)^i)^2+d*(a/(1.5043)^i)^3+e*(a/(1.5043)^i)^4]
     end;
end;
EndProgram;
Data;
1 4837.42
2 3236.92
3 2162.17
4 1442.24
5 960.533

十分感谢!

1stopt非线性拟合求助
1.jpg
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lipenggg

专家顾问 (正式写手)

【答案】应助回帖

引用回帖:
3楼: Originally posted by lipenggg at 2016-04-15 09:09:46
Title "test";
Parameters a,b,c,d,e;
Variable x,y;
StartProgram;
var i:integer;
begin
     for i:= 0 to DataLength -1 do begin
     y:=a*2669/(1.5043^x)*(1+b*(a/(1.5043)^x)+c*(a/(1.50 ...

公式见图片
1stopt非线性拟合求助-1
1111111111111111111111111.png

» 本帖已获得的红花(最新10朵)

一团和气
4楼2016-04-15 09:18:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
danfreda: 金币+10, 非常谢谢,可以运行! 2016-04-15 10:21:41
改成如下试试:

Title "test";
Parameters a,b,c,d,e;
Variable i,y;
Function
     y=a*2669/(1.5043^i)*(1+b*(a/(1.5043)^i)+c*(a/(1.5043)^i)^2+d*(a/(1.5043)^i)^3+e*(a/(1.5043)^i)^4);
Data;
1 4837.42
2 3236.92
3 2162.17
4 1442.24
5 960.533

[ 发自手机版 http://muchong.com/3g ]
2楼2016-04-14 23:30:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lipenggg

专家顾问 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
Title "test";
Parameters a,b,c,d,e;
Variable x,y[Output];
StartProgram;
var i:integer;
begin
     for i:= 0 to DataLength -1 do begin
     y:=a*2669/(1.5043^x)*(1+b*(a/(1.5043)^x)+c*(a/(1.5043)^x)^2+d*(a/(1.5043)^x)^3+e*(a/(1.5043)^x)^4)
     end;
end;
EndProgram;
Data;
1 4837.42
2 3236.92
3 2162.17
4 1442.24
5 960.533
一团和气
3楼2016-04-15 09:09:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

danfreda

铜虫 (小有名气)

送红花一朵
引用回帖:
4楼: Originally posted by lipenggg at 2016-04-15 09:18:08
公式见图片

1111111111111111111111111.png
...

非常感谢,不知道是不是我的版本太低(1.5)程序没能运行,但是这个套路我是模仿说明书里写的啊,刚开始用好多还不懂,谢谢谢谢!
5楼2016-04-15 10:24:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +4 a089 2026-08-26 5/250 2026-08-27 09:24 by Howard28
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +4 苏东坡二世 2026-08-23 4/200 2026-08-27 09:18 by 冥王篇
[基金申请] 基金不中,共勉 +7 eulota 2026-08-26 7/350 2026-08-27 09:02 by hsding2026
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 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 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 怎么看青基中了没有啊 +4 叶九微 2026-08-26 4/200 2026-08-26 10:52 by xiacongxu
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见