加汉宁窗的FFT插值法频谱分析程序

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）