模拟信号求解相位差（1）_simulink测量相位差-CSDN博客

本文链接：https://blog.csdn.net/abcxa1986/article/details/78180608

本文基于《正弦相位的取样数值测量——高庆》在MATLAB上进行仿真实验，也可看成对此方法的验证。

1.理论基础

假设输入的两个模拟信号的表达式分别为：
$\mu_1(t)=\boldsymbol{U_{1m}}\sin(\omega t+\theta_1)$
$\mu_2(t)=\boldsymbol{U_{2m}}\sin(\omega t+\theta_2)$
将模拟信号转化为数字信号的采样序列：
$\mu_1(n)=\boldsymbol{U_{1m}}\sin(\cfrac{2\pi n}{N}+\theta_1)$
$\mu_2(n)=\boldsymbol{U_{2m}}\sin(\cfrac{2\pi n}{N}+\theta_2)$
其中 $N=\cfrac{T}{T_s}$ 为模拟信号一个周期内的采样点数， $T_s$ 是采样器的取样周期
从数学的角度可以得到以下推导：
$\sum_{n=1}^{N}U_1(n)U_2(n)=2\sum_{n=1}^{N}U_m\sin(\cfrac{2\pi n}{N}+\theta_1)\sin(\cfrac{2\pi n}{N}+\theta_2)\\=\sum_{n=1}^{N}U_m\cos(\theta_1-\theta_2)-\sum_{n=1}^{N}U_m\cos(\cfrac{4\pi n}{N}+\theta_1+\theta_2)$
其中用到
$\sin\alpha\sin\beta=\cfrac{1}{2}\cos(\alpha-\beta)-\cfrac{1}{2}\cos(\alpha+\beta)$
$U_m=\cfrac{1}{2}U_{1m}U_{2m}$ ，其中 $U_{1m}=Max{(U_1(n))},U_{2m}=Max{(U_2(n))}$ ，这样就可以通过对采样序列的查询来得到 $U_{1m}U_{2m}$ 的值。
从数学方面可以知道，余弦信号在一个周期内整数倍取样和为0，所以上式可以变成 $上式=\sum_{n=1}^{N}U_m\cos(\theta_1-\theta_2)=NU_m\cos\phi$
其中 $\phi=\theta_1-\theta_2$ 就是所求两个信号的相位差
求两个采样序列乘积的平均值为
$\overline U=\cfrac{1}{N}\sum_{n=1}^NU_1(n)U_2(n)=U_m\cos\phi$
这表明满足采样定理，且 $\cfrac{T}{T_s}=N$ 时，一个周期内的正弦信号瞬间取样值的乘积的平均值正比于它们之间的相位差 $\phi$ 的余弦值，即 $\cos\phi=\cfrac{\overline U}{U_m}$
又根据泰勒公式可以得到：
$\cos\phi=1-\cfrac{\phi^2}{2!}+\cfrac{\phi^4}{4!}-\cfrac{\phi^6}{6!}+\cdots+(-1)^n\cfrac{\phi^{2n}}{(2n)!}\\+(-1)^{n+1}\cfrac{\phi^{2n+2}}{(2n+2)!}\cos\zeta\phi,0<\zeta<1$
可以近似表示为 $\cos\phi\approx1-\cfrac{\phi^2}{2!}\cdots$
我们可以得出被测相位差 $\phi$ ，
$\phi=\sqrt{2-2\cos\phi}=\sqrt{2-2\cfrac{\overline U}{U_m}}(取正值)$
但这样得出的相位差会有较大误差，我们需要进行校正，所采用的方法是直接迭代法，取泰勒展开式的前四项：
$\cos\phi\approx1-\cfrac{\phi^2}{2}+\cfrac{\phi^4}{24}-\cfrac{\phi^6}{720}$
令 $y=\phi^2$ 得
$\cos\phi=1-\cfrac{y}{2}+\cfrac{y^2}{24}-\cfrac{y^3}{720}$
得到迭代公式
$y_i=2(1-\cfrac{\overline U}{U_m})+\cfrac{y_{i-1}^2}{24}-\cfrac{y_{i-1}^3}{720}$
当 $y_i-y_{i-1}<a(a为所给定的精确系数)$ 时，则被测相位差为 $\phi=\sqrt{y_i}(取正值)$

2.流程图

3.MATLAB仿真实现

clear all;

%N = 200;                                                                  %采样数量
F = 1;                                                                     %被采样信号频率1Hz
Ts = 0.01;                                                                 %采样间隔 
Fs = 1/Ts;                                                                 %采样频率
%T = N*Ts;                                                                 %采样时长

%n = 0:N-1;                                                                %时域采样序列(N个采样)
NP = floor( (1/F)/(Ts) );                                                  %1个周期采样点数,floor()向负无穷大方向取整
n = 0:NP-1;

nTs = n*Ts;                                                                %时域采样时间序列
%对下列信号进行研究
g1=inline('2*sin(t)');
g2=inline('3*sin(t+pi/4)');

%采样序列
f1=g1(2*pi*nTs);
f2=g2(2*pi*nTs);
%f2转置
f2=f2.';

U1m=max(max(f1));                                                          %求出f1中的最大值
U2m=max(max(f2));                                                          %求出f2中的最大值
Um=U1m*U2m/2;                                                              %计算Um
Uv=(f1*f2)/NP;                                                             %计算平均值

angle=sqrt(2-2*(Uv/Um));                                                       %获得初始相位差
%下面通过迭代提高相位差的精度
Y0=angle*angle;
a=0.00001;                                                                 %要求

for i=1:10000000000
    Y1=2*(1-Uv/Um)+Y0*Y0/12-Y0*Y0*Y0/360;                                  %迭代公式
    if(Y1-Y0)<a                                                            %若两次迭代数值小于所要求精度，则停止迭代
        break
    else
        Y0=Y1;
    end
end
angle=sqrt(Y1)