C#,数值计算,矩阵的乔莱斯基分解(Cholesky decomposition)算法与源代码

news2024/10/7 8:23:29

一、安德烈·路易斯·乔尔斯基

安德烈·路易斯·乔尔斯基出生于法国波尔多以北的查伦特斯海域的蒙古扬。他在波尔多参加了Lycée e,并于1892年11月14日获得学士学位的第一部分,于1893年7月24日获得第二部分。1895年10月15日,乔尔斯基进入莱科尔理工学院,在当年223名入学学生中排名第88位。他在莱科尔理工学院的教授包括卡米尔·乔丹和发现放射性的著名物理学家亨利·贝克勒尔。在成功的两年后,他于1897年参加了莱科尔理工学院的期末考试。在222名学生中,他提高了自己的地位,在这些考试中排名第38位。随后,他加入军队,成为少尉,并从1897年10月开始在炮兵学校学习。他在1899年完成了学业,并保持了稳步的进步,因为现在他在那一年获得资格的86名学生中排名第五。

突尼斯当时仍然是奥斯曼帝国的正式一部分,但在1883年突尼斯和法国政府签署的《马尔萨公约》之后,它一直是法国的保护国。法国人稳定了经济,建立了现代通讯。从1902年1月到6月,乔尔斯基在突尼斯执行了一项任务,然后从1902年11月到1903年5月,他被派去执行第二项任务。1903年12月31日,乔尔斯基开始在阿尔及利亚服役。那个国家从法国得到的待遇比突尼斯少得多。法国人建立了对该国及其原住民的统治,但发展了新的通讯、医院和医疗服务。他于1904年6月6日离开阿尔及利亚。

乔尔斯基职业生涯中最重要的一步是他于1905年6月进入陆军地理服务部测地部门,据报道他在那里:-

... 他有敏锐的智力和出色的数学能力,有探究精神和独创性的想法。

德拉姆布雷在1798年春天完成了基线测量,这是他对定义米的贡献的一部分。1882年,法国专家回到了德拉姆布雷的基线,但没有重新测量它,而是倾向于通过三角测量间接检查德拉姆布雷的计算。20世纪第一个十年,巴黎子午线修订后,计划对法国进行新的三角测量。调整网格的问题让地理服务部门的官员非常担忧,他们迫切希望找到一种简单、快速、准确的方法。为了用最小二乘法求解条件方程,Cholesky发明了一个非常巧妙的计算程序,该程序立即被证明非常有用:它现在被称为Cholesky方法,我们将在下面描述它。1905年9月26日,他成为中尉,两年后于1907年5月10日结婚。乔尔斯基和他的妻子有三个孩子;一个儿子和两个女儿。

克里特岛是奥斯曼帝国的一部分已经有200年了,但随着帝国的衰落,欧洲的主要国家打算参与克里特岛的未来。1896年,克里特岛发生了一场反对土耳其人的革命,之后希腊试图取得控制权,列强强行达成了和解。法国人负责西蒂亚区,意大利人负责伊拉佩特拉区,英国人负责伊拉克利安区,俄罗斯人负责雷瑟姆农区,而四人共同负责查尼亚区。法国军队在克里特岛的司令部鲁班斯基上校曾是测地学家,希望在岛上进行测绘工作。1906年3月和4月,测地科长布尔乔亚中校进行了初步侦察。乔尔斯基于1907年11月7日被派往克里特岛,作为执行任务的三名军官之一。他在卡沃西平原测量了基线,并使用天文测量来确定基线南端的精确位置。在他的同事被派去执行其他任务后,乔尔斯基进行了三角测量,并对法国和英国的区域进行了地形测量。他在极其困难的条件下成功地进行了三角测量,因为克里特岛冬季高山上降雪量很大。然而,政治事件阻止了土地调查的进行,乔莱斯基于1908年6月25日离开克里特岛。

