wysqh

木蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 4010
帖子: 250
在線: 28.4小時(shí)
蟲號(hào): 884780

[交流] 【求助】matlab數(shù)據(jù)導(dǎo)入Fortran

目前有一個(gè)用matlab產(chǎn)生的一個(gè)128*128特定矩陣存為data.dat文件現(xiàn)在想要將這個(gè)矩陣導(dǎo)入到Fortran中存為eta(i,j) 行列要一一對(duì)應(yīng) 該如何實(shí)現(xiàn)呢？

回復(fù)此樓

» 猜你喜歡

267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑已經(jīng)有7人回復(fù)
一志愿華中科技大學(xué)，080502，354分求調(diào)劑已經(jīng)有3人回復(fù)
一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑已經(jīng)有4人回復(fù)
一志愿武理材料305分求調(diào)劑已經(jīng)有6人回復(fù)
材料專碩英一數(shù)二306 已經(jīng)有3人回復(fù)
085600材料與化工調(diào)劑 324分已經(jīng)有7人回復(fù)
311求調(diào)劑已經(jīng)有3人回復(fù)
26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線已經(jīng)有7人回復(fù)
295求調(diào)劑已經(jīng)有5人回復(fù)
354求調(diào)劑已經(jīng)有5人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2011-04-11 10:31:09

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

snoopyzhao

至尊木蟲 (職業(yè)作家)

★ ★
微塵、夢想(金幣+2): 謝謝回貼…… 2011-04-11 11:30:55

你的 data.dat 是個(gè)文本文件？

直接讀入不就可以了嗎？也不是啥難事兒……

行列對(duì)應(yīng)無非就是循環(huán)的下標(biāo)排列的問題。

給個(gè)列子說明要如何一一對(duì)應(yīng)……

贊一下(1人)

回復(fù)此樓

2樓2011-04-11 10:57:13

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

snoopyzhao

至尊木蟲 (職業(yè)作家)

如果沒有記錯(cuò)，matlab 與 Fortran 的矩陣存貯格式是一致的。

贊一下

回復(fù)此樓

3樓2011-04-11 10:57:41

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wysqh

木蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 4010
帖子: 250
在線: 28.4小時(shí)
蟲號(hào): 884780

引用回帖:

Originally posted by snoopyzhao at 2011-04-11 10:57:13:
你的 data.dat 是個(gè)文本文件？

直接讀入不就可以了嗎？也不是啥難事兒……

行列對(duì)應(yīng)無非就是循環(huán)的下標(biāo)排列的問題。

給個(gè)列子說明要如何一一對(duì)應(yīng)……

data.dat文件是按照一行一行的存的如附件
就想實(shí)現(xiàn) eta(i,j) 每一行與文件里面的每一行相對(duì)應(yīng) 也是128*128排列
是采用read語句實(shí)現(xiàn)嗎？具體程序該怎么寫啊。。

[ Last edited by wysqh on 2011-4-11 at 11:34 ]

贊一下

回復(fù)此樓

4樓2011-04-11 11:33:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingxb

金蟲 (正式寫手)

應(yīng)助: 12 (小學(xué)生)
金幣: 1674.5
帖子: 482
在線: 59.5小時(shí)
蟲號(hào): 34085

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

基本結(jié)構(gòu)如此，如果需要指定格式的話，替換read語句中的*為格式描述。

贊一下

回復(fù)此樓

5樓2011-04-11 14:17:55

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingxb

金蟲 (正式寫手)

應(yīng)助: 12 (小學(xué)生)
金幣: 1674.5
帖子: 482
在線: 59.5小時(shí)
蟲號(hào): 34085

為了確保讀的正確，一般我都是輸出進(jìn)行比較。

贊一下

回復(fù)此樓

6樓2011-04-11 14:18:53

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

snoopyzhao

至尊木蟲 (職業(yè)作家)

引用回帖:

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

基本結(jié)構(gòu)如此，如果需要指定格式的話，替換read語句中的*為格式描述。

其實(shí)一句話就夠了，呵呵……

CODE:

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

贊一下

回復(fù)此樓

