samp10_4.m

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

%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')