maxwell_电机气隙磁密与用matlab进行fft谐波分析(11页)

目录（1）Matlab6.5以上版本软件； ..........................................................错误!未定义书签。

绪论. (1)1 公式分析及计算 (2)1.1傅里叶变换的原理 (2)1.2傅里叶变换的证明 (3)1.3 周期信号的分解 (3)1.4 方波的分解 (5)2 建模与仿真 (7)2.1建模 (7)2.2仿真 (8)3 仿真结果分析 (10)4 小结 (11)参考文献 (13)绪论方波是一种非正弦曲线的波形，通常会于电子和讯号处理时出现。

由于一般电子零件只有“高(1)”和“低(0)”两个值，方波就自然产生，所以理想方波只有“高”和“低”这两个值。

电流的波形为矩形的电流即为方波电流。

不论时间轴上下是不是对称的，只要是矩形就可叫方波，必要时，可加“对称”，“不对称”加以说明。

而在现实世界，方波只有有限的带宽。

因为方波可以快速从一个值转至另一个(即0→1或1→0)，所以方波就用作时钟讯号来准确地触发同步电路。

但是如果用频率定义域来表示方波，就会出然一连串的谐波。

所以方波可用相应频率的基波及其奇次谐波合成。

在电路信号系统的分析中，随着电路规模的加大，微分方程的阶数以及联立后所得的方程的个数也随之加大，加上电器元件的多样化，这些都给解题运算分析电路系统带来了一定的困难。

传统的计算机编程语言，如FORTRAN、C语言等，虽然都可以帮助计算，但在处理高阶微分方程和大规模的联立方程组的问题时大量的时间和精力都花在矩阵处理和图形的生成分析等繁琐易错的细节上。

而MATLAB凭借其强大的矩阵运算能力、简便的绘图功能、可视化的仿真环境以及丰富的算法工具箱，已成为科研和工程技术人员的有力开发工具。

利用MATLAB不仅可以简单快速的求解电路方程，同时，MAYLAB提供的Simulink工具还可以直接建立电路模型，随意改变模型的参数，并且还可以快速得到仿真拟结果，进一步省去了编程的步骤。

电机气隙磁密空间谐波消除方法说实话电机气隙磁密空间谐波消除方法这事，我一开始也是瞎摸索。

我就知道这谐波不消除，电机运行起来那可麻烦了，效率低不说，还容易产生各种问题。

我最先尝试的方法是改变电机的绕组结构。

我想啊，这就好比重新规划一个城市的道路走向似的，如果道路布局不合理，交通就会混乱，电机里的绕组结构不合理，那磁密分布肯定也乱。

我那时候是按照一些传统的书本理论来改变绕组的连接方式，但是效果特别不明显。

我当时就感觉挺挫败的，这理论和实际咋就对不上号呢。

这一失败就让我意识到，光靠理论知识不行，得实际深入探究每个环节的影响。

后来我开始关注电机的铁心材料。

我觉得这铁心就像电机的骨架一样，要是骨架有问题，那整体结构肯定受影响。

我试了好几种不同的铁心材料，想着不同材料的磁导率不一样，会不会对气隙磁密空间谐波有改善呢。

结果发现有一些材料确实有效果，能在一定程度上减少一些谐波，但是没有从根本上消除。

而且有时候呢，一种材料在某个频段的谐波有改善，但在其他频段又不行了，真的是很让人头疼。

然后我又想到了调整气隙长度。

我把这气隙想象成两个物体之间的缝隙，这缝隙大小肯定会影响二者之间的作用力或者某种交流。

当我微调整气隙长度的时候，确实发现磁密空间谐波有了不同的变化。

不过这个调整相当不好把控，因为我发现气隙长度的变化对电机其他性能也有影响，如果只为了消除谐波去改变气隙长度，可能就顾此失彼了。

我还尝试过优化电机的磁极形状呢。

我记得那时候画了好多奇怪的磁极形状，就想找到一种最完美的形状可以让磁密分布均匀，消除谐波。

