谐波分析高精度FFT

基于matlab谐波抑制的仿真研究(毕设)

电力系统谐波抑制的仿真研究 目 录 1 绪论…………………………………………………………………………… 1.1 课题背景及目的………………………………………………………… 1.2国内外研究现状和进展………………………………………………… 1.2.1国外研究现状 …………………………………………………… 1.2.1国内研究现状 …………………………………………………… 1.3 本文的主要内容…………………………………………………………… 2 有源电力滤波器及其谐波源研究……………………………………………… 2.1 谐波的基本概念………………………………………………………… 2.1.1 谐波的定义……………………………………………………… 2.1.2谐波的数学表达………………………………………………… 2.1.3电力系统谐波标准………………………………………………… 2.2 谐波的产生……………………………………………………………… 2.3 谐波的危害和影响……………………………………………………… 2.4 谐波的基本防治方法…………………………………………………… 2.5无源电力滤波器简述…………………………………………………… 2.6 有源电力滤波器介绍…………………………………………………… 2.6.1 有源滤波器的基本原理.……………………………………… 2.6.2 有源电力滤波器的分类.……………………………… 2.7并联型有源电力滤波器的补偿特性…………………………………… 2.7.1谐波源………………………………………………………… 2.7.2有源电力滤波器补偿特性的基本要求…………………………… 2.7.3影响有源电力滤波器补偿特性的因素…………………………… 2.7.4并联型有源电力滤波器补偿特性……………………………… 2.8 谐波源的数学模型的研究……………………………………………… 2.8.1 单相桥式整流电路非线性负荷………………………………… 2.8.2 三相桥式整流电路非线性负荷.………………………………… 3 基于瞬时无功功率的谐波检测方法…………………………………………… 3.1谐波检测的几种方法比较…………………………………………… 3.2三相电路瞬时无功功率理论…………………………………………… 3.2.1瞬时有功功率和瞬时无功功率……………………………………… 3.2.2瞬时有功电流和瞬时无功电流……………………………………… 3.3 基于瞬时无功功率理论的p q -谐波检测算法.…………………… 3.4基于瞬时无功功率理论的p q i i -谐波检测法.…………………… 4并联有源电力滤波器的控制策略…………………………………………… 4.1并联型有源电力滤波器系统构成及其工作原理………………………… 4.2并联有源电力滤波器的控制研究.……………………………… 4.2.1并联有源电力滤波器直流侧电压控制…………………… 4.2.2有源电力滤波器电流跟踪控制技术…………………………… 4.2.2.1 PWM 控制原理………………………………………… 4.2.2.2滞环比较控制方式………………………………………… 4.2.2.3三角波比较方式………………………… 4.3有源电力滤波器的主电路设计 …………………………………………

基于MATLAB的电力谐波分析

目录 摘要 (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计算出的谐波系数的关系。实例证明：准确测量各次谐波参数，对电力系统谐波分析和抑制具有很大意义，可确保系统安全、可靠、经济地运行。同时实验结果表明，该法对设备要求不高，易于实现。 关键字：MA TLAB电力谐波分析 Harmonic Analysis of Electric Power System Based On Matlab Student: Teacher: Electrical and Information Engineering Abstract:The harmonic problem of electric power system has caused the attention of people in1920s and 1930s.Until 1950s,owing to the development of high voltage direct current transportation electricity technology,people published a large number of theses about the electricity power system harmonic problem,which caused by the current transform device.Since 1970s,because of the speedly development of eletricity power electronics technology,the various electric power electronics devices were applied extensively in the electric power system,industry,traffic and family,but the harm which the harmonic creates was serious more and more.Many country of the world all pay attention to the harmonic problem. Summary and Prospects of the first domestic and international power harmonics detection and analysis methods, and power harmonics of the basic concepts of the nature and characteristic parameters of a detailed analysis, given a harmonic suppression measures. Obtained based on the

matlab信号仿真谐波

综合训练① 实验内容：利用matlab绘制频率自定的正弦信号（连续时间和离散时间），复指数信号（连续时间），并举例实际中哪些物理现象可以用正弦信号，复指数信号来表示。绘制成谐波关系的正弦信号（连续时间和离散时间），分析其周期性和频率之间的关系。实验步骤： 一、绘制谐波关系的正弦信号 分析：由于正弦信号可以表示成两个共轭的复指数信号相减，然后再除去两倍的单位虚数得到，故，我们将正弦信号设置为 X=exp(j*pi*n/4)-exp(-j*pi*n/4))/(2*j) 此信号就相当于 x=sin(pi*n/4) 设计程序如下： n=[0:32]; %设置n的取值 x=(exp(j*pi*n/4)-exp(-j*pi*n/4))/(2*j); %限定离散正弦信号 stem(n,x) %绘制该离散正弦信号 通过Matlab所得图形如下：

