逆矩阵、伪逆矩阵:数据的压缩和复原

2024-04-30 17:58

本文主要是介绍逆矩阵、伪逆矩阵:数据的压缩和复原,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

逆矩阵、伪逆矩阵、数据的压缩和复原:这一块知识虽然很简单,但在光学各种实验情况下经常用到,特此总结。

 

矩阵的乘法 (观测矩阵):C = A* B

C:M*1

A:M*N

B:N*1

A可以理解为一个转化矩阵,或者说 观测矩阵。

矩阵B在观测矩阵A上,观察到的效果是矩阵C

这个思想在物理思想上非常重要,尤其是光学。

实际生活中,往往我们需要通过观测矩阵(A)和观察到的结果(C),来求出原矩阵(B)

也就是解非线性齐次方程,M个方程,N个未知数。

 

求逆矩阵当矩阵A必须是方阵的时候,A才是可逆矩阵,这时M=N,有唯一解,B = inv(A) * C,matlab中测试如下。

%测试逆矩阵
A = rand(10,10);
B = rand(10,1);
C = A*B;
BB = inv(A)*C;
B'
BB'结果:ans =0.0986    0.1420    0.1683    0.1962    0.3175    0.3164    0.2176    0.2510    0.8929    0.7032ans =0.0986    0.1420    0.1683    0.1962    0.3175    0.3164    0.2176    0.2510    0.8929    0.7032

求伪逆矩阵:伪逆矩阵是逆矩阵的广义形式。由于奇异矩阵或非方阵的矩阵不存在逆矩阵,但在matlab里可以用函数pinv(A)求其伪逆矩阵。理解其实就是求不存在逆矩阵的逆矩阵。

M>N的时候,也就是说我们探测到的数据量大于原数据量,此时方程数量大于未知数数量,可以直接通过求伪逆求出矩阵B,B = pinv(A) * C

%测试伪逆,M > N
A = rand(10,8);
B = rand(8,1);
C = A*B;
BB = pinv(A)*C;
B'
BB'结果:
ans =0.1375    0.3900    0.9274    0.9175    0.7136    0.6183    0.3433    0.9360ans =0.1375    0.3900    0.9274    0.9175    0.7136    0.6183    0.3433    0.9360

