Корреляционный анализ

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