麻烦高手看看我的校正程序问题出现在什么地方？
f1=6.3;fs=400;
det=1/400;
t=0:det:1023*det;
N=length(t);
y1=sin(2*pi*f1*t);%+rand(1,N);%+200*sin(2*pi*6.3*t);
f=abs(fft(y1));
%plot(f);
=jiaozheng(y1,fs)
以上是信号模型
function =jiaozheng(s,fs)
index=0;
N=800;
L=200;
x=s(1:800);
xf=fft(x);
A=abs(xf);
A(1)=0;
=max(A);
A0=Amax*2/N;
phmax=angle(xf(index));
f=(index-1)*fs/N;
x1=s(201:1000);
=ft(x1,index-1);
phmax1=atan2(Ik,Rk);
dph=phmax1-phmax-2*pi*L*index/N;
dph=mod(dph,2*pi);
ifdph < (-1)*pi
dph=dph+2*pi;
end
ifdph > pi
dph=dph-2*pi;
end
df=(dph)/(2*pi*L/N);
fz=(index+df)*fs/N;
=ft1(df);
Az=Amax*2/Wk1/N/N;
上面是校正函数
function =ft1(f)
Rk1=0.0;
Ik1=0.0;
Num=800;
win=hanning(Num);
win1=win/sum(win);
for i=1:Num
Rk1=Rk1+win1(i)*cos(2*pi*(i-1)*(f));
Ik1=Ik1+win1(i)*sin(2*pi*(i-1)*(f));
end
Ik1=-Ik1;
Wk1=+sqrt(Rk1*Rk1+Ik1*Ik1)/Num;
plot(win1);
上面是ft1函数
function =ft(x,k)
Rk=0.0;
Ik=0.0;
Num=length(x);
for i=1:Num
Rk=Rk+x(i)*cos(2*pi*(i-1)*k/Num);
Ik=Ik+x(i)*sin(2*pi*(i-1)*k/Num);
end
Ik=-Ik;
Wk=sqrt(Rk*Rk+Ik*Ik)/Num;
上面是单点ft函数
为什么校正的幅值总是和实际相差很大呢？
我是参考《离散频谱分析校正理论与技术》这本书的4.5.2做的
如果频率是整数，则幅值为一 如果不是的话就相差很大！！
页:
[1]