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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 如何读入以下的矩阵
81/1返回列表
查看: 1593  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

kathy2008

木虫 (正式写手)

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

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

» 猜你喜欢

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

kathy2008

木虫 (正式写手)
仔细一看，是一个下三角矩阵，只是按5个一组来输出，请高手指点。
赞一下 回复此楼
2楼2011-05-02 10:58:51
已阅   回复此楼   关注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的回帖

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的回帖

snoopyzhao

至尊木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
微尘、梦想(金币+5): 谢谢回复！ 2011-05-02 20:30:56
刚才那个里面关于 lfm 的赋值有点儿问题，尽管不影响使用，呵呵……下面这个是改正后的。
CODE:
program overlap

real, dimension(88,88) :: ol
character(len=80) :: line,lfm
integer :: i, j, k, l, 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
      l=0
      do j=i,88
         if (l < 5) l=l+1
         read(idat,'(a)') line
         write(lfm, '(a,i0,a)') '(tr1,',l,'(tr2,f9.6))'
!        print *, lfm, min(i+4,j)
         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
赞一下(1人) 回复此楼
6楼2011-05-02 15:42:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kathy2008

木虫 (正式写手)
引用回帖:
Originally posted by snoopyzhao at 2011-05-02 15:42:17:
刚才那个里面关于 lfm 的赋值有点儿问题，尽管不影响使用，呵呵……下面这个是改正后的。

[code]
program overlap

real, dimension(88,88) :: ol
character(len=80) :: line,lfm
integer :: i, j, k, l, ...

继续请教: idat = 12，是怎么样
确定把12付给idat这个变量的呢？也就是说为什么要选择12，而不是其他数字，谢谢。
赞一下 回复此楼
7楼2011-05-02 15:57:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)
引用回帖:
Originally posted by kathy2008 at 2011-05-02 15:57:54:
继续请教: idat = 12，是怎么样
确定把12付给idat这个变量的呢？也就是说为什么要选择12，而不是其他数字，谢谢。

其它数字也可以啊，这只是 unit 的一个编号，呵呵……一般大于 10 会比较安全，至于写成 idat，是方便修改……
赞一下 回复此楼
8楼2011-05-02 16:00:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kathy2008 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