% calculateIndicatorCorrelation.m
function [correlationDegree] = calculateIndicatorCorrelation(matterElementMatrix, evaluatedMatterElement)
numIndicators = size(matterElementMatrix, 1);
numLevels = size(matterElementMatrix, 2);
correlationDegree = zeros(1, numLevels);
for i = 1:numLevels
correlationSum = 0;
for j = 1:numIndicators
classicalDomain = matterElementMatrix{j, i}(1, :);
sectionDomain = matterElementMatrix{j, i}(2, :);
minValues = min(evaluatedMatterElement{j}(2, :), sectionDomain, [], 2);
correlationSum = correlationSum + sum(minValues) / sum(classicalDomain - sectionDomain);
end
correlationDegree(i) = correlationSum / numIndicators;
end
end