当前位置：首页 > news >正文

2025-04-26-利用奇异值重构矩阵-美团

news 来源：原创 2025/7/12 23:48:49

2025-04-26-利用奇异值重构矩阵-美团

题目内容

在一家致力于图像处理的科技公司，你被分配到一个新项目，目标是开发一种图像压缩算法，以减少存储空间并加速传输。团队决定使用奇异值分解（ $S V D$ ）对图像进行降维处理，以达到压缩的目的。

现在，你需要编写一个程序，对给定的灰度图像矩阵进行奇异值分解，并重构出近似的低秩矩阵。

请你帮助团队实现一个使用 $N u m P y$ 库的程序，对给定的矩阵进行奇异值分解，共利用前 $(k)$ 个奇异值重构矩阵。具体要求如下：

读取输入矩阵，为一个二维列表，表示灰度图像的像素值矩阵。
读取整数 $(k)$ ，表示使用前 $(k)$ 个奇异值进行矩阵重构。
对矩阵进行奇异值分解( $S V D$ )，获取左奇异矩阵 $(U)$ 、奇异值对角矩阵 $(\Sigma)$ 和右奇异矩阵 $V^T)$ 。
利用前 $(k)$ 个奇异值和对应的奇异向量重构矩阵。
输出重构后的矩阵，每个元素保留两位小数（使用 $ro u n d (x, 2)$ ）。

输入描述

第一行包含两个整数 $(m)$ 和 $(n)$ 表示矩阵的行数和列数。
接下来的 $(m)$ 行，每行包含 $(n)$ 个整数，表示矩阵的元素，元素之间用空格分隔。
最后一行包含一个整数 $(k)$ ，表示用前 $(k)$ 个奇异值进行矩阵重构。

输出描述

输出 $(m)$ 行，每行包含 $(n)$ 个浮点数，表示重构后的矩阵元素，元素之间用空格分隔，结果均保留两位小数，使用 $ro u n d (x, 2)$ 。

补充说明

奇异值分解（ $S V D$ ），对于一个矩阵 $(A)$ 其奇异值分解为：
$\Sigma V^T$
其中：4 5
- $(U)$ 是 $\times m)$ 的正交矩阵，其列为 $AA^T)$ 的特征向量。
- $(\Sigma)$ 是 $\times n)$ 的对角矩阵，对角线上为非负的奇异值，按降序排列。
- $V^T)$ 是 $\times n)$ 的正交矩阵的转置，其列为 $A^TA)$ 的特征向量。
矩阵重构

使用前 $(k)$ 个奇异值和对应的奇异向量，可以近似地重构原矩阵：
$A_k = U_k \Sigma_k V_k^T$
其中：
- $U_k)$ 为 $(U)$ 的前 $(k)$ 列。
- $(\Sigma_k)$ 为 $(\Sigma)$ 中前 $(k)$ 个奇异值构成的对角矩阵。
- $V_k^T)$ 为 $V^T)$ 的前 $(k)$ 行。

样例

输入:
4 5
52 55 61 66 70
63 59 55 90 109
85 104 117 123 119
105 122 145 160 159
2输出:
47.3 54.1 61.41 69.2 70.35
63,02 58.35 55.4 90.03 109.11
84.08 101.15 118.72 123.8 119.52
107.77 125.0 143.24 157.95 158.38

import numpy as np

n,m = map(int, input().split())
mat = []
for _ in range(n):mat.append(list(map(int, input().split())))
k = int(input())

A = np.array(mat)
A

array([[ 52,  55,  61,  66,  70],[ 63,  59,  55,  90, 109],[ 85, 104, 117, 123, 119],[105, 122, 145, 160, 159]])

eig_values, U = np.linalg.eig(A@(A.T))
U

array([[ 0.30050372, -0.02694092,  0.69020128,  0.6577187 ],[ 0.37778194,  0.90391529,  0.05467   , -0.19294855],[ 0.54241183, -0.38739585,  0.36703081, -0.64884685],[ 0.68758106, -0.17926459, -0.62122602,  0.33041599]])

eig_values = sorted(np.sqrt(eig_values),reverse=True)
eig_values

