fzzyfft.m

学matlab还有比这个更好的吗?我千兴万苦才找到的,这是一个学习matlab的视频,功能强大

function [A,B,C,fn,t,w]=fzzyfft(T,M,Nf)
% 利用FFT，计算[0,T}区间上定义的时间波形的Fourier级数展开系数A,B和频谱C,fn 。
%	T		时间波形周期
%	M		用作2的幂次
%	Nf		输出谐波的阶次，决定A,B的长度为(Nf+1)。Nf不要超2^(M-1)。
%	A,B		分别是Fourier级数中cos,sin展开项的系数。A(1)是直流量。
%	C		是定义在[-fs/2,fs/2]上的频谱
%	t,w		是原时间波形数据对
if (nargin<2 | isempty(M));M=8;end										%		<9>
if nargin<3;Nf=6;end														%		<10>
N=2^M;
f=1/T;
w0=2*pi*f; 
dt=T/N;
n=0:1:(N-1);
t=n*dt;
w=time_fun(t,T);															%		<17>
%
W=fft(w);
cn=W/N;
%
z_cn=find(abs(cn)<1.0e-10);
cn(z_cn)=zeros(length(z_cn),1); 										%		<23>
cn_SH=fftshift(cn);														%		<24>
C=[cn_SH cn_SH(1)];
A(1)=C(N/2+1);
A(2:N/2+1)=2*real(C((N/2+2):end));
B(2:N/2+1)=-2*imag(C((N/2+2):end));
if Nf>N/2;error(['第三输入宗量 Nf 应小于 ' int2str(N/2-1)]);end
A(Nf+2:end)=[];
B(Nf+2:end)=[];
n1=-N/2:1:N/2; 
fn=n1*f;