v7 中 GroebnerBasis 的奇怪行为

2023-12-21

我在使用时遇到了一些奇怪的行为GroebnerBasis. In m1下面,我使用希腊字母作为变量,并在m2,我用了一个拉丁字母。他们俩都没有与之相关的规则。为什么根据我选择的变量,我会得到截然不同的答案?

Image:

可复制代码:

Clear["Global`*"]
g = Module[{x}, 
    x /. Solve[
      z - x (1 - b - 
           b x ( (a (3 - 2 a (1 + x)))/(1 - 3 a x + 2 a^2 x^2))) == 0,
       x]][[3]];
m1 = First@GroebnerBasis[\[Kappa] - g, z]
m2 = First@GroebnerBasis[k - g, z]

EDIT:

正如贝利萨留所指出的,我的用法GroebnerBasis并不完全正确,因为它需要多项式输入,而我的则不需要。这个错误是由复制面食引入的,直到现在才被注意到,因为当我使用以下代码完成其余代码时,我得到了预期的答案m1从上面。但是,我并不完全相信这是一种不合理的用法。考虑下面的例子:

x = (-b+Sqrt[b^2-4 a c])/2a;
p = First@GroebnerBasis[k - x,{a,b,c}]; (*get relation or cover for Riemann surface*)
q = First@GroebnerBasis[{D[p,k] == 0, p == 0},{a,b,c},k,
    MonomialOrder -> EliminationOrder]; 

Solve[q==0, b] (*get condition on b for double root or branch point*) 

{{b -> -2 Sqrt[a] Sqrt[c]}, {b -> 2 Sqrt[a] Sqrt[c]}}

哪个是对的。所以我的解释是OK to use GroebnerBasis在这种情况下,但我不太熟悉其背后的深层理论,所以我在这里可能完全错误。

P.S. I heard that if you mention GroebnerBasis three times in your post, Daniel Lichtblau will answer your question :)


这些示例所显示的错误将在版本 9 中修复。我不知道如何在版本 8 及更早版本中规避它。如果我没记错的话,这是由某些代码中的中间数字溢出引起的,该代码正在检查符号多项式系数是否可能为零。

出于某些目的,可能适合指定更多变量以及可能的非默认术语顺序。此外,至少在有效的情况下,清除分母也会有所帮助。也就是说,我不知道这些策略在这个例子中是否有帮助。

我会更多地查看这段代码,但可能不会在不久的将来。

丹尼尔·利希布劳

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

