`interpft`

Reamostra um sinal periódico usando a transformada de Fourier.

Sintaxe

y = interpft(x, n)

Retorna n amostras igualmente espaçadas do sinal periódico cujas amostras são x.

y = interpft(x, n, dim)

Reamostra ao longo da dimensão dim.

Descrição

O sinal é transformado com a FFT, preenchido com zeros (ou truncado) no domínio da frequência até o comprimento n, e transformado de volta. É uma interpolação trigonométrica exata para sinais de banda limitada e periódicos no intervalo amostrado.

Quando n é maior que o comprimento da entrada, o resultado interpola entre os pontos originais, que continuam presentes; a saída é escalada de modo a preservar as amplitudes. Escolher n menor que o comprimento da entrada reduz a taxa de amostragem e descarta as altas frequências.

Argumentos de entrada

x: Vetor ou array de amostras igualmente espaçadas ao longo de um período.
n: Número de amostras na saída.
dim: Dimensão de trabalho. Por padrão, a primeira dimensão não singleton.

Argumentos de saída

y: Sinal reamostrado com n pontos na dimensão escolhida.

Exemplo

Dobrando as amostras de um período de seno

x = [0 1 0 -1];
y = interpft(x, 8);
disp(round(y*1000)/1000)

Saída

         0    0.7070    1.0000    0.7070         0   -0.7070   -1.0000   -0.7070

Veja também

fft ifft interp1 resample