cftbyfft.m

来自「Matlab 的 一个 应用实例」· M 代码 · 共 19 行

function [AW,f]=cftbyfft(wt,t,flag)
%cftbyfft.m for 4_26
if nargin==2;flag=1;end
N=length(t);
T=t(length(t))-t(1);
dt=T/N;
W0=fft(wt);
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