24小时热门版块排行榜    

查看: 1559  |  回复: 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

专家顾问 (正式写手)

引用回帖:
8楼: Originally posted by lipenggg at 2017-02-18 08:08:32
有lib文件夹的1.5版本可以计算共享参数拟合和编程代码,没有的不能。
...

1.5版本有两种,一种是6.3M的,没有lib文件夹。一种是10M,有lib文件夹。

[ 发自手机版 http://muchong.com/3g ]
一团和气
9楼2017-02-18 08:11:16
已阅   回复此楼   关注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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金不中,共勉 +4 eulota 2026-08-26 4/200 2026-08-26 21:29 by 木风9012
[基金申请] 国自然面上复盘~欢迎讨论 +4 晴天加油 2026-08-26 5/250 2026-08-26 21:27 by 木风9012
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 只有每年这种时候来逛逛小木虫 +26 yaoyewhu2008 2026-08-20 28/1400 2026-08-26 13:49 by zzuzxg
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 建议基金发布提前给出明确的时间点 +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
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
信息提示
请填处理意见