15度 有限差分 叠前偏移程序

#include
#include

#define N 128
#define M 256

#define  V 2000//传播速度
#define T 0.002//时间间隔
#define T1 0.004//深度域转换为时间域间隔
#define  X 7//x轴间隔



void main（）
{
	printf（“正确\n“）;
	float lvboyinzi（float *yz）;
	void zheji（float*yzfloat*jieguofloat*z1）;//Z1[]即为Fi
	FILE*out*out1;
	fopen_s（&out“数据.txt““w+“）;
	fopen_s（&out1“shuju.txt““w+“）;
	//out=fopen（“数据.txt““w+“）;//打开可读写文件
	//out1=fopen（“shuju.txt““w+“）;
	int ijk;
	struct shuju
	{
		float z[M+1][N];//M=256N=128
	};
	struct shuju shu[2];

	float yz[3]rjieguo[N]z1[N];
	float s[N]; 

	for（i=0;i<=M;i++）//i 256 j 128
	{
		for（j=0;j		{
			shu[0].z[i][j]=1;
			shu[1].z[i][j]=0;
		}
	}
	shu[0].z[200][64]=1;
	for（i=0;i	{
		for（j=0;j
	}

	//////////////////////////////////////////////////////
	//              求取滤波因子            //

	r=lvboyinzi（yz）;
	printf（“%f\n“r）;
	//////////////////////////////////////////////////////

	///褶积运算//////

	printf（“正确\n“）;

	for（i=0;i	{
		for（j=M-1;j>i;j--）
		{
			for（k=0;k			for（k=0;k			zheji（yzjieguoz1）;

			for（k=0;k			for（k=1;k			s[N-1]=0;//给定SN的值，便于求Si
			for（k=N-2;k>=0;k--）s[k]=shu[