2022年9月7日星期三 BPP的球面N个点的matlab仿真

matlab代码

clc;
close all;
clear all;
fclose('all');
format shortG;
 
nums_points = 1000; % 点数量
radius = 10;
loc = zeros(nums_points, 3);
for ii = 1 : nums_points
    phi = acos(-1.0 + (2.0 * ii - 1.0) / nums_points);
    theta = sqrt(nums_points * pi) * phi;
    loc(ii, 1) = radius * cos(theta) * sin(phi);
    loc(ii, 2) = radius * sin(theta) * sin(phi);
    loc(ii, 3) = radius * cos(phi);
end
 
figure
plot3(loc(:, 1), loc(:, 2), loc(:, 3), '.');hold on;% 画结果
[xx, yy, zz]=sphere(50); 
h2=surf(xx * radius, yy * radius, zz * radius); % 画一个球面做参考
set(h2,'edgecolor','none','facecolor','r','facealpha',0.7);
axis equal;
axis([-radius, radius, -radius, radius, -radius, radius]);
hold off;

参考：在球面上生成均匀分布点方法总结_八步赶蝉的博客-CSDN博客_球面均匀分布