版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

北京石油化工学院2026年研究生招生接收调剂公告

返回列表

kimileegdut

捐助贵宾 (小有名气)

应助: 1 (幼儿园)
帖子: 56
在线: 24小时
虫号: 3526305
注册: 2014-11-08

[求助] 奇异值分解法求广义逆代码求助已有1人参与

以下是从《FORTRAN常用算法程序集》里面找到的求广义逆的奇异值分解法子程序代码，但是用例题测试之后发现计算有误，想请教下大神是不是子程序有误？
SUBROUTINE BGINV(M,N,A,AA,L,EPS,U,V,KA,S,E,WORK)
DIMENSION A(M,N),U(M,M),V(N,N),AA(N,M)
DIMENSION S(KA),E(KA),WORK(KA)
DOUBLE PRECISION A,U,V,AA,S,E,WORK
CALL BMUAV(A,M,N,U,V,L,EPS,KA,S,E,WORK)
IF(L.EQ.0) THEN
   K=1
10 IF(A(K,K).NE.0.0) THEN
      K=K+1
      IF(K.LE.MIN(M,N)) GO TO 10
   END IF
   K=K-1
   IF(K.NE.0) THEN
      DO 40 I=1,N
      DO 40 J=1,M
      AA(I,J)=0.0
      DO 30 II=1,K
30       AA(I,J)=A(I,J)+V(II,I)*U(J,II)/A(II,II)
40    CONTINUE
   END IF
END IF
RETURN
END

回复此楼

» 猜你喜欢

化工调剂303分，过四级已经有5人回复
071000生物学调剂已经有3人回复
一志愿华东理工大学，080500学硕，317分，求调剂已经有10人回复
327求调剂已经有9人回复
0856初试324分求调剂已经有6人回复
08开头看过来！！！已经有6人回复
材料调剂已经有7人回复
一志愿北京科技大学085601材料工程英一数二初试总分335求调剂已经有8人回复
085900土木水利336分求调剂已经有6人回复
319求调剂已经有13人回复

1楼 2015-07-21 22:10:45

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็

应助: 237 (大学生)
金币: 6101.9
散金: 1114
红花: 39
帖子: 1380
在线: 553.8小时
虫号: 506699
注册: 2008-02-18
专业: 模式识别

【答案】应助回帖

这里的源代码必然适合你

http://download.csdn.net/detail/lcfactorization/9022013

赞一下

回复此楼

ส็็็็็็็็็็็็็็็็็็็็

2楼2015-08-19 10:02:30