24小时热门版块排行榜    

查看: 2432  |  回复: 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出错,不知道为什么。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

洞洞啊

新虫 (小有名气)

你好,5.5下载链接能发一下么

发自小木虫Android客户端
3楼2021-05-08 15:14:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

引用回帖:
3楼: Originally posted by 洞洞啊 at 2021-05-08 15:14:31
你好,5.5下载链接能发一下么

Auto2Fit5.5  ,试用版网上是有资源的,找一下。
数值计算
4楼2021-05-11 20:05:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 a1993b522 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 17:01 by BCdvBAEpt3zm
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 16:52 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 18/900 2026-11-01 16:48 by BCdvBAEpt3zm
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 X5bHxOULxFqb 2026-10-30 22/1100 2026-11-01 16:35 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 X5bHxOULxFqb 2026-10-30 14/700 2026-11-01 16:32 by BCdvBAEpt3zm
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 15:04 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 xvJS1BYDxvBS 2026-10-31 6/300 2026-11-01 15:01 by BCdvBAEpt3zm
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 10/500 2026-11-01 14:50 by BCdvBAEpt3zm
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 14:48 by BCdvBAEpt3zm
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +9 zuesli 2026-10-27 15/750 2026-11-01 12:59 by Equinoxhua
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 3/150 2026-11-01 12:53 by BCdvBAEpt3zm
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 7/350 2026-11-01 12:27 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 17/850 2026-11-01 12:24 by BCdvBAEpt3zm
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 12:14 by BCdvBAEpt3zm
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 12:14 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 X5bHxOULxFqb 2026-10-30 18/900 2026-11-01 12:11 by BCdvBAEpt3zm
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 4/200 2026-11-01 11:01 by BCdvBAEpt3zm
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-01 10:43 by BCdvBAEpt3zm
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 03:59 by YXM1FDCSE8Bn
[考博] 申博求助 20+3 Katrinae 2026-10-26 4/200 2026-10-27 09:56 by 北京莱茵编辑
信息提示
请填处理意见