7樓2011-04-11 17:50:29

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wysqh

木蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 4010
帖子: 250
在線: 28.4小時(shí)
蟲號(hào): 884780

引用回帖:

我開始是這樣做的但是運(yùn)行之后會(huì)出現(xiàn)下面的錯(cuò)誤

贊一下

回復(fù)此樓

8樓2011-04-11 22:10:37

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wysqh

木蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 4010
帖子: 250
在線: 28.4小時(shí)
蟲號(hào): 884780

引用回帖:

Originally posted by snoopyzhao at 2011-04-11 17:50:29:
其實(shí)一句話就夠了，呵呵……

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)
是不是少了一個(gè)括號(hào)？

贊一下

回復(fù)此樓

9樓2011-04-11 22:14:10

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

snoopyzhao

至尊木蟲 (職業(yè)作家)

wysqh(金幣+7): 謝謝你的建議~因?yàn)椴皇菍ｉT搞這個(gè)的很多東西不懂 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)
是不是少了一個(gè)括號(hào)？

數(shù)括號(hào)我還是會(huì)的，我不知道你咋看出來少了一個(gè)括號(hào)，呵呵……

下面是我測試過的代碼：

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

贊一下

回復(fù)此樓

10樓2011-04-11 23:16:05

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 354求調(diào)劑 +3	Tyoumou 2026-03-18	5/250	2026-03-18 16:41 by Tyoumou
[考研] 266求調(diào)劑 +5	陽陽哇塞 2026-03-14	9/450	2026-03-18 15:05 by stone_128
[考研] 311求調(diào)劑 +11	冬十三 2026-03-15	12/600	2026-03-18 14:36 by 星空星月
[考博] 環(huán)境領(lǐng)域全國重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3	QGZDSYS 2026-03-13	5/250	2026-03-18 11:13 by QGZDSYS
[考研] 材料與化工求調(diào)劑 +6	為學(xué)666 2026-03-16	6/300	2026-03-17 20:15 by peike
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[考研] 302求調(diào)劑 +9	負(fù)心者當(dāng)誅 2026-03-11	9/450	2026-03-17 17:13 by ruiyingmiao
[考研] 【0856】化學(xué)工程（085602）313 分，本科學(xué)科評(píng)估A類院�；瘜W(xué)工程與工藝，誠求調(diào)劑 +7	小劉快快上岸 2026-03-11	8/400	2026-03-17 16:57 by ruiyingmiao
[考研] 梁成偉老師課題組歡迎你的加入 +8	一鴨鴨喲 2026-03-14	10/500	2026-03-17 15:07 by 一鴨鴨喲
[考研] 藥學(xué)383 求調(diào)劑 +3	藥學(xué)chy 2026-03-15	4/200	2026-03-16 20:51 by 元子^0^
[基金申請(qǐng)] 今年的國基金是打分制嗎？ 50+3	zhanghaozhu 2026-03-14	3/150	2026-03-16 17:07 by 北京萊茵潤色
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染
[考研] 中科大材料專碩319求調(diào)劑 +3	孟鑫材料 2026-03-13	3/150	2026-03-14 18:10 by houyaoxu
[考研] 復(fù)試調(diào)劑 +4	z1z2z3879 2026-03-14	5/250	2026-03-14 16:30 by JourneyLucky
[考研] 招收0805（材料）調(diào)劑 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、
[考研] 材料專碩350 求調(diào)劑 +4	王金科 2026-03-12	4/200	2026-03-13 16:02 by ruiyingmiao
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8	努力奮斗112 2026-03-12	9/450	2026-03-13 15:41 by JourneyLucky
[考研] 求調(diào)劑 +3	程雨杭 2026-03-12	3/150	2026-03-13 15:06 by JourneyLucky
[考研] 085600材料與化工 309分請(qǐng)求調(diào)劑 +7	dtdxzxx 2026-03-12	8/400	2026-03-13 14:43 by jxchenghu
[考研] 283求調(diào)劑，材料、化工皆可 +8	蘇打水7777 2026-03-11	10/500	2026-03-13 09:06 by Linda Hu