马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
我有一笔加速规撷取数据单位为电压值经FFT转成功率频谱后, 要再化成1/3 octave band与 加速度RMS(dB)的频谱下面该如何作?
下面是目前程序
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
close all;
x=xlsread('data.xls','a1:a8192');
fs=1/0.002; N=length(x); freqStep=fs/N; freq=freqStep*(-(N-1)/2:(N-1)/2);
y=fft(y); y1=abs(y); magx1=fftshift(y1); T=16;
fl=[0.891,1.122,1.412,1.778,2.238,2.817,3.547,4.465,5.621,7.077,8.909,11.216,14.120,17.776,22.378,28.173,35.467,44.651,56.212,70.767,89.090,112.158,141.198,177.758,223.784];
fu=[1.122,1.413,1.779,2.240,2.819,3.550,4.469,5.626,7.082,8.916,11.225,14.131,17.790,22.396,28.195,35.495,44.686,56.256,70.823,89.160,112.246,141.310,177.898,223.961,281.950]; Fc=[1,1.25,1.6,2,2.5,3.15,4,5,6.3,8,10,12.5,16,20,25,31.5,40,50,63,80,100,125,160,200,250];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Fl为频宽下限, Fu为频宽上线, Fc为中心频率
?????????????
谢谢:@D
[ 本帖最后由 ChaChing 于 2010-4-20 13:48 编辑 ] |