矩阵方程

对于矩阵 A(n,n) 和 B(n,m) 组成的矩阵方程

[A][X] = [B]

记 X(n,m) 的第i列向量为 Xi(i = 1,2...m), 矩阵B的第i列向量为 Bi(i = 1,2...m), 则上述方程等价为

但在实际计算时不应该分别求解,如果是这样的话就造成计算机资源极大浪费,而应该是对所有向量一次选主元消去,然后分别回代。即可以得到方程的解矩阵X。具体做法是将矩阵A(n,n)和B(n,m)组成增广矩阵[AB],通过选主元消去将AB的第1列至第n列变成上三角矩阵,用解上(下)三角方程组的回带方法解方程组 [Aup][Xi] = [Bi] (i = 1,2...m) 。再把Xi拼接起来就得到 [X] 了。

以下是模块代码:

用以下的矩阵方程来验证

输出结果为

原文发布于微信公众号 - 数值分析与有限元编程(program_fem)

原文发表时间:2017-04-11

本文参与腾讯云自媒体分享计划,欢迎正在阅读的你也加入,一起分享。

发表于

我来说两句

0 条评论
登录 后参与评论

相关文章

来自专栏深度学习思考者

深度学习目标检测算法——Faster-Rcnn

Faster-Rcnn代码下载地址:https://github.com/ShaoqingRen/faster_rcnn 一 前言   Faster rcnn是...

3115
来自专栏闪电gogogo的专栏

《统计学习方法》笔记七(1) 支持向量机——线性可分支持向量机

应用拉格朗日对偶性,通过求解对偶问题得到原始问题的最优解,一是因为对偶问题往往更容易求解,二是自然引入核函数,进而推广到非线性分类的问题。

802
来自专栏小鹏的专栏

机器学习模型训练时候tricks

1. 当训练集的效果(准确率)上不去,和贝叶斯估计(人的表现)存在一定差距的时候:(1)增加模型的复杂度。 或者直接多加几层。(2)现在目前效果非常好的一些模型...

2097
来自专栏用户画像

统计学习方法的步骤

版权声明:本文为博主-姜兴琪原创文章,未经博主允许不得转载。

851
来自专栏Petrichor的专栏

深度学习: smooth L1 loss 计算

RPN的目标函数是分类和回归损失的和,分类采用 交叉熵,回归采用稳定的 Smooth L1,

3194
来自专栏机器学习算法工程师

深入浅出解读卷积神经网络

作者:石文华 编辑:田 旭 卷积神经网络 ? 图1 全连接神经网络结构图 ? 图2 卷积神经网络结构图 卷积神经网络和全连接的神经网络结构上的差异还是比较大的,...

2784
来自专栏用户2442861的专栏

VGG论文笔记/小结

版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/haluoluo211/article/d...

3252
来自专栏瓜大三哥

DeepFace

检测→对齐→表示→分类 在文章中,通过3d模型改进了人脸对齐的方法。然后,通过基于4million人脸图像(4000个个体)训练的一个9层的人工神经网络来进行人...

2486
来自专栏大数据风控

R分类算法-KNN算法

监督学习(Supervised Learning): 从给定标注的训练集中学习出一个函数,根据这个函数为新数据进行标注。 无监督学习(Unsupervi...

3457
来自专栏AI研习社

用Kaggle经典案例教你用CNN做图像分类!

前言 在上一篇专栏《利用卷积自编码器对图片进行降噪》中,我们利用卷积自编码器对 MNIST 数据进行了实验,这周我们来看一个 Kaggle 上比较经典的一...

4036

扫码关注云+社区

领取腾讯云代金券