|
|
马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。
您需要 登录 才可以下载或查看，没有账号？我要加入
x
搜索了论坛发现很多事求时域的1/3倍频程。请问下怎么求频域的倍频程？
我有一个频域-声压曲线。横坐标HZ纵坐标db
可否直接用下面程序（MATLAB在振动信号处理中的应用）转成1/3倍频啊~做出来的数值不对~~能不能指点下啊
f=[1.00 1.25 1.60 2.00 2.50 3.15 4.00 5.00 6.30 8.00];
fc=[f,10*f,100*f,1000*f,10000*f];
%中心频率与下限频率的比值
oc6=2^(1/6);
nc=length(fc);
n=length(x);
nfft=2^nextpow2(n);
a=fft(x,nfft);
for j=1:nc
fl=fc(j)/oc6;
fu=fc(j)*oc6;
nl=round(fl*nfft/sf+1);
nu=round(fu*nfft/sf+1);
if fu>sf/2
m=j-1;break
end
b=zeros(1,nfft);
b(nl:nu)=a(nl:nu);
b(nfft-nu+1:nfft-nl+1)=a(nfft-nu+1:nfft-nl+1);
c=ifft(b,nfft);
yc(j)=sqrt(var(real(c(1:n))));
end
plot(fc(1:m),yc(1:m));
|
|
18 点体能
发表于 2015-11-11 09:26
提升卡
变色卡
显身卡
楼主