gdttest.m

来自「数字信号处理的Matlab实现的随书光盘」· M 代码 · 共 7 行

function [x,y,z]=gdttest(a,t0,si)
%wavelet transform Mexscro Hat
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;   %一阶导数