svfft.m

来自「有关matlab的电子书籍有一定的帮助希望有用」· M 代码 · 共 30 行

function [errstr,P,f] = svfft(x,Fs,valueArray)
%SVFFT spectview Wrapper for FFT.
%  [errstr,P,f] = svfft(x,Fs,valueArray) computes the power spectrum P
%  at frequencies f using the parameters passed in via valueArray:
%
%   valueArray entry     Description
%    ------------         ----------
%          1                Nfft

% Copyright (c) 1988-98 by The MathWorks, Inc.
% $Revision: 1.2 $ $Date: 1997/11/26 20:14:47 $

errstr = '';
P = [];
f = [];

err = 0;
nfft = valueArray{1};
evalStr = 'P = (abs(fft(x,nfft)).^2)/nfft; f=(0:nfft-1)''*Fs/nfft;';

eval(evalStr,'err = 1;')
if err
    errstr = {'Sorry, couldn''t evaluate fft; error message:'
               lasterr };
    return
end

[P,f] = sbswitch('svextrap',P,f,nfft);