matlab中IE333节点求THD

节点电压的总谐波失真率（THD）是计算电力系统中非线性负载引起的电压畸变的一种常用指标。在MATLAB中，可以通过以下步骤计算IE333节点的THD：

1. 首先，使用MATLAB中的FFT函数对节点电压进行傅里叶变换，得到频域中各个谐波分量的幅值和相位角。

2. 然后，根据公式THD = sqrt(（V2^2 + V3^2 + ... + Vn^2）/ V1^2) * 100% 计算总谐波失真率，其中V1表示基波电压幅值，V2、V3、...、Vn表示各个谐波分量的电压幅值。

下面是一个简单的MATLAB代码示例，用于计算IE333节点的THD。

% 节点电压数据
v_ie333 = [230.00 -2.8980+228.5864i -2.8980-228.5864i -0.3576+0.2295i -0.3576-0.2295i -0.1661+0.0059i -0.1661-0.0059i 0.0976+0.0371i 0.0976-0.0371i -0.0338+0.0134i -0.0338-0.0134i];

% 计算基波电压幅值
v1 = abs(v_ie333(1));

% 计算各个谐波分量的电压幅值
v2 = abs(v_ie333(2));
v3 = abs(v_ie333(3));
v4 = abs(v_ie333(4));
v5 = abs(v_ie333(5));
v6 = abs(v_ie333(6));
v7 = abs(v_ie333(7));
v8 = abs(v_ie333(8));
v9 = abs(v_ie333(9));
v10 = abs(v_ie333(10));
v11 = abs(v_ie333(11));

% 计算总谐波失真率
THD = sqrt((v2^2 + v3^2 + v4^2 + v5^2 + v6^2 + v7^2 + v8^2 + v9^2 + v10^2 + v11^2) / v1^2) * 100;
disp(['IE333节点的THD为：' num2str(THD) '%']);

运行上述代码，即可得到IE333节点的THD值。需要注意的是，该代码中采用了硬编码的方式获取各个谐波分量的电压幅值，实际应用时需要根据具体情况进行修改。