仿真轴承振动故障信号

yummyking

纵观整个论坛，就没发现有现成的编好的仿真轴承故障信号的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，为什么我用频谱分析缺求不出来呢？