仿真轴承振动故障信号
纵观整个论坛,就没发现有现成的编好的仿真轴承故障信号的matlab程序,而且某些国内学报类EI仿真的信号,编写的公式和画的图也不对应。于是,我自己画了一个,不知道对不对,希望大家批评指正!扣扣 462427140A=0.5; %冲击幅值
B=0.7; %阻尼系数
fm=0; %调制频率
fn=1300; %轴承固有频率
fs=10000; %采样频率
N=10000; %采集点数
Te=N/fs; %采样时间
T=0.02; %时间间隔
M=(Te)/T; %M值
a=0.0001; %时间间隔
t=0:a:(Te-a); %时间总点数
tt=0:a:T-a; %时间序列
b=(M)*length(tt); %一共的点数
x=zeros(size(t));
for i=1:M
for j=1:length(tt)
S(i,j)=(A*cos(2*pi*fm*(tt(j)))).*exp(-B*(tt(j)-i*T)).*cos(2*pi*fn*(tt(j)-i*T));
end
% Y(1,j)=Y(1,:)+S(i,:);
end
x=reshape(S,1,b);
%绘制时域波形
plot(t,x);xlabel('时间 t/s');ylabel('幅值');title('仿真信号时域波形');
这个故障频率是50hz,为什么我用频谱分析缺求不出来呢?
页:
[1]