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,一经查实,立即删除!

相关文章

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. 主机…

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

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

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

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

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;抽离我们的意识…

[ai笔记11] 论ai韭菜的自我修养

欢迎来到文思源想的ai空间&#xff0c;这是技术老兵学习ai以及观点分享的第11篇内容&#xff01; 上班之后时间确实少了许多&#xff0c;但是最近也没闲着&#xff0c;关于ai的学习一直在探索两个部分&#xff0c;一个是看那本有名的书《这就是ChatGPT》&#xff0c;另外一个则…

Gmail邮箱群发邮件的技巧?邮箱怎么绑定?

Gmail邮箱注册教程指南&#xff1f;如何注册新的Gmail邮箱帐户&#xff1f; Gmail邮箱作为谷歌推出的邮件服务&#xff0c;以其高效、稳定和便捷的特性受到广大用户的喜爱。然而&#xff0c;如何在Gmail中进行有效的群发邮件&#xff0c;接下来&#xff0c;蜂邮将给大家介绍一…

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

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

【深入理解设计模式】 工厂设计模式

工厂设计模式 工厂设计模式是一种创建型设计模式&#xff0c;它提供了一种在不指定具体类的情况下创建对象的接口。在工厂设计模式中&#xff0c;我们定义一个创建对象的接口&#xff0c;让子类决定实例化哪一个类。工厂方法使一个类的实例化延迟到其子类。 工厂设计模式的目…

part2 思维导图(微分部分)

1. 导数和微分 2. 一元微分学 2.1 四则运算 2.2 口诀 2.3 导数计算 2.3.1 计算1 2.3.1 计算2

win32汇编定时器继续学习2

前面学习了win32汇编定时器&#xff0c;还非常不熟悉&#xff0c;继续熟悉&#xff0c;稍微增加一点功能&#xff1b; .386.model flat,stdcalloption casemap:noneinclude windows.inc include user32.inc includelib user32.lib include kernel32.inc includelib kernel32…

2024生物发酵展创新盛会-兰格恒流泵

​ 参展企业介绍 兰格恒流泵有限公司成立于1997年&#xff0c;是英国豪迈国际有限公司&#xff08;Halma&#xff09;的全资子公司。我们始终坚持以客户需求为导向&#xff0c;提供安全、高精度及高品质的流体传输设备和解决方案。兰格产品系列超…

前端架构: 脚手架之Chalk和Chalk-CLI使用教程

Chalk Chalk 是粉笔的意思, 它想表达的是&#xff0c;给我们的命令行中的文本添加颜色类似彩色粉笔的功能 在官方文档当中&#xff0c;它的 Highlights 核心特性 Expressive API Highly performant No dependencies Ability to nest styles 256/Truecolor color support Auto-…

提高工作效率,体验ONLYOFFICE办公软件

ONLYOFFICE办公软件 一、前言二、特点完整办公套件协作过程更容易 三、访问地址 一、前言 随着数字化办公的普及&#xff0c;办公软件在我们的工作中扮演着越来越重要的角色。为了提高工作效率&#xff0c;我们需要一个功能强大、易于使用的办公软件。ONLYOFFICE作为一款全功能…

分布式系统一致性与共识算法

分布式系统的一致性是指从系统外部读取系统内部的数据时&#xff0c;在一定约束条件下相同&#xff0c;即数据&#xff08;元数据&#xff0c;日志数据等等&#xff09;变动在系统内部各节点应该是一致的。 一致性模型分为如下几种&#xff1a; ① 强一致性 所有用户在任意时…