| 查看: 1371 | 回复: 3 | ||
| 【悬赏金币】回答本帖问题，作者a1993b522将赠送您 10 个金币 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
a1993b522银虫 (小有名气)
|
[求助]
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出错，不知道为什么。 |
回复此楼» 猜你喜欢
金属氧化物薄膜空气退火和真空退火对表面粗糙度的影响
已经有0人回复
-
Apex2单晶测试
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有125人回复
-
结构建模
已经有0人回复
-
声子谱计算问题
已经有2人回复
-
已结束
已经有23人回复
-
Error in routine lambda (100): wrong # or too many modes
已经有0人回复
-
祈福国自然面上
已经有55人回复
-
请问castep计算出的能带图横轴数据缩放因子怎么计算
已经有1人回复
-
Result in Physics 从见刊开始 一般多久可以被SCI收录
已经有3人回复
独孤神宇
版主 (知名作家)
- 应助: 489 (硕士)
- 贵宾: 0.008
- 金币: 29683.8
- 散金: 802
- 红花: 121
- 沙发: 1
- 帖子: 5370
- 在线: 850.6小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
4楼2021-05-11 20:05:31
独孤神宇
版主 (知名作家)
- 应助: 489 (硕士)
- 贵宾: 0.008
- 金币: 29683.8
- 散金: 802
- 红花: 121
- 沙发: 1
- 帖子: 5370
- 在线: 850.6小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
a1993b522: 金币+10, ★★★很有帮助, 谢谢，今天研究了一下，发现是由于0<a*1e5*x+b<1的原因，删除了这个约束条件就能够正常运行 2020-10-19 23:18:10
感谢参与，应助指数 +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; |
2楼2020-10-19 21:32:48
洞洞啊
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 110.3
- 散金: 20
- 帖子: 148
- 在线: 10.6小时
- 虫号: 21254898
- 注册: 2020-02-26
- 专业: 化工热力学和基础数据
3楼2021-05-08 15:14:31
10