矩阵病态问题

2023-10-12 06:04
文章标签 问题 矩阵 病态

本文主要是介绍矩阵病态问题,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

  • 病态(ill-conditioned)矩阵
  • 病态(ill-posed)方程

病态矩阵

  • 矩阵的frobenius范数:
    • 对每一个元素平方,求和再开平方

A=\begin{bmatrix} a &b \\ c& d \end{bmatrix}

\begin{Vmatrix} A \end{Vmatrix}=\sqrt{|a|^2+|b|^2+|c|^2+|d|^2}

  • 矩阵的条件数:矩阵的frobenius范数与该矩阵的逆的frobenius范数的乘积

cond(A)=\begin{Vmatrix} A \end{Vmatrix}\cdot\begin{Vmatrix} A^{-1} \end{Vmatrix}

  • 对于同阶矩阵,矩阵的条件数越大,他的病态程度越高 
    • 一般而言,条件数小于100,矩阵健康
    • 介于100-1000之间,有一定的病态趋势
    • 大于1000,病态

病态矩阵的实质

  • 矩阵向量(行或列)之间存在共线性

举例

\begin{bmatrix} 1 & 2\\ 3.00001 & 6 \end{bmatrix}\begin{bmatrix} x_1\\ x_2 \end{bmatrix}=\begin{bmatrix} 6\\ 18 \end{bmatrix}

import numpy as npA = np.array([[1, 2], [3.00001, 6]])
b = np.array([6, 18])
def cond(A):a = np.sum(A**2)**0.5b = np.sum(np.linalg.inv(A)**2)**0.5return a*bx = np.linalg.solve(A, b)
print(x, "  ", cond(A))
[0. 3.]    2500003.000017089


 

\begin{bmatrix} 1 & 2\\ 2 & 6 \end{bmatrix}\begin{bmatrix} x_1\\ x_2 \end{bmatrix}=\begin{bmatrix} 6\\ 18 \end{bmatrix}

import numpy as npA = np.array([[1, 2], [2, 6]])
b = np.array([6, 18])
def cond(A):a = np.sum(A**2)**0.5b = np.sum(np.linalg.inv(A)**2)**0.5return a*bx = np.linalg.solve(A, b)
print(x, "  ", cond(A))
[0. 3.]    22.500000000000004


 

这篇关于矩阵病态问题的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/193790

相关文章

怎样通过分析GC日志来定位Java进程的内存问题

《怎样通过分析GC日志来定位Java进程的内存问题》:本文主要介绍怎样通过分析GC日志来定位Java进程的内存问题,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录一、GC 日志基础配置1. 启用详细 GC 日志2. 不同收集器的日志格式二、关键指标与分析维度1.

Java 线程安全与 volatile与单例模式问题及解决方案

《Java线程安全与volatile与单例模式问题及解决方案》文章主要讲解线程安全问题的五个成因(调度随机、变量修改、非原子操作、内存可见性、指令重排序)及解决方案,强调使用volatile关键字... 目录什么是线程安全线程安全问题的产生与解决方案线程的调度是随机的多个线程对同一个变量进行修改线程的修改操

Redis出现中文乱码的问题及解决

《Redis出现中文乱码的问题及解决》:本文主要介绍Redis出现中文乱码的问题及解决,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录1. 问题的产生2China编程. 问题的解决redihttp://www.chinasem.cns数据进制问题的解决中文乱码问题解决总结

全面解析MySQL索引长度限制问题与解决方案

《全面解析MySQL索引长度限制问题与解决方案》MySQL对索引长度设限是为了保持高效的数据检索性能,这个限制不是MySQL的缺陷,而是数据库设计中的权衡结果,下面我们就来看看如何解决这一问题吧... 目录引言:为什么会有索引键长度问题?一、问题根源深度解析mysql索引长度限制原理实际场景示例二、五大解决

Springboot如何正确使用AOP问题

《Springboot如何正确使用AOP问题》:本文主要介绍Springboot如何正确使用AOP问题,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录​一、AOP概念二、切点表达式​execution表达式案例三、AOP通知四、springboot中使用AOP导出

Python中Tensorflow无法调用GPU问题的解决方法

《Python中Tensorflow无法调用GPU问题的解决方法》文章详解如何解决TensorFlow在Windows无法识别GPU的问题,需降级至2.10版本,安装匹配CUDA11.2和cuDNN... 当用以下代码查看GPU数量时,gpuspython返回的是一个空列表,说明tensorflow没有找到

解决未解析的依赖项:‘net.sf.json-lib:json-lib:jar:2.4‘问题

《解决未解析的依赖项:‘net.sf.json-lib:json-lib:jar:2.4‘问题》:本文主要介绍解决未解析的依赖项:‘net.sf.json-lib:json-lib:jar:2.4... 目录未解析的依赖项:‘net.sf.json-lib:json-lib:jar:2.4‘打开pom.XM

IDEA Maven提示:未解析的依赖项的问题及解决

《IDEAMaven提示:未解析的依赖项的问题及解决》:本文主要介绍IDEAMaven提示:未解析的依赖项的问题及解决,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝... 目录IDEA Maven提示:未解析的依编程赖项例如总结IDEA Maven提示:未解析的依赖项例如

Redis分片集群、数据读写规则问题小结

《Redis分片集群、数据读写规则问题小结》本文介绍了Redis分片集群的原理,通过数据分片和哈希槽机制解决单机内存限制与写瓶颈问题,实现分布式存储和高并发处理,但存在通信开销大、维护复杂及对事务支持... 目录一、分片集群解android决的问题二、分片集群图解 分片集群特征如何解决的上述问题?(与哨兵模

SpringBoot+Redis防止接口重复提交问题

《SpringBoot+Redis防止接口重复提交问题》:本文主要介绍SpringBoot+Redis防止接口重复提交问题,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不... 目录前言实现思路代码示例测试总结前言在项目的使用使用过程中,经常会出现某些操作在短时间内频繁提交。例