matlab画雷克子波

clear; clc

fm = 400000; %主频
dt = 1e-4; %时间域采样间隔
number = 2000; %采样点数
t = linspace(-0.05,0.05,number)';

R = (1-2*(pi*fm*(t-0.02)*1e-3).^2) .* exp(-(pi*fm*(t-0.02)*1e-3).^2);
plot(t,R,'linewidth', 1.5);

ylim([-1,2])
title('Ricker-wavelet');
xlabel('t(ms)');
ylabel('A');
set(gca, 'LineWidth',1.5,'Fontname', 'times new roman','fontsize',22)