1909年3月25日,他被任命为第二步兵团团长,并于1909年9月至1911年9月随部队执行任务。1911年晚些时候,他回到了陆军地理服务局的测地部门,并被分配到阿尔及利亚和突尼斯工作,七年前他已经在那里积累了经验。从1911年10月27日至1912年4月24日,再从1912年10月23日至1913年4月17日,他在阿尔及利亚和突尼斯进行了精密水准测量工作。这项工作旨在为修建铁路线做准备。突尼斯的电网建设于1913-1914年冬季完成,之后有一段时间进行检查和调整。在摩洛哥也必须进行类似的工作,Cholesky设计了一些方法,使工作能够更快地进行,但仍然保持必要的准确性。在测试了他的程序后,这些程序于1912年7月开始在摩洛哥使用,工作于1913年1月完成。

1913年5月25日,乔莱斯基被分配到外交部,并被任命为突尼斯摄政区地形局局长。他开始在那里工作,但在1914年秋天第一次世界大战爆发后,他被调到其他地方工作。1915年1月,他开始组织火炮射击,2月被派往沃斯的一个部队工作。1916年9月,他被派往罗马尼亚执行军事任务,在那里担任罗马尼亚军队地理服务处处长。他于1917年7月晋升,并继续在罗马尼亚服役至1918年2月。

乔莱斯基于1918年8月31日凌晨5点在法国北部的战场上受伤身亡。在他死后,他的一位同事,贝诺伊特司令,在《自然推断法中线性方程组解析法的起源》一书中发表了乔尔斯基的计算方法,以解决一些最小二乘数据拟合问题。应用méthodeála resolution d'un Systemème Definition d'quations lineaires(首席检察官Cholesky)Ⓣ, 发表于1924年的公报géodesique。

Cholesky分解(或Cholesky分解)采用对称正定矩阵AA,并将其写成a=LL'a=LL′

其中,LL是一个下三角矩阵,具有正对角项(有时称为Cholesky三角形),L'L′

是LL的转置。要解Ax=bsolveAx=b,现在需要解LL'x=bLL′

x=b所以把y=L'xy=L′

x表示Ly=bLy=b,解为yy,然后y=L'xy=L′

对x进行求解,得到解。该方法的优点在于,当MM是三角形矩阵时,求解Mx=bMx=b型方程非常简单。

该方法在1924年出版后很少受到关注,但杰克·托德在二战期间将其纳入了伦敦国王学院的分析课程。1948年,福克斯、赫斯基和威尔金森在一篇论文中对该方法进行了分析,而图灵在同一年发表了一篇关于该方法稳定性的论文。

二、斯基分解(Cholesky decomposition)

矩阵分解 (decomposition, factorization)是将矩阵拆解为数个矩阵的乘积,可分为三角分解、满秩分解、QR分解、Jordan分解和SVD(奇异值)分解等,常见的有三种:
(1)三角分解法 (Triangular Factorization)
(2)QR 分解法 (QR Factorization)
(3)奇异值分解法 (SVD,Singular Value Decomposition)
 

三角分解法亦称因子分解法,由消元法演变而来的解线性方程组的一类方法。设方程组的矩阵形式为Ax=b,三角分解法就是将系数矩阵A分解为一个下三角矩阵L和一个上三角矩阵U之积:A=LU,然后依次解两个三角形方程组Ly=b和Ux=y,而得到原方程组的解,例如,杜利特尔分解法、乔莱斯基分解法等就是三角分解法。

乔莱斯基分解法(Cholesky decomposition method),亦称平方根法,属于三角分解法之一。

楚列斯基分解(Cholesky decomposition)是1993年公布的数学名词。

乔莱斯基乔里斯基分解,同。

using System;
using System.Collections;
using System.Collections.Generic;

namespace Legalsoft.Truffer.Algorithm
{
    /// <summary>
    /// 乔莱斯基分解
    /// </summary>
    public static partial class Algorithm_Gallery
    {
        public static int[,] Cholesky_Decomposition(int[,] matrix)
        {
            int n = matrix.GetLength(0);
            int[,] lower = new int[n, n];

            for (int i = 0; i < n; i++)
            {
                for (int j = 0; j <= i; j++)
                {
                    int sum = 0;
                    if (j == i)
                    {
                        for (int k = 0; k < j; k++)
                        {
                            sum += (int)Math.Pow(lower[j, k], 2);
                        }
                        lower[j, j] = (int)Math.Sqrt(matrix[j, j] - sum);
                    }
                    else
                    {
                        for (int k = 0; k < j; k++)
                        {
                            sum += (lower[i, k] * lower[j, k]);
                        }
                        lower[i, j] = (matrix[i, j] - sum) / lower[j, j];
                    }
                }
            }

            return lower;
        }
    }
}

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

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

