用matlab求解模型 $$\begin{aligned} \min Z &= \sum\limits_{i=1}^n (10h_i+10)\\ \text{s.t. } &(x_i-x_j)^2 + (y_i-y_j)^2 \geq (2r+d)^2,\\ &x_i \pm f_i \leq \frac{L}{2}, y_i\pm f_i \leq \frac{L}{2}, i=1,2,\cdots,n\\ &h_1=h_2=\cdots=h_n,\\ &n\in \mathbb{Z}^+ \end{aligned}$$ 其中 $\mathbb{Z}^+$ 表示正整数集合。 。

这是一个带有非线性约束条件的整数规划问题，可以使用MATLAB的混合整数线性规划（MILP）求解器来解决。下面是MATLAB代码：

% 数据
n = 10; % 点的数量
r = 1; % 半径
d = 1; % 最小距离
L = 10; % 边长
f = 1; % 偏移量

% 构造距离矩阵
dist = zeros(n);
for i = 1:n
    for j = 1:n
        dist(i,j) = (i-j)^2 + (i-j)^2;
    end
end

% 定义 MILP 模型
model.sense = 1; % 最小化
model.obj = [repmat(10,n,1); 0]; % 目标函数
model.vtype = [repmat('C',n,1); 'I']; % 变量类型
model.lb = [-inf(n,1); 0]; % 变量下界
model.ub = [inf(n,1); inf]; % 变量上界

% 距离约束
for i = 1:n
    for j = 1:n
        if i ~= j
            model.A = sparse([model.A; zeros(1,n) 1]);
            model.rhs = [model.rhs; (2*r+d)^2];
            model.sense = [model.sense; '>'];
        end
    end
end

% 边界约束
for i = 1:n
    model.A = sparse([model.A; eye(2*n) zeros(n,1)]);
    model.rhs = [model.rhs; repmat(L/2,4,1); 0];
    model.sense = [model.sense; '<'; '<'; '<'; '<'; '='];
end

% 高度约束
model.A = sparse([model.A; ones(1,n) zeros(1,n) -n]);
model.rhs = [model.rhs; 0];
model.sense = [model.sense; '='];

% 求解
result = gurobi(model);

% 输出结果
disp(result.objval);
disp(result.x);

其中，使用了 Gurobi 作为求解器，需要先安装 Gurobi 并获取 Gurobi 的许可证。在运行上述代码之前，需要提前设置好 Gurobi 的路径和许可证。