声波有限差分法正演模拟c语言程序(共5页).docx

精选优质文档-倾情为你奉上#include#include#define fm 30#define dt 0.001#define PI 3.#define Nt 401#define Nx 200#define Nz 200/-加载震源，雷克子波-void fun(float source)FILE *fp;int it,i;float t1,t2,t0;for(i=0;iNt;i+) sourcei=0.0;t0 = 1.0/fm;printf(t0 = %fn, t0); fp=fopen(ricker.txt,w);for(it=0;itNt;it+)t1=it*dt;t2=t1-t0;sourceit=(1.0-2.0*pow(PI*fm*t1,2.0)*exp(