function y = CLR(x)
% 中心对数比变化
b = geomean(x); % 分母为几何平均数
y = log(x/b);
end
function x_original = inverseCLR(clr_data)
% inverseCLR - 逆中心对数比变换
%
% 输入:
% clr_data - 已经进行CLR变换的数据向量
%
% 输出:
% x_original - 逆变换后的原始数据(未标准化的比例数据)
% 计算CLR数据的反对数(指数)
exp_clr_data = exp(clr_data);
% 计算反对数后数据的归一化比例
x_original = exp_clr_data / sum(exp_clr_data);
end