第十二讲41傅里叶伪谱法12.pptxVIP

1.确认修课名单； 2.作业在下周四前交齐，可以重新提交； 3. 考试时间：11月28日下午3：20-5：20； 4. 答疑地点：科研楼东侧403室 5. 考前布置一个大作业，不用提交。;傅里叶伪谱法;有限差分法：直接从微分方程出发，用差商代替微商， 从而将微分方程直接转化为代数方程。;傅里叶伪谱法将波场函数表示为傅里叶级数的展开形式， 空间域的求导转化为波数域的乘积计算， 时间域的求导用二阶差分格式给出。 （Kosloff et al., 1984）;考虑偏微分方程和边界条件:;一般把偏微分方程形式称为强形式(strong form)，由于以上微分方程在U内每一点都成立，因此总可以找到权函数w，使得下式成立：;谱方法;谱方法; 用配点法求解偏微分方程组;Fourier伪谱方法;快速Fourier变换(Fast Fourier transform，FFT)： 例：乘积运算次数 网格节点数 全矩阵 向量乘积FFT 比值 1D(nx=512) 2.6x105 9.2x103 28.4 1D(nx=2096) 4.4x106 4.6x104 94.98 1D(nx=8384) 7.0x107 2.2x105 312.6 全矩阵-向量相乘的运算次数与FFT运算次数的比值可以看作使用 FFT时的加速比。;李本文;李本文;function df=sderld(f,dx) %SDERLD(f,dx) spectral derivative of vector nx=max(size(f)); % initialize k kmax=pi/dx; dk=kmax/(nx/2); for i=1:nx/2; k(i)=i*dk; k(nx/2+i)=-kmax+i*dk; end k=sqrt(-1)*k; %FFT and IFFT ff=fft(f); ff=k.*ff; df=real(ifft(ff)); ;吴宝年等，2012;吴宝年等，2012;吴宝年等，2012;吴宝年等，2012;吴宝年等，2012;吴宝年等，2012;吴宝年等，2012;吴宝年等，2012;龙桂华等，2009;龙桂华等，2009;在数值计算过程中，介质沿水平和纵向分割成255×255个网格，相邻格点之间的间距为10m．震源位于地表下190m处，且位于模型水平方向的正中心，为主频25的Ricker子波，其中子波的时间采样间隔为1ms.;龙桂华等，2009;龙桂华等，2009;利用伪谱法模拟横向非均匀全球模型中的SH 波场;王彦宾, Hiroshi TAKENAKA, 2012;王彦宾, Hiroshi TAKENAKA, 2012;王彦宾, Hiroshi TAKENAKA, 2012;王彦宾, Hiroshi TAKENAKA, 2012;周期性边界条件（Periodic Boundary Conditions, PBC）是边界条件的一种，反映的是如何利用边界条件替代所选部分（系统）受到周边（环境）的影响。可以看作是如果去掉周边环境，保持该系统不变应该附加的条件，也可以看作是由部分的性质来推广表达全局的性质。