CGS与MGS的矩阵正交化-C语言实现

格拉姆-施密特正交化和改进的格拉姆-施密特正交化

格拉姆-施密特正交化CGS

数学公式

代码实现:

过程版

矩阵运算实现的难点在于每次运算都是一个向量,需要for循环进行,会带来运算时在代码中的复杂,进而难以理解代码的过程

Q矩阵是{e...}  R矩阵是上三角矩阵 
Q[0] = A[0];   // -> b1 = a1
R[0][0] = norm(Q[0], m); //-> ||b1||
Q[0] =A[0] / R[0][0];    //-> b1/||b1||

先求 a_{2}^{T}e_{1}      R[1][2] = dotProduct(A[i], Q[j], m);

再用 Q[2] = A[2]; //这样可以使用迭代去求b2

R[j][i] = dotProduct(A[i], Q[j], m); // -> a2Te1Q[i][k] = A[i][k];   //方便下一般实现循环递归减
Q[i][k] -=  R[j][i]*Q[j][k];   //用于递归实现 bn = an - anTen-1en-1-...R[i][i] = norm(Q[i], m);// -> ||b1||
Q[i][k] /= R[i][i];    //-> b1/||b1||

纯净版

void cgs(double** A, int m, int n, double** Q, double** R) {if (n>=1){R[0][0] = norm(A[0], m);for (int k = 0; k < m; k++) {Q[0][k] =A[0][k]/ R[0][0];}}for (int i = 1; i < n; i++) {for (int j = 0; j < i; j++) {R[j][i] = dotProduct(A[i], Q[j], m);}for (int k = 0; k < m; k++){Q[i][k] = A[i][k];}for (int j = 0; j < i; j++) {for (int k = 0; k < m; k++) {Q[i][k] -=  R[j][i]*Q[j][k];}}R[i][i] = norm(Q[i], m);for (int k = 0; k < m; k++) {Q[i][k] /= R[i][i];}}
}

改进的格拉姆-施密特正交化MGS

数学公式

在于先求e 每次都对全部b进行运算

代码实现

过程版

for (int i = 0; i < m; i++)

{

        Q[i]= A[i];

}

for (int i = 0; i < m; i++)
{for (int j = 0; j < n; j++){Q[i][j] = A[i][j];}
}

R[i][i] = norm(Q[i], m);
for (int k = 0; k < m; k++) {Q[i][k] /= R[i][i];
}

先求 b_{2}^{T}e_{1}    

 R[i][j] = dotProduct(Q[i], Q[j], m);

对每一个b仅需运算 每一个循环都会少计算b1,b2,b3..bn所以j = i+1

 for (int j = i + 1; j < n; j++) {for (int k = 0; k < m; k++) {Q[j][k] -= Q[i][k] * R[i][j];}}

纯净版

