基于matlab的地震波检测数值模拟仿真,含仿真操作录像

clc; close all; clear all; load('olah_jan_2012.mat'); rata2 = mean(kmp_geomag); pj_kmp_geomag = length(kmp_geomag); n = 0:pj_kmp_geomag-1; n = n/(1440-1); hari = 31; point = 64; h = ones(1,point-1)/(point-1); conv_kmp_geomag = conv(kmp_geomag,h, 'valid'); pj_conv_kmp_geomag = length(conv_kmp_geomag); awl = 1; akh = 10; pj_magnit = length(magnit); plot(n,kmp_geomag); hold on; plot(n((point/2)-1:(pj_conv_kmp_geomag+((point/2)-2))),conv_kmp_geomag,'r'); for mm=awl:akh text(gempa(mm)/(1440-1),kmp_geomag(gempa(mm)),[' ',num2str(mm),... '. M=',num2str(magnit(mm))],'EdgeColor','k'); plot(gempa(mm)/(1440-1),kmp_geomag(gempa(mm)),'ko',... 'MarkerEdgeColor','k',... 'MarkerFaceColor','g',... 'MarkerSize',6); end for nn = 1:hari plot([nn nn], [min(kmp_geomag) max(kmp_geomag)],':k'); hold on; end plot([0 31],[rata2 rata2],':k'); axis('tight'); ylabel('magnetic flux density (nT)'); xlabel('date'); title('magnetic flux density and the earthquake in January 2012');