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

ddddaaaaa

程序如下：function resultCorrect=spectrumcorrectenergymethodshare(inputDate,correctNum,fs)
%功能：离散频谱校正能量重心法，适用单频点信号校正，只采用了加汉宁窗的结果
%注意:信号的模型为Acos(2*pi*f*t+pha),注意t从0开始，correctNum为采用校正的正点数，汉宁窗通常采用两点就可以获得很高的精度
%输入：inputDate待分析数据，数据长度为偶数，统一为行向量；fs采样频率
1 t& M- {6 s2 Y: N% e! @%输出：resultCorrect校正后的频率，幅值，相位结果
resultCorrect=zeros(1,3);   
0 _2 ?8 h. l* h% z7 m' ?N=length(inputDate);    %数据长度
0 L: g5 X) e9 N) E4 V0 `w=hann(N,'periodic');   %生成汉宁窗
, y5 J, I( A0 I  F1 rfftDate=fft(inputDate.*w');  
/ P+ {1 ?  O" _k=2.667;                %汉宁窗恢复系数
fftDate=fftDate(1:N/2)/N*2;      %单边复数谱  
5 W# z: v8 J$ o! n; g' x' L2 UfftDateMag=abs(fftDate);         %单边幅值谱
0 {9 `) u% s9 m6 @fftDatePower=fftDateMag.^2;      %单边功率谱
[~,maxIndex]=max(fftDatePower);     %功率最大值对应位置
0 p6 X' m5 i% p% @maxAngle=angle(fftDate(maxIndex));  %最大值处对应的相位  
dn=-correctNum:correctNum;
f=sum((maxIndex+dn).*fftDatePower(maxIndex+dn))/sum(fftDatePower(maxIndex+dn));   %归一化校正频率      
' o9 Y5 m+ |. Y* {resultCorrect(1,1)=(f-1)*fs/N;       %频率校正结果,注意matlab下标是从1开始的
resultCorrect(1,2)=sqrt(k*sum(fftDatePower(maxIndex+dn)));    %校正幅值结果            
  J* b! D2 m" J7 R; a' v, \resultCorrect(1,3)=maxAngle+pi*(maxIndex-f);                  %校正相位结果
resultCorrect(1,3)=mod(resultCorrect(1,3),2*pi);
resultCorrect(1,3)=resultCorrect(1,3)-(resultCorrect(1,3)>pi)*2*pi;    %象限定在（-pi,pi]  
end
0 {: W; L8 s* F$ T. S/ X- j2 u; k$ @& n可仿真看下效果，误差还是很小的
t=0:0.01:1-0.01;
* L  G+ T" w; L0 v2 yx=4.2*cos(2*pi*5.4*t+0.4);
8 P- o$ P* }$ d# e+ r9 MresultCorrect=spectrumcorrectenergymethodshare(x,2,100)
1 {, Y  H/ q0 n4 m8 _7 EresultCorrect =

" y( f3 C: K7 x% W  m( k     5.3993    4.1995    0.4021
* c4 h1 D9 n* j7 X
9 V7 }5 c3 q: N
7 Y+ |. V7 T) u+ s0 q) \
' D3 {2 w; e5 w4 ?6 G0 a2 w7 R

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

! B$ }  }9 A, L  S

多言数穷

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