分析：同样的连续型的正弦信号同样也可以用类似方式绘制. x=sym('(exp(j*pi*t/T)+exp(-j*pi*t/T))/2');%函数表示正弦信号 x5=subs(x,5,'T'); %设置周期大小ezplot(x5,[0,10]) %绘制图形 所得结果如下：

二、绘制复指数信号 分析：由于复指数信号有实数部分和虚数部分，所以绘制其图形，我们采取了分别绘制的方法，将实数和虚数分别画出。 实验程序如下： t=[0:.01:10]; %产生时间轴的等差点 y=exp((1+j*10)*t); %设置复指数信号 subplot(211),plot(t,real(y)); %绘制实数信号图形 grid subplot(212),plot(t,imag(y)); %绘制虚数部分图形 grid 实验所得结果如下：

用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之后某点n用复数a+bi表示，那么这个复数的模就是An=根号a*a+b*b，相位就是Pn=atan2(b,a)。根据以上的结果，就可以计算出n点（n≠1，且n<=N/2）对应的信号的表达式为：An/(N/2)*cos(2*pi*Fn*t+Pn)，即2*An/N*cos(2*pi*Fn*t+Pn)。对于n=1点的信号，是直流分量，幅度即为A1/N。 由于FFT结果的对称性，通常我们只使用前半部分的结果，即小于采样频率一半的结果。 好了，说了半天，看着公式也晕，下面以一个实际的信号来做说明。 假设我们有一个信号，它含有2V的直流分量，频率为50Hz、相位为-30度、幅度为3V 的交流信号，以及一个频率为75Hz、相位为90度、幅度为1.5V的交流信号。用数学表达式就是如下： S=2+3*cos(2*pi*50*t-pi*30/180)+1.5*cos(2*pi*75*t+pi*90/180) 式中cos参数为弧度，所以-30度和90度要分别换算成弧度。我们以256Hz的采样率对这个信号进行采样，总共采样256点。按照我们上面的分析，Fn=(n-1)*Fs/N，我们可以知道，每两个点之间的间距就是1Hz，第n个点的频率就是n-1。我们的信号有3个频率：0Hz、50Hz、75Hz，应该分别在第1个点、第51个点、第76个点上出现峰值，其它各点应该接近0。实际情况如何呢？我们来看看FFT的结果的模值如图所示。

基于MATLAB的SVPWM高次谐波分析

基于MATLAB的SVPWM高次谐波分析 彭润泽14721429 （上海大学机电工程与自动化学院，上海200072） 摘要：从空间矢量脉宽（SVPWM）调制产生的机理出发,推导了空间矢量电压的公式，简介了SVPWM的实现方式，指出了五段式和七段式两种主要逼近方法下。在MA TLAB/Simulink的软件环境下对两种合成方式搭建了仿真模型，分析了两种方式产生的SVPWM谱的谐波分布和相应逆变器产生电压的谐波分布情况,并总结出了两者的优缺点。得到了七段SVPWM的正弦度高和五段SVPWM开关损耗少的结论。 关键词：空间矢量调制技术；高次谐波分析；傅里叶变换；MA TLAB/Simulink Study on harmonics of SVPWM based on MATLAB Peng Runze (School of Engineering-Mechatronics and automation, Shanghai University, Shanghai 200072, China) Abstract: Starting off with the principle of space vector pulse width modulation(SVPWM), this paper deduces formulas of the voltage space and points out two main approximation methods. In software MATLAB / Simulink , simulation models are built to analyze the harmonics generated by the distribution of two ways SVPWM harmonic spectrum distribution and the corresponding voltage generated by the inverter. Finally, the paper sums up the advantages and disadvantages of two methods and hence draw conclusions that seven-step SVPWM possesses better sine waveform while five-step one reduces switching loss. Key words: SVPWM; Fourier transform; harmonic; MATLAB/Simulink 0 引言 在全数字高性能交流调速系统中,通常采用数字脉宽调制来代替传统的模拟脉宽调制。空间矢量脉宽调制(SVPWM)技术是近年来研究的一个热点。SVPWM是把三相变流器的端部电压状态在复平面上综合为空间电压矢量,并通过不同的开关状态形成八个空间矢量,利用这八个空间矢量来逼近电压圆,从而形成SVPWM波, 它能在较低的开关频率下获得较好的谐波抑制效果和比SPWM高出15%的基波电压[1-3],同时,SVPWM调制技术还有一个优点即易于实现数字和实时控制[4]，其主电路如图1所示。空间矢量脉宽调制(SVPWM)因其与同步旋转坐标系控制量接口方便，直流电压利用率高、总谐波畸变小等特点被广泛应用于PWM 整流器当中[5-6]。电压型PWM逆变器把部分基波能量转换成谐波能量,造成交流输入电压、电流发生畸变,向系统中注入高次谐波,使输入功率因数降低,电能质量下降,对电力系统包括用户的安全、经济运行产生严重的危害和影响,甚至造成电力设备的损坏[7]；]对于三相电机来说，高次谐波会造成转矩脉动，影响电机控制精度、提高控制难度，还会增加损耗，造成电机发热。所以研究空间矢量脉宽调制产生的高次谐波具有重要的意义。 图1 三相电压型逆变器原理图 Fig.1 schematic diagram of the inverter 1 SVPWM工作原理 空间矢量脉宽调制SVPWM(Space Vector Pulse Width Modulation)，实际上是对应于交流感应电机或永磁同步电机中的三