[454.71798467974327, 34.01883573921087, 7.740636522193636, 3.369832336013565]

_, V_t = np.linalg.eig((A.T)@A)
V_t = V_t.T
V_t

array([[ 0.34686878,  0.39389812,  0.44482554,  0.50704668,  0.51919189],[-0.11153546,  0.30307318,  0.68334892, -0.09530755, -0.64780955],[-0.68515863, -0.46097385,  0.26174063,  0.48804619,  0.1065998 ],[ 0.41943452, -0.20543104, -0.25978011,  0.65168322, -0.53823419],[ 0.47103526, -0.70589809,  0.44628536, -0.26634569,  0.09860531]])

U_k = U[:, :k]
sigma = np.diag(eig_values[:k])
V_tk = V_t[:k, :]
A_k = np.round(U_k@sigma@V_tk, 2)
A_k

array([[ 47.5 ,  53.55,  60.16,  69.37,  71.54],[ 56.16,  76.99,  97.43,  84.17,  69.27],[ 87.02,  93.16, 100.71, 126.32, 136.59],[109.13, 121.31, 134.91, 159.11, 166.28]])

基础概念

求奇异值

奇异值可以通过以下步骤计算：

计算矩阵 $A$ 的转置 $A^T$ 。
计算 $A^TA$ 和 $AA^T$ 。
求解 $A^TA$ 和 $AA^T$ 的特征值。
奇异值是 $A^TA$ 和 $AA^T$ 的特征值的平方根。

numpy求矩阵特征值

import numpy as np# 定义一个矩阵
A = np.array([[1, 2], [3, 4]])# 使用 numpy.linalg.eig 计算特征值和特征向量
eigenvalues, eigenvectors = np.linalg.eig(A)# 打印结果
print("矩阵A：")
print(A)
print("\n特征值：")
print(eigenvalues)
print("\n特征向量：")
print(eigenvectors)

矩阵A：
[[1 2][3 4]]特征值：
[-0.37228132  5.37228132]特征向量：
[[-0.82456484 -0.41597356][ 0.56576746 -0.90937671]]

eigenvalues = np.linalg.eigvals(A)
print("特征值：", eigenvalues)

特征值： [-0.37228132  5.37228132]

numpy求奇异值分解

直接用这个算出的答案不太一样

import numpy as np# 创建一个示例矩阵
A = np.array([[1, 2], [3, 4], [5, 6]])# 进行奇异值分解
U, S, Vh = np.linalg.svd(A, full_matrices=True, compute_uv=True)print("U 矩阵 (左奇异向量):")
print(U)
print("\n奇异值 (S):")
print(S)
print("\nVh 矩阵 (右奇异向量的转置):")
print(Vh)

U 矩阵 (左奇异向量):
[[-0.2298477   0.88346102  0.40824829][-0.52474482  0.24078249 -0.81649658][-0.81964194 -0.40189603  0.40824829]]奇异值 (S):
[9.52551809 0.51430058]Vh 矩阵 (右奇异向量的转置):
[[-0.61962948 -0.78489445][-0.78489445  0.61962948]]

2025-04-26-利用奇异值重构矩阵-美团

2025-04-26-利用奇异值重构矩阵-美团题目内容在一家致力于图像处理的科技公司，你被分配到一个新项目，目标是开发一种图像压缩算法，以减少存储空间并加速传输。团队决定使用奇异值分解（ S V D SVD SVD）对图像进行降…...

编程日记 2025/7/12 23:48:49

《解锁SCSS算术运算：构建灵动样式的奥秘》

SCSS作为CSS预处理器，算术运算功能犹如一颗璀璨明珠，赋予我们动态计算样式属性值的强大能力，让网页样式不再是一成不变的刻板呈现，而是能够根据各种条件和需求灵动变化。在SCSS的世界里，算术运算绝非孤立的存在&…...

编程日记 2025/6/29 5:10:07

STM32Cube-FreeRTOS任务管理工具函数-笔记

