P4783 【模板】矩阵求逆 高斯-约旦消元 求矩阵的逆

最新推荐文章于 2024-02-01 21:43:50 发布
最新推荐文章于 2024-02-01 21:43:50 发布
阅读量380 收藏
点赞数
分类专栏： # 数论——高斯消元/矩阵
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/bjfu170203101/article/details/108203034
版权
以上是我们算法的依据下面的是例题。我们用高斯-约旦消元把分块矩阵： A|E 中的A初等变化为E , 整体矩阵变为：输入右块矩阵即可。#include <bits/stdc++.h>using namespace std;typedef long long ll;const int M = 400+7;const int mod =1e9+7;int a[M][M<<1];ll qpow(ll a,ll b){ ll ans=1; wh...
摘要由CSDN通过智能技术生成

以上是我们算法的依据

下面的是例题。

我们用高斯-约旦消元把分块矩阵： A|E  中的A初等变化为E ,  整体矩阵变为： E|A^ {-1}

输入右块矩阵即可。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int M = 400+7;
const int mod =1e9+7;
int a[M][M<<1];
ll qpow(ll a,ll b)
{
	ll ans=1;
	while(b)
	{
		if(b&1)ans=ans*a%mod;
		a=a*a%mod;
		b>>=1;
	}
	return ans;
}
//返回自由元个数，及无穷解个数 
bool Gauss(int n,int m)//对矩阵a高斯消元得到(E|C) C为等式右边的列矩阵 
{
	int row;//当前处理的行
	int col;//当前处理的列 
  	for(row=1,col=1;row<=n&&co
最低0.47元/天 解锁文章
夕林山寸
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
选主元的高斯-约旦（Gauss-Jordan）消元法解线性方程组/求逆矩阵
weixin_30532973的博客
01-05 802
做数据结构课设时候查的资料，主要是看求逆矩阵方面的知识的。 选主元的高斯-约当（Gauss-Jordan）消元法在很多地方都会用到，例如求一个矩阵的逆矩阵、解线性方程组（插一句：LM算法求解的一个步骤），等等。它的速度不是最快的，但是它非常稳定（来自网上的定义：一个计算方法，如果在使用此方法的计算过程中，舍入误差得到控制，对计算结果影响较小，称此方法为数值稳定的），同时它的求解过程也比较清晰...
（高斯-约旦法）矩阵求逆算法
11-27
（高斯-约旦法）矩阵求逆算法，比较适合编程
高斯约旦消元法求n阶矩阵的逆
weixin_45081640的博客
11-10 633
高斯约旦消元法求n阶矩阵的逆 本代码以四阶矩阵为例，需要求更高次的，只需要将代码中的4修改下即可 // gaussJordan.cpp : 此文件包含 "main" 函数。程序执行将在此处开始并结束。 // #include <iostream> #include <opencv2/opencv.hpp> #include <opencv2/imgproc/imgproc.hpp> #include <opencv2/highgui/highgui.hpp>
矩阵求逆 —— 初等变换法（高斯-约旦消元）
热门推荐
一只萌新
08-21 1万+
P4783 【模板】矩阵求逆 题目描述 求一个N×NN×NN×N的矩阵的逆矩阵。答案对109+710^9+7109+7取模。 1.逆矩阵的定义 假设 AAA 是一个方阵，如果存在一个矩阵 A−1A^{-1}A−1，使得 A−1A=IA^{-1}A=IA−1A=I 并且 AA−1=IAA^{-1}=IAA−1=I 那么，矩阵 A 就是可逆的，A−1A^{-1}A−1 称为 A 的逆矩阵 2.逆矩阵求...
高斯消元法，高斯约旦消元法，线性方程组，矩阵求逆
wanglianda2007的博客
03-22 366
详解高斯消元法，高斯约旦消元法。求解线性方程组，矩阵求逆。
Java求逆矩阵（高斯约旦法）
04-02
本篇文章将深入探讨如何使用高斯约旦法在Java中求解逆矩阵。 高斯约ordan法是一种基于行变换将矩阵化简为阶梯形矩阵，进而得到单位矩阵的方法。逆矩阵的求解可以通过扩展矩阵的方式来实现，即将原始矩阵A与单位矩阵...
全选主元高斯-约旦法矩阵求逆
03-01
全选主元高斯约旦矩阵求逆算法，matlab语言编写，程序中对比了算法与matlab函数求逆结果，结果一致。
Gauss–Jordan elimination:用高斯-约旦消除法求矩阵求逆-开源
04-29
高斯-约旦消元法（Gauss-Jordan Elimination）是线性代数中解决线性方程组的一种重要方法，它同时也被用来求解矩阵的逆。这个方法是基于矩阵运算，结合了高斯消元法和矩阵的行变换，通过一系列的行操作将系数矩阵...
C#矩阵求逆计算-自己封装函数
最新发布
03-24
计算逆矩阵可以使用高斯-约旦消元法或者伴随矩阵法。这里我们将重点介绍伴随矩阵法，因为这种方法更适合编程实现。 1. **伴随矩阵法**： - 首先，构造A的伴随矩阵Adjugate(A)，它的(i, j)元素是A的余子矩阵的行列...
Gauss-Jordan法求逆矩阵
qwezhaohaihong
08-10 1万+
转载来源于： http://student.zjzk.cn/course_ware/web-gcsx/gcsx/chapter3/chapter3.2.htm http://student.zjzk.cn/course_ware/web-gcsx/gcsx/chapter1/chapter1.2.htm#21 先回顾一下高斯消元法： §1.2 消元法与矩阵的初等
高斯约化法求矩阵逆
05-25
这是本人结合测量学高斯约化法求法方程的解用C/C++编写的。
4.6 高斯约当消元法
jhshanvip的博客
04-08 1124
4.6 高斯约当消元法 高斯消元法把矩阵变换为上三角阵，上三角阵还可以继续变换为对角阵。例如上面增广矩阵 [A,b][A, b][A,b] 变换为上三角阵 [24−2201140048] \left[ \begin{matrix} 2 & 4 & -2 & 2\\ 0 & 1 & 1 & 4\\ 0 & 0 & 4 &...
高斯-约旦消元法 理解
weixin_30583563的博客
10-30 343
高斯消元是一种解方程的很巧妙的方法，核心是把方程转换成矩阵形式，然后再通过加减消元，求出值后再回带，就解出了这个方程，这里我就不赘述了。 我一般用高斯-约旦消元法，这种方法是直接转换成单位矩阵求解，减少回带次数，提高精确度，实现方式如下： 下方是一个方程 把它转换成矩阵形式就是: 我们可以这样对其进行变换: 在第一行找到第一个元素，向下寻找有没有绝对值更大的，如果有就交换。...
待定系数法、高斯-约旦消元法、LU分解法及Cholesky分解法实现矩阵求逆
zy4213的博客
02-01 1762
文章主要介绍了常见的矩阵求逆的方法并给出了具体实例及计算过程，主要包含：待定系数法、伴随矩阵求逆法、高斯-约旦消元法、LU分解法和Cholesky分解法。
（Gauss-Jordan）高斯消元法求逆矩阵（含C/C++实现代码）
小谢的博客
05-28 5887
文章目录一、问题引出二、 一、问题引出 给定一个 n×nn \times nn×n 的矩阵 AAA ，我们想求得一个矩阵 BBB 使得 ∣A×B∣|A \times B|∣A×B∣ 即 AAA 矩阵和 BBB 矩阵的积矩阵的行列式为 111 ，那么这个 BBB 矩阵就是 AAA 矩阵的逆矩阵，或者说 A×BA \times BA×B 得到单位矩阵 EEE 二、 ...
高斯约当法矩阵求逆（Matlab实现）
weiming0的博客
03-23 2462
数学模型到具体算法仍有一段距离，需要考虑计算机执行速度和“四舍五入”等问题，这里主要利用矩阵的初等变换进行矩阵求逆，也就是高斯约当法，代码如下： 主函数：main.m %测试 A=[1 -2 3; 2 3 4; 3 5 7]; b=[-10;5;7]; %内置函数 x1=A\b %高斯约当法 x=inverse(A,b) 函数：reverse.m %求可逆矩阵（AX=b） %高斯约当法 function x=inverse(A,b) row=size(A,1); %C用于存储A的逆矩阵 C=ey
矩阵求逆的快速算法
01-18 5593
龚敏敏　
高斯约当法求逆矩阵的算法实现（C++）
u014558968的专栏
06-20 4106
#include"iostream.h" #include"math.h" void main() { float a[10][10],A[10][10],b[10],c[10][10],d=0,f=0; int i=0,j=0,k=1,l=0,m=0,n=0; //------------------------------------------------------------------
C/C++语言实现矩阵求逆运算—高斯约化/消元法
luofl的专栏
07-27 1万+
在网站上搜索了一下，发现C语言代码可能是很原始的版本，编译的时候会报出诸多错误。 不过大致还是能看懂，就看你能不能转过弯来咯。 原文的代码是要存成矩阵形式，而实际代码写的着则是一维数组形式。 可以直接使用的代码如下： #include #include #include // 对矩阵求逆，结果保存在a中 int brinv(double a[], int n) {
请帮我写一个高阶矩阵求逆程序，用c语言，使用给定数组和随机矩阵，使用高斯-约旦消元法来求解逆矩阵
05-31
好的，以下是一个基于高斯-约旦消元法的矩阵求逆的C语言程序。程序中使用了随机矩阵和给定的数组来测试求逆算法的正确性。 ```c #include <stdio.h> #include <stdlib.h> #include <time.h> #define N 5 // 矩阵的大小 // 打印矩阵 void print_matrix(double **A, int n) { for (int i = 0; i < n; i++) { for (int j = 0; j < n; j++) { printf("%f ", A[i][j]); } printf("\n"); } } // 高斯-约旦消元法求解逆矩阵 void invert(double **A, double **B, int n) { // 初始化单位矩阵 for (int i = 0; i < n; i++) { for (int j = 0; j < n; j++) { B[i][j] = (i == j) ? 1.0 : 0.0; } } // 逐列消元 for (int k = 0; k < n; k++) { double pivot = A[k][k]; for (int j = 0; j < n; j++) { A[k][j] /= pivot; B[k][j] /= pivot; } for (int i = 0; i < n; i++) { if (i != k) { double factor = A[i][k]; for (int j = 0; j < n; j++) { A[i][j] -= factor * A[k][j]; B[i][j] -= factor * B[k][j]; } } } } } int main() { double **A = malloc(N * sizeof(double *)); double **B = malloc(N * sizeof(double *)); double *arr = malloc(N * N * sizeof(double)); // 用于存储随机矩阵 srand(time(NULL)); // 用系统时间初始化随机数生成器 // 初始化随机矩阵和数组 for (int i = 0; i < N; i++) { A[i] = malloc(N * sizeof(double)); B[i] = malloc(N * sizeof(double)); for (int j = 0; j < N; j++) { arr[i * N + j] = rand() % 10 + 1; // 生成1~10之间的随机数 A[i][j] = arr[i * N + j]; } } printf("随机矩阵为：\n"); print_matrix(A, N); invert(A, B, N); printf("逆矩阵为：\n"); print_matrix(B, N); // 检查逆矩阵是否正确 for (int i = 0; i < N; i++) { for (int j = 0; j < N; j++) { double sum = 0.0; for (int k = 0; k < N; k++) { sum += A[i][k] * B[k][j]; } if (i == j && abs(sum - 1.0) > 1e-6) { printf("逆矩阵有误！\n"); break; } if (i != j && abs(sum) > 1e-6) { printf("逆矩阵有误！\n"); break; } } } // 释放内存 free(arr); for (int i = 0; i < N; i++) { free(A[i]); free(B[i]); } free(A); free(B); return 0; } ``` 注意：由于本程序使用了动态内存分配，因此必须在程序结束时手动释放内存，否则会导致内存泄漏。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值