From 03de95f41f2f6b733cbf2f7b0671a566c3712da0 Mon Sep 17 00:00:00 2001
From: L_DelOff <51275636+LDelOff@users.noreply.github.com>
Date: Wed, 16 Jun 2021 23:29:35 +0300
Subject: [PATCH] =?UTF-8?q?=D0=9A=D0=BE=D1=80=D1=80=D0=B5=D0=BB=D1=8F?=
 =?UTF-8?q?=D1=86=D0=B8=D0=BE=D0=BD=D0=BD=D1=8B=D0=B9=20=D0=B0=D0=BD=D0=B0?=
 =?UTF-8?q?=D0=BB=D0=B8=D0=B7?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

---
 analysis_corr.m | 52 +++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 52 insertions(+)
 create mode 100644 analysis_corr.m

diff --git a/analysis_corr.m b/analysis_corr.m
new file mode 100644
index 0000000..35009c8
--- /dev/null
+++ b/analysis_corr.m
@@ -0,0 +1,52 @@
+function param=analysis_corr(param)
+%% By L_DelOff
+% считает корреляционные матрицы по всем измерениям по одному катаклизму
+load(param.filename_report,'report');
+R=[1000 10^(0.43*param.M)];
+for i=1:length(R)
+    index=[];
+    for j=1:length(report.data)
+        if report.data(j).R==R(i)
+            index=[index j];
+        end
+    end
+    for j=1:length(index)
+        corr_table(i).RPA_HHeden{j+1,1}=['F-',num2str(report.data(index(j)).N_sat)];
+        corr_table(i).RPA_HHeden{1,j+1}=corr_table(i).RPA_HHeden{j+1,1};
+        corr_table(i).RPA_O2den=corr_table(i).RPA_HHeden;
+        corr_table(i).EP_Eden=corr_table(i).RPA_HHeden;
+        corr_table(i).EP_Etemp=corr_table(i).RPA_HHeden;
+        corr_table(i).R=R(i);
+    end    
+    for j=1:length(index)
+        for k=j:length(index)
+            a=report.data(index(j)).RPA_HHeden;
+            b=report.data(index(k)).RPA_HHeden;
+            Mc=corrcoef(a,b,'Row','pairwise');
+            corr_table(i).RPA_HHeden{k+1,j+1}=num2str(Mc(1,2));
+            corr_table(i).RPA_HHeden{j+1,k+1}=num2str(Mc(2,1));
+            
+            a=report.data(index(j)).RPA_O2den;
+            b=report.data(index(k)).RPA_O2den;
+            Mc=corrcoef(a,b,'Row','pairwise');
+            corr_table(i).RPA_O2den{k+1,j+1}=num2str(Mc(1,2));
+            corr_table(i).RPA_O2den{j+1,k+1}=num2str(Mc(2,1));
+            
+            a=report.data(index(j)).EP_Eden;
+            b=report.data(index(k)).EP_Eden;
+            Mc=corrcoef(a,b,'Row','pairwise');
+            corr_table(i).EP_Eden{k+1,j+1}=num2str(Mc(1,2));
+            corr_table(i).EP_Eden{j+1,k+1}=num2str(Mc(2,1));
+            
+            a=report.data(index(j)).EP_Etemp;
+            b=report.data(index(k)).EP_Etemp;
+            Mc=corrcoef(a,b,'Row','pairwise');
+            corr_table(i).EP_Etemp{k+1,j+1}=num2str(Mc(1,2));
+            corr_table(i).EP_Etemp{j+1,k+1}=num2str(Mc(2,1));
+        end
+    end
+end
+report.corr_table=corr_table;
+param.report=report;
+save(param.filename_report,'report');
+end
\ No newline at end of file