#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[