ssa

1
%% Compare reconstruction and original time series
2
clear; close all; clc;
3

4
    %% Covariance matrix - Picture 1
5
    % Ковариационная матрица
6
    C=loadData(1,"C");
7
    figure();
8
    set(gcf,'name','Covariance matrix');
9
    clf;
10
    imagesc(C);
11
    axis square;
12
    set(gca,'clim',[-1 1]);
13
    colorbar;
14

15
    %% Eigenvalues - Picture 2
16
    % Собственный вектор
17
    LBD=loadData(2,"LBD");
18
    figure();
19
    set(gcf,'name','Eigenvalues')
20
    clf;
21
    plot(LBD,'o-');
22

23
    %% Original time series and reconstruction - Picture 3
24
    % Первоначальный временной ряд и его реконструкция для sET12
25
    seg = 100; % номер сегмента pw для визуализации
26
    win = 1024;
27
    spw=loadData(3,"spw");
28
    tim=loadData(3,"tim")';
29
    sET12=loadData(3,"sET12");
30
    figure();
31
    set(gcf,'name','Original time series and reconstruction'); clf;
32
    plot(tim(1:win),spw(:,seg),'b-',tim(1:win),sET12(:,seg),'r-');
33
    legend('Original','sET12'); xlabel("t,s",'interp','none'); ylabel("sET",'interp','none');
34
    
35
    %% Original time series and reconstruction - Picture 4
36
    % Первоначальный временной ряд и его реконструкция для sET34
37
    seg = 100; % номер сегмента pw для визуализации
38
    win = 1024;
39
    spw=loadData(4,"spw");
40
    tim=loadData(4,"tim")';
41
    sET34=loadData(4,"sET34");
42
    figure();
43
    set(gcf,'name','Original time series and reconstruction'); clf;
44
    plot(tim(1:win),spw(:,seg),'b-',tim(1:win),sET34(:,seg),'r-');
45
    legend('Original','sET34'); xlabel("t,s",'interp','none'); ylabel("sET",'interp','none');
46
    
47
    %% Визуализация АКФ сингулярных троек для сегментов pw - Picture 5
48
    lag  = floor(win/10); % наибольший лаг АКФ <= win/10
49
    lagS = 2*lag;
50
    ns=loadData(5,"ns");
51
    Time=loadData(5,"Time");
52
    Acf_sET12=loadData(5,"Acf_sET12");
53
    figure();
54
    set(gcf,'name','АКФ сингулярных троек sET12 сегментов pw'); clf;
55
    % mesh(ns,lgl,Acf_sET12(1:lag,:),'FaceAlpha',0.5,'FaceColor','flat'); colorbar;
56
    % xlabel("ns",'interp','none'); ylabel("lag",'interp','none');
57
    mesh(ns,Time,Acf_sET12(1:lag,:),'FaceAlpha',0.5,'FaceColor','flat'); colorbar;
58
    xlabel("ns",'interp','none'); ylabel("lag,s",'interp','none');
59
    zlabel("Acf",'interp','none'); grid on;
60
    
61
    %% Огибающие АКФ сингулярных троек sET12 сегментов pw - Picture 6
62
    ns=loadData(6,"ns");
63
    Time=loadData(6,"Time")';
64
    EnvAcf_sET12=loadData(6,"EnvAcf_sET12");
65
    figure();
66
    set(gcf,'name','Огибающие АКФ сингулярных троек sET12 сегментов pw');
67
    clf;
68
    % mesh(ns,lgl,EnvAcf_sET12(1:lag,:),'FaceAlpha',0.5,'FaceColor','flat'); colorbar;
69
    % xlabel("ns",'interp','none'); ylabel("lag",'interp','none');
70
    mesh(ns,Time,EnvAcf_sET12(1:lag,:),'FaceAlpha',0.5,'FaceColor','flat'); colorbar;
71
    xlabel("ns",'interp','none'); ylabel("lag,s",'interp','none');
72
    zlabel("Env_Acf",'interp','none'); grid on;
73

74
    %% Нормированные АКФ сингулярных троек sET12 сегментов pw - Picture 7
75
    ns=loadData(7,"ns");
76
    Time=loadData(7,"Time")';
77
    AcfNrm_sET12=loadData(7,"AcfNrm_sET12");
78
    figure();
79
    set(gcf,'name','Нормированные АКФ сингулярных троек sET12 сегментов pw');
80
    clf;
81
    % mesh(ns,lgl,AcfNrm_sET12(1:lag,:),'FaceAlpha',0.5,'FaceColor','flat'); colorbar;
82
    % xlabel("ns",'interp','none'); ylabel("lag",'interp','none');
83
    mesh(ns,Time,AcfNrm_sET12(1:lag,:),'FaceAlpha',0.5,'FaceColor','flat'); colorbar;
84
    xlabel("ns",'interp','none'); ylabel("lag,s",'interp','none');
85
    zlabel("Acf_Nrm",'interp','none'); grid on;
86
    
87
    %% Частоты нормир-ой АКФ сингуляр-х троек сегментов pw - Picture 8
88
    ns=loadData(8,"ns");
89
    insFrc_AcfNrm=loadData(8,"insFrc_AcfNrm");
90
    smo_insFrc_AcfNrm=loadData(8,"smo_insFrc_AcfNrm")';
91
    figure();
92
    set(gcf,'name','Частоты нормир-ой АКФ сингуляр-х троек сегментов pw');
93
    clf;
94
    p1 = plot(ns,insFrc_AcfNrm,'b','LineWidth',0.8); hold on;
95
    plot(ns,smo_insFrc_AcfNrm,'r','LineWidth',0.4)
96
    % plot(ns,smo_insFrc_AcfNrm,'r','LineWidth',0.8); grid on; % smo_insFrc_AcfNrm
97
    xlabel("ns",'interp','none'); ylabel("insFrc_AcfNrm,Hz",'interp','none');
98
    legend(p1,'sET12'); 
99
    
100
    %% Периодограмма Томсона sET12 сегментов pw - Picture 9
101
    figure();
102
    iGmin=loadVar(9,"iGmin");
103
    iGmax=loadVar(9,"iGmax");
104
    ns=loadData(9,"ns");
105
    pto_sET12=loadData(9,"pto_sET12");
106
    fG=loadData(9,"fG");
107
    set(gcf,'name','Периодограмма Томсона sET12 сегментов pw');
108
    clf;
109
    mesh(ns,fG(iGmin:iGmax),pto_sET12(iGmin:iGmax,:),'FaceAlpha',0.5,'FaceColor','flat');
110
    colorbar; grid on;
111
    xlabel("ns",'interp','none'); ylabel("f,Hz",'interp','none');
112
    zlabel("P(f)",'interp','none');
113
    
114
    %% Агрегирование сегментов очищенной пульсовой волны cpw - Picture 10
115
    figure();
116
    ns=loadData(10,"ns");
117
    pto_fMAX12=loadData(10,"pto_fMAX12");
118
    set(gcf,'name','Частоты основного тона sET сегментов pw');
119
    clf;
120
    plot(ns,pto_fMAX12,'b'); hold on; grid on;
121
    % plot(ns,smo_pto_fMAX12,'r','LineWidth',0.8); grid on;
122
    xlabel("ns",'interp','none'); ylabel("fMAX,Hz",'interp','none');
123

124