我拿以前的老电机做实验，把磁极形状切割打磨，重新组装，然后测试磁密。

可是很多形状都不理想，有些形状下电机甚至都不能正常启动了，想想也是挺好玩的，也是自己瞎鼓捣吧。

现在呢，我觉得要比较彻底地消除电机气隙磁密空间谐波，可能需要综合的方法。

像是在合适的绕组结构基础上，选择较好的铁心材料，再搭配一个经过合理微调的气隙长度以及优化的磁极形状。

气隙磁密的fft分解-概述说明以及解释1.引言1.1 概述概述部分：在电机设计中，气隙磁密是一个重要的参数，它影响到电机的性能和效率。

对气隙磁密进行准确的测量和分析对于电机设计和优化至关重要。

传统的方法往往依赖于频域分析技术，而FFT（快速傅里叶变换）作为一种高效的信号处理工具，在气隙磁密的分析中也发挥着重要作用。

本文将探讨气隙磁密的FFT分解方法，介绍其优势和应用领域，并展望未来的研究方向。

通过本文的研究，有望为电机设计提供更准确、高效的分析手段，推动电机技术的发展和应用。

1.2 文章结构:本文主要分为三个部分，分别是引言、正文和结论。

在引言部分，将介绍本文研究的背景和意义，以及文章的结构安排。

在正文部分，将从气隙磁密的概念及重要性、FFT在信号处理中的应用，以及气隙磁密的FFT 分解方法这三个方面展开深入讨论。

最后在结论部分，将总结气隙磁密的FFT分解的优势，展望未来研究方向，并给出结论。

通过这种结构安排，读者能够清晰地了解本文的内容和逻辑发展。

1.3 目的本文旨在研究气隙磁密的FFT分解方法，探讨其在电磁学领域中的应用和意义。

通过对气隙磁密的概念及重要性进行分析，结合FFT在信号处理中的优势，探讨如何将FFT技术应用于气隙磁密的分解过程中。

通过实验验证和理论分析，以期提高气隙磁密分析的效率和精度，为电磁设备的设计和研发提供理论支持和技术参考。

同时，通过本文的研究，展望未来气隙磁密分析领域的发展方向，为相关研究提供思路和启示。

最终，通过总结研究成果和结论，为读者提供对气隙磁密的FFT分解方法有一个全面的认识，为相关领域的研究工作提供参考和借鉴。

2.正文2.1 气隙磁密的概念及重要性气隙磁密指的是在电机等电磁器件中存在的磁场能量在气隙中的分布情况。

在电机中，气隙磁密的大小和分布直接影响到电机的性能和效率。

通常情况下，我们希望气隙磁密能够尽可能均匀地分布在气隙中，以获得更高的磁场强度和更高的转矩输出。

maxwell有限元模型获得气隙磁密波形以Maxwell有限元模型获得气隙磁密波形为标题的文章引言：在电磁学领域中，研究电磁场分布和磁场行为是非常重要的。

而Maxwell有限元模型是一种常用的数值计算方法，可以用来模拟和分析电磁场问题。

本文将通过Maxwell有限元模型来研究气隙磁密波形，探讨其相关特性和应用。

一、Maxwell有限元模型简介Maxwell有限元模型是一种基于Maxwell方程组的数值计算方法，通过将问题区域离散化为有限数量的单元，建立数学模型来近似求解电磁场分布。

这种模型能够更好地描述电磁场在复杂几何结构中的行为，并得到精确的解析结果。

二、气隙磁密波形的研究意义气隙磁密波形是指磁场在存在气隙的介质中的分布形态。

研究气隙磁密波形可以帮助我们了解电磁场在具有气隙的材料中的传播特性，对于电磁设备的设计和优化具有重要意义。

通过Maxwell有限元模型，我们可以得到气隙磁密波形的定量描述，从而为相关研究提供理论依据和实验验证。

