L_DelOff
/
Earthquake
1
0
Fork 0
Code Issues 1 Pull Requests Packages Projects Releases Wiki Activity

Добавлен анализ величины отклонения параметра

master
L_DelOff 2021-06-16 23:27:33 +03:00
parent 960cc6480f
commit f954b43594
1 changed files with 25 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

25
analysis.m
View File

@ -76,10 +76,35 @@ for i=1:length(param.dates)
report.data(num).EP_Etemp(i)=mymean(Etemp,1e10); report.data(num).EP_Etemp(i)=mymean(Etemp,1e10);
end end
report.data(num).RPA_HHeden_lim=mylim(report.data(num).RPA_HHeden,param);
report.data(num).RPA_O2den_lim=mylim(report.data(num).RPA_O2den,param);
report.data(num).EP_Eden_lim=mylim(report.data(num).EP_Eden,param);
report.data(num).EP_Etemp_lim=mylim(report.data(num).EP_Etemp,param);
param.report=report; param.report=report;
save(fname_rep,'report'); save(fname_rep,'report');
end end
function out=mylim(y,param)
Q1=quantile(y,0.25);
Q2=quantile(y,0.5);
Q3=quantile(y,0.75);
L1=Q2-param.k*(Q3-Q1);
L2=Q2+param.k*(Q3-Q1);
out=[];
for i=1:length(y)
if (y(i)<=L1)||(y(i)>=L2)&&(isnan(y(i)))
out(i)=(y(i)-Q2)*100/(L1-Q2);% отклонение в процентах
%out(i)=1;
else
out(i)=0;
end
end
end
function num=give_number(report,f,R) function num=give_number(report,f,R)
flag=0; % Проверяет, есть ли в отчете данные со спутника f flag=0; % Проверяет, есть ли в отчете данные со спутника f
[~,b]=size(report.data); [~,b]=size(report.data);