计算声学第一课.ppt

计 算 声 学 马 俊 吉林大学物理学院声学与微波物理系 第二章? 纯解析算法----全波场的计算? §2.1 单相准弹性介质地层下单极点源激发的井孔声波全波场的计算? 本节学习在单相（准）弹性介质流体中点源（单极源）激发与辐射声场的数值模拟方法----离散波数法。 声源脉冲的类型 声源脉冲的类型 PF3余弦包络 时域函数 1 声源脉冲的类型 PF3 声源脉冲的类型 PF2变形瑞克子波 时域函数 2 声源脉冲的类型 PF2 快速傅立叶变换子程序 FFT 快速傅立叶变换子程序 SUBROUTINE COMPFR(A,B,N1,M1,INV) A－－－－－－变换量实部数组 B－－－－－－变换量虚部数组 N1,M1－－－－变换量数组元素量及相关量 INV－－－－逻辑变量（.TRUE. OR .FALSE.） 正变换 反变换 A, B既是输入变量又是输出变量 快速傅立叶变换子程序组成及调用方式 ★ SUBROUTINE COMPFR(XR,XI,N2,M2,.INV.) SUBROUTINE FFT(A,B,N1,M1,KS) SUBROUTINE REORD(A,B,