matlab中fft插值,在加窗插值FFT分析中，如何用matlab求得双谱线插值的修正...

下面是我在之前论坛里找到的程序里进行改编的，但还是弄不来，希望懂的人能指导下

f0=50;

fs=50*256;

N=2048;%2048/256=8

rad=180/pi;

n=[0:1:N-1];

xb=[220,0.5,25,0.4,6,0.3,4,0.2,2];

Q=[pi/3,pi/3,pi/6,pi/3,pi/3,pi/3,pi/3,pi/3,pi/3];

s=zeros(1,N);

%s=220*sin(2*pi*f0*n./fs+pi/6)+0.5*sin(4*pi*f0*n./fs+pi/5)+0.2*sin(8*pi*f0*n./fs+pi/7)+0.07*sin(16*pi*f0*n./fs+pi/5)+0.5*sin(32*pi*f0*n./fs+pi/5);

for i=1:9

s=s+xb(i)*sin(2*pi*f0*i*n./fs+Q(i));

end

w=0.338946-0.481973*cos(2*pi*n./N)+0.161054*cos(4*pi*n./N)-0.018027*cos(6*pi*n./N);

r=s.*w;

v=fft(r,N);

u=abs(v);

vf=abs(v)/1024;

%为了求出谐波，那么必须找出基波，在35-65中频谱最大值所对应的就是离基波频率最近的那个点

y=0;

for i=fix(35/(fs/N)):fix(65/(fs/N))

if u(i&#