【实例简介】
【实例截图】
【核心代码】
samplingrate = 0.00005;
[s,t,f]=st(amplitude,0,0.5/samplingrate,samplingrate,1);
T=max(t);
subplot(3,1,1),plot(t,amplitude,'color','k','linewidth',1)
xlabel('时间/s');
ylabel('幅值/V');
xlim([0 T]);
% title('微震波形','FontName','Times New Roman','FontSize',14);
set(gca,'FontName','Times New Roman','FontSize',14);
set(get(gca,'XLabel'),'FontName','Times New Roman','FontSize',14);
set(get(gca,'YLabel'),'FontName','Times New Roman','FontSize',14);
subplot(3,1,2:3),contour(t,f,abs(s),'fill','on');
xlabel('时间/s');
ylabel('频率/Hz');
xlim([0 T]);
ylim([0 3000]);
% title('时频谱','FontName','Times New Roman','FontSize',14);
set(gca,'FontName','Times New Roman','FontSize',14);
set(get(gca,'XLabel'),'FontName','Times New Roman','FontSize',14);
set(get(gca,'YLabel'),'FontName','Times New Roman','FontSize',14);