//離散フーリエ変換

clear
n=32;
xd=(1:n);

for k = 1:n,
  for l =1:n, 
      a(k,l) = exp(-%i*(k-1)*(l-1)*2*%pi/n);
  end;
end;
a;

//入力ベクトル（１/８デューティの方形波）
x=[ones(1,n/8),zeros(1,7*n/8)]';

//入力ベクトル（のこぎり波）
//x=[1:1:n]';

//DFT変換
y=a*x;

xbasc();
//bar(xd,x,0.5);
bar(xd,abs(y),0.5);