HMM 前向算法中的下溢

2023-11-24

我正在实现 HMM 的前向算法来计算给定 HMM 发出给定观察序列的概率。我希望我的算法对于下溢具有鲁棒性。我无法在对数空间中工作,因为前向算法需要概率的乘法和加法。避免下溢的最佳方法是什么?

我已经阅读了一些关于此的资料,但我得到的最好的建议是缩放每个时间步的概率第 6 节在这里。在算法结束时,您将不会得到您想要的(观察序列的)确切概率。另外,除非我弄错了,如果按照上述参考文献中的建议在每个时间步缩放概率,则无法对来自两个不同 HMM 的给定观察序列的概率进行有意义的比较(以找出哪个)更有可能输出序列)。有什么建议么?


在参考文献末尾的方程 32 中,将每个概率值 alpha_t(i) 乘以 C_t。所以最后你将最终概率乘以所有 C_t 的乘积。您可以通过跟踪 log(C_t) 的总和来跟踪所有这些。然后最后你可以算出 log(alpha_t(i)) - SUM_(j

本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

HMM 前向算法中的下溢 的相关文章

  • 如何在 dijkstra 算法中以 O(log n ) 时间更新优先级队列中的键?

    过去一周我一直在研究 dijkstra 算法 我在 java 中有正确的运行代码 它使用数组来计算标准 findMin 函数 该函数为您提供距离最小的顶点 显然它是 O n 现在我希望使用优先级队列 最小堆 来实现它 我的思考过程是 whi
  • 以一定角度遍历二维数组

    通常我们按行或列遍历数组 但这里我想以角度遍历它 我会尝试解释我的意思 因此 假设角度是 45 度 那么它会搜索为 0 0 then 0 1 1 0 then 0 2 1 1 2 0 等等 抱歉 无法上传图像 因为我是新用户 不允许这样做
  • 4 x 3 锁图案

    我遇到了这个 它要求计算在 4x3 网格中可以制作特定长度的锁定图案的方式数 并遵循规则 可能有些点不能包含在路径中 有效的模式具有以下属性 图案可以使用第一次接触的点序列来表示 与绘制图案的顺序相同 从 1 1 到 2 2 的图案与图案不
  • Bellman-Ford 算法检测什么?负重还是负循环?

    如果给定一个图 现在我们要从源头计算最短路径 现在 如果一条边具有负权重 但在到达目的地时有边到后边返回到该边 我的意思是如果没有循环 那么我们就没有负循环 但是here http en wikipedia org wiki Bellman
  • 如何设置K-means openCV c++的初始中心

    我正在尝试使用 OpenCv 和 Kmeans 对图像进行分割 我刚刚实现的代码如下 include opencv2 objdetect objdetect hpp include opencv2 highgui highgui hpp i
  • 加密成本高,解密成本低

    我希望该用户 攻击者加密数据并发送给服务器 现在我想要一种与标准算法完全相反的算法 使用快 难以解密 即很难使用服务器发送的密钥来加密密码等数据 以防止随机攻击 但很容易解密这样服务器在验证用户时消耗的时间非常少 但是对于攻击者来说 每次使
  • 反转二进制网络

    如何反转二元方程 以便找到哪些输入将产生给定的输出 Example Inputs i0 through i8 Outputs o0 through o8 Operators XOR AND 二元方程 1 i0 1 i1 0 i2 1 i3
  • 按百分比减少多边形面积

    我有一个由点 x y 组成的多边形 我想做的是将其减少一个百分比 请记住 我不想只是扩大规模 多边形应该有一种内部边界 其宽度取决于百分比 该内部边界被多边形切断 谁知道可以实现这一目标的算法 输入 点数组 百分比 输出 点数组 你所寻求的
  • 如何在给定目标索引数组的情况下对数组进行就地排序?

    你如何对给定的数组进行排序arr in place给定目标索引数组ind 例如 var arr A B C D E F var ind 4 0 5 2 1 3 rearrange arr ind console log arr gt B E
  • 如何高效地在屏幕上精确绘制N个点?

    这听起来是一个简单的问题 但我发现要获得良好的性能是非常棘手的 我提出的第一个算法是随机绘制点 从一组中检查是否已绘制 否则绘制 如果我们只绘制几个点 那么这种方法效果很好 但当我们接近填满屏幕时 速度会灾难性地减慢 我想出的最好的方法是构
  • 使用主方法求解 T(n) = 2T(n/2) + n/log n 和 T(n) = 4T(n/2) + n/log n 之间的差异

    我最近偶然发现了一个资源 其中 2T n 2 n log ntypeMM 宣布复发无法解决 我接受它作为一个引理 直到今天 另一种资源被证明是矛盾的 在某种意义上 根据资源 下面的链接 其中的 Q7 和 Q18 是建议 分别在问题中的1和2
  • 通过分布式数据库聚合作业优化网络带宽

    我有一个分布式 联合数据库 结构如下 数据库分布在三个地理位置 节点 每个节点集群有多个数据库 关系数据库是 PostgreSQL MySQL Oracle 和 MS SQL Server 的混合体 非关系数据库是 MongoDB 或 Ca
  • 如何从一组重叠的圆计算多边形集?

    这个问题是一些计算细节的扩展这个问题 https stackoverflow com questions 1667310 combined area of overlapping circles 假设有一组 可能重叠的 圆 并且希望计算这组
  • 寻找距离原点最近的 100 颗恒星的算法

    首先让我提出正确的问题 问 有一个文件包含超过一百万个点 x y 每个点代表一颗星星 a b 处有一颗行星地球 现在 任务是构建一种算法 返回距离地球最近的 100 颗恒星 您的算法的时间和空间复杂度是多少 这个问题在各种采访中被问过很多次
  • 如何确定算法函数的复杂度?

    您如何知道算法函数对于特定操作是否需要线性 常数 对数时间 它取决于CPU周期吗 您可以通过三种方式 至少 做到这一点 在网上查找算法 看看它是如何描述其时间复杂度的 根据输入大小 自己检查算法 查看嵌套循环和递归条件等内容 以及每个循环运
  • 解释一下从 N 个给定集合中的每一个中给出第 K 个最大数字的示例?

    今天我尝试解决一个Facebook 编程挑战赛 https facebook interviewstreet com recruit challenges 我遇到的挑战是 酒吧问题 可以找到here https github com alo
  • 寻找将集合映射到整数的双射函数

    对于任意两个序列 a b 其中 a a1 a2 an 且 b b1 b2 bn 0a b具有相同的元素 而不关心它们的顺序 例如 如果 a 1 1 2 3 b 2 1 3 1 c 3 2 1 3 则 f a f b f a f b 我知道有
  • 数组中连续元素的最大乘积

    我在现场面试的时候被问到了这个算法问题 由于没有要求我签署保密协议 我将其发布在这里寻求答案 给定一个数组REAL不包含 0 的数字 找到产生最大乘积的连续元素 该算法应在线性时间内运行 我考虑过以下方法 使用两个数组 第一个是利用DP思想
  • 生成所有多集大小为 n 的分区的算法

    我一直在试图找出一种方法来生成多重集的所有不同的大小为 n 的分区 但到目前为止却空手而归 首先让我展示一下我想要实现的目标 假设我们有一个输入向量uint32 t std vector
  • 在常数空间中创建 1..N 的随机排列

    我正在寻找枚举固定空间中数字 1 N 的随机排列 这意味着我无法将所有数字存储在列表中 原因是 N 可能非常大 超过可用内存 我仍然希望能够一次遍历这样一个数字的排列 只访问每个数字一次 我知道对于某些 N 可以这样做 许多随机数生成器随机

随机推荐

  • 我需要导入什么才能访问我的模型?

    我想运行一个脚本来填充我的数据库 我想通过 Django 数据库 API 访问它 唯一的问题是我不知道需要导入什么才能访问它 如何才能实现这一目标 也导入您的设置模块 import os os environ DJANGO SETTINGS
  • jQuery 模式对话框和 jqGrid

    如何将 Jquery 模态确认与 jqGrid 一起使用 假设当我提交条目时 它将弹出一个模式对话框并显示名称以及发送到服务器的消息 我的方法 dialog confirm dialog autoOpen false resizable f
  • 是否有一个类似 PropertyPlaceholderConfigurer 的类可以与 Spring 一起使用并接受 XML?

    Spring 有一个非常方便的类 称为属性占位符配置器 它采用标准 properties 文件并将其中的值注入到您的 bean xml 配置中 有谁知道有一个类可以做完全相同的事情 并以相同的方式与 Spring 集成 但接受配置的 XML
  • 加速 Python 中的 MSD 计算

    这是向社区发出的呼吁 看看是否有人有提高 MSD 计算实现速度的想法 它主要基于这篇博客文章的实现 http damcb com mean square disp html 目前 对于 5000 个点的 2D 轨迹 当前的实现大约需要 9
  • Flutter:StatelessWidget.build 被多次调用[重复]

    这个问题在这里已经有答案了 我总是将代码放在我的问题中 但这一次这是不可能的 因为错误可能出现在一千行代码中的任何地方 然而 我注意到我的主屏幕的构建方法 StatelessWidget 它是 a 的后代MaterialApp home 属
  • JavaScript 按 id 合并对象[重复]

    这个问题在这里已经有答案了 在 Javascript 中合并两个数组的正确方法是什么 我有两个数组 例如 var a1 id 1 name test id 2 name test2 var a2 id 1 count 1 id 2 coun
  • RequestFactory 理论:为什么 Locator<>.find() 被如此频繁地调用?

    我是 RequestFactory 的新手 但得到了慷慨的帮助托马斯 布罗耶在查看了下面的文档之后 情况变得好多了 RequestFactory 入门 请求工厂移动部件 GWT 2 4 中的 RequestFactory 更改 但你能解释一
  • 如何在 Python 中重置 TCP 套接字?

    我有一个用 Python 编写的套接字代理 当它从一对通信对等点接收到 RST 时 它将通过让套接字被垃圾收集来关闭与两个对等点的连接 这会导致其他对等方看到 FIN 而不是 RST 这意味着代理有效地将 RST 转换为 FIN 我认为这并
  • Anaconda“无法创建进程”

    我是 Python 新手 我刚刚安装了 Python anaconda python 2 7 在启动 Anaconda 时提示发生了 无法创建进程 的情况 那么这里有人可以帮我吗 我感谢每一个帮助 我遇到了完全相同的错误 因为我的用户名包含
  • SVN:将存储库主干移动到另一个分支(带有历史记录)

    我正在使用带有大量存储库的 SVN 设置 我试图通过将一个的主干移动到另一个的分支来合并一些 旧的是新的主题版本 减去我稍后将应用的一些代码修复 所以这对我来说很有意义 简短版本 我想从 RepositoryA trunk 转到 Repos
  • VB“Financial.Pmt”在 C# 中等效吗?

    Microsoft VisualBasic 程序集中有一个内置函数 我可以在 VB 中这样使用它 Financial Pmt dAPR 100 12 iNumberOfPayments dLoanAmount 1 我当前的项目是用C 编写的
  • 如何使用 jQuery 计算 ASP.NET 中 gridview 的行数

    有谁知道如何使用 jQuery 计算 asp GridView 中的行数 如果没有找到行那么我想做一些事情 A GridView只是呈现为标准 HTML 表格 因此只需计算trGridView 下的元素 var totalRows tr l
  • 在 ListBoxFor 中选择值的挑战

    最近在开发我的第一个 ASP Net MVC2 Web 应用程序时 当我需要在列表框中选择多个值时 我遇到了一些问题 我用一些 jQuery 解决了这个问题 但继续编写了一些非常简单的代码来演示 我使用 EF 作为模型 有两个实体 Cust
  • 如何从我的 Android 应用程序中删除 QUERY_ALL_PACKAGES 权限?

    由于 Google 的反馈 我的 Google Play 更新版本已被拒绝 3 次 应使用不太广泛的应用程序可见性方法 我们无法批准您的应用使用 QUERY ALL PACKAGES 权限 因为声明的任务可以使用不太广泛的应用可见性方法来完
  • stl容器是否使用隐式共享?

    众所周知 Qt 小部件使用隐式共享 所以我对 stl 容器感兴趣std vector std string也使用隐式共享 如果没有 为什么 因为它非常有用 如果答案是肯定的 我们如何确定呢 我需要简单的 C stl 程序 该程序显示 stl
  • 将 html 结果保存到 txt 或 html 文件

    我有一个带有 html 代码的变量 以下是 R 控制台中代码变量的输出 h1 My First Heading h1 p My first paragraph p 我尝试将内容保存到txt文件中 write table code file
  • 什么是卷积神经网络中的“线性投影”[关闭]

    Closed 这个问题不符合堆栈溢出指南 目前不接受答案 我正在阅读剩余学习 我有一个问题 3 2中提到的 线性投影 是什么 一旦得到这个看起来很简单 但无法理解 有人可以提供简单的例子吗 首先 重要的是要了解什么x y and F以及为什
  • Google BigQuery 的 JDBC 驱动程序?

    有谁知道 Google BigQuery 的 JDBC 接口或驱动程序吗 请只使用 Java 我已经不再使用 Python 库了 有一个第 3 方 JDBC 驱动程序 可以从以下位置获取 http code google com p sta
  • 不完整类型的 new 在包含在模板中时编译

    考虑这段代码 有一个明显的编译错误 1 struct A struct B B new A error allocation of incomplete type A Using a unique ptr也无济于事 2 struct A s
  • HMM 前向算法中的下溢

    我正在实现 HMM 的前向算法来计算给定 HMM 发出给定观察序列的概率 我希望我的算法对于下溢具有鲁棒性 我无法在对数空间中工作 因为前向算法需要概率的乘法和加法 避免下溢的最佳方法是什么 我已经阅读了一些关于此的资料 但我得到的最好的建