v7 中 GroebnerBasis 的奇怪行为 的相关文章

  • 在 Mathematica 中写入文件时格式化数字

    这是一个延续this https stackoverflow com q 4943706 181759关于数字格式的问题 与我之前的相关question https stackoverflow com q 1687705 181759关于获
  • Mathematica 中的条件数据操作

    我正在努力准备高效数据分析的最佳工具在数学中 我有大约 300 列和 100 000 行 最好的技巧是什么 删除 提取 或简单地 考虑 数据结构的部分 用于绘制例如 我能想到的最棘手的例子之一是 给定一个数据结构 对于第 2 列中的值等于
  • 从稀疏定义列表中挑选无模式下值的算法

    我有以下问题 我正在开发一个随机模拟器 它随机采样系统的配置 并存储每个配置在特定时间实例被访问次数的统计数据 代码大致是这样的 f Integer Integer 0 someplace later in the code e g ind
  • 了解内核-前端通信——为什么我的前端冻结?

    EDIT 只需确认您是否可以重现此内容就会很有用 只需一台计算机即可尝试此操作 无需远程连接 Update似乎其他人无法在 Mac 或 Win7 上重现此问题 因此它要么是 WinXP 特定的 要么是我的机器特定的 此时我要放弃了 最好有一
  • 对函数体评估感到困惑

    我对以下行为感到困惑Function In 1 InlineCellInMessage Function expr DisplayForm Cell BoxData MakeBoxes expr StandardForm Input Hol
  • Mathematica 下调 Lhs 值

    有谁知道 Mathematica 中是否有一个内置函数可以获取 downvalue 规则的 lhs 没有任何保留 我知道如何编写代码来做到这一点 但这对于内置程序来说似乎已经足够基本了 例如 a 1 2 a 2 3 BuiltInIDoNo
  • 如何用矩阵的一些不连续的行和列形成子矩阵

    我有一个 10 x 10 矩阵 我想从这个主矩阵形成一个子矩阵 使用除第一 第二和第八列和行之外的所有行和列 我知道 Part 可以用来形成子矩阵 但这些示例主要是关于仅使用连续的行和列形成子矩阵 如果这是你的矩阵 tst RandomIn
  • 将符号暴露给 $ContextPath

    有多种Internal 有用的上下文函数 例如InheritedBlock https stackoverflow com questions 4198961 what is in your mathematica tool bag 568
  • 在 Mathematica 中查找先前定义的消息

    Mathematica 默认定义了许多有用的消息来表示常见错误 例如使用错误数量的参数调用函数或未找到文件 一般来说 我更喜欢尽可能使用现有的 已定义的消息 因为这样可以更轻松地通过诸如Check Quiet and On Off 然而 我
  • 解决电力塔

    a 2 Power 10 6 10 9 3 Power 4 9 7 5 TwoTower n Nest 2 1 n 最小的是什么n这样TwoTower n gt a This question http www quora com How
  • 调整 Mathematica 中现有 Graphics 对象的样式/属性

    Mathematica 的优势之一是其一致的对象底层表示 因此 要更改绘图的属性而不重做用于生成它的计算 我可以做类似的事情 Replace myplot Graphics x List y List gt Graphics x Flatt
  • NIntegrate - 为什么在这种情况下 Mathematica 8 的速度要慢得多?

    我有一个 Mathematica 代码 我必须在数值上评估数千个与此类似的积分 NIntegrate Pi Cos Pi 2 x y 1 y 1 y Sin 2 Pi x 1 y Sin Pi 2 x y 1 y E x 1 y x 0 1
  • Mathematica 中的输入/输出单元标签是由什么生成的以及如何向它们添加自动计时?

    当 Mathematica 计算单元格时 它会为输入单元格和输出单元格提供细胞标签 http reference wolfram com mathematica ref CellLabel htmls In Line and Out Lin
  • R 中从右到左的运算符结合性可能吗?

    我是 R 新手 我刚刚发现我患有支架恐惧症 https mathematica stackexchange com a 17315 2266 请参阅链接中的评论 我喜欢这种方式magrittr符号 gt 有效 因为它在某些情况下避免了嵌套括
  • 根据 Mathematica 中的另一个列表值拆分列表

    在 Mathematica 中我有一个点坐标列表 size 50 points Table RandomInteger 0 size RandomInteger 0 size i 1 n 以及这些点所属的聚类索引列表 clusterIndi
  • 零件和跨度:有什么理由这“不应该”起作用吗?

    这是一个设计问题 而不是现有功能的问题 我想使用 1 2 3 4 5 1 3 2 5 我预计 1 2 3 2 3 4 5 但它是无效的 During evaluation of In 1 Part pspec Part specificat
  • 图像分析-光纤识别

    我是图像分析新手 您知道如何以仅获取纤维的方式对该图像进行二值化吗 我尝试过不同的阈值技术等 但没有成功 我不介意应该使用什么工具 但我更喜欢 NET or Matlab PS 我不知道该把答案放在哪里 所以我把它放在StackOverfl
  • Mathematica 中的数据缓存

    有一个非常耗时的操作在我的包中生成数据集 我想保存此数据集 并仅在我手动删除缓存文件时才让包重建它 这是我的方法作为包的一部分 myDataset Module fname data fname cached data mx If File
  • Mathematica 二维热方程动画

    我正在努力绘制二维温度梯度图 但遇到了很多麻烦 我当前的方法是定义一个插值函数 然后尝试多次绘制它的图表 然后为该图表表设置动画 这是我到目前为止所拥有的 RT 388 726919 R 1 FUNC NDSolve D T x y t t
  • Mathematica 圆柱分解的计算复杂度是多少

    数学 圆柱分解 http reference wolfram com mathematica ref CylindricalDecomposition html实现一种称为圆柱代数分解的算法 Wolfram MathWorld 的文章圆柱代

