基于MATLAB实现的对机械类振动信号处理，幅值谱，相位谱，自功率谱，自回归谱等+使用说明文档.rar

function Amplitude() clear clc clf conna=database('mydb','',''); %连接数据库 curs=exec(conna,'select * from X118_BA_time$'); %建立并打开游标 curs=fetch(curs); %把数据库中的数据读取到MATLAB中 y=cell2mat(curs.Data); %把数据库读取到的数据赋值给y close(curs); %关闭游标 close(conna); %关闭数据库连接 len = length(y); for i = 1:20 if (len < 2^i) % 如果2^i大于y的长度 N = 2^(i-1); % 则截取y的前2^(i-1)个数据参与计算，后面的部分忽略掉 break; end end fs=1000; %fs的频率为1000 unit='mm'; %单位为mm Y = fft(y, N); % 快速傅里叶变换Y = Y(k)，k = 1,2,...N f = fs/N*(0:N/2); % 频率步长Δf = 1/T = 1/((1/fs)*N) = fs/N Ay = abs(Y)/N; % 双边幅值谱 % 单边幅值谱 A = 2*Ay(1:N/2+1); hold on plot(f, A) title('幅值谱') xlabel('频率 f /Hz') ylabel(['幅值 A /',unit]) hold off grid on