- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
求解涉及 digamma 函数的方程组的最有效方法是什么?
我有一个向量 v,我想求解一个向量 w,使得对于所有 i:
digamma(sum(w)) - digamma(w_i) = v_i
和
w_i > 0
我找到了 gsl 函数 gsl_sf_psi,它是 digamma 函数(使用某种级数计算)。是否有一个恒等式可以用来简化方程?我最好的选择是使用求解器吗?我正在使用C++0x;哪种求解器最容易使用且速度最快?
<小时/>根据我的初步研究,digamma 不容易可逆(搜索逆 digamma 给出了通过二分搜索工作的算法),因此整个系统没有简化是有道理的。
因此,使用求解器现在会留下两个问题:处理 digamma 计算速度非常慢的事实,以及处理 w_i > 0 的限制,否则 digamma(w_i) 将在 w_i = 0 时崩溃。
对于第一个问题,我想也许我应该为最近计算的 digamma 值实现一个缓存 - 我认为这是一个好主意,但不太了解寻根算法的工作原理。
我的想法是解决第二个问题是找到w'_i = log(w_i)。这样,w'_i就在整条线上。我想知道这是否是一个好主意。可能没有直接找到 digamma(exp(w')) 的函数?此外,该算法可能在 w' 空间中采取步骤,但不会改进事情,因为从 w'->w 的映射会失去一些精度,因此 w' 的两个元素可能映射到相同的 w。
仍然存在找到一个好的、快速的寻根算法的问题。我想我可以在一个单独的问题中问这个问题。
谢谢...
最佳答案
我建议使用求解器是最好的主意,主要是因为考虑各种方程的各种稳定性和收敛区域可能很棘手,而且重新发明轮子是没有用的。虽然我从未真正解决过像您提到的那样的系统,但我认为以下库之一很可能有您想要的解决方案:
此外,如果这些都没有完全满足您的要求,您可以了解 GNU Octave 或类似的东西如何解决该系统,然后阅读他们关于用于实现解决该系统所需函数的算法的文档。从这里开始,更多的是弄清楚如何使用算法,如何实现它,以及在哪些情况下值得这样做(Octave、Matlab、Mathematica 的文档非常全面,并列出了在大多数情况下定义了该算法的出版物) ,如果您正在寻找开源/免费替代品,还有 Scilab 和 SageMath,并且有多种方法可以在 C++ 中使用这些例程(但我不确定这有多容易或多困难)
希望有帮助。
关于math - 求解包含二 Gamma 函数的方程组的最有效方法是什么?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/2615924/
我有一个类似于以下的结构。 class A { string title; List bItem; } class B { int pric
本地流 和 远程流 两者都是“媒体流列表 ”。 本地流 包含“本地媒体流 ” 对象 但是,远程流 包含“媒体流 ” 对象 为什么差别这么大? 当我使用“本地流 “- 这个对我有用: localVide
我正在尝试将 8 列虚拟变量转换为 8 级排名的一列。 我试图用这个公式来做到这一点: =IF(OR(A1="1");"1";IF(OR(B1="1");"2";IF(OR(C1="1");"3";I
我正在使用面向对象编程在 Python 中创建一个有点复杂的棋盘游戏的实现。 我的问题是,许多这些对象应该能够与其他对象交互,即使它们不包含在其中。 例如Game是一个对象,其中包含PointTrac
有没有办法获取与 contains 语句匹配的最深元素? 基本上,如果我有嵌套的 div,我想要最后一个元素而不是父元素: Needle $("div:contains('Needle')")
出于某种原因,我无法在 Google 上找到答案!但是使用 SQL contains 函数我怎么能告诉它从字符串的开头开始,即我正在寻找等同于的全文 喜欢 'some_term%'。 我知道我可以使用
我正在尝试创建一个正则表达式来匹配具有 3 个或更多元音的字符串。 我试过这个: [aeiou]{3,} 但它仅在元音按顺序排列时才有效。有什么建议吗? 例如: 塞缪尔 -> 有效 琼 -> 无效 S
嘿所以我遇到了这样的情况,我从数据库中拉回一个客户,并通过包含的方式包含所有案例研究 return (from c in db.Clients.Include("CaseStudies")
如果关键字是子字符串,我无法弄清楚为什么这个函数不返回结果。 const string = 'cake'; const substring = 'cak'; console.log(string.in
我正在尝试将包含特定文本字符串的任何元素更改为红色。在我的示例中,我可以将子元素变为蓝色,但是我编写“替换我”行的方式有些不正确;红色不会发生变化。我注意到“contains”方法通常写为 :cont
我想问一下我是否可以要求/包含一个语法错误的文件,如果不能,则require/include返回一个值,这样我就知道所需/包含的文件存在语法错误并且不能被要求/包含? file.php语法错误 inc
我想为所有包含youtube链接的链接添加一个rel。 这就是我正在使用的东西-但它没有用。有任何想法吗? $('a [href:contains(“youtube.com”)]')。attr('re
我正在尝试在 Elasticsearch 中查询。除搜索中出现“/”外,此功能均正常运行。查询如下所示 GET styling_rules/product_line_filters/_search {
我正在开发名为eBookRepository的ASP.NET MVC应用程序,其中包含在线图书。 电子书具有自己的标题,作者等。因此,现在我正在尝试实现搜索机制。我必须使用Elasticsearch作
我已阅读Firebase Documentation并且不明白什么是 .contains()。 以下是文档中 Firebase 数据库的示例规则: { "rules": { "rooms"
我的问题是我可以给出条件[ 'BookTitleMaster.id' => $xtitid, ] 如下所示 $bbookinfs = $this->BookStockin->BookIssue->fi
我需要能够使用 | 检查模式在他们中。例如,对于像“dtest|test”这样的字符串,像 d*|*t 这样的表达式应该返回 true。 我不是正则表达式英雄,所以我只是尝试了一些事情,例如: Reg
我想创建一个正则表达式来不匹配某些单词... 我的字符:var test = "é123rr;and;ià456;or;456543" 我的正则表达式:test.match(\((?!and)(?!o
我在 XSLT 中有一个名为 variable_name 的变量,如果相关产品具有名称为 A 或 B 或两者均为 A & 的属性,我将尝试将其设置为 1 B.
您好,我想让接待员和经理能够查看工作类型和费率并随后进行更新。但是技术人员只能查看不能更新。该图是否有效? 我读到扩展用例是由发起基本用例的参与者发起的。我应该如何区分技术人员只能启动基本案例而不能启
我是一名优秀的程序员,十分优秀!