24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知(长期有效)
查看: 1570  |  回复: 9

superzhuxun

金虫 (小有名气)

[求助] 1stOpt代码有问题,还请大神帮忙看看,很着急,解决了又追加,谢谢大家。已有1人参与

版本1stOpt 1.5,已知函数表达式,B,N为待拟合常量,r为变量,s为函数积分中间量,如上传图111。
自己编写的代码也上传,和上面的图一起打包,请大家了1stOpt代码有问题,还请大神帮忙看看,很着急,解决了又追加,谢谢大家。
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 1111.rar
  • 2015-07-29 17:31:38, 30.44 K

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
该问题不应该是拟合问题吧?一组r值求出一组f(r,B,N),该值必须满足范围[-10^(-6), 10^(-6))]之间,是这个意思吧?
试了一下,无法满足范围[-10^(-6), 10^(-6))],最多可满足[-10^(-2.1069846), 10^(-2.1069846))],此时:

b: 1.68190766417634
n: -0.136717008660948
2楼2015-07-29 22:53:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

superzhuxun

金虫 (小有名气)

引用回帖:
2楼: Originally posted by dingd at 2015-07-29 22:53:56
该问题不应该是拟合问题吧?一组r值求出一组f(r,B,N),该值必须满足范围之间,是这个意思吧?
试了一下,无法满足范围,最多可满足,此时:

b: 1.68190766417634
n: -0.136717008660948...

大哥,我的程序是正确的吗?我再看看,那估计是我自己的解有问题。之后有什么情况我在跟大哥交流啊,哈哈。

[ 发自小木虫客户端 ]
3楼2015-07-29 23:01:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

superzhuxun

金虫 (小有名气)

引用回帖:
2楼: Originally posted by dingd at 2015-07-29 22:53:56
该问题不应该是拟合问题吧?一组r值求出一组f(r,B,N),该值必须满足范围之间,是这个意思吧?
试了一下,无法满足范围,最多可满足,此时:

b: 1.68190766417634
n: -0.136717008660948...

大哥,我在我的1.5的版本里面把那个取值范围改大了,还是现实函数表达式有误。。。。

[ 发自小木虫客户端 ]
4楼2015-07-29 23:05:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

CODE:
Title "B,N";
Parameters B,N ;
Constant r=[4.5,8.5,17,23,30,40,58,74,104,141,194,241];
ConstStr y=((exp(10*N))/(20*B))*int(s*0.015*(ln(1.418*(10^(-6))*s^2)-7.964*(10^(-7))*s^2)*bessi0(s*r/(20*B))*exp(-(r^2+s^2)/(40*B)),s=0,100)+0.0039+5.47*(10^5)*exp(-0.5*((r+1316.87)/230.58)^2);
MinFunction 0;
         For(r)(y<=10^(-2.1069846));
         For(r)(y>=-10^(-2.1069846));

5楼2015-07-29 23:21:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

superzhuxun

金虫 (小有名气)

引用回帖:
5楼: Originally posted by dingd at 2015-07-29 23:21:00
Title "B,N";
Parameters B,N ;
Constant r=;
ConstStr y=((exp(10*N))/(20*B))*int(s*0.015*(ln(1.418*(10^(-6))*s^2)-7.964*(10^(-7))*s^2)*bessi0(s*r/(20*B))*exp(-(r^2+s^2)/(40*B)),s=0,100)+ ...

大哥,我刚才把你给的程序放到我的1.5版本里,但说是常数变量r数值定义有误,你能再帮我看看嘛?麻烦大哥了啊。
1stOpt代码有问题,还请大神帮忙看看,很着急,解决了又追加,谢谢大家。-1
大哥,请你看看.jpg

6楼2015-07-29 23:33:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

引用回帖:
6楼: Originally posted by superzhuxun at 2015-07-29 23:33:02
大哥,我刚才把你给的程序放到我的1.5版本里,但说是常数变量r数值定义有误,你能再帮我看看嘛?麻烦大哥了啊。

大哥,请你看看.jpg
...

版本问题。
7楼2015-07-30 09:17:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

superzhuxun

金虫 (小有名气)

引用回帖:
7楼: Originally posted by dingd at 2015-07-30 09:17:33
版本问题。...

好的,大哥。最后想麻烦一下大哥,这个东西小弟急用,数据老是变,老是请大哥算也不是个事情,小弟学生,买正版的也买不起。所以还请大哥看看,能不能帮小弟把这个程序改成1.5版本的,小弟感激不尽。

[ 发自小木虫客户端 ]
8楼2015-07-30 10:32:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

引用回帖:
8楼: Originally posted by superzhuxun at 2015-07-30 10:32:01
好的,大哥。最后想麻烦一下大哥,这个东西小弟急用,数据老是变,老是请大哥算也不是个事情,小弟学生,买正版的也买不起。所以还请大哥看看,能不能帮小弟把这个程序改成1.5版本的,小弟感激不尽。
...

没用过1.5版,多少年前的版本了?不知道怎么改。
9楼2015-07-30 13:56:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

superzhuxun

金虫 (小有名气)

引用回帖:
9楼: Originally posted by dingd at 2015-07-30 13:56:27
没用过1.5版,多少年前的版本了?不知道怎么改。...

恩恩,还是麻烦大哥啦,哈哈
10楼2015-07-30 14:24:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 superzhuxun 的主题更新
信息提示
请填处理意见