clear
//ＤＣＴ変換
n=4;
xd=(1:n);

for i = 1:n,
  for j = 1:n, 
      if i==1 then c=1/sqrt(2); else c=1;end;
      a(i,j) = sqrt(2/n)*c*cos((i-1)*(j-1/2)*%pi/n);
  end;
end;
//入力ベクトル（のこぎり波）
x=[1:1:n]',
//DCT変換
y=a*x,
xbasc();
bar(xd,y,0.5);