How to use pmtm

[EDIT: 20110608 10:11 CDT - reformat - WDR]

Hi I need to make the power spectrum of some neuron samples and I have to use the function pmtm, but I don't know exactly how.

My inputs are voltage (mv) and time (s) and I have one data each 0,00005 seconds. This is all i have, I send you my script, I hope you could help me. Thank you so much.

%Este script obtiene los power spectrum de cada canal.

[FileName,PathName]=uigetfile('*.*','Select DATA files');

%uigetfile me pregunta por un archivo ('*.*' cualquier archivo) 

b=input('numero de filas de texto del archivo ');%pregunto 

%Ahora para identificar qué es que en A me dice qué es dato y que es texto, del fichero.

A=importdata([PathName,FileName],'\t',b); %'\t'=separa columnas,b=lineas de texto del archivo
%A=importadata(filename, delimiter, nheaderlines)

%ahora voy a extraer solo los datos, los numeritos, D son la parte de datos de A
D=A.data;

[n,canales]=size(D); %n es el numero de filas de los datos, y canales es el numero de columnas de los datos

t=D(:,1); %la primera columna de datos es el tiempo

for i=2:canales %el resto de columnas son los datos de potenciales
    if i==10
        V=(1/10)*D(:,i);%El canal 10 es 10 veces mayor... hay que cambiarlo.
    elseif i==2
        V=D(:,i+1);
    else
        V=D(:,i);
    end

      x=V;
      nw=4;
      nfft=[];
      fs=20000; %es la frecuencia de muestreado fs=1/tmuestreado, como t muestreado es 0,00005

      [P,f]=pmtm(x,nw,nfft,fs);

      %x es el vector de datos
      %nw es el producto de tiempo de ancho de banda
      %nfft es la frecuencia en la que evaluamos P
      %P es power spectrum calculado por la via multitaper.
      %f es el vector frcuencias.

      figure
      plot(f,P)
      xlabel('frecuencia (Hz)');
      ylabel('PSD');
  end

[Wayne King]

So what exactly is your question about the multitaper method? You've posted code that we cannot execute because we do not have your data.