Simulink三相正弦序分量计算异常及实时幅值相位提取需求

针对你遇到的Simulink实时处理实信号时，无法用abs()和phase()直接提取幅值相位的问题，以下是三种实用的解决方案，适配实时处理场景：

1. 离散傅里叶变换(DFT)提取法 #

适合稳态正弦信号，通过缓存一段采样数据做FFT，精准提取基波分量的幅值与相位：

• Simulink配置：用Buffer模块缓存N个采样点（建议取2的幂次，如128，平衡精度与实时性），将缓存后的数据输入MATLAB Function模块。
• 核心代码：

function [mag, phase_rad] = fft_extract(u, fs, f0)
    N = length(u);
    % 加汉宁窗抑制频谱泄漏
    u_windowed = u .* hanning(N)';
    % 计算FFT并归一化
    Y = fft(u_windowed) / N;
    % 定位基波对应的FFT索引
    k = round(f0 * N / fs);
    % 实信号幅值为FFT分量的2倍
    mag = 2 * abs(Y(k+1));
    % 提取相位（弧度制）
    phase_rad = angle(Y(k+1));
end

• 注意事项：需提前传入采样频率fs与基波频率f0，缓存长度N需根据实时性要求调整。

2. 锁相环(PLL)跟踪法 #

适合动态或频率小幅波动的场景，实时跟踪信号的相位与幅值：

• Simulink配置：直接将实时采集的单相信号输入MATLAB Function模块，无需缓存。
• 核心代码：

function [mag, phase_rad] = pll_extract(u, fs, f0)
    persistent integrator pi_kp pi_ki theta_prev mag_prev
    % 初始化参数
    if isempty(integrator)
        integrator = 0;
        pi_kp = 2*pi*f0*0.1;  % 比例系数
        pi_ki = (2*pi*f0*0.1)^2/4;  % 积分系数
        theta_prev = 0;
        mag_prev = 0;
    end
    % 低通滤波提取幅值
    alpha = 2*pi*10/fs;  % 低通截止频率设为10Hz
    mag = (1-alpha)*mag_prev + alpha*abs(u);
    mag_prev = mag;
    
    % 鉴相+PI调节跟踪相位
    error = u * cos(theta_prev) - (u/mag)*sin(theta_prev);
    integrator = integrator + error/fs;
    theta = theta_prev + (pi_kp*error + pi_ki*integrator)/fs;
    phase_rad = theta;
    theta_prev = theta;
end

• 优势：无需缓存数据，响应速度快，能实时跟踪信号相位变化。

3. 同步参考坐标系(SRF)变换法 #

专为三相平衡信号设计，精度与实时性兼顾：

• Simulink配置：将三相实信号Va、Vb、Vc直接输入MATLAB Function模块。
• 核心代码：

function [mag, phase_rad] = srf_extract(Va, Vb, Vc, fs)
    persistent integrator pi_kp pi_ki theta_prev
    % 初始化参数（假设基波频率50Hz）
    if isempty(integrator)
        integrator = 0;
        pi_kp = 2*pi*50*0.2;
        pi_ki = (2*pi*50*0.2)^2/4;
        theta_prev = 0;
    end
    % Clark变换：三相→αβ坐标系
    alpha = (2/3)*(Va - 0.5*Vb - 0.5*Vc);
    beta = (2/3)*(sqrt(3)/2*Vb - sqrt(3)/2*Vc);
    % Park变换：αβ→dq坐标系
    d = alpha*cos(theta_prev) + beta*sin(theta_prev);
    q = -alpha*sin(theta_prev) + beta*cos(theta_prev);
    % PI调节跟踪相位（平衡时q分量为0）
    error = -q;
    integrator = integrator + error/fs;
    theta = theta_prev + (pi_kp*error + pi_ki*integrator)/fs;
    mag = d;  % 平衡时d轴分量即为相电压幅值
    phase_rad = theta;
    theta_prev = theta;
end

• 优势：利用三相平衡特性，无需额外幅值提取步骤，直接输出准确的幅值与相位。

补充说明 #

你之前用abs()和phase()失效的原因：Simulink中处理的是瞬时实信号，abs()仅能获取当前时刻的绝对值，而非信号的稳态幅值；phase()仅对复信号有效，直接用于实信号会返回无意义的固定值（0或π），无法反映随时间变化的相位信息。

内容的提问来源于stack exchange，提问作者Trickk