function [x,y,z]=gdttest(a,t0,si)
%wavelet transform Mexscro Hat
t=(1:1:si);    %給出時(shí)間序列
x=(1-((t-t0)/a).^2).*exp(-((t-t0)/a).^2/2)/a;   %二階導(dǎo)數(shù)
y=exp(-((t-t0)/a).^2/2)/a;     %高斯函數(shù)
z=(t-t0).*exp(-((t-t0)/a).^2/2)/a;   %一階導(dǎo)數(shù)