相关文章

Java Web(六)--XML

介绍 官网&#xff1a;XML 教程 为什么需要&#xff1a; 需求 1 : 两个程序间进行数据通信&#xff1f;需求 2 : 给一台服务器&#xff0c;做一个配置文件&#xff0c;当服务器程序启动时&#xff0c;去读取它应当监听的端口号、还有连接数据库的用户名和密码。spring 中的…

git使用过的命令记录

目录 git add .git commit --amendgit push -f origin HEAD:mastergit checkout .git stash想把某个pr的修改应用到本地git 撤销&#xff0c;放弃本地修改参考文档 git add . 将本地修改提交到暂存区 git commit --amend 将修改的东西 追加到上一个pr git push -f origin H…

SQL数据库基础语法-查询语句

SQL数据库基础语法-查询语句 Group By #对数据进行分组 >select name,count(id) from student group by name; #查询name字段人数&#xff0c;cont函数进行计数 >select * from users group by users; >select * from users where id1 group by 2; >select * from …

线阵相机参数介绍之轴编码器控制

1.1 功能介绍 编码器是将检测对象的运动与相机拍摄取图相匹配的设备&#xff0c;也即检测对象运动一定距离&#xff0c;相机就拍摄一定行高的图像。 编码器会将检测对象的实际位移转换为固定数量电信号。例如&#xff1a;编码器的精度是2000p/r,该参数的含义是编码器每转一圈输…

webpack打包速度优化思维导图

webpack打包速度优化思维导图 前言附件 前言 去年的时候公司一个项目体积过大&#xff0c;我是m1芯片的macpro&#xff0c;光启动就要1分钟&#xff0c;配置差点都电脑&#xff0c;启动就要3分钟&#xff0c;自然打包速度也会慢很多&#xff0c;我们是gitlab设置成了自动打包的…

数据分析 — 电商用户分析和用户 RFM 模型

目录 一、电商用户分析1、数据字段信息2、数据读取3、数据清洗4、可视化分析1、每年销售额的增长情况2、各个地区分店的销售额3、每个分店每一年的销售额4、销售淡旺季5、新增用户 二、RFM 模型1、RFM 模型的三个维度2、RFM 的客户类型标签3、RFM 模型的二分法思想4、代码 一、…

【算法】动态规划1,最小花费爬楼梯,解码方法

一、动态规划简介 动态规划 , 英文名称 Dynamic Programming , 简称 DP , 不是具体的某种算法 , 是一种算法思想 ; 动态规划 , 没有具体的步骤 , 只有一个核心思想 ; 动态规划 的 核心思想 是 由大化小 , 大规模问题 使用 小规模问题 计算结果 解决 , 类似于 分治算法 ; 二、…

ipad作为扩展屏的最简单方式(仅需在同一局域网下,无需数据线)

ipad和win都下载安装toDesk&#xff0c;并且都处于同一局域网下 连接ipad&#xff0c;在ipad中输入win设备的设备密码和临时密码&#xff0c;连接上后可以看到ipad会是win屏幕的镜像&#xff0c;此时退出连接&#xff0c;准备以扩展模式再次连接。 注意&#xff0c;如果直接从…

SOPHON算能科技新版SDK环境配置以及C++ demo使用过程

目录 1 SDK大包下载 2 获取SDK中的库文件和头文件 2.1 注意事项 2.2 交叉编译环境搭建 2.2.1 首先安装工具链 2.2.2 解压sophon-img包里的libsophon_soc__aarch64.tar.gz&#xff0c;将lib和include的所有内容拷贝到soc-sdk文件夹 2.2.3 解压sophon-mw包里的sophon-mw-s…

pyqt如何实现拖拽打开文件(通过windows的快捷方式打开文件)

桌面端的开发中如何通过windows的快捷方式打开文件&#xff0c;那么如何将需要打开的数据传递给qt程序呢&#xff1f; 研究了一下发现很简单 通过sys.argv可以轻松的实现传参 sys.argv import sys print(sys.argv)这个方法可以获取系统传递给程序的参数&#xff0c;默认是个列…

