目录

• 0 分类
• 1、PLL
• • 1.1 系统结构
• 2、SOGI-based PLL 未完待续
• • 2.1 结构理解
  • 2.2 C语言实现
  • • 2.2.1 通过控制框图直接实现
    • 2.2.2 通过离散SOGI的传递函数来实现
    • • 2.2.2.1 传递函数
      • 2.2.2.2 离散化
      • 2.2.2.3 SOGI的代码实现
  • 2.3 仿真验证
  • 2.3.1 代码验证
  • 2.3.2 simulink验证-模块搭建
  • 2.3.3 simulink验证-传递函数
  • 2.3.4 simulink验证-C语言

0 分类

闭环锁相环技术的主要方案有：

https://blog.csdn.net/weixin_31054839/article/details/113707862

1、PLL

• 输入为：电网交流电压
• 输出为：同相交流电压、正交信号

锁相环(phase locked loop)是一种利用相位同步产生的电压，去调谐压控振荡器以产生目标频率的负反馈控制系统。学过自动控制原理的人都知道，这是一种典型的反馈控制电路，利用外部输入的参考信号控制环路内部振荡信号的频率和相位，实现输出信号频率对输入信号频率的自动跟踪，一般用于闭环跟踪电路。是无线电发射中使频率较为稳定的一种方法，主要有VCO（压控振荡器）和PLL IC（锁相环集成电路），压控振荡器给出一个信号，一部分作为输出，另一部分通过分频与PLL IC所产生的本振信号作相位比较，为了保持频率不变，就要求相位差不发生改变，如果有相位差的变化，则PLL IC的电压输出端的电压发生变化，去控制VCO，直到相位差恢复，达到锁相的目的。能使受控振荡器的频率和相位均与输入信号保持确定关系的闭环电子电路。

1.1 系统结构

由三部分构成：我想知道最后反馈的是什么！！！！！
答：最后反馈的是相位，通过使输出相位控制在目标值wt（也就是输出的wt和输入的wt一直保持一致，那么两者的相位差一直保持为初始相位差），最后达到频率相等，那么最后。

• 鉴相器（PD,Phase Detector）
• 环路滤波器（LF,Loop Filter）
• 压控振荡器（VCO,Voltage Controlled Oscillator）
  有的资料在反馈回路中还包括分频器1/N¹

环节	输入	输出
鉴相器	交流电压	和
环路滤波器(也就是PI环节)		未知
VCO(也就是积分环节)	未知

PLL只是输出到，其他的都是其他环节处理的了！
最后将变为PWM波是其他的环节。

总的来说，PLL就是调节输入输出的频差为0的，鉴相器就是通过相乘求得相位差**（还是频差呢？）**；环路滤波器就是积分器，将输入的交流信号积分成直流信号，最后输出VSO的控制电压信号；VCO最后控制频差。
当输入输出相位差恒定时（也就是频率相等），

• 模拟乘法器
  模拟乘法器的输出电压为：
  

• 低通滤波器LF
  低通滤波器LF将上式中的和频分量滤掉，剩下的差频分量作为压控振荡器的输入控制电压。即为：
  

• 对两边求微分，可得频差的关系式为：
  

• VCO
  压控振荡器的压控特性如图8-4-3所示，该特性说明压控振荡器的振荡频率ωu以ω0为中心，随输入信号电压uc（t）的变化而变化。该特性的表达式为
  
  

https://www.electricaltechnology.org/2019/02/phase-locked-loop-pll.html
https://www.cirmall.com/articles/31424/
https://www.eet-china.com/mp/a119921.html

2、SOGI-based PLL 未完待续

未完待续
[1]杨勇,李世华,朱彬彬,陶雪慧,赵方平,王志冰.基于双阶广义积分的单相光伏并网逆变器灵活功率控制[J].电网技术,2014,38(04):952-958.DOI:10.13335/j.1000-3673.pst.2014.04.020
[1]李永刚,杨雅菲,姜玉霞.弱电网下SOGI-PLL参数对LCL型并网逆变器稳定性影响及参数优化[J].河南理工大学学报(自然科学版),2021,40(04):126-134.DOI:10.16186/j.cnki.1673-9787.2020030094 .
https://liquidsdr.org/blog/pll-howto/
https://imperix.com/doc/implementation/sogi-pll#Simulink-implementation-of-a-SOGI-PLL

2.1 结构理解

首先SOGI只是上文PLL中鉴相器的一部分。
鉴相器：交流电压 => 和
SOGI：交流电压 => 和
因此，SOGI + 【 => 变换】= 鉴相器

下图中阴影部分为SOGI(Second Oder Generalized Integrator)，剩下的部分为滤波器部分。

