N=64;
t=0:N-1;
fs=64;
f0=4.33;a=0.87;ph=45;%f=[0:fs/N/f0:N-1];
y=a*cos（2*pi*t*f0/fs+ph*pi/180）;
y1=y（1:N）;?%N個數據
%%%%獲取對稱窗￥￥￥￥￥￥￥￥￥￥￥￥￥￥￥￥￥
win=hanning（N‘periodic‘）‘;
%$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

win1=win/sum（win）;?%窗歸一
y11=y1.*win1;

y1_fft=fft（y11N）;
a1=abs（y1_fft）;?%fft振幅譜
a11=a1;
a2=[a11（2:N）?0];%%方便畫圖
p1=angle（y1_fft）*180/pi;%fft相位譜
p2=[p1（2:N）