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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » 1stOpt中关于含有积分式子的参数拟合
361/1返回列表
查看: 698  |  回复: 35
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

savoki大王

银虫 (正式写手)
已领完
1stOpt中关于含有积分式子的参数拟合
领取红包 (小木虫手机app专属红包)

扫一扫，下载小木虫客户端

大神，我现在在参数拟合上遇到了一点问题，我是第一次用1stOpt，所以还请大神耐心解答，这对我很重要~~
我的方程是一个定积分里嵌套一个含有定积分的式子，显示我的方程不对，可是我也没看出来哪里有问题
Title "xue";
Parameters kd,kr;
Variables t,C;
Constants C0=100,u=0.15,D=0.04;z=20,a=37.5,b=2.6667,d=50;
Function C=C0*(2*exp(a))/(sqrt(pi))*(int(exp(-y*y-(a*a)/(4*y*y))*exp(-kr*t-((a*a*b*b)*(kd-kr))/(4*y*y))*(bessI(0,2*(sqrt((a*a*b*b*kr*kd)/(4*y*y)))*sqrt(t-(d*d)/(y*y)))+(kr*exp(kr*(t-(d*d)/(y.*y))))*(int(exp(-kr*x)*bessI(0,2*(sqrt((a.*a.*b*b*kr*kd)/(4*y*y)))*sqrt(x)),x=0, t))),y=z/(2*sqrt(D*t)),100000));
Data;
0        0
82.9015544        0.002855028
155.4404145        0.002802157
238.3419689        0.002741733
310.880829        0.002688862
383.4196891        0.005551443
445.5958549        0.008421577
497.4093264        0.02004562
559.5854922        0.03166211
621.761658        0.049109503
694.3005181        0.069464796
766.8393782        0.095650992
839.3782383        0.127668092
901.5544041        0.162608196
974.0932642        0.203371652
1046.632124        0.24705056
1119.170984        0.293644919
1181.34715        0.343162283
1253.88601        0.392672095
1326.42487        0.445097358
1398.963731        0.494607169
1461.139896        0.544124534
1533.678756        0.590718893
1606.217617        0.634397801
1678.756477        0.678076708
1740.932642        0.715932265
1813.471503        0.753780269
1886.010363        0.785797369
1958.549223        0.814899017
2020.725389        0.841092766
2093.264249        0.86436351
2165.803109        0.884718802
2238.341969        0.902158643
2300.518135        0.919606036
2373.056995        0.931214973
2445.595855        0.94282391
2518.134715        0.954432846
2580.310881        0.960218432
2652.849741        0.968911917
2725.388601        0.97468995
2797.927461        0.980467983
2839.378238        0.983353223
2901.554404        0.986223357
2963.73057        0.989093491
3025.906736        0.989048173
3088.082902        0.991918307
3139.896373        0.994795994
3202.072539        0.994750676
3264.248705        0.99762081
3347.150259        0.997560386

发自小木虫IOS客户端
回复此楼

» 猜你喜欢

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2018-03-25 17:37:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
独孤神宇: 金币+1 2018-03-31 11:58:25
几点错误：
1：”Constants C0=100,u=0.15,D=0.04;z=20,a=37.5,b=2.6667,d=50;“应改成“Constant C0=100,u=0.15,D=0.04,z=20,a=37.5,b=2.6667,d=50;”
2：1stOpt不区分大小写，所以应该把“d”全改成“d1”，否则会和“D”混淆从而出错；
3：积分区间下限“y=z/(2*sqrt(D*t))”，因此t不能为0，否则出现被0除的错误，因此数据第一行应该去除。

改正后可以运行但效果不好，楼主再好好检查下公式是否正确。
赞一下 回复此楼
30楼2018-03-26 17:09:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Breeze灬

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
30楼: Originally posted by dingd at 2018-03-26 17:09:48
几点错误：
1：”Constants C0=100,u=0.15,D=0.04;z=20,a=37.5,b=2.6667,d=50;“应改成“Constant C0=100,u=0.15,D=0.04,z=20,a=37.5,b=2.6667,d=50;”
2：1stOpt不区分大小写，所以应该把“d”全改成“d1”，否则 ...

