写一个完全广义基尼系数的MATLAB程序

以下是一个完全广义基尼系数的MATLAB程序：

function [gini] = full_generalized_gini(x, w, p)
% x: 向量，表示不同个体的收入或财富
% w: 向量，表示不同个体的权重
% p: 参数，控制不平等厌恶程度，若p=1，则为标准基尼系数
% gini: 完全广义基尼系数

n = numel(x); % 总个体数
wx = w.*x; % 考虑权重的收入或财富
wx_sum = sum(wx); % 加权总收入或财富
w_sum = sum(w); % 总权重

% 计算基尼系数的分子
gini_num = 0;
for i = 1:n
    for j = 1:n
        gini_num = gini_num + w(i)*w(j)*(abs(wx(i)-wx(j)))^p;
    end
end

% 计算完全广义基尼系数
gini = gini_num/(2*w_sum*wx_sum^p) - (n-1)/(2*n);
end

该程序的输入为一个向量 x，表示不同个体的收入或财富；一个向量 w，表示不同个体的权重；一个参数 p，控制不平等厌恶程度，若 p=1，则为标准基尼系数。输出为完全广义基尼系数。