STM32Cube-FreeRTOS任务管理工具函数-笔记一、获取任务句柄的函数1. 创建任务并获取句柄2. 获取当前任务句柄3. 获取空闲任务句柄4. 根据任务名称获取句柄二、单个任务操作相关函数1. 程序在运行时可以获取或改变一个任务的优先级3. 获取任务信息4. 获取任务信息5. 获取任务名…...

编程日记 2025/7/4 10:07:32

【第三十四周】多模态大模型调研

多模态大模型调研摘要Abstract引言多模态技术的主要方向视觉-语言大模型（Vision-Language Large Models, VLLMs）语音-语言大模型（Speech-Language Large Models, SLLMs）音乐 - 语言大模型（Music-Language Large Model…...

编程日记 2025/7/7 3:46:44

【2025最新】Baichuan-M1-instruct部署教程

首先机器至少要A100、4090、3090 这里选AutoDL的4090D，运行至少要20G显存。这里镜像选基础镜像11.8【更新！！！！！！！！！24G带不动！显存不够】有时候…...

编程日记 2025/7/10 18:02:19

Unity与Unreal Engine（UE）的深度解析及高级用法

以下是Unity与Unreal Engine（UE）的深度解析及高级用法对比，结合技术特性、行业应用与未来发展进行综合阐述：一、核心差异与适用场景对比 1. 技术架构与编程模式 Unity 语言与脚本：主要使用C#，语法简洁且易于学习，适合快速原型开发和中小型项目。支持可视化脚本工具（如…...

编程日记 2025/6/29 7:26:45

网络：TCP三次握手、四次挥手

目录深刻理解三次握手深刻理解四次挥手深刻理解三次握手三次握手时，如果最后一个ACK包，服务器没有收到，此时： 客户端：认为已经建立链接服务器：认为没有建立链接，还在超时等待。而此时…...

编程日记 2025/6/29 8:20:24

spdlog自定义formatter

用了之后发现，spdlog的默认日志记录格式为： [2014-10-31 23:46:59.678] [my_loggername] [info] Some message 但是这个格式不是我想要的，怎么办，这个也简单，上面的内容也就是几个标签的组合而已，spdlog自定…...

编程日记 2025/7/2 15:21:19

Spring AI 实战：第四章、Spring AI多模态之看图说话

引言：从"码农"到"多媒体魔术师" “曾经，我们的代码核心擅长处理文本，就像餐厅里只会做炒饭的厨师。现在有了Spring AI多模态支持，我们突然拥有满汉全席的烹饪技巧！” 作为一名常年与String打交道的开发者，当第一次看到Spring AI可以同时处理图片、音…...

编程日记 2025/6/29 6:09:53

ES6入门---第二单元模块五：模块化

js不支持模块化注意： 需要放到服务器环境 1、如何定义模块？ export 东西例：1.js文件中 console.log(1模块加载了);//显示是否加载了 export const a 12; export const b 5; export let c 101; const a12; const b5; const c101;ex…...

编程日记 2025/6/29 7:57:39

Python 函数装饰器和闭包(变量作用域规则)

本章内容： Python 如何计算装饰器句法 Python 如何判断变量是不是局部的闭包存在的原因和工作原理 nonlocal 能解决什么问题掌握这些基础知识后，我们可以进一步探讨装饰器： 实现行为良好的装饰器标准库中有用的装饰器实现一个参数化装饰器…...

编程日记 2025/7/12 23:39:41

什么是constexpr？

什么是constexpr？ 简单来说，constexpr就是告诉编译器：“我这个变量或函数的值可以在编译时算出来，请帮我提前算好，运行时直接用结果，不用再算了。” • **传统const**只表示变量不可修改，但不…...

编程日记 2025/6/29 5:42:42

如何在 PowerEdge 服务器上设置 NIC 分组

以下文章提供了有关 Windows、VMware 和 Linux 中的 NIC 分组的信息。什么是网络适配器分组？设置 NIC 分组 Windows设置 NIC 分组 VMware设置 NIC 分组 Linux 什么是网络适配器分组（绑定）？ 网络适配器分组是一个术语&#xff0…...

编程日记 2025/6/29 5:22:15

ES6入门---第三单元模块四：Set和WeakSet