N>M的时候,也就是说我们探测到的数据量大于原数据量,此时方程数量小于未知数数量,方程没有唯一解,不能通过B = pinv(A) * C 求出B(压缩感知情况另外考虑 :https://blog.csdn.net/tyfwin/article/details/88902091)

%测试伪逆 M<N
A = rand(8,10);
B = rand(10,1);
C = A*B;
BB = pinv(A)*C;
B'
BB'结果:
ans =0.7360    0.7947    0.5449    0.6862    0.8936    0.0548    0.3037    0.0462    0.1955    0.7202ans =0.7045    0.7584    0.6073    0.5426    1.0107    0.0920    0.2619    0.1060    0.2999    0.6264

 

这篇关于逆矩阵、伪逆矩阵:数据的压缩和复原的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



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

相关文章

大模型研发全揭秘:客服工单数据标注的完整攻略

在人工智能(AI)领域,数据标注是模型训练过程中至关重要的一步。无论你是新手还是有经验的从业者,掌握数据标注的技术细节和常见问题的解决方案都能为你的AI项目增添不少价值。在电信运营商的客服系统中,工单数据是客户问题和解决方案的重要记录。通过对这些工单数据进行有效标注,不仅能够帮助提升客服自动化系统的智能化水平,还能优化客户服务流程,提高客户满意度。本文将详细介绍如何在电信运营商客服工单的背景下进行

基于MySQL Binlog的Elasticsearch数据同步实践

一、为什么要做 随着马蜂窝的逐渐发展,我们的业务数据越来越多,单纯使用 MySQL 已经不能满足我们的数据查询需求,例如对于商品、订单等数据的多维度检索。 使用 Elasticsearch 存储业务数据可以很好的解决我们业务中的搜索需求。而数据进行异构存储后,随之而来的就是数据同步的问题。 二、现有方法及问题 对于数据同步,我们目前的解决方案是建立数据中间表。把需要检索的业务数据,统一放到一张M

关于数据埋点,你需要了解这些基本知识

产品汪每天都在和数据打交道,你知道数据来自哪里吗? 移动app端内的用户行为数据大多来自埋点,了解一些埋点知识,能和数据分析师、技术侃大山,参与到前期的数据采集,更重要是让最终的埋点数据能为我所用,否则可怜巴巴等上几个月是常有的事。   埋点类型 根据埋点方式,可以区分为: 手动埋点半自动埋点全自动埋点 秉承“任何事物都有两面性”的道理:自动程度高的,能解决通用统计,便于统一化管理,但个性化定

使用SecondaryNameNode恢复NameNode的数据

1)需求: NameNode进程挂了并且存储的数据也丢失了,如何恢复NameNode 此种方式恢复的数据可能存在小部分数据的丢失。 2)故障模拟 (1)kill -9 NameNode进程 [lytfly@hadoop102 current]$ kill -9 19886 (2)删除NameNode存储的数据(/opt/module/hadoop-3.1.4/data/tmp/dfs/na

异构存储(冷热数据分离)

异构存储主要解决不同的数据,存储在不同类型的硬盘中,达到最佳性能的问题。 异构存储Shell操作 (1)查看当前有哪些存储策略可以用 [lytfly@hadoop102 hadoop-3.1.4]$ hdfs storagepolicies -listPolicies (2)为指定路径(数据存储目录)设置指定的存储策略 hdfs storagepolicies -setStoragePo

Hadoop集群数据均衡之磁盘间数据均衡

生产环境,由于硬盘空间不足,往往需要增加一块硬盘。刚加载的硬盘没有数据时,可以执行磁盘数据均衡命令。(Hadoop3.x新特性) plan后面带的节点的名字必须是已经存在的,并且是需要均衡的节点。 如果节点不存在,会报如下错误: 如果节点只有一个硬盘的话,不会创建均衡计划: (1)生成均衡计划 hdfs diskbalancer -plan hadoop102 (2)执行均衡计划 hd

hdu1043(八数码问题,广搜 + hash(实现状态压缩) )

利用康拓展开将一个排列映射成一个自然数,然后就变成了普通的广搜题。 #include<iostream>#include<algorithm>#include<string>#include<stack>#include<queue>#include<map>#include<stdio.h>#include<stdlib.h>#include<ctype.h>#inclu

hdu1565(状态压缩)

本人第一道ac的状态压缩dp,这题的数据非常水,很容易过 题意:在n*n的矩阵中选数字使得不存在任意两个数字相邻,求最大值 解题思路: 一、因为在1<<20中有很多状态是无效的,所以第一步是选择有效状态,存到cnt[]数组中 二、dp[i][j]表示到第i行的状态cnt[j]所能得到的最大值,状态转移方程dp[i][j] = max(dp[i][j],dp[i-1][k]) ,其中k满足c

【Prometheus】PromQL向量匹配实现不同标签的向量数据进行运算

✨✨ 欢迎大家来到景天科技苑✨✨ 🎈🎈 养成好习惯,先赞后看哦~🎈🎈 🏆 作者简介:景天科技苑 🏆《头衔》:大厂架构师,华为云开发者社区专家博主,阿里云开发者社区专家博主,CSDN全栈领域优质创作者,掘金优秀博主,51CTO博客专家等。 🏆《博客》:Python全栈,前后端开发,小程序开发,人工智能,js逆向,App逆向,网络系统安全,数据分析,Django,fastapi

烟火目标检测数据集 7800张 烟火检测 带标注 voc yolo

一个包含7800张带标注图像的数据集,专门用于烟火目标检测,是一个非常有价值的资源,尤其对于那些致力于公共安全、事件管理和烟花表演监控等领域的人士而言。下面是对此数据集的一个详细介绍: 数据集名称:烟火目标检测数据集 数据集规模: 图片数量:7800张类别:主要包含烟火类目标,可能还包括其他相关类别,如烟火发射装置、背景等。格式:图像文件通常为JPEG或PNG格式;标注文件可能为X