矩阵定理复习记录

矩阵复习

矩阵导数定理

若A是一个如下矩阵:
A = [ a 11 a 12 a 21 a 22 ] A= \begin{bmatrix}a_{11}&a_{12}\\a_{21}&a_{22}\end{bmatrix} A=[a11a21a12a22]
y是一个向量矩阵:
y ⃗ = [ y 1 y 2 ] \vec{y}=\begin{bmatrix}y_1\\y_2\end{bmatrix} y =[y1y2]
则可得 text 定理:
δ A ∗ y ⃗ δ y ⃗ = A T \frac{δA*\vec{y}}{δ\vec{y}} = A^T δy δAy =AT
δ y ⃗ T ∗ A δ y ⃗ = δ A T ∗ y ⃗ δ y ⃗ = A \frac{δ\vec{y}^T*A}{δ\vec{y}} = \frac{δA^T*\vec{y}}{δ\vec{y}} = A δy δy TA=δy δATy =A
也就是对A*y的矩阵,求偏导y,结果为A的转置矩阵;

还可得另一个定理:
δ y ⃗ T ∗ A y ⃗ δ y ⃗ = A y ⃗ + A T y ⃗ \frac{δ\vec{y}^T*A\vec{y}}{δ\vec{y}} = A\vec{y}+ A^T\vec{y} δy δy TAy =Ay +ATy
若A是一个对称矩阵,也就是 A T = A A^T=A AT=A,则上面的还会等于
2 A y ⃗ 2A\vec{y} 2Ay

δ符号表示求导, y ⃗ 表示一个向量 \vec{y}表示一个向量 y 表示一个向量

这部分的推导过程可参考此篇视频

矩阵平方定理

若矩阵A满足相乘原则,则有定理:
A 2 = A T ∗ A A^2 = A^T*A A2=ATA

单位矩阵

是一种恒等矩阵,对角线上全为1,其余全为0,如下:
I = [ 1 0 0 0 1 0 0 0 1 ] I = \begin{bmatrix}1&0&0\\0&1&0\\0&0&1\end{bmatrix} I= 100010001
任何矩阵与单位矩阵相乘等于本身:
A ∗ I = I ∗ A = A A*I = I * A = A AI=IA=A

逆矩阵

注意逆矩阵不是转置矩阵,若两个矩阵A和B,n*n的方阵,且满足:
A ∗ B = I A*B=I AB=I
也就是矩阵相乘等于单位矩阵,也说明A就是B的逆矩阵,A是可逆的,记:B=A^-1

最小二乘法

若输入量为 x 1 , x 2 . . . x n x_1,x_2...x_n x1,x2...xn,输出量为 y 1 , y 2 . . . y n y_1,y_2...y_n y1,y2...yn,为了你和一条函数曲线,是的输入为 x i x_i xi,输出为 y i y_i yi,我们假定它是一个多项式函数如 y i = a x i 2 + b x i + c y_i = ax_i^2 + bx_i + c yi=axi2+bxi+c,x和y都有观察数据,求 a , b , c a,b,c a,b,c,因为数据又多组,带入矩阵中运算:
[ x 1 2 x 1 1 x 2 2 x 2 1 . . . . . . . . . x n 2 x n 2 1 ] [ a b c ] = [ y 1 y 2 . . . y n ] \begin{bmatrix}x_1^2 & x_1&1\\x_2^2&x_2&1\\...&...&...\\x_n^2&x_n^2&1\end{bmatrix}\begin{bmatrix}a\\b\\c\end{bmatrix}= \begin{bmatrix}y_1\\y_2\\...\\y_n\end{bmatrix} x12x22...xn2x1x2...xn211...1 abc = y1y2...yn
进而用X,A,Y替换:
X ∗ A = Y X*A=Y XA=Y
A矩阵就是我们要求取的未知参数,往往信号观察是在由噪声的环境中的,假设噪声为V,且噪声的均值为0,也就是正和负噪声,则推导公式:
X ∗ A = Y + V X*A=Y+V XA=Y+V
为了使误差最小,使用最小二乘法,二乘差值平方,也就是:
( Y − X ∗ A ) 2 = ( Y − X ∗ A ) T ( Y − X ∗ A ) (Y-X*A)^2 = (Y-X*A)^T(Y-X*A) (YXA)2=(YXA)T(YXA)
对上面的式子A求偏导:
δ ( Y − X ∗ A ) T ( Y − X ∗ A ) δ A \frac{δ(Y-X*A)^T(Y-X*A)}{δA} δAδ(YXA)T(YXA)
推导过程可参考此视频最小二乘法讲解,求出后领偏导函数等于0求极值,也就是误差最小值,得到定理:
A = ( X T X ) − 1 X T Y A = (X^TX)^{-1}X^TY A=(XTX)1XTY

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/36102.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

