通过编程来学习线性代数3-行列式的性质

cover

行列式类

在讨论行列式的一些性质之前，我先根据之前讲到的内容写一个行列式的类。

// linera algebra determinant class
// filename: det.js

let Det;

(function(){
    Det = function (array) {
        this.array = array
        this.length = array.length
        // 阶乘，元素个数
        this._itemLength = null
    
        // 保存逆序数的数组
        this.inverseNumberArray = null
    
        // 保存每个项的数组
        this._items = null
        this.itemValues = null

        // 是否需要重新计算
        this.reCalc = false
    }
    
    // 阶乘
    Det.prototype.itemLength = function () {
        if (this._itemLength == null) {
            this._itemLength = factorial(this.length)
        }
        return this._itemLength
    }
    
    // 逆序数
    Det.prototype.inverseNumber = function (index) {
        if (this.inverseNumberArray == null) {
            this.inverseNumberArray = new Array(this.itemLength())
        }
        if (this.inverseNumberArray[index] == undefined) {
            let sum = 0
            let item = this.items()[index]
            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++
                    }
                }
            }
            this.inverseNumberArray[index] = sum
        }
        return this.inverseNumberArray[index]
    }
    
    // 获取保存每个项的数组
    Det.prototype.items = function () {
        if (this._items == null) {
            this._items = []
    
            let standardIndex = []
            for (let i = 0; i < this.length; i++) {
                standardIndex.push(i)
            }
            generate(this.length, standardIndex, this._items)
        }
        return this._items
    }
    
    // 获取单个项的值
    Det.prototype.itemValue = function (index) {
        if (this.itemValues == null) {
            this.itemValues = new Array(this.itemLength())
        }
        if(this.itemValues[index]==undefined || this.reCalc) {
            let inverseCount = this.inverseNumber(index)
            let data = this.array
            let item = this.items()[index]
            let value = (inverseCount % 2 ? -1 : 1)
            for (let j = 0, n = this.length; j < n; j++) {
                value *= data[j][item[j]]
            }
            this.itemValues[index] = value
        }
        return this.itemValues[index]
    }
    
    Det.prototype.calc = function () {
        let sum = 0
        for (let i = 0, len = this.itemLength(); i < len; i++) {
            sum += this.itemValue(i)
        }
    
        console.log(this.array)
        console.log(sum)
    
        return sum
    }
})()

// filename util.js

/**
 * 列举所有@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)
    }
}

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

/**
 * 交互数组中的两个元素
 * 
 * @param {Array} arr
 * @param {Number} i
 * @param {Number} j
 */
function swap(arr, i, j) {
    var temp = arr[i]
    arr[i] = arr[j]
    arr[j] = temp
    return arr
}

基本的东西准备好了，下面我们往行列式类Det上加一些功能，即行列式的一些性质。

性质1：行列式与它的转置行列式相等

行列式的转置行列式即行变成列，列变成行。

转置行列式

在det.js添加：

// 获取转置行列式
Det.prototype.getTransposedDet = function () {
    let len = this.length
    let newArr = new Array(len)
    for (let i = 0; i < len; i++) {
        if(!newArr[i]) {
            newArr[i] = new Array(len)
        }
        for (let j = 0; j < len; j++) {
            newArr[i][j] = this.array[j][i]
        }
    }
    return new Det(newArr)
}

测试：

let det = new Det([
    [2, 1, -5, 1],
    [1, -3, 0, -6],
    [0, 2, -1, 2],
    [1, 4, -7, 6]
])
det.calc()  // 27

let tdet = det.getTransposedDet()
tdet.calc() // 27

这里可以把tdet.array打印出来，看是不是真的转置成功。

性质2：互换行列式的两行（列），行列式变号

添加如下方法

// 互换行列式的两行（列）
Det.prototype.swap = function (n0, n1, isRow=true) {
    let newArr = JSON.parse(JSON.stringify(this.array))

    if(isRow) {
        newArr = swap(newArr, n0, n1)
    } else {
        let len = this.length
        for (let i = 0; i < len; i++) {
            newArr[i] = swap(newArr[i], n0, n1)
        }
    }
    return new Det(newArr)
}

