通过编程来学习线性代数2-计算行列式的值

cover

全排列及其逆序数

为了计算每一项,我们先要了解如何生成每一个相乘的项,了解什么时候是正数什么时候是负数,关于正负问题就需要了解逆序数的定义。

全排列:

将n个不同的元素排成一列


1 > 1

2 > 1,2 | 2,1

3 > 1,2,3 | 1,3,2 | 2,1,3 | 2,3,1 | 3,1,2 | 3,2,1

...

排列组合的知识中可以知道: n个不同的元素, 从中选取一个放到第一位, 有n钟选法, 剩下n-1个.

继续从这n-1各种继续选取, 放到第二位, 有n-1钟选法.

以此类推, 直到选完为止.

n个元素所有排列的种数: n! = n*(n-1)*...*3*2*1

这里我们用 Heap's algorithm 描述的算法来生成每一项:

// Heap's algorithm

/**
 * 列举所有@param A 数组元素的全排列(排列)
 * 
 * @param {Number} n A长度
 * @param {Array} A 元素
 * @param {Array} result 全部结果
 */
function generate(n, A, result) {
    if (n == 1) {
        result.push(A.slice())
    }
    else {
        for (let i = 0; i < n - 1; i++) {
            generate(n - 1, A, result)
            if (n % 2 == 0) {
                swap(A, i, n - 1)
            }
            else {
                swap(A, 0, n - 1)
            }
        }
        generate(n - 1, A, result)
    }
}

可以在浏览器的控制台中试一下, 如generate(3, [1, 2, 3], arr)或者generate(4, ['a', 'b', 'c', 'd'], arr), 字符数字都行,arr需要事先定义好let arr = []

得到每一项后, 就可以进行计算逆序数了. 一个排列的逆序数决定了这一项是正或负数.

逆序:

n个不同的自然数,规定从小到大为标准次序。当某两个元素的先后次序与标准次序不同时,就称这两个元素组成一个逆序

逆序数:

排列中所有逆序的总数称为此排列的逆序数

计算逆序数的方法,我是直接从第一个开始,依次跟剩下的进行对比:

function calcInverseNumber(item) {
    let sum = 0
    for (let i = 0, len = item.length; i < len; i++) {
        // 取出第i个数
        let digit = item[i]
        // 用第i个数与第i位之后的数进行对比
        for (let j = i + 1; j < len; j++) {
            if (digit > item[j]) {
                sum++
            }
        }
    }
    return sum
}

到这里,生成计算中的每一项和计算每一项的逆序数的方法都有了,接下来就需要一个计算方法。这个计算方法需要传入一个行列式,然后通过generate生成相乘的每一项,再通过calcInverseNumber算出逆序数并相加得出结果。

// 计算行列式的值
/**
 * @param {Array} data 行列式数组
 * 
 */
function calcDeterminantV1(data) {
    let n = data.length
    let standardIndex = []
    for (let i = 0; i < n; i++) {
        standardIndex.push(i)
    }

    let indexArr = []
    generate(n, standardIndex, indexArr)

    let sum = 0
    for (let i = 0, len = factorial(n); i < len; i++) {
        let arr = indexArr[i]
        let inverseCount = calcInverseNumber(arr)

        let item = (inverseCount % 2 ? -1 : 1)
        for (let j = 0; j < n; j++) {
            item *= data[j][arr[j]]
        }
        sum += item
    }

    return sum
}

line 7 获取行列式的长度
line 8-11 生成一个标准排列的数组[0, 1, 2, ..., n-1]
line 13-14 生成0~n-1的所有可能的排列,共n!个,indexArrn!个长度为n的数组,这些数组的元素是0~n-1组成的一个排列。
line 16-26 计算行列式data的值,函数factorial(n)为计算n!的值。
line 17-24遍历indexArr数组,计算每个排列的值。
line 19计算排列arr逆序数line 21判断逆序数inverseCount的正负(奇数为负,偶数为正)。
line 22-25 使data中角标(j, arr[j])对应的数进行相乘,得到item,并追加到总数sum中。

函数factorial(n):

/**
 * 
 * @param {Number} n 
 */
function factorial(n) {
    var result = 1
    for (i = 2; i <= n; i++) {
        result *= i
    }
    return result
}

现在可以试试用这个方法来计算行列式的值了,比如:

[[1, 0, 0],
[0, 3, 0],
[0, 0, 3]]
// sum = 9

[[1, 2, -4],
[-2, 2, 1],
[-3, 4, -2]]
// sum = -14

[[2, 1, -5, 1],
[1, -3, 0, -6],
[0, 2, -1, 2],
[1, 4, -7, 6]]
// sum = 27

目前,计算行列式的值已经告一段落了,下一节将实现一些行列式的延伸。

比如行列式按行(列)展开相关知识。

代码

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 219,110评论 6 508
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 93,443评论 3 395
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 165,474评论 0 356
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 58,881评论 1 295
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 67,902评论 6 392
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,698评论 1 305
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,418评论 3 419
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 39,332评论 0 276
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,796评论 1 316
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,968评论 3 337
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 40,110评论 1 351
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,792评论 5 346
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,455评论 3 331
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 32,003评论 0 22
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 33,130评论 1 272
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 48,348评论 3 373
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 45,047评论 2 355

推荐阅读更多精彩内容