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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 1stopt非线性拟合求助
51/1返回列表
查看: 1294  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2016-04-14 22:14:21
已阅   回复此楼   关注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
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

danfreda
一团和气
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 ]
赞一下(1人) 回复此楼
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的回帖
查看全部 10 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