油气勘探新方法报告

《油气勘探新方法报告》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、油气勘探新方法实验报告指导老师：文晓涛班级：勘查4班姓名：徐明学号：200905060412实验一-----时间域希尔伯特变换一实验目的：掌握希尔伯特变换的特点，熟悉其编程，具体了解希尔伯特变换的应用，明白其原理。、二实验原理：对于实信号x(t)，定义与之对应的复数信号h(t)（1）其中，y(t)为复数信号的虚部，它定义为实信号x(t)的希尔伯特变换（HT）（2）其中,为希尔伯特变换滤波器，即实信号x(t)的希尔伯特变换是实信号与HT滤波器的卷积。三实现步骤：1,确定希尔伯特变换算子的长度，设为N，按下式计算指定长度为N的希尔伯特变换算子：（3）2,希尔伯特算子h(n)与原

2、始信号序列x(n)进行时域卷积，得到变换后的结果：（4）四实验程序如下程序运行环境：VC++-----------------------------------------------------------------------------------------#include#definenn200#definem63inti,j;doublex[nn],y[nn],sum,A[nn],Q[nn],W[nn];voidmain(){FILE*fp1,*fp2,*fp3,*fp4;fp1=fopen("希尔伯特.txt","w");fp2=fopen

3、("瞬时振幅.txt","w");fp3=fopen("瞬时相位.txt","w");fp4=fopen("瞬时频率.txt","w");for(i=0;i0)&&(i-(2*(j-(m-1)/2)+1)>0)&&(i+(2*(j-(m-1)/2)+1)<=nn)&&(i-(2*(j-(m-1)/2)+1)<=nn)){sum=sum+x[i-2*(j-(m-1)/

4、2)-1]/(2*(j-(m-1)/2)+1);}}y[i]=2./3.1415926*sum;fprintf(fp1,"%dt%ft%f",i,x[i],y[i]);}///////////////////////////////////for(i=1;i

5、//for(i=1;i

6、----频率域希尔伯特变换一实验目的：掌握频率域希尔伯特变换的特点；熟悉频率域希尔伯特变换的程序的编写；并能熟练掌握频率域希尔伯特变换的原理。二实验程序：---------------------------------------------------------------------------------------程序运行环境为VC++#include#include#include#definepi3.1415926#definexn128voidInverse(float*X,intn);voidfft(

7、float*Xr,float*Xi,intn,intinv);voidmain(){inti;floatdat2[xn];floatdat_real[xn],dat_image[xn],temp[xn];FILE*fp1;//////////////////////////////数组初始化for(i=0;i