Android 12.0 systemui锁屏页面时钟显示样式的定制功能实现

1.前言 在12.0的系统ROM定制化开发中,在进行systemui的相关开发中,当开机完成后在锁屏页面就会显示时间日期的功能,由于 开发产品的需求要求时间显示周几上午下午接下来就需要对锁屏显示时间日期的相关布局进行分析,然后实现相关功能 效果图如图: 2.systemui锁屏页面时钟显…

Linux搭建JavaEE环境

Linux搭建JavaEE环境 需要在 Linux 下进行 JavaEE 的开发&#xff0c;我们需要安装如下软件 JDKtomcatmysqlidea 安装JDK 安装步骤 首先去官网下载jdk的linux版本 下载地址&#xff1a;https://www.oracle.com/cn/java/technologies/downloads/ 1) mkdir /opt/jdk 2) 通过 …

视频推拉流EasyDSS视频直播点播平台授权出现激活码无效并报错400是什么原因?

视频推拉流EasyDSS视频直播点播平台集视频直播、点播、转码、管理、录像、检索、时移回看等功能于一体&#xff0c;可提供音视频采集、视频推拉流、播放H.265编码视频、存储、分发等视频能力服务&#xff0c;在应用场景上&#xff0c;平台可以运用在互联网教育、在线课堂、游戏…

向导式堆栈管理器Dockge

经过申诉&#xff0c;目前博客的几个域名都恢复了&#xff0c;时间也延长到了 2033 年&#xff0c;后面还会不会出问题&#xff0c;老苏就不知道了 什么是 Dockge ? Dockge 是一款时髦的、易于使用的、响应式的、自托管的 docker-compose.yaml 向导式堆栈管理器&#xff0c;可…

flink反压

flink反压&#xff08;backpressure&#xff09;&#xff0c;简单来说就是当接收方的接收速率低于发送方的发送速率&#xff0c;这时如果不做处理就会导致接收方的数据积压越来越多直到内存溢出&#xff0c;所以此时需要一个机制来根据接收方的状态反过来限制发送方的发送速率&…

10、内网安全-横向移动域控提权NetLogonADCSPACKDC永恒之蓝

用途&#xff1a;个人学习笔记&#xff0c;有所借鉴&#xff0c;欢迎指正&#xff01; 背景&#xff1a; 主要针对内网主机中的 域控提权漏洞&#xff0c;包含漏洞探针和漏洞复现利用。 1、横向移动-系统漏洞-CVE-2017-0146&#xff08;ms17-010&#xff0c;永恒之蓝&#xff0…

ts环境下如何解决第三方npm包类型报错的问题

在 TypeScript 环境下&#xff0c;当引入第三方 npm 包时出现类型报错的问题&#xff0c;可以尝试以下几种解决方法&#xff1a; 1、安装 types 包&#xff1a;许多流行的第三方 npm 包都有对应的 TypeScript 类型声明文件&#xff0c;这些声明文件通常以 types/包名 的形式发…

odoo16-API(Controller)带有验证访问的接口

odoo16-API&#xff08;Controller&#xff09;带有验证访问的接口 目前我使用odoo原生的登录token来验证登陆的有效性 废话不多说直接上代码 # 测试获取session_id import requests class GetOdooData(http.Controller):def getOdooToken(self):# http://localhost:8123访问…

蓝桥杯:C++素数、进制转换

素数 定义&#xff1a; 质数又称素数。一个大于1的自然数&#xff0c;除了1和它自身外&#xff0c;不能被其他自然数整除的数叫做质数&#xff1b;否则称为合数&#xff08;规定1既不是质数也不是合数&#xff09;。 代码&#xff1a; #include<bits/stdc.h> using n…

【flutter】环境安装

安装flutter sdk 下载sdk flutter sdk就包含dart&#xff0c;所以我们只用安装flutter sdk就可以了。 我们去清华大学开源软件镜像站下载&#xff0c;flutter开发中&#xff0c;版本对不上基本项目就跑步起来&#xff0c;如果是团队协同开发的话&#xff0c;建议统一下载指定版…