set数据结构: 类似数组，但是里面不能有重复值，如果有，只显示一个 set用法: let setArr new Set([a,b]); setArr.add(a); 往setArr里面添加一项 let setArr new Set().add(a).add(b).add(c); setArr.delete(b); 删除一项 setArr.ha…...

编程日记 2025/6/29 8:01:42

架构进阶：75页架构规划方法课件【附全文阅读】

本文概述了一个关于架构规划方法的目录及其目的，重点介绍了基于联邦企业架构（FEAF）的架构建模方法，并提及了不同层面的架构建模方法以及培训的目的。以下是对该内容的简洁总结： **架构规划方法目录及其目的** 本文旨在…...

编程日记 2025/7/4 19:26:46

前端面经-VUE3篇（三）--vue Router（二)导航守卫、路由元信息、路由懒加载、动态路由

一、导航守卫 vue Router 中的导航守卫（Navigation Guards） 是一个非常重要的功能，用于在路由切换过程中，拦截、控制、检查或延迟页面跳转。你可以理解为： 🔐 “进门前的保安”，控制哪些页面…...

编程日记 2025/6/29 5:09:19

RTX-3090 Qwen3-8B Dify RAG环境搭建

RTX-3090 Qwen3-8B Dify RAG环境搭建一、环境配置二、操作步骤1、创建容器2、下载`Qwen3-8B`和embedding模型3、安装`transformers`4、安装`vllm`5、安装`flash-attention`6、启动兼容OpenAI API的服务1、方案一:启动`vllm`服务【不支持多任务】2、方案二:Flask和PyTorch实现的…...

编程日记 2025/6/29 7:47:42

Circular Plot系列(三)：【视频教程】复现NCS图表之高大上的单细胞UMAP环形图

高端复杂的UMAP复现： 这又是一个高大上且炫酷的单细胞UMAP图，展示的信息很多，有大类细胞和亚群，以及marker基因和cell count信息，还可以增加其他的分组信息等等。没错，看这个图就是circlize一层层画的。我们…...

编程日记 2025/6/29 5:32:25

MCP智能体多Agent协作系统设计（Multi-Agent Cooperation）

目录 🚀 MCP智能体多Agent协作系统设计（Multi-Agent Cooperation） 🌟 为什么需要多Agent协作？ 🧠 多Agent协作系统架构设计 🛠️ 1. 构建基础智能体基类（Agent Base）…...

编程日记 2025/7/9 16:15:02

栈Stack

一栈:先进后出一种特殊的线性表，其只允许在固定的一端进行插入和删除元素操作。进行数据插入和删除操作的一端称为栈顶，另一端称为栈底。栈中的数据元素遵守后进先出LIFO（Last In First Out）的原则。 Stack<Integer> sta…...

编程日记 2025/6/29 7:29:30

HTML01：HTML基本结构

HTML基本结构 <html> <head><meta charset"UTF-8"><title>我的第一个网页</title> </head> <body>我的第一个网页 </body> </html><body、</body等成对的标签，分别叫开发标签和闭合标签单独…...

编程日记 2025/6/29 5:09:23

QT聊天项目DAY07

1.Win配置和使用GRPC 1.1 克隆GRPC库克隆GRPC git clone -b v1.34.0 https://gitee.com/mirrors/grpc-framework.git 查看Git有没有安装没有安装 1.1.1 安装Git https://git-scm.com/ 一路next 添加Git的路径到系统环境变量下我这次没用管理员权限，并且也没…...

编程日记 2025/7/5 13:53:32

2025年PMP 学习三

4.2制定项目管理计划 4. 项目管理计划 - 内容（输出） 项目目标的制定原因： 3个基准（范围基准、进度基准、成本基准） 子管理计划：范围、需求、进度、成本、质量、资源、沟通、风险、采购等管理计划&#xf…...

编程日记 2025/6/29 5:09:22

软考-软件设计师中级备考 10、文件管理、设备管理

一、文件管理 1、文件目录文件控制块（FCB）：是操作系统为管理文件而设置的数据结构，包含了文件的基本信息（如文件名、文件大小、文件类型等）、存取控制信息（如文件所有者的权限、其他用户的权…...

