请根据以下要求，编写matlab环境中的代码：有一个22节点三相系统，节点电压、支路电流、支路功率、每个节点的负荷功率已知。利用这些条件，对单个换相器进行建模，算出某一时刻某个节点的单相负荷应该改接到哪一相上，以使三相不平衡度满足约束且系统总损耗最小。

在MATLAB环境中，我们可以构建一个数学模型来解决这个问题。首先，我们需要导入必要的库，假设我们已经有了`cvxopt`包用于求解优化问题。这里是一个基本框架的示例代码：

% 假设已有的数据
nodes = 22; % 节点总数
phases = 3; % 相数
V_base = 230; % 基波电压
I_max = 150; % 最大电流

% 负荷功率和节点信息 (简化为每个节点只有一个负荷)
loads = [your_load_power_nodes]; % 负荷功率向量
nodes_positions = [your_node_positions]; % 节点相位连接情况

% 支路电流和电压（实际场景需要提供）
currents = [your_currents]; % 已知电流
voltages = [your_voltages]; % 已知电压

% 计算不平衡度和损耗函数
imbalance_limit = 0.1; % 允许的不平衡度
imbalance = @(x) norm(sum(x(nodes_positions,:) .* currents, 2) - loads);
losses = @(x) V_base^2 * sum(x(nodes_positions,:) .* abs(currents));

% 创建优化变量
x = cvx.Variable(nodes, phases); % x(i,j)为1表示接至相j，0表示不接

% 构造优化问题
prob = Problem(Minimize(losses(x)), ...
                [imbalance(x) <= imbalance_limit]);

% 添加约束，保证每个节点只能接入一相
for i = 1:nodes
    prob.Constraints(x(i,:) == 1);
end

% 解决优化问题
try
    solution = solve(prob);
catch ME
    disp('Optimization failed: ' + ME.message);
    return;
end

% 输出结果
selected_phases = find(solution > 0.5); % 查找接线方案
total_loss = fval(prob); % 总损耗
imbalance_value = imbalance(solution);

% 结果分析
disp(['Selected phase for node: ', num2str(selected_phases), ' Loss: ', num2str(total_loss), ' Imbalance: ', num2str(imbalance_value)]);

请替换上述代码中的`[your_load_power_nodes]`、`[your_node_positions]`、`[your_currents]` 和 `[your_voltages]`为实际的数据。请注意，这只是一个基础模型，并未考虑到换相器的具体细节。对于换相器的控制策略，你可能还需要额外的步骤。