测试：

let det = new Det([
    [2, 1, -5, 1],
    [1, -3, 0, -6],
    [0, 2, -1, 2],
    [1, 4, -7, 6]
])
det.calc()  // 27

let tdet = det.swap(3, 1, true) // -27
// let tdet = det.swap(3, 1, false) // -27
tdet.calc()

性质3：行列式的某一行（列）中所有元素都乘以同一个数k，等于用数k乘以此行列式

添加方法：

// 某一行（列）乘以一个数
/**
 * 
 * @param {Number} n 行/列，从0开始
 * @param {Number} k 数 
 * @param {Boolean} isRow 默认行
 */
Det.prototype.multiply = function (n, k, isRow = true) {
    let newArr = JSON.parse(JSON.stringify(this.array)) //deep copy
    let len = this.length

    if (isRow) {
        for (let i = 0; i < len; i++) {
            newArr[n][i] *= k
        }
    } else {
        for (let i = 0; i < len; i++) {
            newArr[i][n] *= k
        }
    }

    return new Det(newArr)
}

通过计算乘之前和之后的两个结果来看：

let det = new Det([
    [2, 1, -5, 1],
    [1, -3, 0, -6],
    [0, 2, -1, 2],
    [1, 4, -7, 6]
])
det.calc() // 27

det.multiply(1, 5).calc()   // 135

也可以通过计算下面两个行列式来验证，第二个行列式的第一列是第一个行列式第一列的两倍，计算结果分别为27和54

let det = new Det([
    [2, 1, -5, 1],
    [1, -3, 0, -6],
    [0, 2, -1, 2],
    [1, 4, -7, 6]
])
det.calc()  // 27

let det2 = new Det([
    [4, 1, -5, 1],
    [2, -3, 0, -6],
    [0, 2, -1, 2],
    [2, 4, -7, 6]
])
det2.calc() // 54

性质4：行列式中如果有两行（列）元素成比例，则此行列式为零

可通过计算下面行列式验证：

let det = new Det([
    [2, 1, -5, 1],  //1
    [1, -3, 0, -6],
    [0, 2, -1, 2],
    [4, 2, -10, 2]  //2*k
])
det.calc()  // 0

let det2 = new Det([
    [2, 1, -5, 1],  //1
    [1, -3, 0, -6],
    [0, 2, -1, 2],
    [2, 1, -5, 1]  //1
])
det2.calc()  // 0

性质5：若行列式的某一行（列）的元素都是两个数之和，如：

let det = new Det([
    [2 + 3, 1, 2, 1],
    [1 + 4, -3, 5, -6],
    [0 - 5, 2, 2, 2],
    [4 + 1, 2, -4, 4]
])
det.calc()  //80

let det2 = new Det([
    [2, 1, 2, 1],
    [1, -3, 5, -6],
    [0, 2, 2, 2],
    [4, 2, -4, 4]
])
det2.calc() //36

let det3 = new Det([
    [3, 1, 2, 1],
    [4, -3, 5, -6],
    [-5, 2, 2, 2],
    [1, 2, -4, 4]
])
det3.calc() //44

性质6：把行列式的某一行（列）的各元素乘以同一个倍数加到另一行（列）对应的元素上去，行列式不变。

性质6

在Det类上添加方法

// 性质6：把行列式的某一行（列）的各元素乘以同一个倍数加到另一行（列）对应的元素上去，行列式不变。
/**
    * 
    * @param {Number} n0 行/列
    * @param {Number} n1 行/列
    * @param {Number} k 行列式的 n0(行/列) + n1(行/列)*k
    * @param {Boolean} isRow
    */
Det.prototype.plusLine = function (n0, n1, k, isRow=true) {
    if(n0==n1) {
        throw('不能加到同一行或列')
    }
    let newArr = JSON.parse(JSON.stringify(this.array))
    let len = this.length

    if (isRow) {
        for (let i = 0; i < len; i++) {
            newArr[n0][i] += newArr[n1][i] * k
        }
    } else {
        for (let i = 0; i < len; i++) {
            newArr[i][n0] += newArr[i][n1] * k
        }
    }

    return new Det(newArr)
}

