yummyking 发表于 2015-6-6 10:52

仿真轴承振动故障信号

纵观整个论坛,就没发现有现成的编好的仿真轴承故障信号的matlab程序,而且某些国内学报类EI仿真的信号,编写的公式和画的图也不对应。于是,我自己画了一个,不知道对不对,希望大家批评指正!扣扣      462427140


A=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]
查看完整版本: 仿真轴承振动故障信号