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

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

安德烈·路易斯·乔尔斯基出生于法国波尔多以北的查伦特斯海域的蒙古扬。他在波尔多参加了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.mzph.cn/news/695431.shtml

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

相关文章

量化粒度是什么???

目录 一、什么是量化粒度二、量化粒度方式三、总结 一、什么是量化粒度 量化粒度是指在模型量化中&#xff0c;对模型参数进行量化的级别。常见的量化粒度包括&#xff1a; Per-tensor 量化: 对每个 tensor 单独进行量化。 Per-layer 量化: 对每个 layer 的所有 tensor 共享…

解决采集时使用selenium被屏蔽的办法

解决采集时使用selenium被屏蔽的办法 实用seleniumbase uc模式 from seleniumbase import Driver driver Driver(ucTrue) # 使用UC模式UC模式是基于undetected-chromedriver 但做了一些优化更新&#xff0c;使用起来更方便 官方例子&#xff1a; from seleniumbase import …

关于音频编码标准AAC,Opus,MP3的概念、原理、优缺点

音频编码标准是用于压缩数字音频数据的技术规范&#xff0c;以减少存储和传输所需的带宽。以下是关于AAC、Opus和MP3这些音频编码标准的概念、原理、优缺点的简要描述&#xff1a; AAC (Advanced Audio Coding) 概念&#xff1a; AAC是一种音频编码标准&#xff0c;由MPEG-2和…

删除数组用delete和Vue.delete有什么区别

删除数组用delete和Vue.delete有什么区别&#xff1f; 在 JavaScript 中&#xff0c;delete 和 Vue.js 中的 Vue.delete 是两个完全不同的概念&#xff0c;它们在删除数组元素时的作用和效果也有所不同。 JavaScript 中的 delete 关键字&#xff1a; 在原生 JavaScript 中&a…

echarts实现表格图例

折线图的表格图例 <!--折线图--> <template><div class"title-container"><span class"title">{{typeof props.series string ? props.series : props.series.name}}</span></div><div v-if"props.series?…

PhotoSweeper X mac版 v4.8.5 相似重复照片清理工具 兼容 M1/M2

PhotoSweeper X for Mac是一款Mac重复照片/相似照片清理工具&#xff01;PhotoSweeper可以帮你进行&#xff1a;重复相似照片/数码相片查找、对比和删除&#xff0c;轻松清理Mac上的重复图片&#xff0c;非常实用。 应用介绍 PhotoSweeper X for Mac是一款Mac重复照片/相似照片…

实战打靶集锦-025-HackInOS

文章目录 1. 主机发现2. 端口扫描3. 服务枚举4. 服务探查5. 提权5.1 枚举系统信息5.2 探索一下passwd5.3 枚举可执行文件5.4 查看capabilities位5.5 目录探索5.6 枚举定时任务5.7 Linpeas提权 靶机地址&#xff1a;https://download.vulnhub.com/hackinos/HackInOS.ova 1. 主机…

机器学习:训练模型后,发生了什么,最后得到什么

当你训练一个模型&#xff0c;尤其是在机器学习或深度学习的上下文中&#xff0c;实际上是在进行一系列计算和优化步骤&#xff0c;以便模型能够从提供的数据中学习。训练模型的过程大致如下&#xff1a; 1. 初始化模型参数 模型的参数&#xff08;如权重和偏置&#xff09;通…

【图片公式识别】图片公式转Word与LaTeX文档:智能识别与转换

前言 嘿&#xff0c;大家好呀&#xff01;&#x1f44b; 谁都知道&#xff0c;写 Word 文档里的公式可不是一件简单的事情&#xff01;你辛辛苦苦在键盘上敲出的数学公式&#xff0c;结果随着 Word 版本的更新&#xff0c;竟然变成了一张图片&#xff01;&#x1f624; 这简直就…

toFixed四舍五入

