ssa
122 строки · 5.5 Кб
1%% Compare reconstruction and original time series
2clear; close all; clc;
3
4%% Covariance matrix - Picture 1
5% Ковариационная матрица
6C=loadData(1,"C");
7figure();
8set(gcf,'name','Covariance matrix');
9clf;
10imagesc(C);
11axis square;
12set(gca,'clim',[-1 1]);
13colorbar;
14
15%% Eigenvalues - Picture 2
16% Собственный вектор
17LBD=loadData(2,"LBD");
18figure();
19set(gcf,'name','Eigenvalues')
20clf;
21plot(LBD,'o-');
22
23%% Original time series and reconstruction - Picture 3
24% Первоначальный временной ряд и его реконструкция для sET12
25seg = 100; % номер сегмента pw для визуализации
26win = 1024;
27spw=loadData(3,"spw");
28tim=loadData(3,"tim")';
29sET12=loadData(3,"sET12");
30figure();
31set(gcf,'name','Original time series and reconstruction'); clf;
32plot(tim(1:win),spw(:,seg),'b-',tim(1:win),sET12(:,seg),'r-');
33legend('Original','sET12'); xlabel("t,s",'interp','none'); ylabel("sET",'interp','none');
34
35%% Original time series and reconstruction - Picture 4
36% Первоначальный временной ряд и его реконструкция для sET34
37seg = 100; % номер сегмента pw для визуализации
38win = 1024;
39spw=loadData(4,"spw");
40tim=loadData(4,"tim")';
41sET34=loadData(4,"sET34");
42figure();
43set(gcf,'name','Original time series and reconstruction'); clf;
44plot(tim(1:win),spw(:,seg),'b-',tim(1:win),sET34(:,seg),'r-');
45legend('Original','sET34'); xlabel("t,s",'interp','none'); ylabel("sET",'interp','none');
46
47%% Визуализация АКФ сингулярных троек для сегментов pw - Picture 5
48lag = floor(win/10); % наибольший лаг АКФ <= win/10
49lagS = 2*lag;
50ns=loadData(5,"ns");
51Time=loadData(5,"Time");
52Acf_sET12=loadData(5,"Acf_sET12");
53figure();
54set(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');
57mesh(ns,Time,Acf_sET12(1:lag,:),'FaceAlpha',0.5,'FaceColor','flat'); colorbar;
58xlabel("ns",'interp','none'); ylabel("lag,s",'interp','none');
59zlabel("Acf",'interp','none'); grid on;
60
61%% Огибающие АКФ сингулярных троек sET12 сегментов pw - Picture 6
62ns=loadData(6,"ns");
63Time=loadData(6,"Time")';
64EnvAcf_sET12=loadData(6,"EnvAcf_sET12");
65figure();
66set(gcf,'name','Огибающие АКФ сингулярных троек sET12 сегментов pw');
67clf;
68% mesh(ns,lgl,EnvAcf_sET12(1:lag,:),'FaceAlpha',0.5,'FaceColor','flat'); colorbar;
69% xlabel("ns",'interp','none'); ylabel("lag",'interp','none');
70mesh(ns,Time,EnvAcf_sET12(1:lag,:),'FaceAlpha',0.5,'FaceColor','flat'); colorbar;
71xlabel("ns",'interp','none'); ylabel("lag,s",'interp','none');
72zlabel("Env_Acf",'interp','none'); grid on;
73
74%% Нормированные АКФ сингулярных троек sET12 сегментов pw - Picture 7
75ns=loadData(7,"ns");
76Time=loadData(7,"Time")';
77AcfNrm_sET12=loadData(7,"AcfNrm_sET12");
78figure();
79set(gcf,'name','Нормированные АКФ сингулярных троек sET12 сегментов pw');
80clf;
81% mesh(ns,lgl,AcfNrm_sET12(1:lag,:),'FaceAlpha',0.5,'FaceColor','flat'); colorbar;
82% xlabel("ns",'interp','none'); ylabel("lag",'interp','none');
83mesh(ns,Time,AcfNrm_sET12(1:lag,:),'FaceAlpha',0.5,'FaceColor','flat'); colorbar;
84xlabel("ns",'interp','none'); ylabel("lag,s",'interp','none');
85zlabel("Acf_Nrm",'interp','none'); grid on;
86
87%% Частоты нормир-ой АКФ сингуляр-х троек сегментов pw - Picture 8
88ns=loadData(8,"ns");
89insFrc_AcfNrm=loadData(8,"insFrc_AcfNrm");
90smo_insFrc_AcfNrm=loadData(8,"smo_insFrc_AcfNrm")';
91figure();
92set(gcf,'name','Частоты нормир-ой АКФ сингуляр-х троек сегментов pw');
93clf;
94p1 = plot(ns,insFrc_AcfNrm,'b','LineWidth',0.8); hold on;
95plot(ns,smo_insFrc_AcfNrm,'r','LineWidth',0.4)
96% plot(ns,smo_insFrc_AcfNrm,'r','LineWidth',0.8); grid on; % smo_insFrc_AcfNrm
97xlabel("ns",'interp','none'); ylabel("insFrc_AcfNrm,Hz",'interp','none');
98legend(p1,'sET12');
99
100%% Периодограмма Томсона sET12 сегментов pw - Picture 9
101figure();
102iGmin=loadVar(9,"iGmin");
103iGmax=loadVar(9,"iGmax");
104ns=loadData(9,"ns");
105pto_sET12=loadData(9,"pto_sET12");
106fG=loadData(9,"fG");
107set(gcf,'name','Периодограмма Томсона sET12 сегментов pw');
108clf;
109mesh(ns,fG(iGmin:iGmax),pto_sET12(iGmin:iGmax,:),'FaceAlpha',0.5,'FaceColor','flat');
110colorbar; grid on;
111xlabel("ns",'interp','none'); ylabel("f,Hz",'interp','none');
112zlabel("P(f)",'interp','none');
113
114%% Агрегирование сегментов очищенной пульсовой волны cpw - Picture 10
115figure();
116ns=loadData(10,"ns");
117pto_fMAX12=loadData(10,"pto_fMAX12");
118set(gcf,'name','Частоты основного тона sET сегментов pw');
119clf;
120plot(ns,pto_fMAX12,'b'); hold on; grid on;
121% plot(ns,smo_pto_fMAX12,'r','LineWidth',0.8); grid on;
122xlabel("ns",'interp','none'); ylabel("fMAX,Hz",'interp','none');
123
124