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

ddddaaaaa

程序如下：function resultCorrect=spectrumcorrectenergymethodshare(inputDate,correctNum,fs)
5 s3 G' e  y8 Y$ ]- `%功能：离散频谱校正能量重心法，适用单频点信号校正，只采用了加汉宁窗的结果
" X( g9 G  ?: g6 Z%注意:信号的模型为Acos(2*pi*f*t+pha),注意t从0开始，correctNum为采用校正的正点数，汉宁窗通常采用两点就可以获得很高的精度
& h- L! U+ r! W9 i! B' `1 f%输入：inputDate待分析数据，数据长度为偶数，统一为行向量；fs采样频率
. Z2 x/ Z! ]+ E% Y6 T' @%输出：resultCorrect校正后的频率，幅值，相位结果
resultCorrect=zeros(1,3);   
N=length(inputDate);    %数据长度
w=hann(N,'periodic');   %生成汉宁窗
2 Q& h% I9 q& ufftDate=fft(inputDate.*w');  
k=2.667;                %汉宁窗恢复系数
. l5 \6 r+ r+ Y4 g, QfftDate=fftDate(1:N/2)/N*2;      %单边复数谱  
. K) v& r. N+ V9 N* afftDateMag=abs(fftDate);         %单边幅值谱
6 f. x& v& Z# offtDatePower=fftDateMag.^2;      %单边功率谱
[~,maxIndex]=max(fftDatePower);     %功率最大值对应位置
maxAngle=angle(fftDate(maxIndex));  %最大值处对应的相位  
- R# e% g, m0 G6 ~" N6 {+ \dn=-correctNum:correctNum;
f=sum((maxIndex+dn).*fftDatePower(maxIndex+dn))/sum(fftDatePower(maxIndex+dn));   %归一化校正频率      
resultCorrect(1,1)=(f-1)*fs/N;       %频率校正结果,注意matlab下标是从1开始的
7 o  ^. |, X1 B; V& h- T8 I' i6 @0 C* ?resultCorrect(1,2)=sqrt(k*sum(fftDatePower(maxIndex+dn)));    %校正幅值结果            
resultCorrect(1,3)=maxAngle+pi*(maxIndex-f);                  %校正相位结果
resultCorrect(1,3)=mod(resultCorrect(1,3),2*pi);
, Z  |$ c" n9 W! TresultCorrect(1,3)=resultCorrect(1,3)-(resultCorrect(1,3)>pi)*2*pi;    %象限定在（-pi,pi]  
end
可仿真看下效果，误差还是很小的
% h( O) H1 I5 H& d, F4 x7 J9 e# `+ et=0:0.01:1-0.01;
* T5 |* h, b9 `1 t- Tx=4.2*cos(2*pi*5.4*t+0.4);
0 R0 B2 @( b) k1 [5 i7 NresultCorrect=spectrumcorrectenergymethodshare(x,2,100)
resultCorrect =
8 [& ~* H8 Y! }3 x/ Y
# q5 O% f4 u' h' z) E# m     5.3993    4.1995    0.4021



9 G; B, R0 ]/ V1 Y

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

多言数穷

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