核心原理都是基于采样和数字信号处理算法，最核心的算法是离散傅里叶变换（DFT） 及其快速算法快速傅里叶变换（FFT）。

---

1. 基础：数据采集与前提

在计算之前，DSP系统需要通过前端电路获取电压和电流信号。

1. 传感器：使用电压互感器（PT）和电流互感器（CT）或霍尔传感器，将高电压、大电流转换为DSP可处理的小电压信号。
2. 信号调理：包括滤波、放大，将信号调整到适合ADC（模数转换器）的范围内。
3. ADC采样：DSP控制ADC以固定的采样频率 ( f_s ) 对电压 ( u(t) ) 和电流 ( i(t) ) 信号进行同步采样。这是最关键的一步，必须保证电压和电流通道的采样时刻完全一致，否则会引入相位误差。

假设我们得到两组离散的序列：

• 电压序列：( u[n] ), ( n = 0, 1, 2, …, N-1 )
• 电流序列：( i[n] ), ( n = 0, 1, 2, …, N-1 )
  其中 ( N ) 是一个周期内的采样点数。

---

2. 单相电路计算

方法一：时域法（基于瞬时功率理论）

这是最直接、最常用的方法，计算量小，适用于嵌入式系统。

1. 瞬时功率：
   ( p[n] = u[n] \times i[n] )

2. 有功功率（P）：
   有功功率是瞬时功率在一个周期内的平均值。
   ( P = \frac{1}{N} \sum_{n=0}^{N-1} p[n] = \frac{1}{N} \sum_{n=0}^{N-1} (u[n] \times i[n]) )
   物理意义：这直接对应了功率的定义，即单位时间内做功的平均值。

3. 视在功率（S）：
   视在功率是电压有效值和电流有效值的乘积。

   • 电压有效值：( U_{rms} = \sqrt{\frac{1}{N} \sum_{n=0}^{N-1} (u[n])^2} )
   • 电流有效值：( I_{rms} = \sqrt{\frac{1}{N} \sum_{n=0}^{N-1} (i[n])^2} )
   • 视在功率：( S = U_{rms} \times I_{rms} )

4. 无功功率（Q）：
   可以通过视在功率和有功功率计算：
   ( Q = \sqrt{S^2 - P^2} )
   注意：这种方法在谐波存在时，会包含谐波产生的畸变功率，计算结果通常是“广义无功功率”。

5. 功率因数（PF）：
   ( PF = \frac{P}{S} )

方法二：频域法（DFT/FFT法）

这种方法更强大，可以同时分析谐波，但计算量更大。

1. DFT分析：
   对电压序列 ( u[n] ) 和电流序列 ( i[n] ) 分别进行DFT（通常用FFT实现），得到它们的频谱。

   • 电压频谱：( U[k] )，包含幅值 ( U_m[k] ) 和相位 ( \phi_u[k] )
   • 电流频谱：( I[k] )，包含幅值 ( I_m[k] ) 和相位 ( \phi_i[k] )
     其中 ( k=0,1,2,…,N-1 ) 代表谐波次数。( k=1 ) 对应基波。

2. 各次谐波的有功功率：
   ( P_k = \frac{1}{2} U_m[k] I_m[k] \cos(\phi_u[k] - \phi_i[k]) = U_{rms}[k] I_{rms}[k] \cos(\theta_k) )
   其中 ( \theta_k = \phi_u[k] - \phi_i[k] ) 是第 ( k ) 次谐波的相位差。

3. 总有功功率（P）：
   ( P = \sum_{k=1}^{M} P_k )
   ( M ) 是考虑的最高次谐波，通常 ( M \leq N/2 )。

4. 各次谐波的无功功率：
   ( Q_k = U_{rms}[k] I_{rms}[k] \sin(\theta_k) )

5. 总无功功率（Q）：
   ( Q = \sum_{k=1}^{M} Q_k )
   注意：这是基波+谐波的总无功功率。

6. 视在功率（S）：
   ( S = U_{rms} \times I_{rms} )
   （这里的 ( U_{rms} ) 和 ( I_{rms} ) 是包含所有谐波的总有效值）

7. 功率因数（PF）：
   ( PF = \frac{P}{S} )

8. 谐波分析：

   • 总谐波畸变率（THD）：
     • 电压THD：( THD_U = \frac{\sqrt{\sum_{k=2}^{M} U_{rms}[k]^2}}{U_{rms}[1]} \times 100% )
     • 电流THD：( THD_I = \frac{\sqrt{\sum_{k=2}^{M} I_{rms}[k]^2}}{I_{rms}[1]} \times 100% )
   • 各次谐波含有率（HR）：
     • 第 ( k ) 次电压谐波含有率：( HR_U[k] = \frac{U_{rms}[k]}{U_{rms}[1]} \times 100% )
     • 第 ( k ) 次电流谐波含有率：( HR_I[k] = \frac{I_{rms}[k]}{I_{rms}[1]} \times 100% )

