AskOverflow.Dev

AskOverflow.Dev Logo AskOverflow.Dev Logo

AskOverflow.Dev Navigation

  • 主页
  • 系统&网络
  • Ubuntu
  • Unix
  • DBA
  • Computer
  • Coding
  • LangChain

Mobile menu

Close
  • 主页
  • 系统&网络
    • 最新
    • 热门
    • 标签
  • Ubuntu
    • 最新
    • 热门
    • 标签
  • Unix
    • 最新
    • 标签
  • DBA
    • 最新
    • 标签
  • Computer
    • 最新
    • 标签
  • Coding
    • 最新
    • 标签
主页 / coding / 问题 / 79160333
Accepted
zajer
zajer
Asked: 2024-11-06 02:37:45 +0800 CST2024-11-06 02:37:45 +0800 CST 2024-11-06 02:37:45 +0800 CST

当存在多个解时,Matlab mldivide 返回 NaN

  • 772

我有两个矩阵:

A = [ -1 0  0;
       1 1 -1;
       0 -1 1 ];
B = [-1; 0; 1];

我想解下列方程:

Ax=B

当我使用mldivide函数时,我得到一个 NaN 矩阵

X = mldivide(A,B)
X =

   NaN
   NaN
   NaN

知道该问题有多种解决方案,我手动测试了其中一种,即[1;0;1]返回B:

A*[1; 0; 1]

并且,正如预期的那样,我确信它是多个解决方案之一。所以我的问题是:为什么mldivide返回错误的解决方案?

matlab
  • 1 1 个回答
  • 22 Views

1 个回答

  • Voted
  1. Best Answer
    James Tursa
    2024-11-06T03:01:11+08:002024-11-06T03:01:11+08:00

    您的矩阵 A 是奇异的。MATLAB 文档指出,在这些情况下,mldivide 不可靠(“... 使用病态矩阵时,可能会产生不可靠的解决方案...”),并建议改用 lsqminnorm( ) 或 pinv( )(参见“提示”)。例如,

    >> A = [ -1 0  0;
        1 1 -1;
        0 -1 1 ];
    >> B = [-1; 0; 1];
    >> A\B
    Warning: Matrix is singular to working precision. 
     
    ans =
       NaN
       NaN
       NaN
    >> A*ans-B
    ans =
       NaN
       NaN
       NaN
    >> lsqminnorm(A,B)
    ans =
        1.0000
       -0.5000
        0.5000
    >> A*ans-B
    ans =
       1.0e-15 *
        0.4441
       -0.3331
       -0.1110
    >> pinv(A)*B
    ans =
        1.0000
       -0.5000
        0.5000
    >> A*ans-B
    ans =
       1.0e-15 *
        0.2220
       -0.4441
        0.2220
    
    • 1

相关问题

  • 具有复数的 matlab accumarray

  • 当散点图符号在图中半透明时如何使Matlab图例显示不透明的绘图符号

  • 使用具有相同形状的一维卷积,因此它可以与 FFT 一起使用吗?

  • 使用文件范围变量进行单元测试,无需类

  • 在 MATLAB 2022 中使用 ode45 时出现未定义的结果 (NaN)

Sidebar

Stats

  • 问题 205573
  • 回答 270741
  • 最佳答案 135370
  • 用户 68524
  • 热门
  • 回答
  • Marko Smith

    Vue 3:创建时出错“预期标识符但发现‘导入’”[重复]

    • 1 个回答
  • Marko Smith

    为什么这个简单而小的 Java 代码在所有 Graal JVM 上的运行速度都快 30 倍,但在任何 Oracle JVM 上却不行?

    • 1 个回答
  • Marko Smith

    具有指定基础类型但没有枚举器的“枚举类”的用途是什么?

    • 1 个回答
  • Marko Smith

    如何修复未手动导入的模块的 MODULE_NOT_FOUND 错误?

    • 6 个回答
  • Marko Smith

    `(表达式,左值) = 右值` 在 C 或 C++ 中是有效的赋值吗?为什么有些编译器会接受/拒绝它?

    • 3 个回答
  • Marko Smith

    何时应使用 std::inplace_vector 而不是 std::vector?

    • 3 个回答
  • Marko Smith

    在 C++ 中,一个不执行任何操作的空程序需要 204KB 的堆,但在 C 中则不需要

    • 1 个回答
  • Marko Smith

    PowerBI 目前与 BigQuery 不兼容:Simba 驱动程序与 Windows 更新有关

    • 2 个回答
  • Marko Smith

    AdMob:MobileAds.initialize() - 对于某些设备,“java.lang.Integer 无法转换为 java.lang.String”

    • 1 个回答
  • Marko Smith

    我正在尝试仅使用海龟随机和数学模块来制作吃豆人游戏

    • 1 个回答
  • Martin Hope
    Aleksandr Dubinsky 为什么 InetAddress 上的 switch 模式匹配会失败,并出现“未涵盖所有可能的输入值”? 2024-12-23 06:56:21 +0800 CST
  • Martin Hope
    Phillip Borge 为什么这个简单而小的 Java 代码在所有 Graal JVM 上的运行速度都快 30 倍,但在任何 Oracle JVM 上却不行? 2024-12-12 20:46:46 +0800 CST
  • Martin Hope
    Oodini 具有指定基础类型但没有枚举器的“枚举类”的用途是什么? 2024-12-12 06:27:11 +0800 CST
  • Martin Hope
    sleeptightAnsiC `(表达式,左值) = 右值` 在 C 或 C++ 中是有效的赋值吗?为什么有些编译器会接受/拒绝它? 2024-11-09 07:18:53 +0800 CST
  • Martin Hope
    The Mad Gamer 何时应使用 std::inplace_vector 而不是 std::vector? 2024-10-29 23:01:00 +0800 CST
  • Martin Hope
    Chad Feller 在 5.2 版中,bash 条件语句中的 [[ .. ]] 中的分号现在是可选的吗? 2024-10-21 05:50:33 +0800 CST
  • Martin Hope
    Wrench 为什么双破折号 (--) 会导致此 MariaDB 子句评估为 true? 2024-05-05 13:37:20 +0800 CST
  • Martin Hope
    Waket Zheng 为什么 `dict(id=1, **{'id': 2})` 有时会引发 `KeyError: 'id'` 而不是 TypeError? 2024-05-04 14:19:19 +0800 CST
  • Martin Hope
    user924 AdMob:MobileAds.initialize() - 对于某些设备,“java.lang.Integer 无法转换为 java.lang.String” 2024-03-20 03:12:31 +0800 CST
  • Martin Hope
    MarkB 为什么 GCC 生成有条件执行 SIMD 实现的代码? 2024-02-17 06:17:14 +0800 CST

热门标签

python javascript c++ c# java typescript sql reactjs html

Explore

  • 主页
  • 问题
    • 最新
    • 热门
  • 标签
  • 帮助

Footer

AskOverflow.Dev

关于我们

  • 关于我们
  • 联系我们

Legal Stuff

  • Privacy Policy

Language

  • Pt
  • Server
  • Unix

© 2023 AskOverflow.DEV All Rights Reserve