请问这个矩阵求解问题怎么用fortran实现？

心语 · 发表于 2014-9-29 16:24:47

您需要登录才可以下载或查看，没有帐号？立即注册

x

已知矩阵关系式y(P0,N)=g(P0,P)*u(P,N)), 其中矩阵y和u已知，P0、P、N为对应矩阵的维数范围，怎么用fortran求g矩阵呢？求大神指点

心语 · 发表于 2014-9-29 16:25:24

额发错平台了

lqouc · 发表于 2014-9-29 16:48:26

do i=1,po
do j=1,n
do k=1,p
y(i,j)=y(i,j)+g(i,k)*u(k,j)
或者直接y=matmul(g,u)
下次这样的问题请自己搜索下，或者看看书，这根本就是一个线性代数的问题而不是编程问题。

心语 · 发表于 2014-9-29 16:56:37

lqouc 发表于 2014-9-29 16:48
do i=1,po
do j=1,n
do k=1,p

谢谢哈不过我不是求y矩阵，我想求g矩阵

lqouc · 发表于 2014-9-29 17:00:24

心语发表于 2014-9-29 16:56
谢谢哈不过我不是求y矩阵，我想求g矩阵

嗯嗯，没看清楚，那就更是一个现代的问题了，你还是可以经过逆矩阵变化什么的最终转化成为一个左边为g的矩阵乘法算式，然后再用上面的方法。至于转置或者求逆这个都是可以用差不多的循环或者fortran内置的函数实现。

心语 · 发表于 2014-9-29 17:04:49

lqouc 发表于 2014-9-29 17:00
嗯嗯，没看清楚，那就更是一个现代的问题了，你还是可以经过逆矩阵变化什么的最终转化成为一个左边为g的 ...

好哒，谢谢

		自动登录	找回密码
密码			立即注册

请问这个矩阵求解问题 怎么用fortran实现？