基于matlab的电力系统谐波抑制仿真研究--任务书

毕业设计（论文）课题任务书 （2013----2014学年） 学院名称：科技学院 课题名称基于matlab的电力系统谐波分析 学生姓名专业电气工程及自动化学号 指导教师任务书下达时间2013年12月16 日 课题概述： 通过本次毕业设计，使学生对电力系统谐波抑制的问题有比较全面的了解。并且学习电力系统谐波对电力系统电能质量方面及继电保护方面的影响,重点研究电力系统谐波的产生机理并利用MATLAB编程进行仿真，得到滤波后较好的电能质量结果。 要求阅读或检索的参考资料及文献（包括指定给学生阅读的外文资料）： 1 唐统一，吴震春，孙树勤.电力系统谐波[M].徐州：中国矿业大学出版社，1991. 2 曲学基，曲敬铠，于明扬等.电力电子滤波技术及其应用[M].北京：电子工业出版社，2008. 3 浅析配电网谐波的危害与治理[J]．科苑观查，32. 4 George J.Wakileh,徐政(译)．电力系统谐波——基本原理、分析方法和滤波器设计[M]．机械工 业出版社，2007． 5 陈志业,尹华丽．电能质量及其治理新技术[J]．电网技术，2002, 26(7)：67-70． 6、电能质量公用电网谐波,GB/T14549-93,中国电力出版社. 7、翁利民，张启春，舒立平等．配电网高次谐波抑制方法的分析[J]．电力自动化设备，2000， 20(3):22-25． 8、Luo Weiping. The Harmonious Detecting and Analyzing of PV Grid Connected Generation. Proceedings of the 27th Chinese Control Conference July 16-18, 2008, Kunming,Yunnan, China． 9、翁利民．电力电子技术与谐波抑制，无功补偿技术研究综述[J]．电力电容器，2004,3:6-10. 10、李承,皱云屏. 有源电力滤波器抑制谐波的机理分析[J]．电力系统自动化，2003,27(20):31-34. 11、李圣清，朱英浩，周有庆等．电网谐波检测方法的综述[J]．高电压技术，2004，30(3)：39-42． 12、王兆安，黄俊．电力电子技术(第4版)EM3．机械工业出版社，2005.07． 13、刘桂英,粟时平. 并联有源电力滤波器滞环电流比较控制研究[J]. 高电压技术, 2005, 3l(5): 58-60． 14、赵金宝，非线性电力负荷的谐波与抑制技术的研究．广西大学硕士学位论文，2006． 15、王兆安，杨君等．谐波抑制和无功功率补偿(第二版)[M]．机械工业出版社，2006.01．

FFT谐波频谱相谱分析Matlab程序二则

FFT谐波频谱相谱分析Matlab程序二则 by ggihhimm 程序一：从Excel表中读取数据，然后分析其频谱、相谱和总谐波畸变率，并画出茎状图（stem），图示如下（测试数据含直流分量）： %该程序从Excel表中读取数据，然后分析其频谱、相谱和总谐波畸变率 %一组数据+采样率显示频谱下面的阀值应适当调大注意I1rms的定义 clear all; close all; clc; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%% y=xlsread('HCwithDI_2.xls'); %一列数据：excel表格的第一列要知道数据的采样率，为 %保证结果的准确性，最好为整数倍周期 L=length(y); mdflag=1; %是否标点,1是，0否 del=1e-1; %删除的阀值，根据数据幅值大小来设置 Fs=50*L; %采样率 f0=50; %信号基波频率%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%% f=0:L-1;

