fftfcnn.m

simulink和matlab的图形显示接口程序

function [ftdata,ftangle]=fftfcnn(data,m,n,npoint,t)  
%w表示为那个网点，是为了计算相位补偿需要；
%此程序是一个网点的间谐波的计算程序，t是取多长的时间数据计算一次fft,也就是基波的f=1/t;
%注意计算间谐波fft的程序时候，如果取t=80mS,则可以分析12.5Hz的基波频率分量；
%千万要注意这个时候选取的要m>＝4,并且是偶数
%计算一个周期fft用的数据点数应该为(T/20)*128点.
%ftdata中每P行就是一个网点的一个通道的FFT的结果

k=t/20;
p=2*m/k-1;%表示给定的数据可以做多少次fft的计算
q=k*npoint;  %每计算一次fft的所需要数据的点数

   
for j=1:1:3*n 
for i=1:1:p   
 
   tp(i,:)=abs(fft(data(j,(1+(i-1)*q/2):(i*q/2+q/2)),q))/(q/2);
   qth(i,:)=angle(fft(data(j,(1+(i-1)*q/2):(i*q/2+q/2)),q)/(q/2));
   
end
ftdata((1+(j-1)*p):j*p,:)=tp;
ftangle((1+(j-1)*p):j*p,:)=qth;
end
ftdata(:,1)=ftdata(:,1)/2;
ftdata(:,2:q)=ftdata(:,2:q)/sqrt(2);