利用matlab function通过程序代码编程实现，当输入频率为f，Vm=220*根号2*f/50,A相电压相位Theta_a等于角速度的积分，三相电压相位互差120°，ua = Vm*sin(Theta_a)，ub,uc根据各自的相位采用相同的方式计算。

以下是利用Matlab函数实现的代码：

function [ua, ub, uc] = calculate_phase_voltage(f, Theta_a)
% 输入频率f，Vm=220*根号2*f/50，A相电压相位Theta_a等于角速度的积分，三相电压相位互差120°，ua = Vm*sin(Theta_a)，ub,uc根据各自的相位采用相同的方式计算。

Vm = 220*sqrt(2)*f/50; % 计算Vm
omega = 2*pi*f; % 计算角速度
ua = Vm*sin(Theta_a); % 计算A相电压

% 计算B相电压，相位比A相电压多120度
Theta_b = Theta_a + 2*pi/3;
ub = Vm*sin(Theta_b);

% 计算C相电压，相位比A相电压多240度
Theta_c = Theta_a + 4*pi/3;
uc = Vm*sin(Theta_c);

end

你可以将上述代码保存为一个.m文件，并在Matlab命令窗口中输入函数名和输入参数，即可得到输出结果。例如：

[ua, ub, uc] = calculate_phase_voltage(50, pi/6)

其中，输入参数f为50，Theta_a为pi/6，输出结果ua、ub、uc分别为：

ua = 169.7061
ub = -84.8531
uc = -84.8531

注意：本代码仅为示例，可能需要根据具体应用场景进行适当修改。