多种频谱校正方法及matlab代码.pdf
《多种频谱校正方法及matlab代码.pdf》由会员分享,可在线阅读,更多相关《多种频谱校正方法及matlab代码.pdf(6页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、多种频谱校正方法采样间隔归一化成1T,采样长度为N.这样 FFT 离散谱线为0,1)iXiN(,相应的频率分辨率2/(1/)NfN.设 FFT 离散谱线局部极高谱线为m(为了数学上简洁,假定从 0开始,注意在 MATLAB环境下数组实际操作的是从1 开始),记频偏量.我们需要使用谱线m和与之相邻一条次高谱线,记这连续两条谱线中左边一条序号为M(当次高谱线在m左侧时1Mm,反之Mm).下面列出若干算法的计算公式1.加矩形窗的精确谱校正1iiiXUjV111()sin()()cos()MMMMoptMMVVMUUMKUU1211cos()sin()cos()sin()optMMoptMMKMZVU
2、MKMZVUM2121cos()cos()()ZMZMMmZZ2.加矩形窗情形,采用解析单频模型的幅值比校正1,211|()|MMMXMmXX3.加汉宁窗情形,采用解析单频模型的幅值比校正1,2112|()|MMMMXXMmXX4.加矩形窗情形,采用解析单频模型的复比值校正311Re()MMMXMmXX5.加汉宁窗情形,采用解析单频模型的复比值校正3112()MMMMXXMmXX6.加矩形窗情形,采用解析单频模型的复合复比值校正311Re()MmMMXMmXX11mRmmXXX,1111mmLmmmmXXXXXX0.5)0.5)mLmR(7.加汉宁窗情形,采用解析单频模型的复合复比值校正311
3、2Re()MMmMMXXMmXX112mmRmmXXXX,1111221mmmmLmmmmXXXXXXXX0.5)0.5)mLmR(8.加矩形窗,Quin 校正 411Re()Re(),Re()Re()mmLRmmXXXX11LRLRLR,,00RRLR当且其它9.加汉宁窗,Quin 校正 411Re()Re(),Re()Re()mmLRmmXXXX212111LRLRLR,,00RRLR当且其它References1.Schoukens,J.,R.Pintelon,H.Van Hamme.The interpolatedfast Fouriertransform:Acomparative s
4、tudy.IEEE Transactions on Instrumentationand Measurement.1992,41(2):文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10
5、S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S
6、5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5
7、P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P
8、10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F
9、8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8
10、G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5226-232.2.谢明,丁康.频谱分析的校正方法.振动工程学报.1994,7(2):172-179.3.陈奎孚,王建立,张森文
11、.频谱校正的复比值法.振动工程学报(已投).2007.4.Quinn,B.G.EstimatingfrequencybyinterpolationusingFouriercoefficients.IEEETransactions on Signal Processing.1994,42(5):1264-1268.%=这是调用调试=DT=1;N=1024;PHI=pi/3;Ampl=1;CiR=11.9;%cycles in recordFreq=CiR/(DT*N);%frequencyTV=0:N-1;DatVec=Ampl*cos(Freq*TV*2*pi+PHI);FV=fft(DatV
12、ec);figuresubplot(2,1,1);plot(TV,DatVec);subplot(2,1,2);plot(abs(FV(1:round(N/2.56);gridonMV,MI=max(abs(FV);%加矩形窗的解析校正-1FreqShift=SpecCorr(FV,MI,N,1);%加矩形窗的解析单频模型校正-2FreqShift=SpecCorr(FV,MI,N,2);%加汉宁窗的解析单频模型校正-3HanDat=DatVec.*hanning(N,periodic);FV=fft(HanDat);FreqShift=SpecCorr(FV,MI,N,3);%加矩形窗的解析
13、单频模型校正+复比值法-4FV=fft(DatVec);FreqShift=SpecCorr(FV,MI,N,4);%加汉宁窗的解析单频模型校正+复比值法-5HanDat=DatVec.*hanning(N,periodic);文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4
14、L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9
15、C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK
16、9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5
17、C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1
18、W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:
19、CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5FV=fft(HanDat);Fre
20、qShift=SpecCorr(FV,MI,N,5);%加矩形窗的解析单频模型校正+复比值法+左右平均-6FV=fft(DatVec);FreqShift=SpecCorr(FV,MI,N,6);%加汉宁窗的解析单频模型校正+复比值法+左右平均-7HanDat=DatVec.*hanning(N,periodic);FV=fft(HanDat);FreqShift=SpecCorr(FV,MI,N,7);%加矩形窗的解析单频模型校正+Quinn 算法-8FV=fft(DatVec);FreqShift=SpecCorr(FV,MI,N,8);%加汉宁窗的解析单频模型校正+Quinn 算法-9H
21、anDat=DatVec.*hanning(N,periodic);FV=fft(HanDat);FreqShift=SpecCorr(FV,MI,N,9);=这是子程序=%spectrum correction assemble%the sampling interval is 1 s(or unitary)%Input:SpecVec-Discrte Fourier Spectrum from FFT%PeakIdx-the peakindex,noting the matrix in MatLab start from 1%TL-the length(or the point numbe
22、r)of the FFT%method-correction method%output:PeakShift-the corrected peak shifting from the peak in discrete%spectrumfunction PeakShift=SpecCorr(SpecVec,PeakIdx,TL,method)%picking up the secondhighest spectrum lineif(abs(SpecVec(PeakIdx-1)abs(SpecVec(PeakIdx+1)IP=PeakIdx-1,PeakIdx;ShiftCorr=-1;%shif
23、t aligning with the PeakIdxelseIP=PeakIdx,PeakIdx+1;ShiftCorr=0;%shift aligning with the PeakIdxendII=IP(1)-1;%noting that the index of a matrix in MATLABstarts from 1,not zeroif(method=1)%an analyitic solution for rectangular windowU=real(SpecVec(IP);文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK
24、9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5
25、C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1W9C5F8S5文档编码:CK9E8G5P5S5 HT5C4L3P10S6 ZY1
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 多种 频谱 校正 方法 matlab 代码
限制150内