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

为了在MATLAB环境中解决这个问题，我们需要首先组织数据结构，然后构建模型并使用优化工具箱（如`fmincon`或`optimization toolbox`）。这里提供一个基础的框架，假设我们有如下数据结构：

- `node_voltage`: 一个22x24矩阵，表示24小时内每小时每个节点的电压。
- `branch_currents`: 一个22x24x3矩阵，表示每小时每个支路的三相电流。
- `branch_power`: 一个22x24x3矩阵，表示每小时每个支路的三相功率。
- `loads`: 一个22x24矩阵，表示每小时每个节点的负荷功率。
- `voltage_constraint`: 节点电压的最大允许偏差。
- `imbalance_threshold`: 不平衡度的阈值。
- `losses_function`: 计算总损耗的函数。
- `switching_cost`: 改变相位的代价（例如，转换开关的费用）。

接下来，我们可以定义一个决策变量 `phase_shifts`，表示每小时每个节点到另一相的切换（0-2pi，对应A-B-C的循环），以及一个辅助函数 `shift_node_to_phase` 来确定如何调整负荷。

% 定义决策变量
num_hours = size(node_voltage, 2);
num_nodes = size(node_voltage, 1);
phase_shifts = zeros(num_nodes, num_hours);

% 辅助函数，确定负荷应连接的相
load_phase = @(node, hour) rem(2*pi + phase_shifts(node, hour), 2*pi); % 将角度映射回A-B-C

% 模型函数 (损失函数)
function [objective, constraint] = loss_and_constraint(model_params)
    % ... 在此处填充计算不平衡度和总损耗的代码 ...
end

% 优化模型
options = optimoptions('fmincon', 'Display', 'iter');
[phase_shifts_optimal, ~, exitflag] = fmincon(@(x) loss_and_constraint(x), phase_shifts, [], [], [], [], [], [], [], options);

请注意，你需要完成`loss_and_constraint`函数的具体内容，包括计算不平衡度（可能涉及支路功率和节点负荷的差分）、总损耗（可能涉及到`branch_power`和`loads`）、以及不平衡度和损耗的约束。