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

ddddaaaaa

程序如下：function resultCorrect=spectrumcorrectenergymethodshare(inputDate,correctNum,fs)
1 p; O. N7 c- J4 I%功能：离散频谱校正能量重心法，适用单频点信号校正，只采用了加汉宁窗的结果
%注意:信号的模型为Acos(2*pi*f*t+pha),注意t从0开始，correctNum为采用校正的正点数，汉宁窗通常采用两点就可以获得很高的精度
%输入：inputDate待分析数据，数据长度为偶数，统一为行向量；fs采样频率
9 |" v3 C0 a/ h. j9 r( }' A* H%输出：resultCorrect校正后的频率，幅值，相位结果
7 O4 J( r- c1 p+ KresultCorrect=zeros(1,3);   
9 D3 ^; H$ E- z3 o6 J* L6 R- h6 IN=length(inputDate);    %数据长度
+ J. M2 E7 A6 W! z) @% jw=hann(N,'periodic');   %生成汉宁窗
fftDate=fft(inputDate.*w');  
3 M& H2 u( p" Pk=2.667;                %汉宁窗恢复系数
fftDate=fftDate(1:N/2)/N*2;      %单边复数谱  
fftDateMag=abs(fftDate);         %单边幅值谱
fftDatePower=fftDateMag.^2;      %单边功率谱
4 R9 g$ B+ K: G/ c$ h; E[~,maxIndex]=max(fftDatePower);     %功率最大值对应位置
- `9 W# \& z+ i- j( s7 L. wmaxAngle=angle(fftDate(maxIndex));  %最大值处对应的相位  
dn=-correctNum:correctNum;
/ }0 g2 I* |" h# K' C, of=sum((maxIndex+dn).*fftDatePower(maxIndex+dn))/sum(fftDatePower(maxIndex+dn));   %归一化校正频率      
resultCorrect(1,1)=(f-1)*fs/N;       %频率校正结果,注意matlab下标是从1开始的
6 {- _; P* V8 a- f, SresultCorrect(1,2)=sqrt(k*sum(fftDatePower(maxIndex+dn)));    %校正幅值结果            
resultCorrect(1,3)=maxAngle+pi*(maxIndex-f);                  %校正相位结果
6 L' Q' o* v2 w5 ]  z+ |resultCorrect(1,3)=mod(resultCorrect(1,3),2*pi);
* ~+ `7 x, U" v) qresultCorrect(1,3)=resultCorrect(1,3)-(resultCorrect(1,3)>pi)*2*pi;    %象限定在（-pi,pi]  
2 }. w# q  Z" [) send
可仿真看下效果，误差还是很小的
" M7 ^9 F) i8 S$ S) vt=0:0.01:1-0.01;
x=4.2*cos(2*pi*5.4*t+0.4);
resultCorrect=spectrumcorrectenergymethodshare(x,2,100)
6 g2 J8 D7 F' s' IresultCorrect =
# i' @- V% X  B2 P4 l9 D( F* U
; f8 {8 r6 B: d8 C     5.3993    4.1995    0.4021
8 b' M# f' v3 E# z% b


4 J; s+ z! m: H4 n. F4 _9 U

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

9 N+ @. i/ ~+ g$ v

多言数穷

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