window.octave

来自「各种工程计算的库函数」· OCTAVE 代码 · 共 25 行

#!/usr/bin/octave
#
# script to read window from file and plot transformed result

wfile = fopen("window.txt", "r", "native");

[win, winsize] = fscanf(wfile, "%g");

#win = win * (1 / (max(win) - min(win)));
winlen = length(win);
winpad = win;
winpad(winlen:(2 * winlen)) = 0;
wint = fft(winpad);
wins = 20 * log10(abs(wint) / (winsize / 2));
wins = wins * (1 / (max(wins) - min(wins))) + 1;

hold on;
plot(win);
plot(wins(1:winsize));

fclose(wfile);

keyboard("Done> ");