在vue项目使用数据可视化 echarts ,柱状图、折线图、饼状图使用示例详解及属性详解

官网地址:Apache ECharts ​一、下载插件并在页面中引入 npm install echarts --save 页面导入: import * as echarts from echarts 全局导入: main.js 中,导入并注册到全局 import echarts from echarts Vue.prototype.$echart…

Clone函数

概述 Clone函数是一种用于复制的计算机函数。在程序编写中,除了自定义一个拷贝构造函数来实现对象复制外,还可以实现一个clone函数。这需要借助编译器实现的一个隐藏拷贝构造函数,这样的做法,更省心。 中文名clone函数外文名clon…

C# 使用FFmpeg.Autogen对byte[]进行编解码

C# 使用FFmpeg.Autogen对byte[]进行编解码,参考:https://github.com/vanjoge/CSharpVideoDemo 入口调用类: using System; using System.IO; using System.Drawing; using System.Runtime.InteropServices; using FFmpeg.AutoGen;namespace F…

C++11异步与通信之 packaged_task

概念简介 packaged_task 用于包装可调用目标(Callable)为一个对象,如lambda,普通函数,小括号重载等,用于异步调用。 其返回值或所抛异常被存储于能通过 std::future 对象访问的共享状态中,和promise类似。 将函数的调用与函数返…

时序预测 | MATLAB实现EEMD-LSTM、LSTM集合经验模态分解结合长短期记忆神经网络时间序列预测对比

时序预测 | MATLAB实现EEMD-LSTM、LSTM集合经验模态分解结合长短期记忆神经网络时间序列预测对比 目录 时序预测 | MATLAB实现EEMD-LSTM、LSTM集合经验模态分解结合长短期记忆神经网络时间序列预测对比效果一览基本介绍模型搭建程序设计参考资料 效果一览 基本介绍 时序预测 | …

小龟带你妙写排序之快速排序

快速排序 一. 快速排序原理二. 题目三. 快速排序的思路分析(图文结合)四.代码 一. 快速排序原理 先从数据序列中选一个元素,并将序列中所有比该元素小的元素都放到它的右边或左边,再对左右两边分别用同样的方法处之直到每一个待处…

runtime error: member access within misaligned address(力扣最常见错误之一)

runtime error: member access within misaligned address(力扣最常见错误之一) 前言原因和解决办法总结 前言 最近博主在刷力扣时,明明代码逻辑都没问题,但总是报下面这个错误: runtime error: member access within…

django实现登录和登录的鉴权

1、创建数据库的管理员表 在models.py 中定义admin表,为了简单,表里只有用户名和密码还有默认加的id 三个字段 from django.db import models# Create your models here.class Admin(models.Model):username models.CharField(verbose_name"用户…

源码框架-三勾软件

javaspringboot微信小程序商城SAAS前后端源码: 三勾商城是开发友好的微信小程序商城,框架支持SAAS,支持发布 iOS Android 公众号 H5 各种小程序(微信/支付宝/百度/头条/QQ/钉钉/淘宝)等多个平台,不可多得的二开神器…

训练用于序列分类任务的 RoBERTa 模型的适配器

