声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

12
返回列表 发新帖
楼主: zhaoyixu

[FFT] 学写程序-比值校正法

[复制链接]
发表于 2010-8-15 17:06 | 显示全部楼层
deltf=1/(1+A(k)/A(k-1));为什么这么算呢
回复 支持 反对
分享到:

使用道具 举报

发表于 2011-5-11 22:05 | 显示全部楼层
Y=y(1:N/2)/N*2;
这句是什么意思啊?
发表于 2011-5-15 09:41 | 显示全部楼层
多谢分享!!!!!!!!!!!
发表于 2011-8-11 21:28 | 显示全部楼层
回复 6 # zhwang554 的帖子

您这段程序看上去没什么问题,是不是可看做采样率是1Hz,实际使用过程中采样率不一样。如果我将此程序加上采样率,现在设置为fs=2Hz,得出的结果也就不对了,为什么呢?是我使用错了吗?
close all;clc;clear all;
N=1024;
fs=2;
w=2*pi;
% t=(0:2*N-2)/fs;
t=(-N+1:N-1)/fs;
y=1.0*exp(j*(w*t*49.1/N+50.0*pi/180))+0.8*exp(j*(w*t*149.2/N+100*pi/180))+0.6*exp(j*(w*t*249.3/N+150*pi/180))+0.4*exp(j*(w*t*349.4/N+200*pi/180))+0.2*exp(j*(w*t*449.5/N+250*pi/180));
y1 = y(N:end);
win =  hanning(N)';;
win1 = win/sum(win);
y11= y1.*win1;
y11_fft = fft(y11,N);
a1 = abs(y11_fft);
p1 = mod(phase(y11_fft)*180/pi,360);
y2 = y(1:2*N-1);
win =  hanning(N)';;
winn =  conv(win,win);%apFFT须要卷积窗
win2 = winn/sum(winn);
y22= y2.*win2;
y222=y22(N:end)+[0 y22(1:N-1)];%构成长N的apFFT输入数据
y2_fft = fft(y222,N);
a2 = abs(y2_fft);
p2=mod( phase(y2_fft)*180/pi,360);
      ee=(p1-p2)/180/(1-1/N);
      aa=(a1.^2)./a2;
subplot(4,1,1),stem(a2,'.');title('apFFT振幅谱');ylim([0,1]);xlim([0 N/2]);grid
subplot(4,1,2),stem(p2,'.');title('初相位校正谱');ylim([0,400]);xlim([0 N/2]);grid
subplot(4,1,3);stem(ee,'.');title('频率校正谱');ylim([-1,1]);xlim([0 N/2]);grid
subplot(4,1,4);stem(aa,'.');title('振幅校正谱');ylim([0,1.5]);xlim([0 N/2]);grid
      disp('相位校正值')
     p2(50:100:450)
      disp('频率校正值')
    ee(50:100:450)+(49:100:450)
      disp('振幅校正值')
      aa(50:100:450)

相位校正值
ans =
   66.6762  161.6313  249.8269  234.3468  245.5667
频率校正值
ans =
   49.2541  149.4448  247.7433  349.6674  447.6453
振幅校正值
ans =
    0.4387    0.3441    0.2137    0.3907    0.6415

发表于 2011-12-6 19:27 | 显示全部楼层
回复 1 # zhaoyixu 的帖子

若被分析信号是谐波信号,不是由单一频率组成,那么你这种方法就不能对多个频率进行校正,那请问你有什么其他好的办法吗?
发表于 2011-12-6 19:29 | 显示全部楼层
还有就是hamming窗又怎么编程?
发表于 2011-12-10 09:31 | 显示全部楼层
这种方法只是对单频率信号而言,若是多频率信号了?
发表于 2011-12-11 10:51 | 显示全部楼层
若是间隔不是很密集的多个频率成分,可忽略掉相互之间的干涉,分别进行校正。
发表于 2012-3-15 10:30 | 显示全部楼层
Thank you very much!!!!
发表于 2012-6-1 23:10 | 显示全部楼层
回复 19 # shajue 的帖子

不是采样频率fs=2Hz,而是fs=1024Hz,此时频率分辨率为1Hz
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-16 07:35 , Processed in 0.054384 second(s), 15 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表