随机推荐

  • 奇怪的重复模板 - 变化

    关于CRP http en wikipedia org wiki Curiously recurring template pattern如果我想实现它的轻微变化 使用模板模板参数 我会收到编译错误 template
  • 自动从我的服务器上传文件到S3的方法

    我有一个项目要完成 但我陷入了困境 我们能否以一定的时间间隔自动将文件从服务器上传到 S3 中 例如 应上传一部分数据 然后每隔一小时上传下一组数据 请按照以下步骤操作 创建一个脚本文件 例如将其命名为 script sh 将其添加到 sc
  • 检测移动设备上隐藏/不可见的蓝牙设备

    我有目标设备的蓝牙 mac id 然而 蓝牙处于隐藏 不可见模式 两款手机都是安卓手机 非根 我做了我的研究 我 知道 如果我知道 mac id 我可以连接到设备 在大多数情况下 但是我只想扫描邻近度而不是连接 即如果目标设备在蓝牙范围内并
  • 了解 VTune 报告

    这是现有线程的后续内容 http stackoverflow com questions 12724887 caching in a high performance financial application 我发现这不是阻碍我的应用程序
  • 使用node js返回所有带有过滤文件的文件

    我想返回文件夹和子文件夹中存在的文件 但需要对扩展名以 html htm 或 aspx 结尾的文件进行过滤 我有一个代码只返回带有扩展名的文件Index html Default htm Index aspx也需要文件的其余部分 但不知道如
  • 如何从 JSPX 输出

    一些 html 标签将给定属性的 任何 值解释为 true gt 选项标签浮现在脑海中 我经常最终做这样的事情
  • 在对象中实现打包/拆包

    我有一个仅包含属性的类 我希望对其进行打包 解包 什么collections abc我应该实现这种行为吗 class Item object def init self name age gender self name name self
  • ObjCMongoDB 坏了?我无法让它工作

    全面披露 我不是一个受过教育的程序员 我的全部编程经验都在 Javascript 和 Objective C 中 所以现在你知道你在处理什么了 小心行事 没有突然的动作 我正在编写一个与 Mac 服务器通信的 iPad 应用程序 我想使用
  • 根据所选的自定义结帐字段值添加自定义电子邮件收件人

    我需要 Woocommerce 根据为 Field Checkout 选择的选项向不同的个人发送自定义电子邮件 从技术上讲 自定义字段是报告他们购买的产品变体的人 但我不确定如何根据购买的产品变体 因此如下 首先 我使用以下代码建立了自定义
  • 从 Gradle 托管设备运行的 Android 模拟器复制数据

    在使用管理的 Android 模拟器上运行 UITests 时Gradle 管理的设备 https developer android com studio test gradle managed devices 如何复制屏幕截图 黄瓜报告
  • 减慢绑定DataGrid的刷新率

    我的 WPF C 应用程序中有一个 DataGrid 它绑定到实体框架集合 每行都有经常更改的绑定列 每秒多次 这会导致该列基本上不可读 因为它经常更改 如何强制 WPF 仅每 0 5 秒或 1 秒显示一个新值 即使该值每 0 1 秒更改一
  • 在 Snow Leopard 中控制 OS X 音量

    这是后续通过 Objective C 控制 Mac OS X 中正在运行的应用程序的音量 https stackoverflow com questions 5109971 controlling volume of running app
  • 使用pyenv安装anaconda,无法配置虚拟环境

    我正在使用 pyenv 来管理 python 安装和虚拟环境 我希望 anaconda 成为这样的安装之一 并且能够使用 anaconda python 创建虚拟环境 使用pyenv install anaconda3 2019 03成功安
  • 如何确定给定类型 (System.Type) 是否继承自特定基类(在 .Net 中)?

    这可能是一个简单的答案 我只是错过了一些东西 但是这里 如果我有一个类型 即 一个实际的 System Type 不是一个实例 我如何判断它继承自另一个特定的基类型 Use the 是子类 http msdn microsoft com e
  • 如何从带有列表的嵌套字典构建多索引 Pandas DataFrame

    我有以下字典 d key1 sub key1 a b c d e key2 sub key2 1 2 3 5 8 9 10 在 的帮助下this https stackoverflow com questions 13575090 cons
  • jQuery - 替换字符串中的所有括号

    我试过这个 mystring mystring replace g replace g replace replace 它适用于所有双引号和单引号 但对于括号 这只替换字符串中的第一个括号 如何使用 JavaScript 替换字符串中的所有
  • 如何正确解析传入的 HTTP 请求

    我使用 WinSck 创建了一个 C 应用程序 它实现了一个小型 仅处理我需要的一些功能 http 服务器 这用于使用 http 请求与外界进行通信 它可以工作 但有时请求无法正确处理 因为解析失败 现在我非常确定请求的格式正确 因为它们是
  • Rails 中无法转储文件错误

    我正在使用 Rails 3 0 3 开发一个应用程序 我在以下代码中收到 无法转储文件 错误 if post update attributes params post redirect to post path post lt ERROR
  • 将 substring() 的 endIndex 减 1 的原因[重复]

    这个问题在这里已经有答案了 可能的重复 为什么 from 索引是包含的 而 end 索引是排除的 https stackoverflow com questions 6339158 why from index is inclusive b
  • v7 中 GroebnerBasis 的奇怪行为

    我在使用时遇到了一些奇怪的行为GroebnerBasis In m1下面 我使用希腊字母作为变量 并在m2 我用了一个拉丁字母 他们俩都没有与之相关的规则 为什么根据我选择的变量 我会得到截然不同的答案 Image 可复制代码 Clear