matlab fft 基波,fft 谐波分析_matlab fft频谱分析_傅里叶变换谐波分析

变频器是工业调速领域中应用最广泛的设备之一,它将工频50HZ变换成各种频率的交流电源,以实现电机的变速运行。变频器完成交-直-交变换,由整流电路、逆变电路、控制电路组成,其中整流输入部分及逆变输出部分主要由非线性电力电子元件组成,在变频器运行时,会在输入输出侧产生大量谐波。

由于变频器输出PWM波含有开关频率整数倍的高次谐波,变频器谐波分析的关键在于要求功率在低基波频率情况下拥有高采样频率、宽带宽、运算能力强、运算速度快等特点。

一变频器谐波的产生

01整流输入侧

变频器的整流输入侧一般采用三相桥式不控整流将交流电整成直流电,因为由二极管组成的整流桥的非线性特性,会产生谐波电流。一般情况,整流输入的电压波形为正弦波,输入电流波形为畸变的“马鞍波”,谐波电流会从用户变频器端通过电网传播,影响公用电网。

输入谐波电流如图1所示。“马鞍波”的波形由变频器整流拓扑结构、参数、负载的不同而不同。

图1 变频器输入电流波形图

02逆变输出侧

变频器的逆变输出侧,由SPWM波控制由IGBT等全控型器件组成的桥式电路,从而形成电压、频率可调的三相交流电。变频器输出的线电压是正弦脉宽、幅值相等的矩形波(PWM波),如图2所示。

图2 变频器输出线电压波形图

二变频器谐波分析

01变频器输入侧谐波电流分析

变频器输入电流的谐波含量由于只是整流部分非线性造成的,并且变频器的输入输出为三个幅值相等、相位相差120°的波形,根据谐波分析理论,不存在3的整数倍谐波(3、6、9、12……)以及偶次谐波(2、4、6、8……)所以一般只是5次(250Hz)、7次(350Hz)、11次(550Hz)、13次(650Hz)、17次(850Hz)……谐波。如图3所

  • 0
    点赞
  • 5
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
好的,这是一个信号处理相关的问题。首先,根据傅里叶级数的原理,一个周期锯齿波可以分解为基波和一系列谐波。基波的频率就是该周期锯齿波的频率,即f=276/200 Hz。谐波的频率是基波频率的整数倍,即2f、3f、4f、5f等。 接下来,我们需要对采样数据进行处理,以获取基波和谐波的频率和幅度比值。根据采样定理,采样频率要大于信号最高频率的两倍,因此我们选择2000Hz的采样频率,采样500个数据点,可以得到采样时间间隔为T=1/2000秒。 首先,我们需要对采样数据进行离散傅里叶变换(DFT),以获取频域信息。MATLAB中可以使用fft函数进行DFT。代码如下: ``` N = 500; % 采样点数 T = 1/2000; % 采样时间间隔 t = (0:N-1)*T; % 采样时间序列 y = sin(2*pi*276/200*t); % 周期锯齿波信号 Y = fft(y); % 离散傅里叶变换 f = (0:N-1)*(1/T)/N; % 频率向量 ``` 其中,N为采样点数,T为采样时间间隔,t为采样时间序列,y为采样数据,Y为DFT结果,f为频率向量。 接下来,我们需要设计一个FIR滤波器,以滤除除基波和2~5次谐波以外的频率分量。由于要求误差小于0.1%,我们需要对滤波器进行精细设计。可以使用firls函数进行多通带滤波器的设计。代码如下: ``` % 设计多通带滤波器 fs = 1/T; % 采样频率 f1 = 2*f; % 通带1,基波的两倍频率 f2 = 3*f; % 通带2,基波的三倍频率 f3 = 4*f; % 通带3,基波的四倍频率 f4 = 5*f; % 通带4,基波的五倍频率 f5 = 6*f; % 阻带1,避免基波的六倍频率被滤除 f6 = fs/2; % 阻带2,避免采样频率的一半被滤除 a = [1 1 0 0 0 0]; % 通带增益 dev = [0.001 0.001 0.01 0.01 0.01 0.01]; % 通带波动和阻带衰减 f = [f1 f2 f3 f4 f5 f6]/(fs/2); % 频率向量 b = firls(100,f,a,dev); % 滤波器系数 ``` 其中,fs为采样频率,f1~f6为通带和阻带的频率,a为通带的增益,dev为通带波动和阻带衰减,b为滤波器系数。由于要求误差小于0.1%,我们选择100个滤波器系数进行滤波器设计。 接下来,我们使用滤波器对采样数据进行滤波。可以使用filter函数进行滤波。代码如下: ``` % 滤波 yf = filter(b,1,y); % 滤波后的数据 Yf = fft(yf); % 滤波后的DFT结果 ``` 其中,yf为滤波后的数据,Yf为滤波后的DFT结果。 最后,我们需要计算基波和2~5次谐波的频率和幅度比值。由于要求误差在1%以内,我们可以通过找到离基波和谐波频率最近的采样频率来计算幅度比值。代码如下: ``` % 计算基波和谐波的频率和幅度比值 f_base = f(abs(f-f(1))<0.001); % 基波频率 f_harm = [2 3 4 5]*f_base; % 谐波频率 A_base = abs(Yf(abs(f-f_base)<0.001)); % 基波幅度 A_harm = abs(Yf(abs(f-f_harm)<0.001)); % 谐波幅度 ratio = A_harm./A_base; % 幅度比值 ``` 其中,f_base为基波频率,f_harm为谐波频率,A_base为基波幅度,A_harm为谐波幅度,ratio为幅度比值。由于要求误差在1%以内,我们选择离基波和谐波频率最近的采样频率进行幅度计算。 完整的代码如下: ``` N = 500; % 采样点数 T = 1/2000; % 采样时间间隔 t = (0:N-1)*T; % 采样时间序列 y = sin(2*pi*276/200*t); % 周期锯齿波信号 Y = fft(y); % 离散傅里叶变换 f = (0:N-1)*(1/T)/N; % 频率向量 % 设计多通带滤波器 fs = 1/T; % 采样频率 f1 = 2*f; % 通带1,基波的两倍频率 f2 = 3*f; % 通带2,基波的三倍频率 f3 = 4*f; % 通带3,基波的四倍频率 f4 = 5*f; % 通带4,基波的五倍频率 f5 = 6*f; % 阻带1,避免基波的六倍频率被滤除 f6 = fs/2; % 阻带2,避免采样频率的一半被滤除 a = [1 1 0 0 0 0]; % 通带增益 dev = [0.001 0.001 0.01 0.01 0.01 0.01]; % 通带波动和阻带衰减 f = [f1 f2 f3 f4 f5 f6]/(fs/2); % 频率向量 b = firls(100,f,a,dev); % 滤波器系数 % 滤波 yf = filter(b,1,y); % 滤波后的数据 Yf = fft(yf); % 滤波后的DFT结果 % 计算基波和谐波的频率和幅度比值 f_base = f(abs(f-f(1))<0.001); % 基波频率 f_harm = [2 3 4 5]*f_base; % 谐波频率 A_base = abs(Yf(abs(f-f_base)<0.001)); % 基波幅度 A_harm = abs(Yf(abs(f-f_harm)<0.001)); % 谐波幅度 ratio = A_harm./A_base; % 幅度比值 ``` 注意,由于采样数据是周期锯齿波,因此我们可以对任意一个周期的数据进行处理,无需对整个数据序列进行处理。另外,由于要求误差精度较高,滤波器设计和幅度计算都需要注意细节。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值