24小时热门版块排行榜    

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

专家顾问 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +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的回帖
查看全部 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

专家顾问 (正式写手)

【答案】应助回帖

引用回帖:
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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 C79jjtjAKjEn 2026-09-15 4/200 2026-09-18 06:44 by manzzl
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 Aj1rhIDL5ixY 2026-09-14 6/300 2026-09-18 01:54 by cITaGg3p5edr
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-18 01:40 by cITaGg3p5edr
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-18 01:04 by cITaGg3p5edr
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 s3fFTmArrBt6 2026-09-13 6/300 2026-09-17 23:28 by cuPZTDXS3VOd
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-13 8/400 2026-09-17 23:28 by cuPZTDXS3VOd
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-17 22:39 by cuPZTDXS3VOd
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +8 QUjhNVAcOSff 2026-09-13 12/600 2026-09-17 22:15 by cuPZTDXS3VOd
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 5BDX0d0WFp7t 2026-09-15 4/200 2026-09-17 16:29 by 0ngsZVSCTS7z
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 0pYnUiPDfdnk 2026-09-14 4/200 2026-09-17 14:15 by ObVzqQQh1rrE
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 40ms4Wlo4umh 2026-09-14 5/250 2026-09-17 14:05 by ObVzqQQh1rrE
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-17 13:15 by ObVzqQQh1rrE
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-17 13:05 by ObVzqQQh1rrE
[找工作] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +8 3n8v2C8RimXI 2026-09-13 8/400 2026-09-17 11:27 by UAIjTRMuUFoj
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 23jxep3nCNZb 2026-09-14 4/200 2026-09-17 04:04 by k06BKNblrGRH
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-16 03:14 by DgNGHc3h5tPl
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:26 by C79jjtjAKjEn
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 07:02 by 5BDX0d0WFp7t
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
请填处理意见