三、Maxwell有限元模型在气隙磁密波形研究中的应用1. 气隙磁密波形的建模：首先，我们需要将问题区域进行网格划分，将气隙和周围材料分别建立为不同的单元。

然后，根据Maxwell方程组和边界条件，建立相关的方程和约束条件。

最后，利用数值方法求解这些方程，得到磁场在气隙中的分布情况。

2. 气隙磁密波形的分析：通过Maxwell有限元模型，我们可以得到气隙磁密波形的定量描述，包括磁场的分布形态和强度变化等。

这些数据可以帮助我们分析气隙对电磁场的影响，如磁场集中、场强衰减等。

3. 气隙磁密波形的优化：通过Maxwell有限元模型的模拟和分析，我们可以对气隙磁密波形进行优化设计。

例如，可以调整气隙的位置、形状和尺寸，以改变磁场的分布和强度，从而达到设计要求。

四、气隙磁密波形的应用案例1. 电机设计：在电机中，气隙磁密波形的分布和强度对电机的性能有着重要影响。

通过Maxwell有限元模型，可以预测和优化电机的磁场分布，提高电机的效率和输出功率。

利用matlab仿真对电力系统谐波治理摘要：随着国民经济和科学技术的蓬勃发展，冶金、化学等现代化大工业和电气化铁路的发展，电网负荷加大，电力系统中的非线性负荷(硅整流设备、电解设备、电力机车)及冲击性、波动性负荷(电弧炉、轧钢机、电力机车运行)使得电网发生波形畸变(谐波)、电压波动、闪变、三相不平衡，非对称性(负序)和负荷波动性日趋严重。

电能质量的下降严重地影响了供用电设备的安全、经济运行，降低了人民的生活质量。

所以在世界各国都十分重视电能质量的管理。

引言新兴负荷的出现对电能质量的要求更高电能质量问题逐渐引起普遍重视，主要原因如下：（1）大量基于计算机的控制设备和电子装置投入使用，其性能对电压质量非常敏感。

（2）调速电机和无功补偿装置，导致系统谐波水平不断上升，从而对电力系统的容量和安全运行产生影响。

（3）电力用户不断增长的电能质量意识迫使电力公司提高供电质量，设法解决诸如电压中断，电压跌落和开关暂态等电能质量问题。

衡量电能质量的主要指标是电网频率和电压质量。

频率质量指标为频率允许偏差；电压质量指标包括允许电压偏差、允许波形畸变率(谐波)、三相电压允许不平衡度以及允许电压波动和闪变。

国家技术监督局已公布了上述电能质量的五个国家标准。

电能质量的具体指标。

1．电网频率我国电力系统的标称频率为50Hz，GB／T15945-1995《电能质量一电力系统频率允许偏差》中规定：电力系统正常频率偏差允许值为±0.2Hz，当系统容量较小时，偏差值可放宽到±0.5Hz，标准中没有说明系统容量大小的界限。

在《全国供用电规则》中规定"供电局供电频率的允许偏差：电网容量在300万千瓦及以上者为±0.2HZ；电网容量在300万千瓦以下者，为±0.5HZ。

实际运行中，从全国各大电力系统运行看都保持在不大于±0.1HZ范围内。

2．电压偏差GBl2325－90《电能质量一供电电压允许偏差》中规定：35kV及以上供电电压正负偏差的绝对值之和不超过额定电压的10％；10kV及以下三相供电电压允许偏差为额定电压的土7％；220V单相供电电压允许偏差为额定电压的7％～10％。

matlab谐波分析总结一基本思路为直观分析显示整流装置的谐波特性，使用matlab的simulink搭建整流电路，利用matlab的fft函数分析其电压与电流波形的谐波特性，并利用matlab的绘图工具，直观的显示谐波的相关参数。

输出详细参数到文件。

包括以下想法：1：用simulink搭建一个由多个不同幅值及相位的正弦波，输出到workspace的simout参数，主要是为了验证算法的正确性。