另外一个理解方式为：
1、这个是SOGI的结构

2.2 C语言实现

2.2.1 通过控制框图直接实现

1、第一种

x_a=(ualfa_0 - SOGI_ualfa.integral_2)*k_sogi;
SOGI_ualfa.integral_2=SOGI_ualfa.integral_2+w0_sogi*(x_a-SOGI_ualfa.integral_3)*T_ADC;
ualfa_1 = SOGI_ualfa.integral_2;
SOGI_ualfa.integral_3 = SOGI_ualfa.integral_3 + w0_sogi*SOGI_ualfa.integral_2*T_ADC;
ualfa_q = SOGI_ualfa.integral_3;
//参考来源 => https://blog.csdn.net/sinat_40597468/article/details/125582569

其中：
ualfa_0对应图xxx中；
SOGI_ualfa.integral_2对应图xxx中的；
SOGI_ualfa.integral_3对应图图xxx中的；
ualfa_1对应图xxx中的；
ualfa_q对应图xxx中的;
T_ADC为采样时间；
这里：k_sogi=1；w0_sogi = 314.15926；

2、第二种。
这个是包含了SOGI + 后续的、转换 + PLL的整个过程。

void PLL(float v_s1)
{
     //SOGI正交变换得到两个正交量v1_s(α)、v2_s(β)
     v0_s = ((v_s1-v1_s)*K-v2_s)*W;
     v1_s += v0_s*Ts;//得到V1(s)(α)
     v3_s += v1_s*Ts;
     v2_s = v3_s*W;//得到V2(s)(β)
     //αβ/dq变换及PI运算
     v_q = sin(phase_angle)*(-1)*v1_s + cos(phase_angle)*v2_s;
     v_d = <u>cos</u>(phase_angle)*v1_s + sin(phase_angle)*v2_s;//如果不需要就删掉
     d_v_q = v_q_ref - v_q;//error
     v_q_temp += (d_v_q-d_v_q1) * Kp_PLL + d_v_q * Ki_PLL ;
     d_v_q1 = d_v_q;
     if(v_q_temp > 20 * pi)
         v_q_temp = 20 * pi;
     else if(v_q_temp < -20*pi)
         v_q_temp = -20 * pi;   //用于对误差参数的限幅，防止误差参数过大造成不能输出正常大小的phase_angle最终不能正常进行锁相
     //由转速得到相角，并取模
     phase_angle += Ts * (100 * pi + v_q_temp);//这里输出的角度是累计误差的角度，同时也是信号当前的角度
     if(phase_angle > 2 * pi)
     {
         phase_angle -= 2 * pi;
     }
     if(phase_angle < 0)
     {
         phase_angle += 2 * pi;
     }
     angles[resultsIndex] = phase_angle;
}
参考来源 => https://zhuanlan.zhihu.com/p/358573305

2.2.2 通过离散SOGI的传递函数来实现

参考来源 => https://www.cnblogs.com/swear/p/12682551.html

2.2.2.1 传递函数

由SOGI的结构图可以得到传递函数：（没有推过，但是好多来源都是这样写的该传递函数）

• 同相传递函数：
  
• 正交信号传递函数：
  

2.2.2.2 离散化

离散化有两种形式：
①利用不同离散化方法的公式来手动离散，从s域 => z域 => 差分方程 来进行离散化；
②利用MATLAB离散。

参考：https://blog.csdn.net/qq_42370406/article/details/127061779?spm=1001.2014.3001.5501

• 手动离散：
  离散方法为双线性变换：
  
  可得，同相传递函数和正交信号传递函数的 z域形式为(过程参考该节后面的来源)：
  
  
  式(7)和式(8)可化简为：
  
  
  根据式(9)和式(10)，得同相传递函数和正交信号传递函数的差分方程形式为(过程参考该节后面的来源)：
  
  

• MATLAB离散：
  离散代码为：

%%定义s为传递函数
s = tf('s');    

%%定义各参数
k = 0.5;
Wn = 100*pi;    %%50Hz
Ts = 1e-4;          %%10kHz

%%写出传递函数
Hd_s = k*Wn*s/(s^2+k*Wn*s+Wn^2);
Hq_s = k*Wn^2/(s^2+k*Wn*s+Wn^2);

Hd_z = c2d(Hd_s,Ts,'tustin')
Hq_z = c2d(Hq_s,Ts,'tustin')

运行结果：

Hd_z =
 
  0.007791 z^2 - 0.007791
  -----------------------
  z^2 - 1.983 z + 0.9844
 
Sample time: 0.0001 seconds
Discrete-time transfer function.