---

3. 三相电路计算

三相电路分为平衡系统和不平衡系统，计算方法有所不同。

情况一：三相三线制系统

常用**“二瓦计法”** 的理论延伸，即计算两个线电压和两个线电流的功率。

1. 采样：采样 ( u_{ab}[n] ), ( u_{cb}[n] ), ( i_a[n] ), ( i_c[n] )。（理论上任意两相线电压和对应的线电流都可以）
2. 瞬时功率：
   ( p_{3\phi}[n] = u_{ab}[n] \times i_a[n] + u_{cb}[n] \times i_c[n] )
3. 总有功功率（P）：
   ( P = \frac{1}{N} \sum_{n=0}^{N-1} p_{3\phi}[n] )

情况二：三相四线制系统（或通用方法）

直接对每一相功率求和，并考虑中线。

1. 采样：采样 ( u_a[n] ), ( u_b[n] ), ( u_c[n] ), ( i_a[n] ), ( i_b[n] ), ( i_c[n] )，如果有中线，还需采样 ( i_n[n] )。
2. 瞬时功率：
   ( p_{3\phi}[n] = u_a[n] \times i_a[n] + u_b[n] \times i_b[n] + u_c[n] \times i_c[n] )
3. 总有功功率（P）：
   ( P = \frac{1}{N} \sum_{n=0}^{N-1} p_{3\phi}[n] )

无功功率和功率因数的计算

对于三相系统，通常先计算各相的有功功率 ( P_a, P_b, P_c ) 和无功功率 ( Q_a, Q_b, Q_c )（使用单相时域法或频域法），然后：

• 总有功功率：( P_{total} = P_a + P_b + P_c )
• 总无功功率：( Q_{total} = Q_a + Q_b + Q_c )
• 总视在功率：( S_{total} = \sqrt{P_{total}^2 + Q_{total}^2} )
  （注意：在三相不平衡系统中，这种定义与 ( S = S_a + S_b + S_c ) 不等效，但前者更常用）
• 总功率因数：( PF_{total} = \frac{P_{total}}{S_{total}} )

频域法（对称分量法/FFT法）

对三相电压电流分别进行FFT，得到各相各次谐波的幅值和相位后，可以：

1. 计算各相、各次谐波的 ( P_k ) 和 ( Q_k )。
2. 将三相的同一谐波次的 ( P_k ) 相加，得到系统该次谐波的总有功功率。
3. 将所有谐波次的有功功率相加，得到系统总有功功率。
4. 无功功率计算同理。
5. 还可以进行正序、负序、零序分量的分析，这对于电能质量评估非常重要。

---

总结：DSP实现的流程

无论是单相还是三相，在DSP中的通用流程可以概括为以下框图：

flowchart TD
    A[模拟信号 u(t), i(t)] --> B[传感器与信号调理];
    B --> C[ADC同步采样];
    C --> D[“得到离散序列<br>u[n], i[n]”];
    D --> E{选择计算方法};
    E -- 时域法 --> F[“计算瞬时功率 p[n] = u[n] × i[n]”];
    F --> G[“平均求得有功功率 P”];
    G --> H[计算电压/电流有效值];
    H --> I[“计算视在功率 S<br>及功率因数 PF”];
    I --> J[“间接计算无功功率 Q”];
    
    E -- 频域法 --> K[对u[n]和i[n]进行FFT];
    K --> L[“得到频谱 Um[k], φu[k]<br>Im[k], φi[k]”];
    L --> M[“计算各次谐波功率<br>Pk, Qk”];
    M --> N[“计算总有功 P<br>总无功 Q”];
    N --> O[计算总视在功率S及PF];
    L --> P[计算THD及各次谐波含有率];
    
    J & O & P --> Q[输出最终结果];

核心要点：

• 时域法：简单、快速、资源占用少，适合对谐波分析要求不高的场合，可以准确计算有功功率和功率因数。
• 频域法：功能强大，能同时得到有功、无功、功率因数以及完整的谐波频谱和THD，是电能质量分析的首选，但计算复杂度高。
• 同步采样：是保证计算精度的生命线。
• 窗函数：在实际应用中，如果采样未整周期同步，会存在频谱泄漏，需要加窗函数（如汉宁窗、海明窗）来改善FFT效果。
• 计算精度：取决于ADC分辨率、采样率、传感器精度和算法稳定性。

希望这份详细的解释能帮助你理解DSP在电力参数计算中的应用。