线性代数之行列式的C#研究实现
最近学习机器学习 才发现以前数学没有学好 开始从线性代数开始学起 读完行列式一章写了些C#的代码学习一下。
直接上C#代码:
using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
using System.Runtime.InteropServices;
using System.IO; namespace LYF.Math
{
/// <summary>
/// 行列式 Determinant
/// </summary>
[SerializableAttribute]
[ComVisibleAttribute(true)]
public class Determinant<T> where T : IComparable, IFormattable, IConvertible, IComparable<T>, IEquatable<T>
{
T[,] tarr = null;
public Determinant(int n)
{
tarr = new T[n, n];
} public Determinant(T[,] arrT)
{
if (arrT == null || arrT.GetLength(0) != arrT.GetLength(1) || arrT.GetLength(0) < 1)
{
throw new MathException("不正确的数组(数组必须行列数相同且大于1)");
}
else
{
tarr=new T[arrT.GetLength(0),arrT.GetLength(0)];
SetItem(arrT);
}
} /// <summary>
/// 获取元素值
/// </summary>
/// <param name="i"></param>
/// <param name="j"></param>
/// <returns></returns>
public T this[int i, int j]
{
//实现索引器的get方法
get
{
return GetItem(i, j);
} //实现索引器的set方法
set
{
SetItem(i, j, value);
}
} /// <summary>
/// 获取元素的余子式
/// </summary>
/// <param name="i"></param>
/// <param name="j"></param>
/// <returns></returns>
public Determinant<T> A(int i, int j)
{
if (N == 1)
{
return null;
}
else if (i>N||j>N)
{
return null;
}
else
{
Determinant<T> a = new Determinant<T>(N - 1);
for (int m = 1; m <= N - 1; m++)
{
for (int n = 1; n <= N - 1; n++)
{
int p = m, q = n;
if (p >= i)
{
p = m + 1;
}
if (q >= j)
{
q = n + 1;
}
a[m, n] = this[p,q];
}
}
return a;
}
} /// <summary>
/// 设置行列式的值
/// </summary>
/// <param name="i">行数(从1开始)</param>
/// <param name="j">列数(从1开始)</param>
/// <param name="value">值</param>
public void SetItem(int i, int j, T value)
{
if (tarr == null)
{
throw new MathException("行列式未正确初始化");
}
else if (i > N || j > N)
{
throw new MathException("超出行列式索引范围");
}
else
{
tarr[i - 1, j - 1] = value;
}
} public void SetItem(T[,] arrT)
{
if (arrT == null || tarr == null)
{
throw new MathException("不能为空");
}
else if (arrT.GetLength(0) != N || arrT.GetLength(1) != N)
{
throw new MathException("传入阶数不同");
}
else
{
for (int m = 0; m <=N-1; m++)
{
for (int n = 0; n <= N- 1; n++)
{
this[m + 1, n + 1] = arrT[m, n];
}
}
}
} /// <summary>
/// 设置行列式的值
/// </summary>
/// <param name="i">行数(从1开始)</param>
/// <param name="j">列数(从1开始)</param>
/// <param name="value">值</param>
public T GetItem(int i, int j)
{
if (tarr == null)
{
throw new MathException("行列式未正确初始化");
}
else if (i > N || j > N)
{
throw new MathException("超出行列式索引范围");
}
else
{
return tarr[i-1, j-1];
}
} /// <summary>
/// 输出行列式信息
/// </summary>
/// <returns></returns>
public override string ToString()
{
StringBuilder sbRs = new StringBuilder();
if(tarr!=null)
{
for (int m = 0; m <= N - 1; m++)
{
for (int n = 0; n <= N - 1; n++)
{
sbRs.Append(string.Format("{0}\t", tarr[m, n]));
}
sbRs.Append("\n");
} }
return sbRs.ToString();
} /// <summary>
/// 获取行列式的阶数
/// </summary>
public int N
{
get{
if (tarr != null)
{
return tarr.GetLength(0);
}
else
{
return 0;
}
} } private string typeName = string.Empty;
private string GetType()
{
if (string.IsNullOrEmpty(typeName))
{
typeName=typeof(T).Name;
File.AppendAllText("E:\\op.txt", typeName);
}
return typeName; } /// <summary>
/// 获取行列式的值
/// </summary>
public T Value
{
get
{
if (N == 1)
{
return tarr[0, 0];
}
else if (N == 2)
{
return Minus(MUL(tarr[0, 0], tarr[1, 1]), MUL(tarr[0, 1], tarr[1, 0]));
}
else
{
T sum = default(T);
for (int i = 1; i <= N; i++)
{
if ((1+i) % 2 == 0)
{
//余子式正值
sum = Add(sum, MUL(this[1, i], this.A(1, i).Value));
}
else
{
//余子式负值
sum = Minus(sum, MUL(this[1, i], this.A(1, i).Value));
}
}
return sum;
} }
} /// <summary>
/// 加法
/// </summary>
/// <param name="left"></param>
/// <param name="right"></param>
/// <returns></returns>
private T Add(T left, T right)
{
switch (GetType())
{
case "Int16":
return ((T)(object)((short)(object)left + (short)(object)right));
case "Int32":
return ((T)(object)((int)(object)left + (int)(object)right));
case "Int64":
return ((T)(object)((long)(object)left + (long)(object)right));
case "Single":
return ((T)(object)((float)(object)left + (float)(object)right));
case "Double":
return ((T)(object)((double)(object)left + (double)(object)right));
case "Decimal":
return ((T)(object)((decimal)(object)left + (decimal)(object)right));
}
throw new MathException("不支持的操作类型");
} /// <summary>
/// 减法
/// </summary>
/// <param name="left"></param>
/// <param name="right"></param>
/// <returns></returns>
private T Minus(T left, T right)
{
switch (GetType())
{
case "Int16":
return ((T)(object)((short)(object)left - (short)(object)right));
case "Int32":
return ((T)(object)((int)(object)left - (int)(object)right));
case "Int64":
return ((T)(object)((long)(object)left - (long)(object)right));
case "Single":
return ((T)(object)((float)(object)left - (float)(object)right));
case "Double":
return ((T)(object)((double)(object)left - (double)(object)right));
case "Decimal":
return ((T)(object)((decimal)(object)left - (decimal)(object)right));
}
throw new MathException("不支持的操作类型");
} /// <summary>
/// 乘法
/// </summary>
/// <param name="left"></param>
/// <param name="right"></param>
/// <returns></returns>
private T MUL(T left, T right)
{
switch (GetType())
{
case "Int16":
return ((T)(object)((short)(object)left * (short)(object)right));
case "Int32":
return ((T)(object)((int)(object)left * (int)(object)right));
case "Int64":
return ((T)(object)((long)(object)left * (long)(object)right));
case "Single":
return ((T)(object)((float)(object)left * (float)(object)right));
case "Double":
return ((T)(object)((double)(object)left * (double)(object)right));
case "Decimal":
return ((T)(object)((decimal)(object)left * (decimal)(object)right));
}
throw new MathException("不支持的操作类型");
} }
}
以上代码就是对行列式的封装 可以求值获得余子式 很基本的东西 求值的话主要用了递归的方式 因为泛型的原因导致计算过程重复拆箱装箱 不过目前好像也没有什么太好的方法了。反正就是学习 所以性能无所谓了。
然后就是调用了直接上调用代码:
int[,] aaa = new int[4, 4]{{1,2,3,6},
{4,5,7,8},
{7,8,9,10},
{3,8,4,3}}; //LYF.Math.Determinant<int> d = new Determinant<int>(4);
LYF.Math.Determinant<int> d = new Determinant<int>(aaa);
d.SetItem(aaa);
Console.WriteLine("当前行列式:");
Console.WriteLine(d.ToString());
Console.WriteLine("余子式M11:");
Console.WriteLine(d.A(1, 1).ToString());
Console.WriteLine("余子式M12:");
Console.WriteLine(d.A(1, 2).ToString());
Console.WriteLine("余子式M22:");
Console.WriteLine(d.A(2, 2).ToString());
Console.WriteLine("N="+d.N);
Console.WriteLine("行列式的值为:"+d.Value.ToString());
Console.Read();
执行结果如下:
线性代数之行列式的C#研究实现的更多相关文章
- CSUOJ 1979 古怪的行列式
Description 这几天,子浩君潜心研究线性代数. 行列式的值定义如下: 其中,τ(j1j2...jn)为排列j1j2...jn的逆序数. 子浩君很厉害的,但是头脑经常短路,所以他会按照行列式值 ...
- 线性代数笔记24——微分方程和exp(At)
原文:https://mp.weixin.qq.com/s/COpYKxQDMhqJRuMK2raMKQ 微分方程指含有未知函数及其导数的关系式,解微分方程就是找出未知函数.未知函数是一元函数的,叫常 ...
- 【线性代数】5-3:克莱姆法则,逆和体积(Cramer's Rule,Inverses,and Volumes)
title: [线性代数]5-3:克莱姆法则,逆和体积(Cramer's Rule,Inverses,and Volumes) categories: Mathematic Linear Algebr ...
- [Swust OJ 643]--行列式的计算(上三角行列式变换)
题目链接:http://acm.swust.edu.cn/problem/643/ Time limit(ms): 1000 Memory limit(kb): 65535 Description ...
- TOJ4537: n阶行列式
4537: n阶行列式 Time Limit(Common/Java):1000MS/3000MS Memory Limit:65536KByteTotal Submit: 28 ...
- Python人工智能之-三大数学难点 !
1. 微积分: 定积分与不定积分.全微分.最小二乘法.二重积分.微分方程与差分方程等... 2. 线性代数: 行列式.矩阵.向量.线性方程组.矩阵的特性和特性向量.二次型等... 3. 概率论和统计学 ...
- NOIP考点
NOIP考点 基础算法 图 树 数论 数据结构 动态规划 搜索 其他算法 省选知识点汇总 图论 数据结构 字符串相关算法及数据结构 数学 计算几何 搜索 动态规划 其他算法 转自:巨佬的博客 加*号是 ...
- [省选]板块(shenben已经AFO!!!)
shenben已经AFO!!! 部分摘抄自网络 同样的,加粗是重点,星号是选学 图论 网络流(dinic,ISAP选一个,费用流写EK就行.*zkw费用流),二分图 点分治,边分治,*动态点分治 树链 ...
- woj1013 Barcelet 字符串 woj1014 Doraemon's Flashlight 几何
title: woj1013 Barcelet 字符串 date: 2020-03-18 18:00:00 categories: acm tags: [acm,字符串,woj] 字符串,字典序. 1 ...
随机推荐
- 使用hive客户端java api读写hive集群上的信息
上文介绍了hdfs集群信息的读取方式,本文说hive 1.先解决依赖 <properties> <hive.version>1.2.1</hive.version> ...
- marked插件在线实时解析markdown的web小工具
访问地址: https://mdrush.herokuapp.com/ github项目: https://github.com/qcer/MDRush 实现简介: 1.动态数据绑定 借助Vuejs, ...
- MySQL(十五)之数据备份中mysqldump详解
前言 其实前面一篇数据备份已经是非常的详细了,这里我想单独的讲解一下mysqldump,相信很多程序员都是用过这个命令的! 一.MySQL数据库的备份与还原 1.1.MySQL数据库备份 1)语法 m ...
- 前端页面卡顿、也许是DOM操作惹的祸?
界面上UI的更改都是通过DOM操作实现的,并不是通过传统的刷新页面实现 的.尽管DOM提供了丰富接口供外部调用,但DOM操作的代价很高,页面前端代码的性能瓶颈也大多集中在DOM操作上,所以前端性能优化 ...
- OpenCV 学习笔记(模板匹配)
OpenCV 学习笔记(模板匹配) 模板匹配是在一幅图像中寻找一个特定目标的方法之一.这种方法的原理非常简单,遍历图像中的每一个可能的位置,比较各处与模板是否"相似",当相似度足够 ...
- 【ASP.NET MVC 学习笔记】- 10 Controller和Action(1)
本文参考:http://www.cnblogs.com/willick/p/3331521.html 1.继承IController接口,示例代码将当前请求的Controller和Action打印到浏 ...
- 白夜追凶 :手 Q 图片的显示和发送逻辑
欢迎大家前往腾讯云社区,获取更多腾讯海量技术实践干货哦~ 作者:陈舜尧 导语: "这张图片在快捷发图栏背景是黑色的,为啥发到AIO(会话窗口)里背景就变成白的了?" 通过一个bug ...
- (转)Spark JAVA RDD API
对API的解释: 1.1 transform l map(func):对调用map的RDD数据集中的每个element都使用func,然后返回一个新的RDD,这个返回的数据集是分布式的数据集 l ...
- JSP异常之org.apache.jasper.JasperException(转)
According to TLD or attribute directive in tag file, attribute items does not accep t any expression ...
- 注销/etc/passwd带来的系统登陆不上
今天在修改虚拟机密码上的时候,将/etc/passwd中root所在的哪行注销掉了,想象是注销了,root登陆时应该不要输入密码,结果是系统进度条走到最后的时候 进入不了系统了. 结果去普及了下/et ...