计算机习题2.doc

計算機習題2 2004.3.24 2.Abel變換。使用計算機習題1得到的33層高度處的折射率N按Abel變換計算出33層的彎角： (1) 然後再做逆變換，決定出33層的折射指數： (2) 在(1)和(2)式中已設處（在標準大氣資料中已設為6471km，即地球平均半徑=6371km加上100km）的折射指數。若這裡的，則(1)和(2)式的上限應用。為方便起見，(1)和(2)式可改寫為 (3) (4) 其中 (5) 至於x現在則是啞變數，(3)到(5)式是相當適合機器計算的形式。 (3)和(4)式所示的積分式可用IMSL副程式DQDAWS或DQDAGS計算。在計算(3)式右邊時，由於並不是以x的連續函數形式出現，故必須先用三次樣條函數進行內插，然後再計算的一次導數。在計算(4)式時，同樣的要將(3)式計算出來的值做三次樣條函數內插。 因此在本題中需用到3個副程式： (a) 一個副程式由標準大氣資料計算出N和a，並送進主程式(這是第10題的內容)。 (b) 一個副程式處理(3)式中的被積項，在這副程式中需決定出，並且計算DQDAWS或DQDAGS所需的函數。 (c) 一個副程式處理(4)式中的被積項，也就是決定和計算DQDAWS或D