资源描述：
《matlab振动分析程序》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、matlab振动分析程序%这是练习60000rpm振动loadyitaizhendongb_6.txt;%加载测试源数据x=yitaizhendongb_6;%赋值N=6400;fs=25600;%源数据中有6400个，我们测试时的采样频率为25600Hzn=0:N-1;t=n/fs;figure(1)%%%%%%%%%%%%%%%%%原始振动时域，频域图，自相关后时域，自相关频域图subplot(2,2,1),plot(t,x);%画出时域图title('振动原始信号'),xlabel('时间/s'),ylabel('幅值/g');axis([00.25-100100]);g
2、ridon;y=fft(x,N);%计算频谱mag=abs(y);f=n*fs/N;subplot(2,2,2),plot(f(1:N/2),mag(1:N/2)*2/N);%奈奎斯特频率之前的频谱及真实幅值title('频域图'),xlabel('频率/Hz'),ylabel('幅值/g');set(gca,'xtick',0:2000:f(N/2));%设置X轴的刻度间距axis([01200006]);gridon;Lag=100;[c,lags]=xcorr(x,Lag,'unbiased');%原始信号自相关subplot(2,2,3),plot(lags/fs,c)
3、;title('自相关后信号'),xlabel('时间/s'),ylabel('幅值/g');ylim([-100,100]);gridonm=length(c)-1;%m值为奇数，会出现警告，故减一处理z=fft(c,m);%原始信号自相关后，求其频谱mag_z=abs(z);ff=(0:m-1)*fs/m;subplot(2,2,4),plot(ff(1:m/2),mag_z(1:m/2)*2/m);title('自谱图'),xlabel('频率/Hz'),ylabel('幅值/g');xlim([012000]),gridon%[u,v]=ginput;%holdon%p
4、lot(u,v,'or');%显示所取的点的位置,'or'表示点击的那个点会用圆圈标记显示%若后面还有plot作图，而又不想放在2,2,4图中，必须使用holdoff，我们这里没有图了，故没有使用holdoff%text(u(1),v(1),['u(1)=',num2str(u(1)),'v(1)=',num2str(v(1))]);%text(u(2),v(2),['u(2)=',num2str(u(2)),'v(2)=',num2str(v(2))]);%text(u(3),v(3),['u(3)=',num2str(u(3)),'v(3)=',num2str(v(3))]
5、);%text(u(4),v(4),['u(4)=',num2str(u(4)),'v(4)=',num2str(v(4))]);%text(u(5),v(5),['u(5)=',num2str(u(5)),'v(5)=',num2str(v(5))]);%text(u(6),v(6),['u(6)=',num2str(u(6)),'v(6)=',num2str(v(6))]);%text(u(7),v(7),['u(7)=',num2str(u(7)),'v(7)=',num2str(v(7))]);%在自谱图取几个特征点，并显示坐标z=cceps(x);%计算倒谱figur
6、e(2)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%倒谱图subplot(2,1,1),plot(t,z);%画出倒谱图title('倒谱'),xlabel('时间/s'),ylabel('幅值/g');xlim([0,0.005]);%可以设置x轴范围，可以看的更清楚点grid[u,v]=ginput;text(u(1),v(1),'1');text(u(2),v(2),'2');text(u(3),v(3),'3');text(u(4),v(4),'4');k=[1,2,3,4];ft=[1/u(1),1/u(2),1/u(
7、3),1/u(4)]subplot(2,1,2),stem(k,ft,'r');title('特征频率'),ylabel('f/Hz'),xlabel('取点数');xlim([14]);grid[u1,v1]=ginput;text(u1(1),v1(1),[num2str(u1(1)),'，',num2str(v1(1))]);text(u1(2),v1(2),[num2str(u1(2)),'，',num2str(v1(2))]);text(u1(3),v1(3),[num2str(u
