基于价格变动的市场分割指数分析

引言

市场分割程度是衡量经济一体化的重要指标之一。本文通过构建相邻省市商品价格的相对变动模型，分析不同地区之间的市场分割情况，并计算市场分割指数的变化趋势。

本文来自对论文的复现：赵儒煜, & 孙宁志. (2019). 市场分割结构及其对区域经济的影响. 南昌大学学报(人文社会科学版),50,(6), 33-44.

1. 参数设置与数据初始化

我们设定以下参数：

• T = 10：分析的年份数
• M = 8：商品种类数
• K = 61：相邻省市对数

数据存储方式：

• P_i(K, M, T) 和 P_j(K, M, T) 分别存储两个相邻地区的商品价格。

2. 随机生成价格数据

为了模拟实际情况，我们假设每年的商品价格在 [10, 20] 之间变化，并随机生成数据。

3. 计算相对价格指数差分

相对价格指数采用对数增长率计算。

4. 计算市场分割指数

市场分割指数通过计算去均值后的相对价格指数方差得到。

5. 相邻省市市场分割排名

对每年的市场分割指数进行降序排序，并输出排名。

6. 绘制市场分割指数趋势

我们绘制市场分割指数随时间的变化趋势。

7. 结果分析

• 市场分割指数越大，说明该年份市场的分割程度较高，地区间商品价格联动性较弱。
• 指数下降，表明市场一体化程度提高，不同地区的价格差异在缩小。

8. 结论

本研究通过分析相邻省市的价格变动情况，构建了一种衡量市场分割程度的方法。研究结果可为区域经济政策、市场一体化研究和商品流通管理提供量化支持。

完整代码：

%% 参数设置
T = 10;      % 年份数
M = 8;       % 商品种类数
K = 61;      % 相邻省市对数

%% 数据初始化
% 用三维数组存储各地区价格数据：P_i 和 P_j 的尺寸均为 K x M x T
P_i = zeros(K, M, T);
P_j = zeros(K, M, T);

%% 随机生成各年份数据
% 模拟生成的数据取值范围为 [10,20]
for t = 1:T
    % 随机生成 61 x 16 的数据矩阵，每个元素均匀分布在 [10,20]
    data = 10 + 10 * rand(K, 16);
    % 前8列为地区 i 的 8 类商品价格，后8列为地区 j 的 8 类商品价格
    P_i(:, :, t) = data(:, 1:8);
    P_j(:, :, t) = data(:, 9:16);
end

%% 计算相对价格指数的差分形式
% 根据公式：
%   ΔQ_{ij,t}^m = ln(P_i(t)/P_i(t-1)) - ln(P_j(t)/P_j(t-1))
% 其中 t 从 2 开始比较
DeltaQ = zeros(K, M, T-1);
for t = 2:T
    % 计算各地区各商品的价格增长率（对数差分）
    growth_i = log(P_i(:,:,t) ./ P_i(:,:,t-1));
    growth_j = log(P_j(:,:,t) ./ P_j(:,:,t-1));
    % 差分形式的相对价格指数
    DeltaQ(:,:,t-1) = growth_i - growth_j;
end

%% 去均值处理及计算市场分割指数
% 对每一对 (i,j) 在每个时间段，先计算8个产品的均值，再去均值后计算方差
meanDeltaQ = mean(DeltaQ, 2);  % 沿产品维度取均值，尺寸：K x 1 x (T-1)
q = DeltaQ - meanDeltaQ;       % 去均值后的相对价格差分数据
var_q = var(q, 0, 2);          % 沿产品维度计算方差，结果尺寸：K x 1 x (T-1)
% 整体市场分割指数（所有省市对的平均值）
marketSegIndex = squeeze(mean(var_q, 1));  % 尺寸：(T-1) x 1

%% 输出每个时间段内相邻省市的分割程度排名
% 对于每个时间段（对应年份 t=2,...,T），输出 61 对省市的分割程度排名
for t = 1:T-1
    % 取出当前时间段各省市对的市场分割指数，结果为 K x 1 向量
    segIndices = squeeze(var_q(:,1,t));
    % 对分割指数降序排序（分割程度越大，排名越靠前）
    [sortedSeg, sortOrder] = sort(segIndices, 'descend');
    
    fprintf('时间段 %d（对应年份 %d）的相邻省市分割程度排名：\n', t+1, t+1);
    fprintf('排名\t省市对编号\t分割程度（方差）\n');
    for r = 1:K
        fprintf('%d\t%d\t\t%.4f\n', r, sortOrder(r), sortedSeg(r));
    end
    fprintf('\n');
end

%% 绘制整体市场分割指数随时间变化趋势图
figure;
plot(2:T, marketSegIndex, '-o', 'LineWidth', 2);
xlabel('年份');
ylabel('市场分割指数（平均价格变动方差）');
title('整体市场分割指数随时间的演变趋势');
grid on;