| 查看: 1073 | 回复: 6 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
想用Fortran语言实现矩阵连乘，并且矩阵里面有一个变量 已有3人参与
|
|||
|
想用Fortran实现对矩阵M的连乘，而且矩阵里面含有一个变量EM，当EM＜=15，EM=6.7. 当EM＞15，EM=8.7，但是用以下程序不能实现。希望大家帮忙解决，不胜感激！ PROGRAM MAIN USE IMSL COMPLEX:: M(2,2),P(2,2) REAL*8 E,TE,TM,EM,TD,H,N COMPLEX*8 M11,M12,M21,M22,A1,A2,A3 TD=0.8D0 TM=1.0D0 EM=7.75D0 N=10 C====================================================== OPEN(UNIT=30,FILE="DNA-DAN01.DAT"  DO E=6.89,9.750,0.001D0 C================================================ cc================================================ P=M DO I=1,N IF (I .LT. 15) THEN EM=6.7 ELSE IF (I .GE. 15) THEN EM=8.7 END IF C==================================================== H=acos((E-EM)/(2*TM)) M11=CMPLX((E-EM)/TD,0) M12=CMPLX(-1,0) M21=CMPLX(1,0) M22=CMPLX(0,0) M=RESHAPE ((/M11,M21,M12,M22/),(/2,2/)) P=MATMUL(M,P) END DO C===================================================== A1=4*SIN(H)*SIN(H) A2=(P(1,2)-P(2,1)+(P(1,1)-P(2,2))*COS(H))**2 A3=(P(1,1)+P(2,2))**2*SIN(H)**2 TE=A1/(A2+A3) Write(30,*) E,TE END DO END |
回复此楼» 猜你喜欢
求个博导看看
已经有19人回复
-
论文终于录用啦！满足毕业条件了
已经有13人回复
-
2025年遐想
已经有5人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- fortran求矩阵本征值出错
已经有6人回复
- fortran中如何创建由二维矩阵为元素的1维矩阵
已经有4人回复
- 急：求用Fortran或者其他语言编写一个程序，处理下面的数据！！！！！！！！！！
已经有37人回复
- Fortran怎么从一个dat文件中读取不同类型数据赋值给不同的变量？
已经有8人回复
- 请问在fortran里面怎么求复数矩阵的行列式！！！
已经有4人回复
- Fortran程序，如何实现采用矩阵保存单链表表头？
已经有6人回复
- fortran中对变量取值
已经有3人回复
- fortran到底怎么弄稀疏矩阵啊？
已经有8人回复
- matlab如何实现矩阵多行数据的循环调用
已经有5人回复
- matlab中如何通过for语句得到已有矩阵的部分数据
已经有3人回复
- fortran如何实现张量相乘（并乘）
已经有9人回复
- 用fortran存储稀疏矩阵
已经有4人回复
- fortran能求带字母矩阵的本征值吗？（矩阵元素带有字母，能得到形式解吗）
已经有5人回复
- fortran程序怎样调用lapack库进行矩阵对角化？
已经有7人回复
- 请问，我想用fortran计算统计中的p值以及95%信度空间
已经有4人回复
- 【求助】fortran求解矩阵
已经有7人回复
- 【求助】fortran中是如何实现将两个变量的值在同一行进行打印出来的
已经有7人回复
- 【求助】矩阵相乘
已经有8人回复
- 【求助】下面这个公式怎么用FORTRAN实现
已经有7人回复
- 【求助】fortran：如何求解 复数矩阵的秩【已完结】
已经有5人回复
- 【求助】求指点，在Fortran里面，怎样实现“数组维度可调”的数组?
已经有5人回复
pippi6
铁杆木虫 (著名写手)
工程和科学数值计算咨询
- 应助: 413 (硕士)
- 贵宾: 0.002
- 金币: 7116.5
- 散金: 15
- 红花: 63
- 帖子: 1639
- 在线: 798.9小时
- 虫号: 2469437
- 注册: 2013-05-14
- 专业: 计算数学与科学工程计算
3楼2014-11-13 20:22:16