В этом примере рассмотрим данные о времени жизни 10 идентичных машин со следующими 6 потенциальными прогностическими $$$$параметрами - постоянными, линейными, квадратичными, кубическими, логарифмическими и периодическими. Набор данных machineDataCellArray.mat содержит C который является 1x10 массив ячеек матриц, где каждый элемент массива ячеек является матрицей, которая содержит данные времени жизни машины. Для каждой матрицы в массиве ячеек первый столбец содержит время, в то время как другие столбцы содержат переменные данных.

Загрузите данные о времени жизни и визуализируйте их по времени.

load('machineDataCellArray.mat','C')
display(C)

C=1×10 cell array
  Columns 1 through 4

    {219x7 double}    {189x7 double}    {202x7 double}    {199x7 double}

  Columns 5 through 8

    {229x7 double}    {184x7 double}    {224x7 double}    {208x7 double}

  Columns 9 through 10

    {181x7 double}    {197x7 double}

for k = 1:length(C)
    plot(C{k}(:,1), C{k}(:,2:end));
    hold on;
end

Просмотрите 6 различных индикаторов состояния - постоянный, линейный, квадратичный, кубический, логарифмический и периодический - для всех 10 машин на графике.

Визуализируйте модность потенциальных прогностических особенностей.

trendability(C)

На графике гистограммы обратите внимание, что элементы Var2 и Var5 имеют значения изменчивости 1. Следовательно, эти особенности более подходят для оставшихся прогнозов срока службы, поскольку они являются лучшими показателями работоспособности машины.

В этом примере рассмотрим данные о времени жизни 10 идентичных машин со следующими 6 потенциальными прогностическими $$$$параметрами - постоянными, линейными, квадратичными, кубическими, логарифмическими и периодическими. Набор данных machineDataTable.mat содержит T, который является 1x10 массив ячеек таблиц, где каждый элемент массива ячеек содержит таблицу данных о времени жизни для машины.

Загрузите и просмотрите данные.

load('machineDataTable.mat','T');
display(T)

T=1×10 cell array
  Columns 1 through 4

    {219x7 table}    {189x7 table}    {202x7 table}    {199x7 table}

  Columns 5 through 8

    {229x7 table}    {184x7 table}    {224x7 table}    {208x7 table}

  Columns 9 through 10

    {181x7 table}    {197x7 table}

head(T{1},2)

ans=2×7 table
    Time    Constant    Linear    Quadratic    Cubic     Logarithmic    Periodic
    ____    ________    ______    _________    ______    ___________    ________

       0     3.2029     11.203     7.7029      3.8829      2.2517        0.2029 
    0.05     2.8135     10.763     7.2637      3.6006      1.8579       0.12251 

Обратите внимание, что каждая таблица в массиве ячеек содержит переменную времени жизни 'Time' и переменные данных 'Constant', 'Linear', 'Quadratic', 'Cubic', 'Logarithmic', и 'Periodic'.

Вычислить модность с помощью Time в качестве переменной времени жизни.

Y = trendability(T,'Time')

Y=1×6 table
    Constant     Linear     Quadratic     Cubic     Logarithmic    Periodic 
    _________    _______    _________    _______    ___________    _________

    0.0035529    0.99984     0.63753     0.92057      0.99582      0.0041995

Из результирующей таблицы значений изменчивости следует, что линейные, кубические и логарифмические элементы имеют значения, близкие к 1. Следовательно, эти три признака более подходят для прогнозирования оставшегося срока службы, поскольку они являются лучшими показателями работоспособности машины.

Рассмотрим данные о сроке службы 4 машин. Каждая машина имеет 4 кода отказа для индикаторов потенциального состояния - $$$$напряжение, ток и мощность. trendabilityEnsemble.zip представляет собой коллекцию из 4 файлов, где каждый файл содержит расписание времени жизни данных для каждой машины - tbl1.mat, tbl2.mat, tbl3.mat и tbl4.mat. Можно также использовать файлы, содержащие данные для нескольких машин. Для каждого графика организация данных выглядит следующим образом:

При выполнении вычислений в массивах TALL MATLAB ® использует либо параллельный пул (по умолчанию при наличии Toolbox™ Parallel Computing), либо локальный сеанс MATLAB. Для выполнения примера с использованием локального сеанса MATLAB измените глобальную среду выполнения с помощью mapreducer функция.

mapreducer(0)

Извлеките сжатые файлы, прочитайте данные в расписаниях и создайте fileEnsembleDatastore с использованием данных расписания. Дополнительные сведения о создании хранилища данных ансамбля файлов см. в разделе fileEnsembleDatastore.

unzip trendabilityEnsemble.zip;
ens = fileEnsembleDatastore(pwd,'.mat');
ens.DataVariables = {'Voltage','Current','Power','FaultCode','Machine'};
% Make sure that the function for reading data is on path
addpath(fullfile(matlabroot,'examples','predmaint','main')) 
ens.ReadFcn = @readtable_data;
ens.SelectedVariables = {'Voltage','Current','Power','FaultCode','Machine'};

Визуализировать модность потенциальных прогностических особенностей с помощью 'Machine' в качестве переменной-члена и группировать данные о сроке службы по 'FaultCode'. Группирование данных о сроке службы гарантирует, что trendability вычисляет метрику для каждого кода неисправности отдельно.

trendability(ens,'MemberVariable','Machine','GroupBy','FaultCode');

Evaluating tall expression using the Local MATLAB Session:
- Pass 1 of 1: Completed in 0.083 sec
Evaluation completed in 0.2 sec
Evaluating tall expression using the Local MATLAB Session:
- Pass 1 of 1: Completed in 0.037 sec
Evaluation completed in 0.14 sec
Evaluating tall expression using the Local MATLAB Session:
- Pass 1 of 1: Completed in 0.13 sec
Evaluation completed in 0.15 sec

trendability возвращает график гистограммы с элементами, ранжированными по значениям их модности. Более высокое значение модности указывает на более подходящий прогностический параметр. Например, элемент-кандидат Current имеет наивысшую степень модности для машин с FaultCode 1.

rmpath(fullfile(matlabroot,'examples','predmaint','main')) % Reset path