２：算出THD%二算法及验证1：Sine叠加输出sine.mdl文件其中含4个Sine Wave，其参数如下表所示。

Sinewave Amplitude bias Frequency(rad/sec)Phase(rad) SampleTime1 2 0.7 2*pi*50 0 -12 0.5 0 2*pi*50*5 Pi/180*90 -13 1 0 2*pi*50*9 pi/180*45 -14 0.3 0 2*pi*50*26 Pi/180/(-135) -1表达的波形为f(t)=2*sin(2*pi*50*t) +0.5*sin(2*pi*50*5*t+pi/2)+1*sin(2*pi*50*9*t+pi/4) +0.3*sin(2*pi*50*26*t-pi*3/4)为不同幅值与相位的50Hz的基波，5次、9次、26谐波的叠加。

含基波、奇次、偶次、高次谐波。

在基波上加了0.7的偏置，模拟直流分量。

示波器输出到workspace的参数名仿真参数10个周波，每周波采样点2048个使用1/50/2048的采样频率，是为了每个周波采2048个点，便于准确的FFT分析。

理论上可以分析1024次以内的谐波。

simulink的scope的输出simulink的workspace的输出ScopeData.signals.values共10*2048个点。

之所以采10个周波，是为了保证可以避开初始的过渡状态，虽然当前的仿真没有过渡状态，但六脉动整流如果负载有电容的话会有。

基于Matlab的加窗FFT电力系统谐波分析目录摘要： (1)1绪论............................................................................................. 错误!未定义书签。

1.1课题背景、研究意义....................................................... 错误!未定义书签。

1.2 谐波的危害与来源.......................................................... 错误!未定义书签。

1.2.1 谐波来源................................................................... 错误!未定义书签。

1.2.2 电力系统谐波的危害 (3)1.3 谐波检测.......................................................................... 错误!未定义书签。

1.4 谐波的标准与指标.......................................................... 错误!未定义书签。

1.5 国内外关于谐波的研究现状 (5)2谐波分析测量............................................................................. 错误!未定义书签。

2.1 傅里叶级数与系数.......................................................... 错误!未定义书签。

2.2 傅里叶级数的复指数形式.............................................. 错误!未定义书签。

目录摘要 (2)Abstract (2)1：绪论 (2)1.1课题背景 (2)1.2谐波的产生 (3)1.3电网中谐波的危害 (5)1.4研究谐波的重要性 (5)2：谐波的限制标准和常用措施 (7)2.1国外谐波的标准和规定 (8)2.1.1谐波电压标准 (8)2.1.2谐波电流的限制 (9)2.2我国谐波的标准和规定 (9)2.2.1谐波电压标准 (10)2.2.2谐波电流的限制 (11)2.3谐波的限制措施 (12)3：谐波的检测与分析 (15)3.1电力系统谐波检测的基本要求 (15)3.2国内外电力谐波检测与分析方法研究现状 (15)3.3谐波的分析 (18)3.3.1电力系统电压（或电流）的傅立叶分析 (19)3.3.2基于连续信号傅立叶级数的谐波分析 (19)4：电力谐波基于FFT的访真 (21)4.1快速傅立叶变换的简要和计算方法 (21)4.1.1快速傅立叶变换的简要 (21)4.1.2快速傅立叶变换的计算方法 (21)4.2 FFT应用举例 (22)5：结论 (28)附录： (28)参考文献： (30)致谢： (30)基于MATLAB的电力谐波分析学生：指导老师：电气信息工程学院摘要：电力系统的谐波问题早在20世纪20年代就引起人们的注意，到了50年代和60年代，由于高压直流输电技术的发展，发表了有关换流器引起电力系统谐波问题的大量论文。

70年代以来，由于电力电子技术的飞速发展，各种电力电子装置在电力系统、工业、交通及家庭中的应用日益广泛，谐波所造成的危害也日趋严重。

世界各国都对谐波问题予以充分的关注。

