用单尺度数学形态学，选择扁平结构元素，处理西储大学轴承故障信号，得到信号的时域图和频域图matlab代码

以下是使用单尺度数学形态学，选择扁平结构元素，处理西储大学轴承故障信号，得到信号的时域图和频域图的 Matlab 代码：

clc;clear;close all;

% 读取数据
load('12k_Drive_End.mat'); % 轴承故障信号

% 绘制原始信号时域图
figure(1);
plot(t, X012_DE_time); % X012_DE_time 为数据中的轴承故障信号
xlabel('Time (s)'); ylabel('Amplitude');
title('Original Signal in Time Domain');

% 执行形态学开运算，选择扁平结构元素
se = strel('rectangle',[7, 1]); % 选择扁平结构元素
X012_DE_time_open = imopen(X012_DE_time, se); % 形态学开运算

% 绘制处理后的信号时域图
figure(2);
plot(t, X012_DE_time_open);
xlabel('Time (s)'); ylabel('Amplitude');
title('Processed Signal in Time Domain');

% 绘制处理后的信号频域图
Fs = 12000; % 采样频率
L = length(X012_DE_time_open);
Y = fft(X012_DE_time_open); % 傅里叶变换
P2 = abs(Y/L); % 双边频谱
P1 = P2(1:L/2+1); % 单边频谱
P1(2:end-1) = 2*P1(2:end-1);
f = Fs*(0:(L/2))/L; % 频率向量

figure(3);
plot(f, P1); 
xlabel('Frequency (Hz)'); ylabel('Amplitude');
title('Processed Signal in Frequency Domain');

其中，数据文件 `12k_Drive_End.mat` 可以在西储大学轴承故障数据集中下载。注意，该代码仅处理了其中的一个轴承故障信号，如果需要处理其他信号，需要对代码进行相应修改。