/*
?????A?Program?of?synthetic?seismogram?generation.
????????????????????Author:?Zhang?Junhua?????Date:?Sep.2006?????*/

#include?“stdio.h“??//?standard?input/output?head?file.
#include?“math.h“??//?math.?head?file.
#include?“malloc.h“??//?memory?allocation?head?file.

#define?fm??35???//?The?main?frequency?of?Ricker?wavelet?you?can?change?it.
#define?dt??0.002??//?The?sample?interval?you?can?also?set?it?to?0.001?or?0.004.
#define?XL??0.060??//?The?length?of?wavelet（s）?you?can?change?it.
#define?HL??0.300??//?The?length?of?h（t）?namely?reflection?coefficient.
????????????????????????//?Here?the?length?is?0.3?second?you?can?change?it.

#define?PI?3.1415926

main（）
{
??FILE?*fp1*fp2*fp3;

??float?*x*h*y;?//?define?the?dynamic?array.
??float?t;
??int?iMN;
??void?conv（）;

??M=XL/dt+1;??//?The?number?of?x（t）.?For?this?example?M=31.
??N=HL/dt+1;??//?The?number?of?h（t）.?For?this?example?N=151.

//?allocate?dynamic?array.
??x=（float?*）calloc（Msizeof（float））;
??h=（float?*）calloc（Nsizeof（float））;
??y=（float?*）calloc（M+N-1sizeof（float））;

//?Generate?Ricker?wavelet.

??fp1=fopen（“Ricker.txt“