本文首先对目前国内外电力谐波检测与分析方法进行了综述与展望，并对电力谐波的基本概念、性质和特征参数进行了详细的分析，给出了谐波抑制的措施。

并得出基于连续信号傅立叶级数的各次谐波系数的计算公式，推导了该计算公式与MATLAB函数FFT计算出的谐波系数的关系。

实例证明：准确测量各次谐波参数，对电力系统谐波分析和抑制具有很大意义，可确保系统安全、可靠、经济地运行。

FFT是离散傅立叶变换的快速算法，可以将一个信号变换到频域。

有些信号在时域上是很难看出什么特征的，但是如果变换到频域之后，就很容易看出特征了。

这就是很多信号分析采用FFT变换的原因。

另外，FFT可以将一个信号的频谱提取出来，这在频谱分析方面也是经常用的。

虽然很多人都知道FFT是什么，可以用来做什么，怎么去做，但是却不知道FFT之后的结果是什意思、如何决定要使用多少点来做FFT。

现在圈圈就根据实际经验来说说FFT结果的具体物理意义。

一个模拟信号，经过ADC 采样之后，就变成了数字信号。

采样定理告诉我们，采样频率要大于信号频率的两倍，这些我就不在此罗嗦了。

采样得到的数字信号，就可以做FFT变换了。

N个采样点，经过FFT之后，就可以得到N个点的FFT结果。

为了方便进行FFT运算，通常N取2的整数次方。

假设采样频率为Fs，信号频率F，采样点数为N。

那么FFT之后结果就是一个为N 点的复数。

每一个点就对应着一个频率点。

这个点的模值，就是该频率值下的幅度特性。

具体跟原始信号的幅度有什么关系呢？假设原始信号的峰值为A，那么FFT的结果的每个点（除了第一个点直流分量之外）的模值就是A的N/2倍。

而第一个点就是直流分量，它的模值就是直流分量的N倍。

而每个点的相位呢，就是在该频率下的信号的相位。

第一个点表示直流分量（即0Hz），而最后一个点N的再下一个点（实际上这个点是不存在的，这里是假设的第N+1个点，也可以看做是将第一个点分做两半分，另一半移到最后）则表示采样频率Fs，这中间被N-1个点平均分成N等份，每个点的频率依次增加。

例如某点n所表示的频率为：Fn=(n-1)*Fs/N。

由上面的公式可以看出，Fn所能分辨到频率为为Fs/N，如果采样频率Fs 为1024Hz，采样点数为1024点，则可以分辨到1Hz。

1024Hz的采样率采样1024点，刚好是1秒，也就是说，采样1秒时间的信号并做FFT，则结果可以分析到1Hz，如果采样2秒时间的信号并做FFT，则结果可以分析到0.5Hz。

气隙磁密的fft分解全文共四篇示例，供读者参考第一篇示例：气隙磁密是电机设计与分析中一个重要的参数，它可以直接影响电机的性能和效率。

在电机中，气隙磁密是指在电机磁场中气隙处的磁感应强度，通常用磁感应强度B表示。

气隙磁密的大小取决于电机的磁路设计、磁场分布和材料特性等因素。

传统的方法中，气隙磁密通常通过有限元分析等数值方法来进行计算和分析。

有限元分析通常需要耗费大量的计算资源和时间。

为了提高计算效率和精度，研究人员提出了使用快速傅里叶变换（FFT）分解技术对气隙磁密进行分析的方法。

FFT是一种将信号从时间域转换到频率域的技术，它可以将复杂的信号分解成若干个简单的正弦和余弦波形。

在电机中，气隙磁密可以看作是一个由不同频率的正弦和余弦波形组成的信号，因此可以利用FFT技术将气隙磁密信号进行频谱分解，得到不同频率成分的磁密波形，从而更好地分析电机的磁场特性。

通过FFT分解气隙磁密信号，可以获得以下几方面的优势：1. 提高计算效率：与传统的有限元分析相比，FFT技术可以在较短的时间内完成气隙磁密的分析，提高计算效率，减少设计周期。

