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/) Powered by OURSOGO.COM