公式里面有贝塞尔函数，软件手册里说里面只能是实数，可是代码里是含有虚拟变量而且还要积分的贝塞尔函数，这个怎么解决呢
赞一下 回复此楼
33楼2018-03-27 13:16:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
盐津梅李2楼
2018-03-25 17:38   回复  
savoki大王(金币+2): 谢谢参与
已获得2个金币 发自小木虫IOS客户端
Jessie666w3楼
2018-03-25 18:06   回复  
savoki大王(金币+1): 谢谢参与
已获得1个金币 发自小木虫IOS客户端
vv01474楼
2018-03-25 18:14   回复  
savoki大王(金币+1): 谢谢参与
已获得1个金币 发自小木虫Android客户端
風VS雲5楼
2018-03-25 18:15   回复  
savoki大王(金币+2): 谢谢参与
已获得2个金币 发自小木虫Android客户端
11931200306楼
2018-03-25 18:21   回复  
savoki大王(金币+2): 谢谢参与
已获得2个金币 发自小木虫Android客户端
KDME7楼
2018-03-25 18:25   回复  
savoki大王(金币+2): 谢谢参与
已获得2个金币 发自小木虫Android客户端
蒋鲁工程师8楼
2018-03-25 18:28   回复  
savoki大王(金币+2): 谢谢参与
已获得2个金币 发自小木虫Android客户端
hbabc9楼
2018-03-25 18:38   回复  
savoki大王(金币+1): 谢谢参与
已获得1个金币 发自小木虫Android客户端
cy042710楼
2018-03-25 19:10   回复  
savoki大王(金币+2): 谢谢参与
已获得2个金币 发自小木虫IOS客户端
独孤神宇11楼
2018-03-25 19:30   回复  
savoki大王(金币+2): 谢谢参与
已获得2个金币 发自小木虫Android客户端
zyren12312楼
2018-03-25 20:54   回复  
savoki大王(金币+1): 谢谢参与
已获得1个金币 发自小木虫Android客户端
墨雨潇潇13楼
2018-03-25 21:09   回复  
savoki大王(金币+1): 谢谢参与
已获得1个金币 发自小木虫Android客户端
润色论文14楼
2018-03-25 21:48   回复  
savoki大王(金币+2): 谢谢参与
已获得2个金币 发自小木虫IOS客户端
yuekong15楼
2018-03-25 22:05   回复  
savoki大王(金币+1): 谢谢参与
已获得1个金币 发自小木虫Android客户端
钢铁年代16楼
2018-03-25 22:26   回复  
savoki大王(金币+2): 谢谢参与
已获得2个金币 发自小木虫Android客户端
化一而万17楼
2018-03-25 23:12   回复  
savoki大王(金币+1): 谢谢参与
已获得1个金币 发自小木虫Android客户端
czw201718楼
2018-03-25 23:30   回复  
savoki大王(金币+2): 谢谢参与
已获得2个金币 发自小木虫Android客户端
哈了沃1919楼
2018-03-25 23:41   回复  
savoki大王(金币+1): 谢谢参与
已获得1个金币 发自小木虫Android客户端
低调滴分割线20楼
2018-03-26 06:22   回复  
savoki大王(金币+1): 谢谢参与
已获得1个金币 发自小木虫Android客户端
bashan21楼
2018-03-26 06:32   回复  
savoki大王(金币+2): 谢谢参与
已获得2个金币 发自小木虫Android客户端
xhmaohan22楼
2018-03-26 07:08   回复  
savoki大王(金币+2): 谢谢参与
已获得2个金币 发自小木虫Android客户端
gongzilu23楼
2018-03-26 07:17   回复  
savoki大王(金币+2): 谢谢参与
已获得2个金币 发自小木虫Android客户端
youngen24楼
2018-03-26 08:39   回复  
savoki大王(金币+1): 谢谢参与
已获得1个金币 发自小木虫Android客户端
M塞北25楼
2018-03-26 09:44   回复  
savoki大王(金币+2): 谢谢参与
已获得2个金币 发自小木虫Android客户端
劣漳坎26楼
2018-03-26 09:49   回复  
savoki大王(金币+1): 谢谢参与
已获得1个金币 发自小木虫Android客户端
虾登琳27楼
2018-03-26 09:59   回复  
savoki大王(金币+2): 谢谢参与
已获得2个金币 发自小木虫Android客户端
joel_scp28楼
2018-03-26 10:10   回复  
savoki大王(金币+1): 谢谢参与
已获得1个金币 发自小木虫Android客户端
Mr__Right29楼
2018-03-26 10:13   回复  
savoki大王(金币+1): 谢谢参与
已获得1个金币 发自小木虫Android客户端
wdxmu31楼
2018-03-26 19:10   回复  
savoki大王(金币+1): 谢谢参与
已获得1个金币 发自小木虫IOS客户端
农加人32楼
2018-03-27 11:09   回复  
savoki大王(金币+2): 谢谢参与
已获得2个金币 发自小木虫Android客户端
库克塔塔34楼
2018-03-27 18:36   回复  
savoki大王(金币+1): 谢谢参与
已获得1个金币 发自小木虫Android客户端
『ljx』35楼
2018-03-27 19:30   回复  
savoki大王(金币+2): 谢谢参与
已获得2个金币 发自小木虫Android客户端
美人老夏36楼
2018-03-27 22:32   回复  
savoki大王(金币+1): 谢谢参与
已获得1个金币 发自小木虫IOS客户端
相关版块跳转 我要订阅楼主 savoki大王 的主题更新
361/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