SOGO論壇

標題: 求助~請問fortran要如何將原來矩陣內插成新矩陣 [列印本頁]

作者: seedhank 時間: 2014-3-7 08:34:27 標題: 求助~請問fortran要如何將原來矩陣內插成新矩陣

請教各位大大,小的現在遇到的問題是我要在一個[60*120]的A矩陣中對兩點進行線性內插,內插點數的方式是在A(1,j)和A(1,j+1)中間內插出8個點存到另一B矩陣,B矩陣大小為[60*1200]
以下是我現在的寫法：
DO i=1,60
DO j=1,120
      DO k=1,1200,10
         xxx=(base2(i,j)-base2(1,j+1))/9
         baseall(i,k)=base2(i,j)
         baseall(i,k+1)=base2(i,j) + xxx
         baseall(i,k+2)=base2(i,j) + (2* xxx)
         baseall(i,k+3)=base2(i,j) + (3* xxx)
         baseall(i,k+4)=base2(i,j) + (4* xxx)
         baseall(i,k+5)=base2(i,j) + (5* xxx)
         baseall(i,k+6)=base2(i,j) + (6* xxx)
         baseall(i,k+7)=base2(i,j) + (7* xxx)
         baseall(i,k+8)=base2(i,j) + (8* xxx)
         baseall(i,k+9)=base2(i,j+1)
      END DO
END DO
END DO

小的的寫法造成一直重複地計算錯誤答案,還請各位大大堤點小的要如何善用不同語法將式子改寫,麻煩各位大大提點!!謝謝各位大大!!

歡迎光臨 SOGO論壇 (https://oursogo.com/)