Hq_z =
 
  0.0001224 z^2 + 0.0002448 z + 0.0001224
  ---------------------------------------
          z^2 - 1.983 z + 0.9844
 
Sample time: 0.0001 seconds
Discrete-time transfer function.

MATLAB离散化的系数为：

2.2.2.3 SOGI的代码实现

无论通过上述两种哪种方式，最后都求得差分方程的系数，因此接下来可以通过代码来进行实现SOGI的结果：

%%
%%函数声明
function [uo,quo]   = Orthogonal_Generator(ui,Ts,w,k)

%%
%%定义各中间变量
persistent x;
persistent y;
persistent temp;
persistent b0;
persistent b2;
persistent a1;
persistent a2;
persistent qb0;
persistent qb1;
persistent qb2;
persistent u0;          %%代表ui(k)
persistent u1;          %%代表ui(k-1)
persistent u2;          %%代表ui(k-2)
persistent osg_u0;  %%代表uo(k)
persistent osg_u1;  %%代表uo(k-1)
persistent osg_u2;  %%代表uo(k-2)
persistent osg_qu0; %%代表uqo(k)
persistent osg_qu1; %%代表uqo(k-1)
persistent osg_qu2; %%代表uqo(k-2)

%%
%%初始化各中间变量
if isempty(x)     x= 0; 
end
if isempty(y)     y= 0; 
end
if isempty(temp)  temp= 0; 
end
if isempty(b0)    b0= 0; 
end
if isempty(b2)    b2= 0; 
end
if isempty(a1)    a1= 0; 
end
if isempty(a2)    a2= 0; 
end
if isempty(qb0)   qb0= 0; 
end
if isempty(qb1)   qb1= 0; 
end
if isempty(qb2)   qb2= 0; 
end
if isempty(u0)    u0= 0; 
end
if isempty(u1)    u1= 0; 
end
if isempty(u2)    u2= 0; 
end
if isempty(osg_u0)   osg_u0= 0; 
end
if isempty(osg_u1)   osg_u1= 0; 
end
if isempty(osg_u2)   osg_u2= 0; 
end
if isempty(osg_qu0)   osg_qu0= 0; 
end
if isempty(osg_qu1)   osg_qu1= 0; 
end
if isempty(osg_qu2)   osg_qu2= 0; 
end

%%
%%各系数赋值
x      = 2*k*w*Ts;
y      = w*Ts*w*Ts;
temp   = 1/(x+y+4.0);
b0     = x*temp;
b2     = (-1.0)*b0;
a1     = (2.0)*(4.0-y)*temp;
a2     = (x-y-4)*temp;
qb0    = (k*y)*temp;
qb1    = qb0*(2.0);
qb2    = qb0;

%%
%%计算过程，对应式11离散序列
u0            = ui;
osg_u0        = (b0*(u0-u2)) + (a1*osg_u1) + (a2*osg_u2);
osg_u2        = osg_u1;
osg_u1        = osg_u0;
%%对应式12离散序列
osg_qu0       = (qb0*u0) + (qb1*u1) + (qb2*u2) + (a1*osg_qu1) + (a2*osg_qu2);
osg_qu2       = osg_qu1;
osg_qu1       = osg_qu0;
%%更新序列值
u2            = u1;
u1            = u0;
%%输出
uo            =osg_u0;
quo           =osg_qu0;

2.2.2节 通过离散传递函数来实现
参考来源 => https://www.cnblogs.com/swear/p/12682551.html

可以参考：https://blog.csdn.net/qq_42370406/article/details/127061779?spm=1001.2014.3001.5501

2.3 仿真验证

2.3.1 代码验证

对于d轴，
从增益中可以看出：1、观察单个曲线，在基频50Hz 处输出信号，相对于输入信号的相移为0。这也就解释了为什么v’ 与v 同幅同相位；偏离50Hz时，增益会立刻衰减，这也就是为什么利用SOGI能实现滤波的原因。 这种算法能够实现，当输出电压AD 采样存在干扰时，也能够实现准确的dq 分量计算。2、观察多个曲线，k是控制宽带的，k越大，宽带越大。

对于q轴的多个图像：k控制的是增益放大或衰减的速度，以及相位滞后的程度。k太小，q轴会有宽带，且k越小，宽带越窄。

以上表格的内容是通过传递函数算出的，下面将通过simulink的仿真对该数据进行验证。

2.3.2 simulink验证-模块搭建

模型：

输入为50Hz，幅值为10的正弦波时，对于幅值：dq轴两个波形都是10；对于相位：d轴和输入波同相位，q轴滞后d轴90°。

2.3.3 simulink验证-传递函数

2.3.4 simulink验证-C语言