|
本帖最后由 jkl1102925008 于 2022-6-23 21:13 编辑
存在待求逆矩阵 a=
1.46989040000000e-09 8.08439760000000e-10 -9.40729890000000e-09
8.08439760000000e-10 4.44641850000000e-10 -5.17401410000000e-09
-9.40729890000000e-09 -5.17401410000000e-09 6.02067090000000e-08
subroutine Kalman_Gain_calc
use lapack95
real(kind=8) :: a(3,3),b(3)
integer :: ipiv(3)
integer :: n,info1,info2
b=0.0
ipiv=0
a=0.0
n=3
a=st_sim_cov
call sgetrf(n,n,a,n,ipiv,info1)
call sgetri(n,a,n,ipiv,b,n,info2)
end subroutine
计算结果 a=
54127.1604943661 2.562316454118586E-315 -9.407298939834163E-009
22340.6659669089 4.446417947376324E-010 -5.174014106046343E-009
0.000000000000000E+000 -5.174014106046343E-009 6.020670895168223E-008
而MATLAB INV(a) 计算结果为
-6.61095895614289e+15 2.60655716871615e+15 -808961396190492
2.60655716871615e+15 1.05500250710316e+17 9.47370245921866e+15
-808961396190492 9.47370245921866e+15 687746102360067
大佬们帮我看看,是我哪里搞错了吗,问题出在哪里,求求了,嘤嘤嘤!
-----------------------------------------------------------------------------------------
使用sgetrf+sgetri计算的结果是错误的。。。
|
|