如何在Matlab/Excel/Python中绘制含4变量的3D柱状图？

核心问题分析 #

你的代码报错大概率是因为未正确处理非网格状的(x,y)样本数据，或是颜色映射逻辑错误。下面给出完整修正代码，实现按第四变量着色的3D柱状图，同时补充Python实现方案（Excel方案局限性较大，仅做简要说明）。

MATLAB 实现代码 #

假设你的CSV数据列名为：Temperature, Concentration, Velocity, FourthVariable（需替换为实际第四变量名）：

% 1. 读取CSV数据（优先用readtable保留列名）
data = readtable('sample_data.csv');

% 提取变量
x = data.Temperature;
y = data.Concentration;
z = data.Velocity;
c = data.FourthVariable; % 替换为你的第四变量列名

% 2. 处理数据：将离散(x,y)转换为网格结构
unique_x = sort(unique(x));
unique_y = sort(unique(y));
[X_grid, Y_grid] = meshgrid(unique_x, unique_y);

% 初始化高度和颜色网格（若有重复(x,y)样本，这里取均值，可按需修改聚合方式）
Z_grid = zeros(size(X_grid));
C_grid = zeros(size(X_grid));
for i = 1:length(unique_x)
    for j = 1:length(unique_y)
        match_idx = (x == unique_x(i)) & (y == unique_y(j));
        if any(match_idx)
            Z_grid(j,i) = mean(z(match_idx));
            C_grid(j,i) = mean(c(match_idx));
        end
    end
end

% 3. 绘制3D柱状图并设置颜色
figure('Position', [100 100 800 600])
bar_obj = bar3(Z_grid);

% 为每个柱子分配对应颜色
colormap jet; % 可替换为parula、viridis等色图
for bar_group = 1:length(bar_obj)
    for face_idx = 1:length(bar_obj(bar_group).Faces)
        % 映射到网格位置获取颜色值
        [grid_y, grid_x] = ind2sub(size(Z_grid), bar_obj(bar_group).Faces(face_idx, 1));
        color_val = interp1(min(C_grid(:)), max(C_grid(:)), C_grid(grid_y, grid_x), colormap);
        bar_obj(bar_group).FaceVertexCData(face_idx, :) = color_val;
    end
    bar_obj(bar_group).EdgeColor = 'none'; % 隐藏柱子边框，提升视觉效果
end

% 4. 添加标签与颜色条
xlabel('Temperature', 'FontSize', 12);
ylabel('Concentration', 'FontSize', 12);
zlabel('Velocity', 'FontSize', 12);
color_bar = colorbar;
color_bar.Label.String = 'Your Fourth Variable Name'; % 替换为第四变量名称
title('3D Bar Chart with Color-Coded Fourth Variable', 'FontSize', 14);

常见报错修正点 #

• 若报错变量未定义：检查CSV列名是否与代码中一致，或改用readmatrix读取后按索引提取列（如x = data(:,1)）
• 若颜色不匹配：确保C_grid的维度与Z_grid完全一致，避免索引越界

Python 实现方案（更灵活易维护） #

用matplotlib的bar3d直接绘制每个离散样本的柱子，颜色映射更直观：

import pandas as pd
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D

# 读取CSV数据
data = pd.read_csv('sample_data.csv')

# 提取变量
x = data['Temperature'].values
y = data['Concentration'].values
z = data['Velocity'].values
c = data['FourthVariable'].values

# 创建3D画布
fig = plt.figure(figsize=(10, 7))
ax = fig.add_subplot(111, projection='3d')

# 设置柱子的宽度/深度（根据数据范围调整）
dx = 0.4
dy = 0.4

# 颜色映射配置
norm = plt.Normalize(c.min(), c.max())
cmap = plt.cm.jet

# 逐个绘制柱子并设置颜色
for xi, yi, zi, ci in zip(x, y, z, c):
    bar_color = cmap(norm(ci))
    ax.bar3d(xi - dx/2, yi - dy/2, 0, dx, dy, zi, color=bar_color, alpha=0.8)

# 配置标签与颜色条
ax.set_xlabel('Temperature', fontsize=12)
ax.set_ylabel('Concentration', fontsize=12)
ax.set_zlabel('Velocity', fontsize=12)
mappable = plt.cm.ScalarMappable(norm=norm, cmap=cmap)
mappable.set_array(c)
plt.colorbar(mappable, label='Your Fourth Variable Name', pad=0.1)
plt.title('3D Bar Chart with Color-Coded Fourth Variable', fontsize=14)

plt.show()

Excel 实现说明 #

Excel原生3D柱状图不支持直接按第四变量着色，需通过以下步骤实现：

1. 用数据透视表将Temperature/Concentration作为行/列，Velocity作为值
2. 插入3D柱状图后，手动逐个设置柱子颜色（效率极低）
3. 或通过VBA批量匹配第四变量值与颜色，但灵活性远不如MATLAB/Python

内容的提问来源于stack exchange，提问作者Kilara