24小时热门版块排行榜    

查看: 1761  |  回复: 9

wysqh

木虫 (小有名气)


[交流] 【求助】matlab数据导入Fortran

目前有一个用matlab产生的一个128*128特定矩阵 存为data.dat文件  现在想要将这个矩阵导入到Fortran中存为eta(i,j) 行列要一一对应 该如何实现呢?
回复此楼

» 猜你喜欢

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★
微尘、梦想(金币+2): 谢谢回贴…… 2011-04-11 11:30:55
你的 data.dat 是个文本文件?

直接读入不就可以了吗?也不是啥难事儿……

行列对应无非就是循环的下标排列的问题。

给个列子说明要如何一一对应……
2楼2011-04-11 10:57:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
如果没有记错,matlab 与 Fortran 的矩阵存贮格式是一致的。
3楼2011-04-11 10:57:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wysqh

木虫 (小有名气)


引用回帖:
Originally posted by snoopyzhao at 2011-04-11 10:57:13:
你的 data.dat 是个文本文件?

直接读入不就可以了吗?也不是啥难事儿……

行列对应无非就是循环的下标排列的问题。

给个列子说明要如何一一对应……

data.dat文件 是按照一行一行的存的 如附件
就想实现 eta(i,j) 每一行与文件里面的每一行相对应 也是128*128排列
是采用read语句实现吗?具体程序该怎么写啊。。

[ Last edited by wysqh on 2011-4-11 at 11:34 ]
4楼2011-04-11 11:33:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingxb

金虫 (正式写手)


wysqh(金币+3): 多谢交流~ 2011-04-12 08:27:27
do i = 1, 128
do j = 1, 128
  read(file, *) eta(i,j)
end do
end do

基本结构如此,如果需要指定格式的话,替换read语句中的*为格式描述。
5楼2011-04-11 14:17:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingxb

金虫 (正式写手)


为了确保读的正确,一般我都是输出进行比较。
6楼2011-04-11 14:18:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
引用回帖:
Originally posted by dingxb at 2011-04-11 14:17:55:
do i = 1, 128
do j = 1, 128
  read(file, *) eta(i,j)
end do
end do

基本结构如此,如果需要指定格式的话,替换read语句中的*为格式描述。

其实一句话就够了,呵呵……
CODE:
read(unit=5, *) ((eta(i,j), j = 1,128),i=1,128)

7楼2011-04-11 17:50:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wysqh

木虫 (小有名气)


引用回帖:
Originally posted by dingxb at 2011-04-11 14:17:55:
do i = 1, 128
do j = 1, 128
  read(file, *) eta(i,j)
end do
end do

基本结构如此,如果需要指定格式的话,替换read语句中的*为格式描述。

我开始是这样做的 但是运行之后会出现下面的错误

8楼2011-04-11 22:10:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wysqh

木虫 (小有名气)


引用回帖:
Originally posted by snoopyzhao at 2011-04-11 17:50:29:
其实一句话就够了,呵呵……
CODE:
read(unit=5, *) ((eta(i,j), j = 1,128),i=1,128)


read(unit=5, *) ((eta(i,j), j = 1,128),i=1,128)
是不是少了一个括号?
9楼2011-04-11 22:14:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
wysqh(金币+7): 谢谢你的建议~因为不是专门搞这个的 很多东西不懂 2011-04-12 08:27:11
引用回帖:
Originally posted by wysqh at 2011-04-11 22:14:10:
read(unit=5, *) ((eta(i,j), j = 1,128),i=1,128)
是不是少了一个括号?

数括号我还是会的,我不知道你咋看出来少了一个括号,呵呵……

下面是我测试过的代码:
CODE:
program readata
real(kind=kind(1.0d0)), dimension(128,128) :: mat
integer :: i,j

read(unit=5,fmt=*) ((mat(i,j), j = 1,128),i=1,128)

!do i=1,128
!   write(*,*) (mat(i,j), j = 1,128)
!end do

end program readata

10楼2011-04-11 23:16:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wysqh 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 070300化学求调剂 +3 苑豆豆 2026-03-20 3/150 2026-03-22 14:48 by ColorlessPI
[考研] 生物学调剂 +5 Surekei 2026-03-21 5/250 2026-03-22 14:39 by tcx007
[考研] 求调剂 +7 Auroracx 2026-03-22 7/350 2026-03-22 12:38 by 素颜倾城1988
[考研] 384求调剂 +3 子系博 2026-03-22 4/200 2026-03-22 11:04 by 搏击518
[考博] 招收博士1-2人 +3 QGZDSYS 2026-03-18 4/200 2026-03-22 10:25 by QGZDSYS
[考研] 生物学一志愿985,分数349求调剂 +4 zxts12 2026-03-21 7/350 2026-03-22 09:57 by zxts12
[考研] 0703化学297求调剂 +3 Daisy☆ 2026-03-20 3/150 2026-03-21 17:45 by ColorlessPI
[考研] 307求调剂 +3 余意卿 2026-03-18 3/150 2026-03-21 17:31 by ColorlessPI
[考研] 0805材料320求调剂 +3 深海物语 2026-03-20 3/150 2026-03-21 15:46 by 无际的草原
[考研] 求调剂 +3 白QF 2026-03-21 3/150 2026-03-21 13:12 by zhukairuo
[考研] 265求调剂 +3 Jack?k?y 2026-03-17 3/150 2026-03-21 03:17 by JourneyLucky
[考研] 初始318分求调剂(有工作经验) +3 1911236844 2026-03-17 3/150 2026-03-21 02:33 by JourneyLucky
[考研] 材料 336 求调剂 +3 An@. 2026-03-18 4/200 2026-03-21 01:39 by JourneyLucky
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 一志愿重庆大学085700资源与环境专硕,总分308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-21 00:39 by JourneyLucky
[考研] 304求调剂 +7 司空. 2026-03-18 7/350 2026-03-20 23:08 by JourneyLucky
[考研] 323求调剂 +3 洼小桶 2026-03-18 3/150 2026-03-20 22:54 by JourneyLucky
[考研] 一志愿武汉理工材料工程专硕调剂 +9 Doleres 2026-03-19 9/450 2026-03-20 22:36 by JourneyLucky
[考研] 材料与化工求调剂 +7 为学666 2026-03-16 7/350 2026-03-19 14:48 by 尽舜尧1
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
信息提示
请填处理意见