2. 提高精度：FFT技术可以将磁密信号分解为不同频率成分，更好地揭示气隙中的磁场变化规律，提高分析精度。

3. 提供更多信息：通过分解不同频率的磁密波形，可以获得气隙中不同频率成分的磁场分布情况，为电机设计和性能优化提供更多信息。

在实际应用中，研究人员可以将FFT技术与有限元分析等方法结合，综合利用它们的优势，更好地分析电机的气隙磁密情况。

还可以通过对气隙磁密信号的频谱分析，进一步优化电机的磁路设计，提高电机的性能和效率。

利用FFT技术对气隙磁密进行分解分析，可以提高计算效率、精度和信息量，为电机设计与性能优化提供更多有益的参考。

随着计算机技术的不断进步和应用软件的完善，预计这一技术在电机行业的应用将越来越广泛，为电机的研发和应用带来更多好处。

第二篇示例：气隙磁密是电力系统中重要的参数之一，用于描述磁场在气隙中的分布情况。

clcclear all;format long;Ns=1000;order=13; I与**********************read the position and flux density************************ fid=fopen(,",' r'); %open the original file fidnew = fopen(* ',' w'); Swrite the new file while feof(fid)==0tline = fgetl(fid); %tline if ^ischar(tline), break, end temp=abs(tline);Nlength=length(tline); isemptyline=O; % (if Nlength==0 isemptyline=l;endallspace=O; % .isspace=0;for i=l:Nlength T=temp(i);if T==32isspace=isspace+l;endif isspace==Nlength allspace=l;breakend<endfindalpha=O; % for j=l:Nlength T=temp(j);if ((T>=65)&(T>=90)) | ((T>=97)&(T>=122)) findalpha=l;break; end end )if (^f indalpha) & (""al 1 space) & (i sempt5r l ine==O) % fprintf(fidnew, tline); fprintf(fidnew, ' \n'); endendfclose(fid); fclose(fidnew);fidl=fopen(* ',' r');■flux_position =fscanf(fidl, > [2, Ns]); fclose(fidl);%********************************read file finish***************************************** flux_position=flux_position,;posl=flux_posit ion(:, 1); (pos_delta=posl(2);pos_length=length(posl);pos_last=posl(pos_length);for i=l:l:pos_length %copy and get another part of position pos2(i)=pos_last+i*pos_delta; ( end posl=posl';fluxl=flux_position(:, 2);flux2=-flux_pos it ion(:, 2);pos=[posl, pos2];%combine and get all part of position > fluxl=fluxf ;flux2=flux2,;flux=[fluxl, flux2];%combine and get all part of flux density value figure;plot(posl, fluxl,' r') ;%plot origional waveform hold on;grid on;fftl=fft(flux, Ns);j=0； amp_har=zeros(l, (order+l)/2);for m=l:2:orderj=j+l;fftl=fft(flux, Ns);fund_ele_front=fftl(m+1);1fund_ele_back=fftl(Ns+l-m);amp_har(j)=(abs(fund_ele_front))/Ns*2;fftl=O*fftl;fftl(m+l)=fund_ele_front;fftl(Ns+l-m)=fund_ele_back;fftl=ifft(fftl,Ns);fftl=real(fftl);plot(posl, fftl);hold on;endYk=(l:2:order);figure;bar(k, amp_har);grid on;peak_b=max(fftl)srms_b=*peak_b clcclear all;format long;Ns=1000;order=7;$**********************read the position and flux density************************ fid=fopenC '," r'); %open the original file fidnew = fopenC ',' w'); %write the new file while feof(fid) =0 Atline = fgetl(fid); Stline if "ischar(tline), break, end temp=abs(tline);Nlength=length(tline); isemptyline=0; % ] if Nlength==0 isemptyline=l; end allspace=0; % * isspace=0;for i=l:Nlength T=temp(i);if T==32 isspace=isspace+l; I end if isspace==Nlength allspace=l;break endendfindalpha=0; % for j=l:Nlength T=temp(j);if ((T>=65)&(T>=90))!((T>=97)&(T>=122)) findalpha=l; break;end endif (^f indalpha) & (""al 1 space) & (i sempt5r l ine==O) %fprintf(fidnew, tline);fprintf(fidnew, ' \n');end(endfclose(fid);fclose(fidnew);fidl=fopen(' ',' r');flux_position =fscanf(fidl, ' %f', [2, Ns]);fclose(fidl);$********************************read filefinish***************************************** flux_position=flux_position,;posl=flux__position(:, 1);)pos_delta=posl (2);pos_length=length(posl);pos_last=posl(pos_length);for i=l:l:pos_length %copy and get another part of positionpos2(i)=pos_last+i*pos_delta;endpos^posl* ;fluxl=flux_position(:, 2);flux2=-flux_position(:, 2);pos=[posl, pos2];%combine and get all part of positionfluxl=fluxf ;flux2=flux2,;flux=[fluxl, flux2];%combine and get all part of flux density value figure; plot(posl, fluxl, ' r');%plot origional waveformhold on;grid on;fftl=fft (flux, Ns);j=0;amp_har=zeros(l, (order+l)/2);for m=l:2:orderfftl=fft(flux, Ns);fund_ele_front=fftl(m+1);f und_e1e_bac k=ff11(Xs+l-m);amp_har(j)=(abs(fund_ele.front))/Ns*2;fftl=O*fftl;fftl(m+l)=fund_e1e_front;fftl(Ns+l-m)=fund_ele_back;fftl=ifft(fftl, Ns);fftl=real(fftl);plot (posl, fftl);hold on;end k=(l:2:order);figure;bar(k, amp_har);grid on;%peak_b=max(fftl) %rms_b=*peak_b。

