cftbyfft.m

书名应为《掌握和精通MATLAB》附有完整的源程序

function [AW,f]=cftbyfft(wt,t,flag)
%cftbyfft.m
%本程序采用FFT计算连续时间Fourier变换。输出幅频谱数据对(f,AW)。
%	输入量(wt,t)为已经窗口化了的时间函数wt(t)，它们分别是长度为N的向量。
%	对于"非平凡"取值时段有限的情况，应使该时段与窗口长度相比足够小。以
%	提高频率分辨率。
%	对于"非平凡"取值时段无限的情况，窗口长度的选取应使窗口外的函数值小
%	到可忽略，以提高近似精度。
%	输入宗量flag控制输出CFT的频率范围。
%		flag取非0时（缺省使用），频率范围在[0,fs）；
%		flag取0时，频率范围在[-fs/2,fs/2)。
if nargin==2;flag=1;end
N=length(t);
T=t(length(t))-t(1);
dt=T/N;
W0=fft(wt);																	%<16>
W=dt*W0;
df=1/T;
n=0:1:(N-1);
%
if flag==0
   n=-N/2:(N/2-1);
   W=fftshift(W);	
end
f=n*df;
AW=abs(W);	
if nargout==0
	plot(f,AW);grid,xlabel('频率f');ylabel('|w(f)|')
end