离散频谱的能量重心校正法程序

ddddaaaaa

程序如下：function resultCorrect=spectrumcorrectenergymethodshare(inputDate,correctNum,fs)
" x6 [5 n8 f8 ~  a( m" l%功能：离散频谱校正能量重心法，适用单频点信号校正，只采用了加汉宁窗的结果
%注意:信号的模型为Acos(2*pi*f*t+pha),注意t从0开始，correctNum为采用校正的正点数，汉宁窗通常采用两点就可以获得很高的精度
9 z  z& v7 L4 Y" s%输入：inputDate待分析数据，数据长度为偶数，统一为行向量；fs采样频率
%输出：resultCorrect校正后的频率，幅值，相位结果
$ P- w) P) S, Y0 p  E4 k0 dresultCorrect=zeros(1,3);   
N=length(inputDate);    %数据长度
: r, p4 x3 g7 Ew=hann(N,'periodic');   %生成汉宁窗
& @2 P$ T! [! A; \% K! _# z+ S. P' QfftDate=fft(inputDate.*w');  
k=2.667;                %汉宁窗恢复系数
6 r$ R/ v/ G2 m0 G0 i6 LfftDate=fftDate(1:N/2)/N*2;      %单边复数谱  
fftDateMag=abs(fftDate);         %单边幅值谱
. L4 p5 _7 k) c+ I+ X- bfftDatePower=fftDateMag.^2;      %单边功率谱
[~,maxIndex]=max(fftDatePower);     %功率最大值对应位置
: x. B: \5 q: f# ImaxAngle=angle(fftDate(maxIndex));  %最大值处对应的相位  
dn=-correctNum:correctNum;
f=sum((maxIndex+dn).*fftDatePower(maxIndex+dn))/sum(fftDatePower(maxIndex+dn));   %归一化校正频率      
resultCorrect(1,1)=(f-1)*fs/N;       %频率校正结果,注意matlab下标是从1开始的
5 t2 ^3 D" {$ E6 R0 W& T* vresultCorrect(1,2)=sqrt(k*sum(fftDatePower(maxIndex+dn)));    %校正幅值结果            
resultCorrect(1,3)=maxAngle+pi*(maxIndex-f);                  %校正相位结果
" t' S; L/ }& e: i, N$ SresultCorrect(1,3)=mod(resultCorrect(1,3),2*pi);
) H# h% V* ^. q9 ?4 AresultCorrect(1,3)=resultCorrect(1,3)-(resultCorrect(1,3)>pi)*2*pi;    %象限定在（-pi,pi]  
end
6 G  H, `% K/ z可仿真看下效果，误差还是很小的
8 Q6 x5 B3 `0 vt=0:0.01:1-0.01;
x=4.2*cos(2*pi*5.4*t+0.4);
+ h% r$ ]8 ^8 L2 d* p( {; j& v; o  TresultCorrect=spectrumcorrectenergymethodshare(x,2,100)
: L& e, |, D; s; m; MresultCorrect =
' a0 ~1 I% b9 P$ V, i- e( L
     5.3993    4.1995    0.4021
/ z& M* f0 G' w+ W3 j

游客，如果您要查看本帖隐藏内容请回复

5 X* r4 T8 `8 q- T8 [$ A: A/ p

多言数穷

离散频谱的能量重心校正法程序