实用MATLABfft程序

实用Matlab快速傅里叶变换程序程序，很有效，直接使用

%  先从示波器数据中取样，因示波器采样点数太多，每3个数据取1个
clear;

% load PHA_i.tmp;
% N=length（PHA_i）;
% x1=zeros（1N）;
% x1=PHA_i（:2）‘;
% x=x1（200:800）;

% 
% load ia1.txt; 
% b=ia1（250:3:55002）;     %取整周期

% load ia2.txt; 
% b=ia2（625:3:56752）;     %取整周期

% 
load u1.txt; 
b=u1（222:3:53442）;     %取整周期

% load emf.txt; 
% b=emf（730:1:61502）;     %取整周期


% 
% load u2.txt; 
% b=u2（450:3:58102）;     %取整周期

% x=v0（2:100001）;

M=length（b）;            %M为2pi区间内的采样点数
M_end=50*M;
 for k=0:M_end/M-1;
     for i=1:M;
        x（k*M+i）=b（i）;     %原始信号
     end
 end
%----------------------------------------------  
n=length（x）; %N=1024; %总采样点
Ts=9.6e-5*2/3;   %Ts=tp/N; %采样周期
tp=n*Ts;     %tp=0.2; %采样总时间
detf=1/tp;   %detf=1/tp; %频域分辨率
t=0:Ts:Ts*（n-1）;      %t=0:Ts:Ts*（N-1）; %时域采样序列
f=0:detf:（n-1）*detf;  %f=0:detf:（N-1）*detf; %频域采样序列
fftx=fft（x）/（n/2）;    %FFT运算
y=abs（fftx）;

figure（2）
axes（‘FontSize‘[18]）;
plot（fy‘LineWidth‘2）;        %半对数图显示
ylabel（‘Vol line to line（V）‘‘fontsize‘[25]）;
% ylabel（‘ia（A）‘‘fontsize‘[25]）
xlabel（‘harmonic（Hz）‘‘fontsize‘[25]）;
grid;
% print -djpeg ‘uasSpetrum.jpeg‘;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1245  2019-01-17 16:35  实用MATLABfft程序\Ifft_fromScope.m

     目录          0  2019-02-25 14:52  实用MATLABfft程序

----------- ---------  ---------- -----  ----

                 1245                    2