编程日记 2025/7/11 17:23:46

Linux环境下的进程创建-fork函数的使用，进程退出exit和_exit的区别，以及进程等待waitpid和status数据的提取方法

目录一、进程创建 1.fork函数 1）进程调用fork函数是如何创建子进程的 2）代码示范 2.写时拷贝二、进程退出 1.退出码 1）什么是退出码？ 2）为什么要有退出码？ 3）退出码是怎么做到的？ …...

编程日记 2025/6/28 22:36:01

【数据结构与算法】常见排序算法详解（C++实现）

目录一、排序的基本概念二、插入排序 2.1 直接插入排序 2.2 折半插入排序 2.3 希尔排序三、交换排序 3.1 冒泡排序 3.2 快速排序四、选择排序 4.1 简单选择排序 4.2 堆排序五、归并排序六、基数排序七、计数排序结语一、排序的基本概念排序就是重新…...

编程日记 2025/6/29 5:09:29

STM32GPIO输入实战-按键key模板及移植

STM32GPIO输入实战-按键key模板及移植一，按键模板展示二，按键模板逻辑1，准备工作：头文件与全局变量2，读取硬件状态：key_read_raw()3,核心处理：key_process_simple() 的四行代码三，…...

编程日记 2025/6/29 5:09:24

LeetCode 1128.等价多米诺骨牌对的数量：计数

【LetMeFly】1128.等价多米诺骨牌对的数量：计数力扣题目链接：https://leetcode.cn/problems/number-of-equivalent-domino-pairs/ 给你一组多米诺骨牌 dominoes 。形式上，dominoes[i] [a, b] 与 dominoes[j] [c, d] 等价当且仅当 (a …...

编程日记 2025/7/11 21:24:42

Spring MVC设计与实现

DispatcherServlet的初始化与请求处理流程初始化阶段 Servlet 生命周期触发：当 Web 容器（如 Tomcat）启动时，根据注解/配置，DispatcherServlet 的 init() 方法被调用。初始化 WebApplicationContext 根 WebApplicat…...

编程日记 2025/6/29 7:54:02

日语学习-日语知识点小记-进阶-JLPT-N1阶段（1）：语法单词

日语学习-日语知识点小记-进阶-JLPT-N1阶段（1）：语法单词 1、前言（1）情况说明（2）工程师的信仰（3）高级语法N1语法和难点一、N1语法学习内容（高级语法&#xff…...

编程日记 2025/7/12 23:42:48

stm32week14

stm32学习十.GPIO 2.基本结构基本结构： F1与其它的的最大区别是上下拉电阻的位置施密特触发器是一种整形电路，可以将非标准方波，整形成方波图中MOS管的输出规则： 3.8中工作模式 ①输入浮空： 上下拉电阻均不工…...

编程日记 2025/6/29 5:37:12

WPF中Binding

绑定ViewModel中的数据添加数据上下文方法一：在XAML中添加 <Window.DataContext><local:MainWindowViewModel /> </Window.DataContext>方法二：在界面层的cs文件中添加 this.DataContext new MainWindowViewModel();绑定 publ…...

编程日记 2025/6/29 7:36:07

Google Agent space时代，浅谈Agent2Agent (A2A) 协议和挑战！

如果说去年Google Cloud大会大家还在数“AI”这个词被提了多少次，那么今年，绝对是“Agent”的主场！开发者主题演讲几乎被它“刷屏”，展区的许多 Demo 也都号称是 Agent 应用。但我得诚实地说，大会现场关于 Agents 的 …...

编程日记 2025/7/8 14:28:25

爬虫的应用

在自然语言处理（NLP）领域，文本数据的预处理是至关重要的基础环节。它如同工匠雕琢璞玉前的打磨工作，直接影响后续模型分析与挖掘的效果。本文将基于 Python，以电商平台的差评和优质评价文本数据为例，详细展…...

编程日记 2025/6/29 8:11:18

