Fortran Coder

查看: 355|回复: 4

[线性代数] 矩阵的运算

[复制链接]

22

帖子

10

主题

0

精华

熟手

F 币
123 元
贡献
72 点
发表于 2019-3-1 09:39:45 | 显示全部楼层 |阅读模式
利用留数定理,用高斯积分公式,被积函数是个矩阵能进行计算吗?
回复

使用道具 举报

120

帖子

8

主题

0

精华

宗师

F 币
2906 元
贡献
509 点
发表于 2019-3-1 10:27:42 | 显示全部楼层
Fortran语法的一大优势就在于对矩阵的便利操作
只要你是对矩阵的每一个元素做同样操作,编程的写法和常量没有任何区别
不知道你具体是怎么计算,感觉语法上写是没问题,但是各矩阵元同时保持精度会有困难

22

帖子

10

主题

0

精华

熟手

F 币
123 元
贡献
72 点
 楼主| 发表于 2019-3-1 20:20:09 | 显示全部楼层
liudy02 发表于 2019-3-1 10:27
Fortran语法的一大优势就在于对矩阵的便利操作
只要你是对矩阵的每一个元素做同样操作,编程的写法和常量没 ...

是的,非常感谢你,我有个问题就是对矩阵积分,但具体操作是对每一个矩阵元积分,那最后的结果是对每一个矩阵元积分的结果怎样处理。

120

帖子

8

主题

0

精华

宗师

F 币
2906 元
贡献
509 点
发表于 2019-3-1 23:17:17 | 显示全部楼层
gaoxing4700 发表于 2019-3-1 20:20
是的,非常感谢你,我有个问题就是对矩阵积分,但具体操作是对每一个矩阵元积分,那最后的结果是对每一个 ...

你既然称是对矩阵积分,说明每个矩阵元是一个函数啊
就按普通积分的办法用取离散点的值加权求和的办法,对Fortran来说和普通积分没有任何区别
不知道你所谓对每一个矩阵元积分结果怎么处理是啥意思
积分完毕之后的结果矩阵的每个矩阵元不就是对应矩阵元函数的积分么
例如A(x)=(1,x;x,x2)这个矩阵你要在-1,1上积分,
用高斯法也无非就是在这个区间上求出某个勒让德多项式的零点xi和对应的权重pi
然后积分结果矩阵就是Sigma(pi*A(xi))么

120

帖子

8

主题

0

精华

宗师

F 币
2906 元
贡献
509 点
发表于 2019-3-1 23:21:12 | 显示全部楼层
gaoxing4700 发表于 2019-3-1 20:20
是的,非常感谢你,我有个问题就是对矩阵积分,但具体操作是对每一个矩阵元积分,那最后的结果是对每一个 ...

不过我猜你所谓的被积函数可能并不是一个简单的矩阵
而是一个含积分参量的,以矩阵为变量的函数,比如exp(A*x)之类的
这个和我说的就完全不是一回事了,那个难点不在于对矩阵积分,而是求矩阵函数
那已经远远远远超出Fortran论坛应该讨论的事情,应该去找数值线性代数的东东……
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

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

GMT+8, 2019-8-21 12:46

Powered by Discuz! X3.2

© 2001-2017 Comsenz Inc.

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