IEC 61000-4-7谐波计算(Matlab代码实现）

程序辅导帮

于 2024-09-17 10:51:14 发布

阅读量245

点赞数 12

文章标签： matlab 开发语言

本文链接：https://blog.csdn.net/2201_75454341/article/details/142312371

版权

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

⛳️赠与读者

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现

⛳️赠与读者

👨‍💻做科研，涉及到一个深在的思想系统，需要科研者逻辑缜密，踏实认真，但是不能只是努力，很多时候借力比努力更重要，然后还要有仰望星空的创新点和启发点。当哲学课上老师问你什么是科学，什么是电的时候，不要觉得这些问题搞笑。哲学是科学之母，哲学就是追究终极问题，寻找那些不言自明只有小孩子会问的但是你却回答不出来的问题。建议读者按目录次序逐一浏览，免得骤然跌入幽暗的迷宫找不到来时的路，它不足为你揭示全部问题的答案，但若能让人胸中升起一朵朵疑云，也未尝不会酿成晚霞斑斓的别一番景致，万一它居然给你带来了一场精神世界的苦雨，那就借机洗刷一下原来存放在那儿的“躺平”上的尘埃吧。

或许，雨过云收，神驰的天地更清朗.......🔎🔎🔎

💥1 概述

此函数将根据以下内容计算谐波频谱值
IEC 61000-4-7标准，该标准规定了所研究系统的10分钟电压和电流测量信号，并计算了测量信号（在这种情况下为50 Hz信号）每10个周期的谐波平均值
此函数的输入为：
1- 测量信号的标称值：频率、电压和电流。例如，f = 50%Hz
Vn = 400%V（电网的相间电压）
Inom=475%A（被测系统的标称电流）
2- 测量的电压和电流信号及其时间
测量信号（U1、U2、U3、I1、I2、I3、U1_时间、U2_时间、U3_时间、I1_时间、I2_时间、I3_时间）
3-数字是1或0的输入，以激活显示在该函数中计算的不同谐波频谱数字的数字。
PS：此功能假设您具有不同测量信号的相同采样率，并且所有测量信号具有相同的长度。
此函数的输出是：
1-电流和电压信号的总谐波失真
2-测量信号的K因子。（在互联网上查找K和THD之间的区别，或者尝试从这个函数中找出答案）
3-三个相位测量电压和电流的不同光谱最大值
4-最后三个输出是光谱顺序的x轴。

📚2 运行结果

部分代码：

function [THD_U,THD_I,K_U,K_I,MaxHarm_I,MaxIHarm_I,MaxHFHarm_I,MaxHarm_U,MaxIHarm_U,MaxHFHarm_U,Order_Harm,Order_IHarm,Order_HFHarm] = IEC_Harmonics_final(f,Vn,Inom,U1,U2,U3,I1,I2,I3,U1_time,U2_time,U3_time,I1_time,I2_time,I3_time,figures)
% This function will calculate the Harmonics spectrum values according to
% the IEC 61000-4-7 Norm which spesify the need for 10 min Voltage and
% current measured signals of the system under study and calulate the
% Harmonics mean values every 10 Periods of the measured signal (In this case 50 Hz signal)
% This Fucntion inputs are:
% 1- The nominal values of the measured signal: Frequency,
% Voltage, and the current. For example f = 50%Hz
% Vn = 400%V (Phase to Phase Voltage of the grid)
% Inom = 475%A (The nominal
% current of the system under measurement)
% 2- The measured Voltages and current signals with there time
% measurement signals (U1,U2,U3,I1,I2,I3,U1_time,U2_time,U3_time,I1_time,I2_time,I3_time)
% 3- figures is an input of 1 or 0 in order to activate showing the figures of the different harmonics spectrums figures calculated in this function.
%PS: This function assume that you have the same sampling rate of the
%different measured signals and that all measured signals has the same
%length.
Fn = f;
Un = Vn;
In = Inom;
Delta_Time_Orginal = I1_time(2)-I1_time(1);
F_Sample = abs(1/Delta_Time_Orginal);
%------------------------------------------------------------------------------------------------------
%---------------------- The Harmonics Spectrum of the complet input signal ----------------------------
%------------------------------------------------------------------------------------------------------
Tp = length(I3_time);
T_S = abs(Tp*Delta_Time_Orginal);
[nyquist_number_I1,Bezugsfrequenz_I1,Signal_amp_I1] = Fast_Harmonics (I1,F_Sample,T_S);
[nyquist_number_I2,Bezugsfrequenz_I2,Signal_amp_I2] = Fast_Harmonics (I2,F_Sample,T_S);
[nyquist_number_I3,Bezugsfrequenz_I3,Signal_amp_I3] = Fast_Harmonics (I3,F_Sample,T_S);
[nyquist_number_U1,Bezugsfrequenz_U1,Signal_amp_U1] = Fast_Harmonics (U1,F_Sample,T_S);
[nyquist_number_U2,Bezugsfrequenz_U2,Signal_amp_U2] = Fast_Harmonics (U2,F_Sample,T_S);
[nyquist_number_U3,Bezugsfrequenz_U3,Signal_amp_U3] = Fast_Harmonics (U3,F_Sample,T_S);
if figures
figure;
subplot(2,1,1);