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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 如何读入以下的矩阵
51/1返回列表
查看: 1600  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

kathy2008

木虫 (正式写手)

[求助] 如何读入以下的矩阵

如题，请指点。不像一般的下三角矩阵。请指点。谢谢
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-05-02 10:40:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)
下载不下来，呵呵……

如果你知道矩阵的行数或列数，自己计算一下其中的数据个数，然后把它们全部读到一个一维数组中，或者直接到二维数组也可以……

类似于这样……

read(*,*) (a(i), i = 1, 100)

或者

(read(*,*) ((a(i,j), i = 1,10), j=i,10)

仅供参考……
赞一下 回复此楼
3楼2011-05-02 11:13:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

kathy2008

木虫 (正式写手)
仔细一看，是一个下三角矩阵，只是按5个一组来输出，请高手指点。
赞一下 回复此楼
2楼2011-05-02 10:58:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kathy2008

木虫 (正式写手)
引用回帖:
Originally posted by snoopyzhao at 2011-05-02 11:13:46:
下载不下来，呵呵……

如果你知道矩阵的行数或列数，自己计算一下其中的数据个数，然后把它们全部读到一个一维数组中，或者直接到二维数组也可以……

类似于这样……

read(*,*) (a(i), i = 1, 100)

...

有88行，有88列，是一个88*88的下三角矩阵。
赞一下 回复此楼
4楼2011-05-02 15:29:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

【答案】应助回帖

kathy2008(金币+10): 神速！！！！ 2011-05-02 15:36:51
比我想象的复杂，下面是一个实现，仅供参考……
CODE:
program overlap

real, dimension(88,88) :: ol
character(len=80) :: line,lfm
integer :: i, j, k, ios
integer :: idat = 12

   open(unit=idat, file='overlap.txt')

   do i=1,88
      do j=1,88
         ol(i,j) = 0.0
      end do
   end do

   do i=1,88,5
      read(idat,'(a)', iostat = ios ) line
      if (ios < 0) exit
      read(idat,'(a)') line
      do j=i,88
         read(idat,'(a)') line
         write(lfm, '(a,i0,a)') '(tr1,',min(i+4,j),'(tr2,f9.6))'
         read(line,lfm) (ol(j,k),k=i,min(i+4,j))
!        write(*,lfm) (ol(j,k),k=i,min(i+4,j))  
      end do
      read(idat,'(a)', iostat = ios) line
      if (ios < 0) exit
   end do

end program overlap
赞一下 回复此楼
5楼2011-05-02 15:30:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