【数字信号】基于matlab轴承故障仿真信号时域波形图+幅度谱图【含Matlab源码 921期】

一、简介

基于matlab轴承故障仿真信号时域波形图+幅度谱图

二、源代码

close all
clear
clc
fs = 20e3;                  % 采样频率
fn = 3e3;                   % 固有频率
y0 = 5;                      % 位移常数
g = 0.1;                     % 阻尼系数
T = 0.01;                   % 重复周期
N = 4096;                  % 采样点数
NT = round(fs*T);      % 单周期采样点数
t = 0:1/fs:(N-1)/fs;      % 采样时刻
t0 = 0:1/fs:(NT-1)/fs;  % 单周期采样时刻
K = ceil(N/NT)+1;       % 重复次数
y = [];
for i = 1:K
    y = [y,y0*exp(-g*2*pi*fn*t0).*sin(2*pi*fn*sqrt(1-g^2)*t0)];
end
y = y(1:N);
Yf = fft(y);                % 频谱
figure(1)
plot(t,y);

三、运行结果

四、备注

版本：2014a