Fortran Coder

查看: 119|回复: 5

[求助] fortran中有关矩阵的计算

[复制链接]

22

帖子

10

主题

0

精华

熟手

F 币
121 元
贡献
71 点
发表于 2019-2-20 16:55:42 | 显示全部楼层 |阅读模式
在fortran中矩阵的行列式怎样表示  det(T)=0
回复

使用道具 举报

638

帖子

0

主题

0

精华

大师

F 币
386 元
贡献
240 点

规矩勋章元老勋章新人勋章水王勋章

QQ
发表于 2019-2-20 23:00:19 | 显示全部楼层
本帖最后由 vvt 于 2019-2-20 23:09 编辑

语法中没有这个函数,需要自己书写,或使用第三方函数库。一般通过LU分解实现,分解后的矩阵对角线上元素的乘积即为原始矩阵行列式的值。

86

帖子

8

主题

0

精华

专家

F 币
542 元
贡献
311 点
发表于 2019-2-21 08:04:33 | 显示全部楼层
vvt 发表于 2019-2-20 23:00
语法中没有这个函数,需要自己书写,或使用第三方函数库。一般通过LU分解实现,分解后的矩阵对角线上元素的 ...

如果不是厄米矩阵的话,感觉应该还是用QR分解计算更合理

22

帖子

10

主题

0

精华

熟手

F 币
121 元
贡献
71 点
 楼主| 发表于 2019-2-21 08:37:10 | 显示全部楼层
vvt 发表于 2019-2-20 23:00
语法中没有这个函数,需要自己书写,或使用第三方函数库。一般通过LU分解实现,分解后的矩阵对角线上元素的 ...

您好,LU分解是怎样实现的啊,能说的详细一点吗?非常感谢

22

帖子

10

主题

0

精华

熟手

F 币
121 元
贡献
71 点
 楼主| 发表于 2019-2-21 08:40:09 | 显示全部楼层
liudy02 发表于 2019-2-21 08:04
如果不是厄米矩阵的话,感觉应该还是用QR分解计算更合理

这个矩阵不是厄米矩阵

1359

帖子

12

主题

5

精华

论坛跑堂

Fcode跑堂伙计

F 币
407 元
贡献
365 点

新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章

发表于 2019-2-21 08:48:27 | 显示全部楼层
用 lapack 或 MKL(内含lapack)库,论坛里应该有不少帖子。
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

QQ|捐赠本站|Archiver|关于我们 About Us|群聊|Fcode

GMT+8, 2019-3-25 07:46

Powered by Discuz! X3.2

© 2001-2017 Comsenz Inc.

快速回复 返回顶部 返回列表