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']);