可以给我一个python或者c++或者matlab的程序来实现吗

以下是一个使用 Python 实现骨架提取的简单示例程序，其中使用了 Floyd 算法计算最短路径，使用了 skimage 库中的 skeletonize 函数进行骨架提取。

import numpy as np
from skimage.morphology import skeletonize

# 计算可达矩阵
def get_reachable_matrix(adj_matrix):
    n = adj_matrix.shape[0]
    reachable_matrix = np.zeros((n, n))
    for k in range(n):
        for i in range(n):
            for j in range(n):
                if adj_matrix[i][k] and adj_matrix[k][j]:
                    reachable_matrix[i][j] = 1
    return reachable_matrix

# 骨架提取
def skeleton_extraction(reachable_matrix):
    skeleton_matrix = skeletonize(reachable_matrix)
    return skeleton_matrix.astype(np.uint8)

# 示例
adj_matrix = np.array([[0, 1, 1, 0],
                       [1, 0, 1, 0],
                       [1, 1, 0, 1],
                       [0, 0, 1, 0]])
reachable_matrix = get_reachable_matrix(adj_matrix)
skeleton_matrix = skeleton_extraction(reachable_matrix)
print(skeleton_matrix)

这里的邻接矩阵 adj_matrix 表示一个简单无向图，输出的结果即为其骨架矩阵。除此之外，你还可以根据自己的需求进行修改，比如使用 Dijkstra 算法计算最短路径，或者使用其他的骨架提取算法等。