介绍 NLP当前的趋势包括下载和微调具有数百万甚至数十亿参数的预训练模型。然而,存储和共享如此大的训练模型非常耗时、缓慢且昂贵。这些限制阻碍了 RoBERTa 模型开发更多用途和适应性更强的 NLP 技术,该模型可以从多个任务中学习并针对多个任务进行学习;在本文中,我们将重…

Kafka:安装和配置

producer:发布消息的对象,称为消息产生者 (Kafka topic producer) topic:Kafka将消息分门别类,每一个消息称为一个主题(topic) consumer:订阅消息并处理发布消息的对象…

模拟 枚举

分享牛客算法基础精选题单题目打卡!!! 目录 字符串的展开 多项式输出 机器翻译 : 铺地毯 : [NOIP2016]回文日期 字符串的展开 原题链接 : 字符串的展开 思路 : 模拟 代码 : #include<iostream> #include<cstring> #include<algorithm> using na…

Java课题笔记~ ServletContext

单个Servlet的配置对象 web.xml <servlet><servlet-name>FirstServlet</servlet-name><servlet-class>com.ambow.test.FirstServlet</servlet-class><init-param><param-name>charset</param-name><param-value>utf-8&…

centos自动同步北京时间

1、安装ntpdate服务 yum -y install ntpdate 2、加入自动任务计划 查找ntpdate的路径&#xff1a; which ntpdate 复制这个路径。 编辑自动任务计划并加入ntpdate&#xff1a; crontab -e # 每小时第30分钟同步AD域控时间 30 * * * * /usr/sbin/ntpdate -u 192.168.2.8 > …

DP——动态规划

DP——动态规划 动态规划算法动态规划的一般步骤特殊DP——背包0-1背包问题完全背包问题 总结 动态规划算法 当涉及到解决具有重叠子问题的优化问题时&#xff0c;动态规划是一种常用的算法技术。它通过将问题分解为一系列重叠子问题&#xff0c;并使用递归或迭代的方式来解决…

Spring Cloud Gateway系例—GatewayFilter 工厂

目录 6.1.AddRequestHeader 6.2.AddRequestHeadersIfNotPresent 6.3.AddRequestParameter 6.4.AddResponseHeader 6.5.CircuitBreaker 6.5.1. 熔断指定的状态码 6.6.CacheRequestBody 6.7.DedupeResponseHeader 6.8.FallbackHeaders 6.9.JsonToGrpc 6.10.LocalRespo…

TypeScript 非空断言

TypeScript 非空断言 发布于 2020-04-08 15:20:15 17.5K0 举报 一、非空断言有啥用 介绍非空断言前&#xff0c;先来看个示例&#xff1a; function sayHello(name: string | undefined) {let sname: string name; // Error } 对于以上代码&#xff0c;TypeScript 编译器…

用户端Web自动化测试-L1

目录&#xff1a; Web自动化测试价值与体系环境安装与使用自动化用例录制自动化测试用例结构分析web浏览器控制常见控件定位方法强制等待与隐式等待常见控件交互方法自动化测试定位策略搜索功能自动化测试用户端Web自动化测试 1.Web自动化测试价值与体系 功能测试场景: UI 自…

IntelliJ Idea 编译时控制台上中文输出乱码

猜测原因是IDEA启动时未指定编码信息&#xff0c;故与系统编码保持一致&#xff08;windows中文系统默认为GBK编码&#xff09;,当以UTF-8编码进行编译在控制台会以GBK编码输出,从而导致乱码 解决方案 指定Idea启动时JVM的默认编码为UTF-8 Help -> Edit Custom Options P…

本地图片的image加密解密-Python 3.10-win10

本地图片的image加密解密- Python 3.10 pyt3int22 -根据1zip下图片批量生成加密的-物体识别.py import ioimport os import base64 import json # 指定图片文件夹 image_dir = "./1zip/" base64code_dir = "./base64code/" base64_to_dir = "./bas…