一、问题描述 3.145.toFixed(2) // 3.15 3.155.toFixed(2) // 3.15 3.1551.toFixed(2) // 3.16 3.1550000000000001.toFixed(2) // 3.16 3.15500000000000001.toFixed(2) // 3.15 二、原因分析&#xff1a…

板块一 Servlet编程:第八节 文件上传下载操作 来自【汤米尼克的JavaEE全套教程专栏】

板块一 Servlet编程&#xff1a;第八节 文件的上传下载操作 一、文件上传&#xff08;1&#xff09;前端内容&#xff08;2&#xff09;后端内容 二、文件下载&#xff08;1&#xff09;前端的超链接下载&#xff08;2&#xff09;后端下载 在之前的内容中我们终于结束了Servle…

Oracle 总结

Oracle常规命令 conn 用户/密码--客户端连接的命令 show user —显示当前登录的用户增删改查 增 insect into 表名(列名列表) values(值列表);删 删除符合条件数据 delete from EMP where ENAME ‘xxx’;删除所有EMP 所有数据 delete from EMP;截断表 truncate table EMP;注…

速盾网络:CDN节点无法读取您网站数据是什么意思

在使用CDN服务时&#xff0c;有时候可能会遇到提示&#xff1a;“CDN节点无法读取您网站数据”的情况&#xff0c;这意味着什么呢&#xff1f;让我们一起来解读速盾网络关于这个问题的解释。 意思解释&#xff1a; 当CDN节点无法读取您网站数据时&#xff0c;通常表示CDN节点…

Web应用程序的自动测试工具WebDriver简介

WebDriver 是一个开源的自动化测试工具&#xff0c;它提供了一套用于Web应用程序的自动化测试接口。通过这些接口&#xff0c;我们可以编写脚本来模拟用户的各种操作&#xff0c;如点击、拖拽、输入等&#xff0c;从而实现Web应用程序的自动化测试。 WebDriver 支持多种…

qtday2作业

思维导图 使用手动连接&#xff0c;将登录框中的取消按钮使用qt4版本的连接到自定义的槽函数中&#xff0c;在自定义的槽函数中调用关闭函数 将登录按钮使用qt5版本的连接到自定义的槽函数中&#xff0c;在槽函数中判断ui界面上输入的账号是否为"admin"&#xff0c;…

Linux搭建FISCO BCOS的第一个区块链网络

一、前言 FISCO BCOS是由金融区块链合作联盟&#xff08;深圳&#xff09;与微众银行共同发起的开源区块链项目&#xff0c;支持多链多账本&#xff0c;满足金融行业复杂业务需求。本文将介绍如何在Ubuntu操作系统上使用Linux命令搭建FISCO BCOS的第一个区块链网络。 目录 一…

K8S部署Java项目(Springboot项目)pod状态:CrashLoopBackOff

天行健&#xff0c;君子以自强不息&#xff1b;地势坤&#xff0c;君子以厚德载物。 每个人都有惰性&#xff0c;但不断学习是好好生活的根本&#xff0c;共勉&#xff01; 文章均为学习整理笔记&#xff0c;分享记录为主&#xff0c;如有错误请指正&#xff0c;共同学习进步。…

昇腾ACL应用开发之模型转换ATC

一.前提条件 在前面的章节中我们已经安装了包含模型转换的套件包CANN-TOOLKIT&#xff0c;默认的安装路径会在/usr/local/Ascend里面&#xff0c;我们将该套件所需要的东西加入到环境变量中以便我们调用&#xff1a; 将source /usr/local/Ascend/ascend-toolkit/set_env.sh加入…

人为物累,心为形役

一、人是什么 你是你&#xff0c;他是他&#xff0c;我是我&#xff0c;有什么区别吗&#xff0c;直到自我发现我与你不同时&#xff0c;不同是什么&#xff0c;身体结构&#xff1f;人生经历&#xff1f;所拥有的一切&#xff1f;暂时搁置这些的话&#xff0c;抽离我们的意识…