Day 41 题目描述做法 /*** Definition for a binary tree node.* public class TreeNode {* int val;* TreeNode left;* TreeNode right;* TreeNode() {}* TreeNode(int val) { this.val val; }* TreeNode(int val, TreeNode left, TreeNode right…...

编程日记 2025/6/29 6:28:38

Java后端开发day40--异常File

（以下内容全部来自上述课程） 异常异常：异常就是代表程序出现的问题 1. 异常的分类 1.1 Error 代表的是系统级别的错误（属于严重问题） 系统一旦出现问题，sun公司会把这些错误封装成Error对象。 Error…...

编程日记 2025/6/29 8:06:24

集成算法学习

集成算法通过组合多个弱学习器提升模型性能，以下是核心内容详解： 一、核心思想 - 组合优势：结合多个简单模型（如决策树），通过“少数服从多数”或“加权平均”等策略，降低方差、偏差或过拟合风险…...

编程日记 2025/6/29 8:06:36

工业认知智能：从数据分析到知识创造

工业认知智能：从数据分析到知识创造引言当前制造业面临的知识管理困境令人震惊：68%的工艺知识存储于老员工头脑中，30%的企业因知识传承断层导致质量事故。麦肯锡研究显示，应用认知智能技术的企业，其工艺创新速度提升3-5倍。本文将系统阐述工业认知智能的"感知-理…...

编程日记 2025/6/29 7:36:20

8.1 Python+Docker+企业微信集成实战：自动化报告生成与CI/CD部署全攻略

Python+Docker+企业微信集成实战：自动化报告生成与CI/CD部署全攻略关键词：PDF报告生成, Word文档自动化, 企业微信集成, Docker容器化, CI/CD流水线 1. 多格式报告生成实战通过扩展报告输出格式，满足不同用户的文档需求。我们使用Python生态的成熟库实现PDF/Word生成，并…...

编程日记 2025/6/29 5:11:13

25.5.4数据结构|哈夫曼树学习笔记

知识点前言一、搞清楚概念 ●权：___________ ●带权路径长度：__________ WPL所有的叶子结点的权值*路径长度之和 ●前缀编码：____________ 二、构造哈夫曼树 n个带权值的结点，构造哈夫曼树算法： 1、转化成n棵树组成的…...

编程日记 2025/6/29 5:11:02

统计学中的p值是什么？怎么使用？

李升伟整理在统计学中，p值（p-value）是帮助研究者判断假设检验结果是否具有统计显著性的重要指标。以下是关于p值的详细解释和使用方法： 1. p值的定义 p值表示在原假设（H0）为真的情况下，观察…...

编程日记 2025/6/29 6:50:45

22：一维码与二维码区别

一维码（条形码） 一维条码即指条码条和空的排列规则，常用的一维码的码制包括：EAN码、39码、交叉25码、UPC码、128码、93码，ISBN码，及Codabar（库德巴码）等。条码是由一组规则排列的条…...

编程日记 2025/6/29 8:33:51

Java学习手册：SQL 优化技巧

一、SQL 查询优化选择合适的索引列 ：索引可以显著提高查询速度，但需要选择合适的列来创建索引。通常，对于频繁作为查询条件的列、连接操作的列以及排序或分组操作的列，应该考虑创建索引。例如，在一个订单表中&#xf…...

编程日记 2025/6/29 3:53:49

《Vue3学习手记8》

vue3中的一些API shallowRef ( ) 和shallowReactive ( ) shallowRef (浅层响应式) 1.作用:创建一个响应式数据，但只对顶层属性进行响应式处理。 2.用法: const originalref(...) const original2shallowRef(original) 3.特点:只跟踪引用值的变化，不关心…...

编程日记 2025/7/12 20:48:54

平衡二叉搜索树模拟实现1-------AVL树(插入,删除,查找)

本章目标 1.AVL树的概念 2.AVL树的模拟实现 1.AVL树的概念 1.AVL树是最先被发明的平衡二叉搜索树,AVL树是一颗空树或者具有以下的性质它的左右子树都是AVL树,并且左右高度差不超过1,AVL树是一颗高度平衡二叉搜索树,通过高度差去控制平衡 2.为什么高度差是1? 当结点个数为8…...

编程日记 2025/7/10 15:35:05