void mgs(double** A, int m, int n, double** Q, double** R) {for (int i = 0; i < m; i++){for (int j = 0; j < n; j++){Q[i][j] = A[i][j];}}for (int i = 0; i < n; i++) {R[i][i] = norm(Q[i], m);for (int k = 0; k < m; k++) {Q[i][k] /= R[i][i];}// 计算 R元素 b2Tei b3Tei...  Q[i]=e[i]for (int j = i + 1; j < n; j++) {R[i][j] = dotProduct(Q[i], Q[j], m);}// 更新 Qfor (int j = i + 1; j < n; j++) {for (int k = 0; k < m; k++) {Q[j][k] -= Q[i][k] * R[i][j];}}}
}

参考博客:

下文有博主有python实现

图解格拉姆-施密特正交化和改进的格拉姆-施密特正交化_matlab实现格拉姆-施密特正交化-CSDN博客

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

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

相关文章

软件测试——白盒测试

目录 学习视频来源&#xff1a;1.01什么是白盒测试(Av809778070,P1)_哔哩哔哩_bilibili 1.语句覆盖 2.判定覆盖 3.条件覆盖 4.判定条件覆盖 5.条件组合覆盖 6.路径覆盖 7.基本路径覆盖 学习视频来源&#xff1a;1.01什么是白盒测试(Av809778070,P1)_哔哩哔哩_bilibili …

正则表达式运用

已经写了表达式&#xff0c;下一步就是匹配字符串得到结果 使用matcher的源码&#xff08;匹配&#xff09;普通方法&#xff0c;find&#xff08;寻找&#xff09;合适的代码&#xff0c;看字符串是否匹配成功 是否可以匹配上 匹配么&#xff0c;匹配就留下&#xff0c;fin…

基于SVD的点云配准

前言 为了更加深入透彻的理解点云配准,所以打算从头实现一次配准的过程,从关键点提取、特征提取、对应关系确定、SVD确定变换矩阵,一步一步完成。 首先介绍一下配准过程中用到的一些函数 1. ISS3D关键点提取 具体来说,它使用ISS(Intrinsic Shape Signatures)关键…

Go 语言基础

1. 理解语句和标记 Go程序由一系列语句组成,而Go语句由一系列标记组成。这里,标记指的是语句中任何有意义的对象,包括如下。 关键字:用于特殊 Go 功能的保留字。操作符和标点符号:用于数学运算或执行特殊操作的符号。标识符:用于标识诸如存储位置或执行操作的函数等事物的…

热线系统自主可控改造:构建安全、高效的服务新纪元

一、引言&#xff1a;热线系统的重要性与挑战 在信息化时代&#xff0c;热线系统作为企业与客户沟通的关键渠道&#xff0c;承载着重要的信息传递和服务功能。然而&#xff0c;随着技术的不断进步和网络安全威胁的日益严峻&#xff0c;传统的热线系统面临着诸多挑战。这些挑战…

性能工具之 JMeter 常用组件介绍(三)

文章目录 一、常用组件介绍二、Sampler&#xff1a;取样器三、Controller:控制器&#xff08;逻辑控制器&#xff09;四、Pre Processor:预处理五、Post Processor:请求之后的处理六、Assertions:断言七、Timer:定时器八、Test Fragment&#xff1a;片段九、Config Element:配置…

STM32 SPI驱动读取LSM6DSRTR

提示&#xff1a;通过SPI驱动读取传感器数据 文章目录 前言一、LSM6DSRTR二、配置步骤1.配置SPI2.引入 LSM驱动库3.结果 总结 前言 制作一个倾角传感器&#xff0c;通过SPI读取LSM6DSRTR的加速度数据转换为角度&#xff0c;不用IIC的原因是考虑IIC通讯的协议过于繁琐&#xff…

springCloudAlibaba之服务熔断组件---sentinel

sentinel组件学习 sentinel学习sentinel容错机制使用代码方式进行QPS流控-流控规则初体验使用SentinelResource注解进行流控使用注解方式 sentinel学习 服务雪崩 服务雪崩效应&#xff1a;因服务提供者的不可用导致服务调用者的不可用&#xff0c;并将不可用逐渐放大的过程。 …

心链10----查询修改加入队伍业务实现

心链 — 伙伴匹配系统 接口设计 查询队伍列表 :::success 分页展示队伍列表&#xff0c;根据名称、最大人数等搜索队伍 P0&#xff0c;信息流中不展示已过期的队伍 从请求参数中取出队伍名称等查询条件&#xff0c;如果存在则作为查询条件不展示已过期的队伍&#xff08;根…

docker实战命令大全

文章目录 1 环境准备1.1 移除旧版本Docker1.2安装工具包1.3配置docker yum源 2 安装最新docker2.1 设置开机自启docker2.2配置加速器 3 实操-镜像3.1搜索镜像3.2下载镜像3.3查看镜像3.4 删除镜像 4 实操-容器4.1运行nginx容器4.2 查看容器4.3启动容器4.5关闭容器4.6查看容器日志…

功能强大的开源数据中台系统 DataCap 2024.03.5 发布

推荐一套基于 SpringBoot 开发的简单、易用的开源权限管理平台&#xff0c;建议下载使用: https://github.com/devlive-community/authx 推荐一套为 Java 开发人员提供方便易用的 SDK 来与目前提供服务的的 Open AI 进行交互组件&#xff1a;https://github.com/devlive-commun…

2个不同node版本的前端项目,用到的一个node版本管理工具nvm,可一键切换node版本

背景&#xff1a; 对于需要在电脑运行2个不同的前端项目&#xff0c;但是使用的node版本不同&#xff0c;为了快捷切换和管理node版本&#xff0c;所以需要用到一个非常方便的工具&#xff0c;那就是下边的nvm&#xff0c;一个命令就可以非常方便的切换你想要的node。 1、下载…

分布式光纤测温DTS与光纤光栅FBG解调仪有什么区别?

分布式光纤测温DTS和光纤光栅FBG解调仪之间存在本质区别。分布式光纤测温DTS是一种完全分布式的温度监测技术&#xff0c;意味着光纤的整个长度都充当传感器&#xff0c;可以感知任何位置的温度变化。然而&#xff0c;由于空间分辨率的限制&#xff0c;目前国内外的大多数技术只…

Three.js中的Raycasting技术:实现3D场景交互事件的Raycaster详解

前言 在Web开发中&#xff0c;Three.js是一个极为强大的库&#xff0c;它让开发者能够轻松地在浏览器中创建和展示3D图形。随着3D技术在网页设计、游戏开发、数据可视化等领域的广泛应用&#xff0c;用户与3D场景的交互变得日益重要。而要实现这种交互&#xff0c;一个核心的技…

探索地产营销管理系统业务架构的设计与应用

随着城市化进程的加速和人们对居住环境需求的不断提升&#xff0c;地产行业正迎来前所未有的发展机遇和挑战。在这个背景下&#xff0c;地产营销管理系统作为地产开发企业营销管理的重要工具&#xff0c;扮演着至关重要的角色。本文将深入探讨地产营销管理系统业务架构的设计与…

java中的面试高频问题----2

一、进程、线程、协程有什么区别&#xff1f; 1.进程&#xff1a;进程是操作系统中独立运行的程序实例&#xff0c;每个进程都有自己的内存空间和系统资源&#xff1b;进程之间相互独立&#xff0c;每个进程有自己的内存地址空间&#xff0c;一个进程无法直接访问另一个进程的…

H3C运维工程师面经

H3C运维工程师面经 1、常用linux命令:top的作用1.1.系统监控与性能分析:1.2.进程管理:1.3.资源使用统计:1.4.自定义视图与交互操作:2、接触过的linux系统、小版本是多少?2.1CentOS的版本CentOS 8系列:CentOS 7系列:CentOS 6系列:2.2Ubuntu版本Ubuntu LTS版本:Ubuntu非…

手机怎么压缩图片?通过三种压缩操作

手机怎么压缩图片&#xff1f;在智能手机日益普及的今天&#xff0c;拍照分享已成为日常生活的一部分。然而&#xff0c;高质量的照片往往占用较大的存储空间&#xff0c;且在网络上传输时速度较慢。那么&#xff0c;如何在手机上压缩图片呢&#xff1f;本文将介绍三种实用的手…

文件无法在当前环境下执行在 x86_64 系统上运行 ARM 可执行文件

目录 遇到的问题是由于"..."文件无法在当前环境下执行。这个错误通常是因为二进制文件的格式不兼容&#xff0c;可能是因为它是为不同的架构编译的。例如&#xff0c;如果二进制文件是为 x86 架构编译的&#xff0c;但你在 ARM 设备上尝试运行它&#xff0c;就会出现…

Oracle创建索引的LOGGING | NOLOGGING区别

在Oracle中&#xff0c;创建索引时的LOGGING和NOLOGGING选项主要影响索引创建过程中产生的重做日志&#xff08;redo log&#xff09;的数量。这两个选项对于性能和数据恢复能力有着显著的影响。以下是关于这两个选项的详细解释和区别&#xff1a; LOGGING 定义&#xff1a;当…