C#计算矩阵的逆矩阵方法实例分析
作者:北风其凉 发布时间:2021-11-06 18:24:29
标签:C#,矩阵
本文实例讲述了C#计算矩阵的逆矩阵方法。分享给大家供大家参考。具体如下:
1.代码思路
1)对矩阵进行合法性检查:矩阵必须为方阵
2)计算矩阵行列式的值(Determinant函数)
3)只有满秩矩阵才有逆矩阵,因此如果行列式的值为0(在代码中以绝对值小于1E-6做判断),则终止函数,报出异常
4)求出伴随矩阵(AdjointMatrix函数)
5)逆矩阵各元素即其伴随矩阵各元素除以矩阵行列式的商
2.函数代码
(注:本段代码只实现了一个思路,可能并不是该问题的最优解)
/// <summary>
/// 求矩阵的逆矩阵
/// </summary>
/// <param name="matrix"></param>
/// <returns></returns>
public static double[][] InverseMatrix(double[][] matrix)
{
//matrix必须为非空
if (matrix == null || matrix.Length == 0)
{
return new double[][] { };
}
//matrix 必须为方阵
int len = matrix.Length;
for (int counter = 0; counter < matrix.Length; counter++)
{
if (matrix[counter].Length != len)
{
throw new Exception("matrix 必须为方阵");
}
}
//计算矩阵行列式的值
double dDeterminant = Determinant(matrix);
if (Math.Abs(dDeterminant) <= 1E-6)
{
throw new Exception("矩阵不可逆");
}
//制作一个伴随矩阵大小的矩阵
double[][] result = AdjointMatrix(matrix);
//矩阵的每项除以矩阵行列式的值,即为所求
for (int i = 0; i < matrix.Length; i++)
{
for (int j = 0; j < matrix.Length; j++)
{
result[i][j] = result[i][j] / dDeterminant;
}
}
return result;
}
/// <summary>
/// 递归计算行列式的值
/// </summary>
/// <param name="matrix">矩阵</param>
/// <returns></returns>
public static double Determinant(double[][] matrix)
{
//二阶及以下行列式直接计算
if (matrix.Length == 0) return 0;
else if (matrix.Length == 1) return matrix[0][0];
else if (matrix.Length == 2)
{
return matrix[0][0] * matrix[1][1] - matrix[0][1] * matrix[1][0];
}
//对第一行使用“加边法”递归计算行列式的值
double dSum = 0, dSign = 1;
for (int i = 0; i < matrix.Length; i++)
{
double[][] matrixTemp = new double[matrix.Length - 1][];
for (int count = 0; count < matrix.Length - 1; count++)
{
matrixTemp[count] = new double[matrix.Length - 1];
}
for (int j = 0; j < matrixTemp.Length; j++)
{
for (int k = 0; k < matrixTemp.Length; k++)
{
matrixTemp[j][k] = matrix[j + 1][k >= i ? k + 1 : k];
}
}
dSum += (matrix[0][i] * dSign * Determinant(matrixTemp));
dSign = dSign * -1;
}
return dSum;
}
/// <summary>
/// 计算方阵的伴随矩阵
/// </summary>
/// <param name="matrix">方阵</param>
/// <returns></returns>
public static double[][] AdjointMatrix(double [][] matrix)
{
//制作一个伴随矩阵大小的矩阵
double[][] result = new double[matrix.Length][];
for (int i = 0; i < result.Length; i++)
{
result[i] = new double[matrix[i].Length];
}
//生成伴随矩阵
for (int i = 0; i < result.Length; i++)
{
for (int j = 0; j < result.Length; j++)
{
//存储代数余子式的矩阵(行、列数都比原矩阵少1)
double[][] temp = new double[result.Length - 1][];
for (int k = 0; k < result.Length - 1; k++)
{
temp[k] = new double[result[k].Length - 1];
}
//生成代数余子式
for (int x = 0; x < temp.Length; x++)
{
for (int y = 0; y < temp.Length; y++)
{
temp[x][y] = matrix[x < i ? x : x + 1][y < j ? y : y + 1];
}
}
//Console.WriteLine("代数余子式:");
//PrintMatrix(temp);
result[j][i] = ((i + j) % 2 == 0 ? 1 : -1) * Determinant(temp);
}
}
//Console.WriteLine("伴随矩阵:");
//PrintMatrix(result);
return result;
}
/// <summary>
/// 打印矩阵
/// </summary>
/// <param name="matrix">待打印矩阵</param>
private static void PrintMatrix(double[][] matrix, string title = "")
{
//1.标题值为空则不显示标题
if (!String.IsNullOrWhiteSpace(title))
{
Console.WriteLine(title);
}
//2.打印矩阵
for (int i = 0; i < matrix.Length; i++)
{
for (int j = 0; j < matrix[i].Length; j++)
{
Console.Write(matrix[i][j] + "\t");
//注意不能写为:Console.Write(matrix[i][j] + '\t');
}
Console.WriteLine();
}
//3.空行
Console.WriteLine();
}
3.Main函数调用
static void Main(string[] args)
{
double[][] matrix = new double[][]
{
new double[] { 1, 2, 3 },
new double[] { 2, 2, 1 },
new double[] { 3, 4, 3 }
};
PrintMatrix(matrix, "原矩阵");
PrintMatrix(AdjointMatrix(matrix), "伴随矩阵");
Console.WriteLine("行列式的值为:" + Determinant(matrix) + '\n');
PrintMatrix(InverseMatrix(matrix), "逆矩阵");
Console.ReadLine();
}
4.执行结果
希望本文所述对大家的C#程序设计有所帮助。
0
投稿
猜你喜欢
- 前言:时间过得真是快,现在已经是2022年了。作为开发来说,时间处理是非常繁琐的。从Java 8开始有了Java 8 T
- 前言 短时间提升自己最快的手段就是背面试题,最近总结了Java常用的面试题,分享给大家,希望大家都能圆梦大厂,加油,我命由我不由天
- 在POI的使用过程中,对大多数API User来说经常面临两个问题,这也是GridExcel致力解决的问题。问题1. 仅使用简单的导入导出功
- Java 线程池ExecutorService1.线程池1.1什么情况下使用线程池单个任务处理的时间比较短.将需处理的任务的数量大.1.2使
- 本篇实例内容是关于C#读取CAD文件的,直接看代码//在不使用任务插件的情况下读取DWG文件的缩略图,以便在没有安装AutoCAD的计算机上
- 程序如下:View Code /* * Hanoi塔游戏 问题描述: * 汉诺塔:汉诺塔(又称河内塔)问
- 本文实例讲述了java实现单词搜索迷宫游戏。分享给大家供大家参考。具体分析如下:我们在杂志上,经常能够看到找单词的小游戏,在一个二维表格中,
- SpringBoot2之PUT请求接收不了参数的解决办法,这个问题,关乎两个Filter过滤器,是spring3和3.5之后提供的,目的就是
- web 容器的设计开发一个web容器涉及很多不同方面不同层面的技术,例如通信层的知识,程序语言层面的知识等等,且一个可用的web容器是一个比
- SpringMvc中普通类注入Service为null场景:使用Quartz定时器时,普通的java类需要注入spring的service类
- 一、使用@Profile1.1、@Profile修饰类开发环境package com.example.demo.config;import
- 原因:feign传值出错无法接收到传值由于是POST所以添加@RequestBody进行尝试解决:错误原因是未添加@RequestBody尝
- 尝试了各种防止中文乱码的方式,但是还是乱码;最后还是细节问题导致;解决方式:以及俩种方式是百度的,我的问题不是这俩块1.在requestMa
- 添加MyBatis的代码并修改以下部分:1.添加MyBatisConfigpackage myshop.config;import java
- LeetCode54. 螺旋矩阵 java实现题目难度 中给定一个包含 m x n 个元素的矩阵(m 行, n 列),请按照顺时针螺旋顺序,
- 一、项目需求二、项目思路1、菜单制作2、中奖号码生成 getNumber (随机数 Math.random)3、购买号码和中奖号码比对 生成
- 算法文章,总是带给我们无穷的思考和兴趣,一个问题,多种解决方法,看你如何去思考它,对于标题所引出的问题,我觉得,使用递归是比较有效的方法,当
- 图片上传功能是我们web里面经常用到的,获得的方式也有很多种,这里我用的是request.getInputStream()获取文件流的方式。
- 条件:1、android:ellipsize=”marquee”2、TextView必须单行显示,即内容必须超出TextView
- 本文实例讲述了Spring实战之属性占位符配置器用法。分享给大家供大家参考,具体如下:一 配置文件<?xml version=&quo