f=f*Fs/L/f0; %输出图像横轴为谐波次数 yfft=fft(y); yfft(L/2+1:L)=[]; f(L/2+1:L)=[]; mag=abs(yfft)/L*2; %各次谐波幅值 ang=angle(yfft)*180/pi; mag(1)=mag(1)*0.5; %直流分量幅值 mag1=mag; mag1(mag<=del)=[]; %去掉偶数次谐波 ang1=ang; ang1(mag<=del)=[]; %去掉相应相角 f1=f; f1(mag<=del)=[]; %去掉相应序号 L1=length(f1); subplot(211); hdl1=stem(f1(1:L1),mag1(1:L1),'b:o'); set(hdl1,'MarkerFaceColor','red','Marker','square'); ylabel('各次谐波幅值/A'); tem=get(gca,'XLim'); xwidth=tem(end)-tem(1); x=zeros(1,L1); zhi=zeros(2,L1); dotstr=cell(2,L1); xright=1.5; %右调因子x(1)=f1(1)+xwidth/100*xright; %调整位置 if f1(L1)

matlab谐波分析程序

clc clear all; format long; Ns=1000; order=13; ! %**********************read the position and flux density************************ fid=fopen('','r'); %open the original file fidnew = fopen('','w'); %write the new file while feof(fid)==0 > tline = fgetl(fid); %tline if ~ischar(tline), break, end temp=abs(tline); Nlength=length(tline); isemptyline=0; % { if Nlength==0 isemptyline=1; end allspace=0; % 、 isspace=0; for i=1:Nlength T=temp(i); if T==32 isspace=isspace+1; % end if isspace==Nlength allspace=1; break end < end findalpha=0; % for j=1:Nlength T=temp(j); ! if ((T>=65)&(T>=90))|((T>=97)&(T>=122)) findalpha=1;

break; end end ） if (~findalpha)&(~allspace)&(isemptyline==0) % fprintf(fidnew,tline); fprintf(fidnew,'\n'); end - end fclose(fid); fclose(fidnew); fid1=fopen('','r'); · flux_position =fscanf(fid1,'%f',[2,Ns]); fclose(fid1); %********************************read file finish***************************************** 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 position pos2(i)=pos_last+i*pos_delta; （ end pos1=pos1'; flux1=flux_position(:,2); flux2=-flux_position(:,2); pos=[pos1,pos2];%combine and get all part of position > flux1=flux1'; flux2=flux2'; flux=[flux1,flux2];%combine and get all part of flux density value figure; plot(pos1,flux1,'r');%plot origional waveform " hold on; grid on; fft1=fft(flux,Ns);

基于MATLAB的谐波分析FFT

目录 （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工具还可以直接建立电路模型，随意改变模型的参数，并且还可以快速得到仿真拟结果，进一步省去了编程的步骤。MATLAB具有数值计算功能；图形处理及可视化功能；可视化建模及动态仿真功能等等。它给用户带来的是最直观，最简洁的程序开发环境。它的语言简洁紧凑，使用方便灵活，程序书写形式自由，利用起丰富的库函数避开繁杂的子程序编程任务，压缩了一切不必要的编程工作。同时，它的运算符也很丰富。由于MATLAB是用C 语言编写的，MATLAB提供了和C语言几乎一样多的运算符，灵活使用MATLAB的运算符将使程序变得极为简短。它的程序的可移植性，基本上不做修改就可以在各种型号的计算机和操作系统上运行。 本文应用MATLAB来验证定理：方波可用相应频率的基波及其奇次谐波合成。

使用Matlab的FFT分析六脉动整流的谐波

matlab谐波分析总结 一基本思路 为直观分析显示整流装置的谐波特性，使用matlab的simulink搭建整流电路，利用matlab的fft函数分析其电压与电流波形的谐波特性，并利用matlab的绘图工具，直观的显示谐波的相关参数。输出详细参数到文件。 包括以下想法： 1：用simulink搭建一个由多个不同幅值及相位的正弦波，输出到workspace的simout参数，主要是为了验证算法的正确性。 ２：算出THD% 二算法及验证 1：Sine叠加输出 sine.mdl文件 其中含4个Sine Wave，其参数如下表所示。

Sinewave Amplitud e bias Frequency(rad/sec ) Phase(rad) Sample Time 1 2 0.7 2*pi*50 0 -1 2 0.5 0 2*pi*50*5 Pi/180*90 -1 3 1 0 2*pi*50*9 pi/180*45 -1 4 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次以内的谐波。

maxwell_电机气隙磁密与用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之后某点n用复数a+bi表示，那么这个复数的模就是An=根号a*a+b*b，相位就是Pn=atan2(b,a)。根据以上的结果， 就可以计算出n点（n≠1，且n<=N/2）对应的信号的表达式为：An/(N/2)*cos(2*pi*Fn*t+Pn)，即2*An/N*cos(2*pi*Fn*t+Pn)。对于n=1点的信号，是直流分量，幅度即为A1/N。 由于FFT结果的对称性，通常我们只使用前半部分的结果，即小于采样频率一半的结果。