楼主: poptide

6128 4

[其他] 里昂惕夫逆矩阵计算求助 [推广有奖]

已卖：121份资源

副教授

21%

还不是VIP/贵宾

楼主

poptide 发表于 2016-9-3 20:15:10 |AI写论文

20论坛币

A1 = cell2mat(A);
E = eye(1435);
L = inv(E - A1);
我在运用投入产出表计算里昂惕夫时遇到一个问题，主要是在求里昂惕夫逆矩阵时矩阵A中有元素是0，根源在于投入产出表中有一部分为0.理论上来说应该是允许为0，但是在逆矩阵运算该如何解决？（高精度也没成功）

把分母output那一列的0改成一个很小的数即可。关键问题不是A中有0，而是如果分母是0，那么有些系数就是nan，matlab遇到nan就没法计算逆矩阵。

沙发

wisdom101126 发表于 2016-9-3 20:15:11

把分母output那一列的0改成一个很小的数即可。
关键问题不是A中有0，而是如果分母是0，那么有些系数就是nan，matlab遇到nan就没法计算逆矩阵。

已有 1 人评分	经验	收起理由
giresse	+ 60	精彩帖子

总评分: 经验 + 60 查看全部评分

藤椅

dgy203 发表于 2016-9-3 20:37:35

矩阵中有元素为0与矩阵可逆没有什么关系，只要行列式不为0就行。如果矩阵不可逆，也可以用广义逆求之，pinv函数

板凳

wwwseeking 发表于 2016-9-4 14:45:39

用pseudo inverse, 或者借用线性方程组的思想求解inverse，比如用lsqr，查下matlab help 很容易的

报纸

熊洪 发表于 2016-9-8 13:21:25

用pseudo inverse, 或者借用线性方程组的思想求解inverse，比如用lsqr，查下matlab help 很容易的

jg-xs1
拉您进交流群