| 查看: 716 | 回复: 2 | |||
ldoop
|
[交流]
【求助】请大家帮忙看看这个程序【已完结】 已有2人参与
|
|
大家好，请帮忙看看我的这个小程序： 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值是复数，应如何处理？或者避免？ 谢谢 [ Last edited by nono2009 on 2010-9-30 at 08:43 ] |
回复此楼» 猜你喜欢
出来了
已经有9人回复
-
有多少人是今天查系统知道结果的？
已经有14人回复
-
只有每年这种时候来逛逛小木虫
已经有28人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有5人回复
-
为什么国自然不能直接公布
已经有4人回复
-
系统查不到
已经有9人回复
-
科研孤儿太难了
已经有22人回复
-
怎么查啊
已经有4人回复
-
国合里面能看到了
已经有7人回复
-
国际合作可查了，中了面上
已经有18人回复
mshwangg
至尊木虫 (正式写手)
- 程序强帖: 5
- 应助: 206 (大学生)
- 金币: 10702.8
- 散金: 100
- 红花: 19
- 帖子: 597
- 在线: 195.4小时
- 虫号: 576702
- 注册: 2008-06-21
- 专业: 物理学I
2楼2010-09-29 22:16:35
ldoop 
铁杆木虫 (著名写手)
- 应助: 7 (幼儿园)
- 金币: 3485.8
- 散金: 3151
- 红花: 16
- 沙发: 1
- 帖子: 2225
- 在线: 285.5小时
- 虫号: 533053
- 注册: 2008-03-25
- 专业: 群落生态学
3楼2010-09-30 00:27:34
28