pujiaozheng.txt

谱校正的MATLAB程序

%========================这是调用调试==================
DT=1;
N=1024;
PHI=pi/3;
Ampl=1;
CiR=11.9; %cycles in record
Freq=CiR/(DT*N); %frequency

TV=[0:N-1];

DatVec=Ampl*cos(Freq*TV*2*pi+PHI);

FV=fft(DatVec);
figure
subplot(2,1,1);plot(TV,DatVec);
subplot(2,1,2);plot(abs(FV(1:round(N/2.56))));grid on

[MV,MI]=max(abs(FV));

%加矩形窗的解析校正--1
FreqShift=SpecCorr(FV,MI,N,1);

%加矩形窗的解析单频模型校正--2
FreqShift=SpecCorr(FV,MI,N,2);

%加汉宁窗的解析单频模型校正--3
HanDat=DatVec.*hanning(N,'periodic')';
FV=fft(HanDat);
FreqShift=SpecCorr(FV,MI,N,3);

%加矩形窗的解析单频模型校正+复比值法--4
FV=fft(DatVec);
FreqShift=SpecCorr(FV,MI,N,4);

%加汉宁窗的解析单频模型校正+复比值法--5
HanDat=DatVec.*hanning(N,'periodic')';
FV=fft(HanDat);
FreqShift=SpecCorr(FV,MI,N,5);