帮我写个无人机蜂群干扰，不同分布方式产生不同干扰效果的极坐标方向图的matlab代码

以下是一个简单的无人机蜂群干扰的极坐标方向图的 Matlab 代码，其中不同分布方式通过修改变量 `mode` 来实现：

clc;
clear;

% 随机生成100个无人机坐标
x = rand(1, 100) * 100;
y = rand(1, 100) * 100;

% 极坐标转换
theta = atan2(y, x);
r = sqrt(x.^2 + y.^2);

% 干扰模式（1-随机，2-集中，3-散开）
mode = 1;

% 根据不同的干扰模式生成不同的随机扰动
if mode == 1
    % 产生随机扰动
    theta = theta + randn(1, 100) * 0.2;
    r = r + randn(1, 100) * 5;
elseif mode == 2
    % 将所有无人机集中到坐标原点
    theta = theta - pi/2;
    r = r / 5;
elseif mode == 3
    % 将所有无人机散开
    theta = theta + rand(1, 100) * pi;
    r = r * 2;
end

% 处理 theta 的范围
theta(theta < 0) = theta(theta < 0) + 2*pi;
theta(theta > 2*pi) = theta(theta > 2*pi) - 2*pi;

% 统计每个方向上的无人机数量
N = 360;
count = zeros(1, N);
for i = 1:N
    idx = find(theta >= (i-1)/N*2*pi & theta < i/N*2*pi);
    count(i) = length(idx);
end

% 绘制极坐标方向图
theta = linspace(0, 2*pi, N+1);
theta(end) = [];
polarplot(theta, [count, count(1)], '-r');
title(sprintf('Interference Pattern (Mode %d)', mode));

这里使用了随机干扰、集中干扰和散开干扰三种模式。你可以根据需要修改干扰模式和相应的参数，例如扰动的强度、集中的程度等，以产生不同的干扰效果。