1/3倍频程怎么做，着急！【matlab吧】_百度贴吧
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&签到排名：今日本吧第个签到，本吧因你更精彩，明天继续来努力！
本吧签到人数：0成为超级会员，使用一键签到本月漏签0次！成为超级会员，赠送8张补签卡连续签到：天&&累计签到：天超级会员单次开通12个月以上，赠送连续签到卡3张
关注：178,292贴子：
1/3倍频程怎么做，着急！收藏
我现在有用tes1350A采集的数据，想做频谱分析，着急！ 哥哥姐姐们谁能帮帮我
江苏盘锦天燃气模温机，一台省得让你怀疑人生的燃气锅炉
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%s = xlsread('ay.xls');%输入时程数据sf=256;
%采样频率x=s(:,2); %定义三分之一倍频程的中心频率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);%大于并接近n的2的幂次方长度nfft=2^nextpow2(n);%FFT变换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
%以每个中心频率段为通带进行带通频率滤波
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(b(1:n))));end%绘制输入时程曲线图形subplot(2,1,1);t=0:1/sf:(n-1)/plot(t,x);xlabel('时间(s)');ylabel('加速度(g)');%绘制三分之一倍频程有效值图形subplot(2,1,2);plot(fc(1:m),yc(1:m));xlabel('频率(Hz)');ylabel('有效值');
回复 bin2049 :你帮我看看看这个图对吗
就是用上面的程序做出来的
能不能做成直方图啊
围观、、、、、、、、、、这里是一个小尾巴 纯手打
终于看到一个跟噪声有关的
登录百度帐号推荐应用