clc
clear
fs = 20e3; % 采样频率
fn = 3e3; % 固有频率
y0 = 3; % 位移常数
g = 0.1; % 阻尼系数
T = 0.02; % 重复周期
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);
tfrgabor(y',128,32);
end
%% 滚动轴承故障稳态信号
clc
clear
fs=20e3; %采样频率
fn=3e3; %固有频率
y0=3; %位移常数
g=0.1; %阻尼系数
T=0.02; %重复周期
N=14200; %采样点数
NT=round(fs*T); %单周期采样点数
t=0:1/fs:(N-1)/fs; %采样时刻
t0=0:1/fs:(NT-1)/fs; %单周期采样时刻
K=ceil(N/NT)+1; %重复次数
tau=0.02*T; %随机波动
y=[];
for i=1:K
y=[y,y0*exp(-g*2*pi*fn*(t0-tau)).*sin(2*pi*fn*(t0-tau)).*stepfun((t0-tau),0)]; %故障信号
end
y=y(1:N);
y1=awgn(y,-10); %加噪声
subplot 211 %输出冲击波形
plot(t,y);
title('原始信号'),xlabel('时间(t/s)'),ylabel('幅值');
subplot 212 %输出带有噪声的冲击波形
plot(t,y1);
title('原始信号+噪声'),xlabel('时间(t/s)'),ylabel('幅值');
%% 滚动轴承故障稳态信号
clc
clear
fs=20e3; %采样频率
fn=3e3; %固有频率
y0=3; %位移常数
g=0.1; %阻尼系数
T=0.02; %重复周期
N=14200; %采样点数
NT=round(fs*T); %单周期采样点数
t=0:1/fs:(N-1)/fs; %采样时刻
t0=0:1/fs:(NT-1)/fs; %单周期采样时刻
K=ceil(N/NT)+1; %重复次数
tau=0.02*T; %随机波动
y=[];
for i=1:K
y=[y,y0*exp(-g*2*pi*fn*(t0-tau)).*sin(2*pi*fn*(t0-tau)).*stepfun((t0-tau),0)]; %故障信号
end
y=y(1:N);
y1=awgn(y,-10); %加噪声
subplot 211 %输出冲击波形
plot(t,y);
title('原始信号'),xlabel('时间(t/s)'),ylabel('幅值');
subplot 212 %输出带有噪声的冲击波形
plot(t,y1);
title('原始信号+噪声'),xlabel('时间(t/s)'),ylabel('幅值');
%% 滚动轴承故障稳态信号
clc
clear
fs=20e3; %采样频率
fn=3e3; %固有频率
y0=3; %位移常数
g=0.1; %阻尼系数
T=0.02; %重复周期
N=14200; %采样点数
NT=round(fs*T); %单周期采样点数
t=0:1/fs:(N-1)/fs; %采样时刻
t0=0:1/fs:(NT-1)/fs; %单周期采样时刻
K=ceil(N/NT)+1; %重复次数
tau=0.02*T; %随机波动
y=[];
for i=1:K
y=[y,y0*exp(-g*2*pi*fn*(t0-tau)).*sin(2*pi*fn*(t0-tau)).*stepfun((t0-tau),0)]; %故障信号
end
y=y(1:N);
y1=awgn(y,-10); %加噪声
subplot 211 %输出冲击波形
plot(t,y);
title('原始信号'),xlabel('时间(t/s)'),ylabel('幅值');
subplot 212 %输出带有噪声的冲击波形
plot(t,y1);
title('原始信号+噪声'),xlabel('时间(t/s)'),ylabel('幅值');
%% 滚动轴承故障稳态信号
clc
clear
fs=20e3; %采样频率
fn=3e3; %固有频率
y0=3; %位移常数
g=0.1; %阻尼系数
T=0.02; %重复周期
N=14200; %采样点数
NT=round(fs*T); %单周期采样点数
t=0:1/fs:(N-1)/fs; %采样时刻
t0=0:1/fs:(NT-1)/fs; %单周期采样时刻
K=ceil(N/NT)+1; %重复次数
tau=0.02*T; %随机波动
y=[];
for i=1:K
y=[y,y0*exp(-g*2*pi*fn*(t0-tau)).*sin(2*pi*fn*(t0-tau)).*stepfun((t0-tau),0)]; %故障信号
end
y=y(1:N);
y1=awgn(y,-10); %加噪声
subplot 211 %输出冲击波形
plot(t,y);
title('原始信号'),xlabel('时间(t/s)'),ylabel('幅值');
subplot 212 %输出带有噪声的冲击波形
plot(t,y1);
title('原始信号+噪声'),xlabel('时间(t/s)'),ylabel('幅值');
%% 滚动轴承故障稳态信号
clc
clear
fs=20e3; %采样频率
fn=3e3; %固有频率
y0=3; %位移常数
g=0.1; %阻尼系数
T=0.02; %重复周期
N=14200; %采样点数
NT=round(fs*T); %单周期采样点数
t=0:1/fs:(N-1)/fs; %采样时刻
t0=0:1/fs:(NT-1)/fs; %单周期采样时刻
K=ceil(N/NT)+1; %重复次数
tau=0.02*T; %随机波动
y=[];
for i=1:K
y=[y,y0*exp(-g*2*pi*fn*(t0-tau)).*sin(2*pi*fn*(t0-tau)).*stepfun((t0-tau),0)]; %故障信号
end
y=y(1:N);
y1=awgn(y,-10); %加噪声
subplot 211 %输出冲击波形
plot(t,y);
title('原始信号'),xlabel('时间(t/s)'),ylabel('幅值');
subplot 212 %输出带有噪声的冲击波形
plot(t,y1);
title('原始信号+噪声'),xlabel('时间(t/s)'),ylabel('幅值');
%% 滚动轴承故障稳态信号
clc
clear
fs=20e3; %采样频率
fn=3e3; %固有频率
y0=3; %位移常数
g=0.1; %阻尼系数
T=0.02; %重复周期
N=14200; %采样点数
NT=round(fs*T); %单周期采样点数
t=0:1/fs:(N-1)/fs; %采样时刻
t0=0:1/fs:(NT-1)/fs; %单周期采样时刻
K=ceil(N/NT)+1; %重复次数
tau=0.02*T; %随机波动
y=[];
for i=1:K
y=[y,y0*exp(-g*2*pi*fn*(t0-tau)).*sin(2*pi*fn*(t0-tau)).*stepfun((t0-tau),0)]; %故障信号
end
y=y(1:N);
y1=awgn(y,-10); %加噪声
subplot 211 %输出冲击波形
plot(t,y);
title('原始信号'),xlabel('时间(t/s)'),ylabel('幅值');
subplot 212 %输出带有噪声的冲击波形
plot(t,y1);
title('原始信号+噪声'),xlabel('时间(t/s)'),ylabel('幅值');
%% 滚动轴承故障稳态信号
clc
clear
fs=20e3; %采样频率
fn=3e3; %固有频率
y0=3; %位移常数
g=0.1; %阻尼系数
T=0.02; %重复周期
N=14200; %采样点数
NT=round(fs*T); %单周期采样点数
t=0:1/fs:(N-1)/fs; %采样时刻
t0=0:1/fs:(NT-1)/fs; %单周期采样时刻
K=ceil(N/NT)+1; %重复次数
tau=0.02*T; %随机波动
y=[];
for i=1:K
y=[y,y0*exp(-g*2*pi*fn*(t0-tau)).*sin(2*pi*fn*(t0-tau)).*stepfun((t0-tau),0)]; %故障信号
end
y=y(1:N);
y1=awgn(y,-10); %加噪声
subplot 211 %输出冲击波形
plot(t,y);
title('原始信号'),xlabel('时间(t/s)'),ylabel('幅值');
subplot 212 %输出带有噪声的冲击波形
plot(t,y1);
title('原始信号+噪声'),xlabel('时间(t/s)'),ylabel('幅值');
%% 滚动轴承故障稳态信号
clc
clear
fs=20e3; %采样频率
fn=3e3; %固有频率
y0=3; %位移常数
g=0.1; %阻尼系数
T=0.02; %重复周期
N=14200; %采样点数
NT=round(fs*T); %单周期采样点数
t=0:1/fs:(N-1)/fs; %采样时刻
t0=0:1/fs:(NT-1)/fs; %单周期采样时刻
K=ceil(N/NT)+1; %重复次数
tau=0.02*T; %随机波动
y=[];
for i=1:K
y=[y,y0*exp(-g*2*pi*fn*(t0-tau)).*sin(2*pi*fn*(t0-tau)).*stepfun((t0-tau),0)]; %故障信号
end
y=y(1:N);
y1=awgn(y,-10); %加噪声
subplot 211 %输出冲击波形
plot(t,y);
title('原始信号'),xlabel('时间(t/s)'),ylabel('幅值');
subplot 212 %输出带有噪声的冲击波形
plot(t,y1);
title('原始信号+噪声'),xlabel('时间(t/s)'),ylabel('幅值');
%% 滚动轴承故障稳态信号
clc
clear
fs=20e3; %采样频率
fn=3e3; %固有频率
y0=3; %位移常数
g=0.1; %阻尼系数
T=0.02; %重复周期
N=14200; %采样点数
NT=round(fs*T); %单周期采样点数
t=0:1/fs:(N-1)/fs; %采样时刻
t0=0:1/fs:(NT-1)/fs; %单周期采样时刻
K=ceil(N/NT)+1; %重复次数
tau=0.02*T; %随机波动
y=[];
for i=1:K
y=[y,y0*exp(-g*2*pi*fn*(t0-tau)).*sin(2*pi*fn*(t0-tau)).*stepfun((t0-tau),0)]; %故障信号
end
y=y(1:N);
y1=awgn(y,-10); %加噪声
subplot 211 %输出冲击波形
plot(t,y);
title('原始信号'),xlabel('时间(t/s)'),ylabel('幅值');
subplot 212 %输出带有噪声的冲击波形
plot(t,y1);
title('原始信号+噪声'),xlabel('时间(t/s)'),ylabel('幅值');
%% 滚动轴承故障稳态信号
clc
clear
fs=20e3; %采样频率
fn=3e3; %固有频率
y0=3; %位移常数
g=0.1; %阻尼系数
T=0.02; %重复周期
N=14200; %采样点数
NT=round(fs*T); %单周期采样点数
t=0:1/fs:(N-1)/fs; %采样时刻
t0=0:1/fs:(NT-1)/fs; %单周期采样时刻
K=ceil(N/NT)+1; %重复次数
tau=0.02*T; %随机波动
y=[];
for i=1:K
y=[y,y0*exp(-g*2*pi*fn*(t0-tau)).*sin(2*pi*fn*(t0-tau)).*stepfun((t0-tau),0)]; %故障信号
end
y=y(1:N);
y1=awgn(y,-10); %加噪声
subplot 211 %输出冲击波形
plot(t,y);
title('原始信号'),xlabel('时间(t/s)'),ylabel('幅值');
subplot 212 %输出带有噪声的冲击波形
plot(t,y1);
title('原始信号+噪声'),xlabel('时间(t/s)'),ylabel('幅值');
%% 滚动轴承故障稳态信号
clc
clear
fs=20e3; %采样频率
fn=3e3; %固有频率
y0=3; %位移常数
g=0.1; %阻尼系数
T=0.02; %重复周期
N=14200; %采样点数
NT=round(fs*T); %单周期采样点数
t=0:1/fs:(N-1)/fs; %采样时刻
t0=0:1/fs:(NT-1)/fs; %单周期采样时刻
K=ceil(N/NT)+1; %重复次数
tau=0.02*T; %随机波动
y=[];
for i=1:K
y=[y,y0*exp(-g*2*pi*fn*(t0-tau)).*sin(2*pi*fn*(t0-tau)).*stepfun((t0-tau),0)]; %故障信号
end
y=y(1:N);
y1=awgn(y,-10); %加噪声
subplot 211 %输出冲击波形
plot(t,y);
title('原始信号'),xlabel('时间(t/s)'),ylabel('幅值');
subplot 212 %输出带有噪声的冲击波形
plot(t,y1);
example2_11.zip_冲击仿真信号_冲击信号仿真_周期性冲击_振动冲击信号_轴承仿真信号
版权申诉
5星 · 超过95%的资源 160 浏览量
2022-07-14
04:16:52
上传
评论 1
收藏 1KB ZIP 举报
四散
- 粉丝: 54
- 资源: 1万+
评论1