验证：

// 性质6
let det = new Det([
    [1, 1, 2, 1],
    [1, -3, 5, 3],
    [0, 2, 2, 2],
    [1, 2, -4, 4]
])
det.calc()  //-112

// 第1列每行对应元素加上第4列每行对应元素乘3
let det1 = det.plusLine(0, 3, 3, false)
det1.calc()  //-112

代码

参考

最后编辑于 ：
©著作权归作者所有,转载或内容合作请联系作者
  • 人面猴
    序言：七十年代末，一起剥皮案震惊了整个滨河市，随后出现的几起案子，更是在滨河造成了极大的恐慌，老刑警刘岩，带你破解...
    沈念sama阅读 216,692评论 6赞 501
  • 死咒
    序言：滨河连续发生了三起死亡事件，死亡现场离奇诡异，居然都是意外死亡，警方通过查阅死者的电脑和手机，发现死者居然都...
    沈念sama阅读 92,482评论 3赞 392
  • 救了他两次的神仙让他今天三更去死
    文/潘晓璐 我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
    开封第一讲书人阅读 162,995评论 0赞 353
  • 道士缉凶录：失踪的卖姜人
    文/不坏的土叔 我叫张陵，是天一观的道长。 经常有香客问我，道长，这世上最难降的妖魔是什么？ 我笑而不...
    开封第一讲书人阅读 58,223评论 1赞 292
  • 港岛之恋（遗憾婚礼）
    正文 为了忘掉前任，我火速办了婚礼，结果婚礼上，老公的妹妹穿的比我还像新娘。我一直安慰自己，他们只是感情好，可当我...
    茶点故事阅读 67,245评论 6赞 388
  • 恶毒庶女顶嫁案：这布局不是一般人想出来的
    文/花漫 我一把揭开白布。 她就那样静静地躺着，像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上，一...
    开封第一讲书人阅读 51,208评论 1赞 299
  • 城市分裂传说
    那天，我揣着相机与录音，去河边找鬼。 笑死，一个胖子当着我的面吹牛，可吹牛的内容都是我干的。 我是一名探鬼主播，决...
    沈念sama阅读 40,091评论 3赞 418
  • 双鸳鸯连环套：你想象不到人心有多黑
    文/苍兰香墨 我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
    开封第一讲书人阅读 38,929评论 0赞 274
  • 万荣杀人案实录
    序言：老挝万荣一对情侣失踪，失踪者是张志新（化名）和其女友刘颖，没想到半个月后，有当地人在树林里发现了一具尸体，经...
    沈念sama阅读 45,346评论 1赞 311
  • 护林员之死
    正文 独居荒郊野岭守林人离奇死亡，尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,570评论 2赞 333
  • 白月光启示录
    正文 我和宋清朗相恋三年，在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,739评论 1赞 348
  • 活死人
    序言：一个原本活蹦乱跳的男人离奇死亡，死状恐怖，灵堂内的尸体忽然破棺而出，到底是诈尸还是另有隐情，我是刑警宁泽，带...
    沈念sama阅读 35,437评论 5赞 344
  • 日本核电站爆炸内幕
    正文 年R本政府宣布，位于F岛的核电站，受9级特大地震影响，放射性物质发生泄漏。R本人自食恶果不足惜，却给世界环境...
    茶点故事阅读 41,037评论 3赞 326
  • 男人毒药：我在死后第九天来索命
    文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹，春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,677评论 0赞 22
  • 一桩弑父案，背后竟有这般阴谋
    文/苍兰香墨 我抬头看了看天上的太阳。三九已至，却和暖如春，着一层夹袄步出监牢的瞬间，已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,833评论 1赞 269
  • 情欲美人皮
    我被黑心中介骗来泰国打工， 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留，地道东北人。 一个月前我还...
    沈念sama阅读 47,760评论 2赞 369
  • 代替公主和亲
    正文 我出身青楼，却偏偏与公主长得像，于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子，可洞房花烛夜当晚...
    茶点故事阅读 44,647评论 2赞 354