line 22, in <module> ax.bar(x + i*width, data[i], width, label=str(2022-i))

这个错误提示的意思是在程序运行时出现了错误，具体是在第22行，即绘制柱状图时出现了错误。

经过查看代码和数据，我发现错误原因可能是数据的维度不匹配。由于你所提供的数据是10年的数据，而在第22行代码中，绘图时使用了一个循环，循环次数为数据的行数，也就是10次。因此，我猜测你的数据应该是一个10行1列的数据，而不是一个1行10列的数据。

为了解决这个问题，你可以将数据进行转置，使得数据变成10行1列。具体修改方法如下：

1. 将数据转置，即用 `data.T` 替换原来的 `data`。
2. 在循环中，将 `data[i]` 替换为 `data.T[i]`。

修改后的代码如下所示：

import matplotlib.pyplot as plt
import numpy as np

# 设置中文字体
plt.rcParams['font.sans-serif']=['SimHei']
plt.rcParams['axes.unicode_minus'] = False

# 数据准备
indicators = ['年末总人口(万人)', '男性人口(万人)', '女性人口(万人)', '城镇人口(万人)', '乡村人口(万人)']
data = np.array([[141175, 141260, 141212, 141008, 140541, 140011, 139232, 138326, 137646, 136726],
                 [72206, 72311, 72357, 72039, 71864, 71650, 71307, 70857, 70522, 70063],
                 [68969, 68949, 68855, 68969, 68677, 68361, 67925, 67469, 67124, 66663],
                 [92071, 91425, 90220, 88426, 86433, 84343, 81924, 79302, 76738, 74502],
                 [49104, 49835, 50992, 52582, 54108, 55668, 57308, 59024, 60908, 62224]])

# 绘图
fig, ax = plt.subplots(figsize=(10,6))
x = np.arange(len(indicators))
width = 0.15

for i in range(data.T.shape[0]):
    ax.bar(x + i*width, data.T[i], width, label=str(2022-i))

ax.set_title('近10年中国人口统计图-直方图')
ax.set_xticks(x + 2*width)
ax.set_xticklabels(indicators, rotation=30)
ax.legend()

plt.show()

运行修改后的代码，你应该可以得到一个包含近十年中国人口统计图-直方图的图表。如果还有问题，欢迎继续向我提问。