% 定义线段首尾两个端点的坐标 start_point = [0, 0, 0]; end_point = [0, 0, 1000]; % 定义圆柱半径和每段长度 radius = 10; segment_length = 10; % 计算线段向量和长度 line_vector = end_point - start_point; line_length = norm(line_vector); % 计算每段的方向向量 direction_vector = line_vector / line_length * segment_length; % 初始化质心坐标和质量 centroid = [0, 0, 0]; total_mass = 0; % 循环处理每段 for i = 1:floor(line_length / segment_length) % 计算当前段的起点和终点坐标 segment_start = start_point + (i-1) * direction_vector; segment_end = segment_start + direction_vector; % 计算圆柱体积

cylinder_volume = pi * radius^2 * segment_length;

% 计算质量
segment_mass = cylinder_volume * density; % density是圆柱体密度
total_mass = total_mass + segment_mass;

% 计算当前段的质心坐标
segment_centroid = (segment_start + segment_end) / 2;
centroid = centroid + segment_centroid * segment_mass;

end

% 计算整个物体的质心坐标
centroid = centroid / total_mass;