24小时热门版块排行榜

返回列表

【悬赏金币】回答本帖问题，作者a1993b522将赠送您 10 个金币

a1993b522

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 5.4
散金: 104
帖子: 124
在线: 192.7小时
虫号: 2414043
注册: 2013-04-12
性别: GG
专业: 工程热物理相关交叉领域

[求助] 1stOpt中非线性拟合怎样限制参数的范围已有1人参与

下载了5.5试用版，按照操作说明编写了一个小程序用于曲线拟合，其中需要求解a,b,c三个参数，但c的范围为[20,30],采用下列程序计算后，发现最终结果没有c的输出（不限制c的范围时，可以输出a，b，c三个参数的值），不知道是不是写错了，另外就是在计算过程中需要确保0<a*1e5*x+b<1，但貌似没有任何作用，求大神给个提示，网上实在是找不到相关资料。
Algorithm=UGO[100];
Variables x,y;
Parameters a,b,c[20,30];
Constant h=6.6262e-34,v=3e8,kb=1.381e-23;
Function y=(a*1e5*x+b)*2*pi*h*v^2/x^5/exp(h*v/kb/x/100/c);
0<a*1e5*x+b<1;
另外就是在操作的过程中，当Function中的函数“y=(a*1e5*x+b)*2*pi*h*v^2/x^5/exp(h*v/kb/x/100/c)”写成“y=(a*1e5*x+b)*2*pi*h*v^2/x^5/exp(h*v/kb/x/c)，程序运行之后就提示Function出错，不知道为什么。

回复此楼

» 猜你喜欢

新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有89人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复
求助NH4V4O10晶体的CIF文件已经有0人回复
英国全奖博士招聘-深度学习与量子物理已经有0人回复
间接带隙半导体有效质量求助已经有0人回复
投稿chemical physical letters不送审？已经有2人回复

1楼 2020-10-19 08:38:36

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

独孤神宇

版主 (知名作家)

应助: 490 (硕士)
贵宾: 0.008
金币: 31014.8
散金: 802
红花: 122
沙发: 1
帖子: 5600
在线: 855.5小时
虫号: 3522474
注册: 2014-11-06
性别: GG
专业: 机械动力学
管辖: 计算模拟

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
a1993b522: 金币+10, ★★★很有帮助, 谢谢，今天研究了一下，发现是由于0<a*1e5*x+b<1的原因，删除了这个约束条件就能够正常运行 2020-10-19 23:18:10

Variables x,y;
Parameters a,b,c=[20,30];
Constant h=6.6262e-34,v=3e8,kb=1.381e-23;
Function y=(a*1e5*x+b)*2*pi*h*v^2/x^5/exp(h*v/kb/x/100/c);
0<a*1e5*x+b<1;
data;
//x,y;