www.gusucode.com > 股票波动分析源码程序 > 股票波动分析源码程序/GeomagEarthQuakebyBUdis/EartQGeomag.m
clc; close all; clear all; nm_sta = 'kak'; thn = '2012'; bln = '01'; jml_hri = [31 29 31 30 31 30 31 31 30 31 30 31]; max_hri = jml_hri(str2num(bln)); nmfl_sv = [nm_sta thn '_' bln '.mat']; jdl = [nm_sta ' ' thn '-' bln]; load (nmfl_sv); n = 1:length(x); hri = n/(24*60); h = sqrt(x.^2 + y.^2); point = 91; b = ones(1,point)/point; a = 1; min_x = min(x); max_x = max(x); min_y = min(y); max_y = max(y); min_z = min(z); max_z = max(z); min_h = min(h); max_h = max(h); nmfl_xls = ['Isi_gmp_' bln '_' thn '_M6toM10.xls']; [num, txt, raw] = xlsread(nmfl_xls); nmfl_xls_5 = ['Isi_gmp_' bln '_' thn '_M5toM6.xls']; [num_5, txt_5, raw_5] = xlsread(nmfl_xls_5); N_num_tbl = length(num(:,1)); dt_num_tbl = (num(:,14)-1440); N_num_tbl_5 = length(num_5(:,1)); dt_num_tbl_5 = (num_5(:,14)-1440); subplot(211); plot(hri,h); axis tight; grid on; hold on; for p = 1:max_hri plot([p p], [min_h max_h],'--r'); hold on; end for pp = 3:N_num_tbl plot(dt_num_tbl(pp)/1440,h(dt_num_tbl(pp)),'o',... 'MarkerEdgeColor','k',... 'MarkerFaceColor','r',... 'MarkerSize',8); hold on; text(dt_num_tbl(pp)/1440,h(dt_num_tbl(pp)),[' ',num2str(num(pp,1))],... 'EdgeColor','k'); end for pp = 3:N_num_tbl_5 plot(dt_num_tbl_5(pp)/1440,h(dt_num_tbl_5(pp)),'o',... 'MarkerEdgeColor','k',... 'MarkerFaceColor','g',... 'MarkerSize',6); hold on; end title(['Komponen H ( ' jdl ' )']); xlabel(['Date']); ylabel(['magnitude']); subplot(212); plot(hri,z); axis tight; grid on; hold on; for p = 1:max_hri plot([p p], [min_z max_z],'--r'); hold on; end for pp = 3:N_num_tbl plot(dt_num_tbl(pp)/1440,z(dt_num_tbl(pp)),'o',... 'MarkerEdgeColor','k',... 'MarkerFaceColor','r',... 'MarkerSize',8); hold on; text(dt_num_tbl(pp)/1440,z(dt_num_tbl(pp)),[' ',num2str(num(pp,1))],... 'EdgeColor','k'); end for pp = 3:N_num_tbl_5 plot(dt_num_tbl_5(pp)/1440,z(dt_num_tbl_5(pp)),'o',... 'MarkerEdgeColor','k',... 'MarkerFaceColor','g',... 'MarkerSize',6); hold on; end title(['Komponen Z ( ' jdl ' )']); xlabel(['Date']); ylabel(['magnitude']);