资源描述：
《使用matlab的fft分析六脉动整流的谐波》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、matlab谐波分析总结一基本思路为直观分析显示整流装置的谐波特性，使用matlab的simulink搭建整流电路，利用matlab的fft函数分析其电压与电流波形的谐波特性，并利用matlab的绘图工具，直观的显示谐波的相关参数。输出详细参数到文件。包括以下想法：1：用simulink搭建一个由多个不同幅值及相位的正弦波，输出到workspace的simout参数，主要是为了验证算法的正确性。２：算出THD%二算法及验证1：Sine叠加输出sine.mdl文件其中含4个SineWave，其参数如下表所示。locatedintheTomb,DongShenJiabang,deferthe
2、nextdayfocusedontheassassination.Linping,Zhejiang,1ofwhichliquorwinemasters(WuzhensaidinformationisCarpenter),whogotAfewbayonets,duetomissedfatal,whennightcameSinewaveAmplitudebiasFrequency(rad/sec)Phase(rad)SampleTime120.72*pi*500-120.502*pi*50*5Pi/180*90-13102*pi*50*9pi/180*45-140.302*pi*50*26P
3、i/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次以内的谐波。locatedint
4、heTomb,DongShenJiabang,deferthenextdayfocusedontheassassination.Linping,Zhejiang,1ofwhichliquorwinemasters(WuzhensaidinformationisCarpenter),whogotAfewbayonets,duetomissedfatal,whennightcamesimulink的scope的输出simulink的workspace的输出ScopeData.signals.values共10*2048个点。之所以采10个周波，是为了保证可以避开初始的过渡状态，虽然当前的仿真
5、没有过渡状态，但六脉动整流如果负载有电容的话会有。实际进行FFT分析时，只抽取一个稳定状态周期的2048个点。2：代码程序名分fenxi_FFT.m%基波频率为50HzN=2048*4;%每周波采2048个点fp=30;%拟显示的倍频次数，需小于N/2用于使前面的参数显示清晰wq=5;%稳定周期个数，为保证下行程序数组不出错，workspce的数据周期必须大于wq值kh=0.01;%幅值比例系数，输出到文本时太小的谐波参数不输出fid=fopen('output.txt','wt');%准备写入文件，文件在当前目录下fprintf(fid,'每周波采样点数=%d',N);y=Scop
6、eData.signals.values;%读取workspce数据,第一个点前面为示波器输出到workspce参数名t=0:0.02/N:0.02/N*(size(y)-1);%每个周期为0.02秒，每周期采Ｎ个点，共size(y)-1个点locatedintheTomb,DongShenJiabang,deferthenextdayfocusedontheassassination.Linping,Zhejiang,1ofwhichliquorwinemasters(WuzhensaidinformationisCarpenter),whogotAfewbayonets,duetom
7、issedfatal,whennightcamesubplot(2,1,1);plot(t,y);xlabel('时间/s');ylabel('幅值');title('原始信号');gridon;subplot(2,1,2);pn=wq*N;%取数据的基准点，主要目的是等系统稳定plot(t(pn:pn+N-1),y(pn:pn+N-1));%数据周期必须小于wq，将出错xlabel('时间/s');ylabel('幅值');tit
