有MATLAB的代码吗

当然可以，以下是基于 MATLAB 实现的室内定位泰勒展开迭代算法示例代码：

function [x, y] = indoor_positioning_taylor(wifi_data, wifi_info, pos_init, max_iter, eps)
% wifi_data: 接收到的WiFi信号强度数据
% wifi_info: 参考WiFi信号强度数据和对应位置信息
% pos_init: 初始位置
% max_iter: 最大迭代次数
% eps: 迭代精度

% 定义非线性函数
function f = fun(x, y, r, rss_ref)
    f = (x - r(1))^2 + (y - r(2))^2 - (10^((rss_ref - 30) / 10))^2;
end

% 定义雅可比矩阵
function J = jac(x, y, r, rss_ref)
    J = [
        2 * (x - r(1)), 2 * (y - r(2)), -2 * 10^((rss_ref - 30) / 10) / log(10) * 10^((rss_ref - 30) / 10) * (x - r(1));
        2 * (x - r(1)), 2 * (y - r(2)), -2 * 10^((rss_ref - 30) / 10) / log(10) * 10^((rss_ref - 30) / 10) * (y - r(2))
    ];
end

% 初始化
pos = pos_init;
A = zeros(length(wifi_info), 3);
b = zeros(length(wifi_info), 1);

% 迭代
for i = 1 : max_iter
    % 构建线性方程组
    for j = 1 : length(wifi_info)
        rss_ref = wifi_info(j, 1);
        pos_ref = wifi_info(j, 2 : 3);
        A(j, :) = jac(pos(1), pos(2), pos_ref, rss_ref)(1 : 2, :);
        b(j, :) = -jac(pos(1), pos(2), pos_ref, rss_ref)(3, :);
    end

    % 加权最小二乘法求解
    x = (A' * diag(wifi_data) * A) \ (A' * diag(wifi_data) * b);

    % 更新位置
    pos = pos + x';

    % 检查迭代精度
    if norm(x) < eps
        break
    end
end

% 返回定位结果
x = pos(1);
y = pos(2);
end

该算法和 Python 版本的实现基本相同，只是语法有些不同。需要注意的是，MATLAB 中的数组默认是列向量，因此需要对矩阵进行转置操作。具体使用时，只需将该代码保存为一个 .m 文件，然后在 MATLAB 中调用即可。