%Samp10_5
x=zeros(1,500);   %設置空矢量
%給出檢測數據
for t=1:1:500  
    if (t<200)
     x(t)=50.*exp(t/300).*sin(2*pi*0.01*t);
 elseif (t>=200)&(t<300)
   x(t)=50.*exp(t/300).*sin(2*pi*0.01*200)+30;
elseif ((t>=300)&(t<400))
    x(t)=50.*exp(t/300).*sin(2*pi*0.01*200)+10;
else
    x(t)=50.*exp(t/300).*sin(2*pi*0.01*200)+10+250*sin(2*pi*0.003*t);
end
end
%結束檢測數據
x=x';   %將檢測數據轉置以備后面應用
x=x-mean(x);   %去掉平均值
    a=1;   %給出小波尺度為1
for n=(1:1:500)
    t0=n*1;
    [g,fai,v]=gdttest(a,t0,500);   %調用小波函數的計算
    xx(n)=sum(x.*g');    %一階導數的計算（墨西哥草帽波）
    yyy(n)=sum(x.*(fai)');%高斯函數的計算
    zzz(n)=sum(x.*v');   %二階導數的計算
end
subplot(3,1,1),plot(x);   %繪出原始數據波形
ylabel('檢測數據')
subplot(3,1,2),plot(yyy);   %繪出高斯函數為基函數的小波變換
ylabel('高斯函數結果')
subplot(3,1,3),plot(zzz); %繪出高斯函數的一階導數為基函數的小波變換
ylabel('一階導數結果')
xlabel('時間/s')