th=(0:127)/128*2*pi;
x=cos(th);
y=sin(th);   %計算復平面上的單位圓
f=(abs(fft(ones(10,1),128)))';  % 計算一步頻率響應的幅值
stem3(x,y,f,'*');    %繪制三維枝干圖
xlabel('實部');
ylabel('虛部');
zlabel('幅值');
title('頻率響應');    %進行標注