MATLAB: Function definitions are not permitted…. trying to make an FFT thing

fftfunctionMATLAB

Function definitions are not permitted…. trying to make an FFT thing

I'm trying to get this working

close all
clear all
function[X,freq]=positiveFFT(x,Fs);
N=length(x);
k=0:N-1;
T=N/Fs;
freq=k/T; 
X=fft(x)/N; 
cutOff = ceil(N/2);
X = X(1/cutOff);
freq = freq(1:cutOff);

Any ideas?

Best Answer

I think it might be right. Just discovered that the function line has to be the top line of code otherwise it doesn't work.

Related Solutions

MATLAB: Quick question about elementwise

Hi Tom, I don't think you need the dot multiplication on q, just

   x=exp(-(t/t0).^2).*cos(w0*(1-q*t).*t);

On whether it's right or not, it's not readily apparent to me what you're trying to do with (cos(w0*(1-q*t).*t)

MATLAB: Suggestions for quantifying spectral leakage

"But then I realised this wouldn't work as if there was a 1MHz bin, there wouldn't be any spectral leakage."

The above statement is not correct. The spacing of the DFT bins is determined by the number of samples and the sampling frequency.

Leakage is a result of the convolution of the signal's spectrum with that of the window -- they are separate things.

For example:

      Fs = 4e6;
      t = 0:1/Fs:0.2-1/Fs;
      x = cos(2*pi*1e6*t);
      xdft = fft(x'.*hamming(length(t)));
      fundbin = (1e6/5)+1;
      xdft([fundbin-4:fundbin+4])

The frequency spacing is 5 Hz, you see that there is power not only at 1 MHz, but also at +/- 5 Hz

Best Answer

Related Solutions

MATLAB: Quick question about elementwise

MATLAB: Suggestions for quantifying spectral leakage

Related Question