樹形図で重心法を用いるやり方

樹形図を作成する際に重心法を用いたいのですがどうすればいいのでしょうか。

以下に示すのはChatGPTで作成したサンプルのコードです。

% サンプルデータの生成
X = rand(100, 2); % 100個の2次元データ

% クラスター数
k = 3;

% k-means クラスタリングの実行
[idx, centroids] = kmeans(X, k);

% クラスタリング結果から距離行列を計算
D = pdist(centroids);

% 距離行列を階層的クラスタリングに変換
Z = linkage(D, 'average'); % ここでは平均法を使用

% 樹形図のプロット
dendrogram(Z);
title('重心法によるクラスタリングの樹形図');