基于Matlab的电网谐波电流实时检测方法
宋婷;刘崇新
【期刊名称】《胜利油田职工大学学报》
【年(卷),期】2007(0)5
【摘要】随着电力电子等非线性用电设备的大量应用,电网当中的谐波和无功问题日益严重.谐波和无功的存在不仅降低了电力系统的供电效率,而且严重影响了电能质量,危害公共电网的安全运行,干扰临近用电设备的正常运行.讨论了基于瞬时无功功率理论的谐波电流检测方法,并应用matlab仿真平台进行了仿真研究,给出了仿真波形,验证了该方法的正确性.
【总页数】2页(P73-74)
【作者】宋婷;刘崇新
【作者单位】胜利油田职工大学,山东,东营,257000;西安交通大学电气学院,陕西,西安,710049
【正文语种】中文
【中图分类】TM761+.12
【相关文献】
1.基于瞬时无功功率理论的改进型煤矿电网谐波电流检测方法
2.一种基于瞬时无功功率理论的电网谐波电流检测方法
3.基于瞬时无功功率理论的电网谐波电流检测方法仿真研究
4.基于高通滤波器的电网谐波电流检测方法研究
5.基于MATLAB的谐波电流检测方法的建模与仿真
因版权原因，仅展示原文概要，查看原文内容请购买。

