24小时热门版块排行榜

>论坛更新日志 (3615)
>虫友互识 (502)
>导师招生 (270)
>考博 (229)
>文献求助 (199)
>基金申请 (163)
>休闲灌水 (110)
>文学芳草园 (50)
>仿真模拟 (45)
>硕博家园 (45)
>博后之家 (41)
>论文投稿 (37)
>考研 (34)
>招聘信息布告栏 (24)
>找工作 (22)
>催化 (20)

返回列表

ldoop

铁杆木虫 (著名写手)

应助: 7 (幼儿园)
金币: 3485.8
散金: 3151
红花: 16
沙发: 1
帖子: 2225
在线: 285.5小时
虫号: 533053
注册: 2008-03-25
专业: 群落生态学

[交流] 【求助】mathematica 求助

大家好，请帮忙看看我的这个小程序：

n = Table[{140, 80, 57, 41, 31, 24, 20, 16, 13, 11, 10, 8, 7, 6, 5, 4,
   4, 4, 3, 3, 2, 2, 2, 1, 1, 1, 1, 1, 1, 1, 1}];
f[j_, r_, x_, a_] :=
  n[[j]] Log[r x] + LogGamma[a + 1] + LogGamma[n[[j]]] +
Log[Hypergeometric2F1[n[[j]], n[[j]], n[[j]] + a, x - r x]] -
Log[x] - Log[n[[j]]] - LogGamma[n[[j]] + a] -
Log[HypergeometricPFQ[{1, 1, 1}, {2, 1 + a},
   x] - (1 - r) HypergeometricPFQ[{1, 1, 1}, {2, 1 + a}, x - r x]];
L[r_, x_, a_] := Sum[N[f[j, r, x, a]], {j, 1, Length[n]}];
NMaximize[{L[x, r, a], 0 < x, x < 1, 0 < r, r < 0.2, a > 0}, {x, r, a}]

为什么计算不出来值呢？
如果限制条件是r<1，可以得到值，再改成其他的就得到SystemException["MemoryAllocationFailure", 请大家帮忙看看原因？

另外为什么当n中含有29，30，会产生复数，也就是n=29 or 30 计算的L值是复数，应如何处理？或者避免？
谢谢

回复此楼

» 猜你喜欢

1楼 2010-09-29 19:36:29

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

相关版块跳转我要订阅楼主 ldoop 的主题更新

返回列表