%Samp10_4
si=500;t0=250;a=10;
t=(1:1:si);    %給出時間序列
x=(1-((t-t0)/a).^2).*exp(-((t-t0)/a).^2/2)/a;   %二階導數
y=exp(-((t-t0)/a).^2/2)/a;     %高斯函數
z=(t-t0).*exp(-((t-t0)/a).^2/2)/a;   %一階導數
subplot(3,1,1),plot(t,x)
title('高斯小波基')
subplot(3,1,2),plot(t,y)
title('小波基的一階導數')
subplot(3,1,3),plot(t,z)
title('小波基的二階導數')
xlabel('時間/s')