clcclear all;format long;Ns=1000;order=13;%**********************read the position and fluxdensity************************fid=fopen('B.dat','r'); %open the original filefidnew = fopen('b1.dat','w'); %write the new filewhile feof(fid)==0tline = fgetl(fid); %tline?if ~ischar(tline), break, endtemp=abs(tline);Nlength=length(tline);isemptyline=0; %if Nlength==0isemptyline=1;endallspace=0; %isspace=0;for i=1:NlengthT=temp(i);if T==32isspace=isspace+1;endif isspace==Nlengthallspace=1;breakendendfindalpha=0; %for j=1:NlengthT=temp(j);if ((T>=65)&(T>=90))|((T>=97)&(T>=122))findalpha=1;break;endendif(~findalpha)&(~allspace)&(isemptyline==0) % fprintf(fidnew,tline);fprintf(fidnew,'\n');endendfclose(fid);fclose(fidnew);fid1=fopen('b1.dat','r');flux_position =fscanf(fid1,'%f',[2,Ns]);fclose(fid1);%********************************read filefinish*****************************************flux_position=flux_position';pos1=flux_position(:,1);pos_delta=pos1(2);pos_length=length(pos1);pos_last=pos1(pos_length);for i=1:1:pos_length %copy and get another part of positionpos2(i)=pos_last+i*pos_delta;endpos1=pos1';flux1=flux_position(:,2);flux2=-flux_position(:,2);pos=[pos1,pos2];%combine and get all part of positionflux1=flux1';flux2=flux2';flux=[flux1,flux2];%combine and get all part of flux density value figure;plot(pos1,flux1,'r');%plot origional waveformhold on;grid on;fft1=fft(flux,Ns);j=0;amp_har=zeros(1,(order+1)/2);for m=1:2:orderj=j+1;fft1=fft(flux,Ns);fund_ele_front=fft1(m+1);fund_ele_back=fft1(Ns+1-m);amp_har(j)=(abs(fund_ele_front))/Ns*2;fft1=0*fft1;fft1(m+1)=fund_ele_front;fft1(Ns+1-m)=fund_ele_back;fft1=ifft(fft1,Ns);fft1=real(fft1);plot(pos1,fft1);hold on;endk=(1:2:order);figure;bar(k,amp_har);grid on;peak_b=max(fft1)rms_b=0.707*peak_bclcclear all;format long;Ns=1000;order=7;%**********************read the position and flux density************************fid=fopen('B.dat','r'); %open the original file fidnew = fopen('b1.dat','w'); %write the new file while feof(fid)==0tline = fgetl(fid); %tline?if ~ischar(tline), break, endtemp=abs(tline);Nlength=length(tline);isemptyline=0; %if Nlength==0isemptyline=1;endallspace=0; %isspace=0;for i=1:NlengthT=temp(i);if T==32isspace=isspace+1;endif isspace==Nlengthallspace=1;breakendendfindalpha=0; %for j=1:NlengthT=temp(j);if ((T>=65)&(T>=90))|((T>=97)&(T>=122))findalpha=1;break;endendif(~findalpha)&(~allspace)&(isemptyline==0) % fprintf(fidnew,tline);fprintf(fidnew,'\n');endendfclose(fid);fclose(fidnew);fid1=fopen('b1.dat','r');flux_position =fscanf(fid1,'%f',[2,Ns]);fclose(fid1);%********************************read filefinish*****************************************flux_position=flux_position';pos1=flux_position(:,1);pos_delta=pos1(2);pos_length=length(pos1);pos_last=pos1(pos_length);for i=1:1:pos_length %copy and get another part of positionpos2(i)=pos_last+i*pos_delta;endpos1=pos1';flux1=flux_position(:,2);flux2=-flux_position(:,2);pos=[pos1,pos2];%combine and get all part of positionflux1=flux1';flux2=flux2';flux=[flux1,flux2];%combine and get all part of flux density value figure;plot(pos1,flux1,'r');%plot origional waveformhold on;grid on;fft1=fft(flux,Ns);j=0;amp_har=zeros(1,(order+1)/2);for m=1:2:orderj=j+1;fft1=fft(flux,Ns);fund_ele_front=fft1(m+1);fund_ele_back=fft1(Ns+1-m);amp_har(j)=(abs(fund_ele_front))/Ns*2;fft1=0*fft1;fft1(m+1)=fund_ele_front;fft1(Ns+1-m)=fund_ele_back;fft1=ifft(fft1,Ns);fft1=real(fft1);plot(pos1,fft1);hold on;endk=(1:2:order);figure;bar(k,amp_har);grid on;%peak_b=max(fft1)%rms_b=0.707*peak_b。