写出matlab左上半圆生成随机点的代码

以下是Matlab代码，用于在左上半圆内生成指定数量的随机点：

% 设置圆的半径和圆心坐标
radius = 1;
center_x = 0;
center_y = 1;

% 设置生成点的数量
num_points = 1000;

% 生成随机点
points = zeros(num_points, 2);
count = 0;
while count < num_points
    % 生成随机点
    x = -1 + 2*rand();
    y = 0 + 1*rand();
    % 判断是否在左上半圆内
    if x^2 + (y-center_y)^2 < radius^2
        count = count + 1;
        points(count, :) = [x, y];
    end
end

% 绘制散点图
scatter(points(:,1), points(:,2), 'filled');
axis([-1 1 0 1]);

这段代码利用了圆的基本方程 x^2 + y^2 = r^2 的性质，通过生成坐标在 [-1,1] 范围内的随机点，然后筛选出在左上半圆内的点，最终用散点图进行可视化。