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

ddddaaaaa

程序如下：function resultCorrect=spectrumcorrectenergymethodshare(inputDate,correctNum,fs)
; i: X4 i$ z3 f( f: |4 V2 S" J%功能：离散频谱校正能量重心法，适用单频点信号校正，只采用了加汉宁窗的结果
%注意:信号的模型为Acos(2*pi*f*t+pha),注意t从0开始，correctNum为采用校正的正点数，汉宁窗通常采用两点就可以获得很高的精度
%输入：inputDate待分析数据，数据长度为偶数，统一为行向量；fs采样频率
% e# f+ z8 ]8 M& e%输出：resultCorrect校正后的频率，幅值，相位结果
resultCorrect=zeros(1,3);   
- q+ b* m: L7 qN=length(inputDate);    %数据长度
6 }9 o& d( ~; X; R: I; r$ o/ rw=hann(N,'periodic');   %生成汉宁窗
fftDate=fft(inputDate.*w');  
- F. t8 G  G4 I/ ~* ~$ _% ]( f/ dk=2.667;                %汉宁窗恢复系数
0 _0 p9 \5 F! G0 q; ~fftDate=fftDate(1:N/2)/N*2;      %单边复数谱  
( @$ ]5 C8 W, L( Y1 ~) offtDateMag=abs(fftDate);         %单边幅值谱
fftDatePower=fftDateMag.^2;      %单边功率谱
6 _5 ]. l" w: B[~,maxIndex]=max(fftDatePower);     %功率最大值对应位置
maxAngle=angle(fftDate(maxIndex));  %最大值处对应的相位  
4 a; @1 L. L. Y- idn=-correctNum:correctNum;
f=sum((maxIndex+dn).*fftDatePower(maxIndex+dn))/sum(fftDatePower(maxIndex+dn));   %归一化校正频率      
resultCorrect(1,1)=(f-1)*fs/N;       %频率校正结果,注意matlab下标是从1开始的
resultCorrect(1,2)=sqrt(k*sum(fftDatePower(maxIndex+dn)));    %校正幅值结果            
& j9 d; c# {3 V4 t, zresultCorrect(1,3)=maxAngle+pi*(maxIndex-f);                  %校正相位结果
* X* J3 d% B  a' \resultCorrect(1,3)=mod(resultCorrect(1,3),2*pi);
resultCorrect(1,3)=resultCorrect(1,3)-(resultCorrect(1,3)>pi)*2*pi;    %象限定在（-pi,pi]  
; [; ]$ w. i5 d9 A. nend
可仿真看下效果，误差还是很小的
2 s# N3 b. I: `, o( s1 z, @t=0:0.01:1-0.01;
x=4.2*cos(2*pi*5.4*t+0.4);
0 y; a5 w$ k! e& Z4 t( m2 g0 J9 `resultCorrect=spectrumcorrectenergymethodshare(x,2,100)
1 Y% S* \2 W  X# }resultCorrect =

! m. ~# r7 _" ]! G9 j6 h2 w. w/ ~  r, z     5.3993    4.1995    0.4021

3 c' k/ V. Z5 }$ Y- Z; x# R5 l
9 e0 k, C9 f9 z, s2 K  W6 Y
7 j3 `1 N0 V$ B/ K' O& ^# O

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

, ?; V7 X4 @3 \  }1 o+ M

多言数穷

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