地震子波波形显示及一维地震合成记录制作

最新推荐文章于 2023-10-08 16:14:06 发布

张三風儿

最新推荐文章于 2023-10-08 16:14:06 发布

阅读量2k

点赞数 5

文章标签：地震勘探多次波地震成像压制和消除

本文链接：https://blog.csdn.net/CSDN_fairy/article/details/107095149

版权

//最小相位子波
#include <stdio.h>
#include <math.h>
#include <malloc.h>
#define PI 3.1415962

void waveletfun(int nt_wavelet,float dt,float *wavelet,float fm,float m1m2);
void reflectfun(int nt_record,float *Vp,float *des,float *reflect);
void convolutionfun(int nt_record, int nt_wavelet, float *record, float *reflect, float *wavelet);

// 主程序部分
void main()
{
   
	// 变量的定义
	int nt_wavelet=150;
	float *wavelet;
	float dt=0.002;
	int it;
	int j;
	float fm=10;
	float m1m2=2.0;
	int nt_record = 600;
	float *reflect;
	float *des;
	float *Vp;
	float *record=0;

	// 定义文件
	FILE *waveletfp;
	waveletfp = fopen ("minphase_wavelet.txt","w");
	FILE *reflectfp;
	reflectfp = fopen ("minphase_reflect.txt","w");
	FILE *recordfp;
	recordfp = fopen ("minphase_record.txt&