In Diagnostic Feature Designer, you explore features interactively, using tools for signal processing, feature generation, and ranking. Once you select the set of the features that perform best, you can generate a MATLAB® function that reproduces the calculations for those features. You can apply this function directly to a larger set of measurement data and increase the number of members in your feature set. You can also modify the function to suit your application, and incorporate part or all of the function into other code.
This example shows how to generate a MATLAB function to compute a set of features, and how to validate that function with the original data set.
The example assumes that you are familiar with ensemble data concepts and with basic operations in the app, such as data import, signal processing, and feature generation. For more information on these concepts and operations, see Identify Condition Indicators for Predictive Maintenance Algorithm Design.
This example uses ensemble data generated from a transmission system model in Using Simulink to Generate Fault Data. Outputs of the model include:
Vibration measurements from a sensor monitoring casing vibrations
Tachometer sensor, which issues a pulse every time the shaft completes a rotation
Fault code indicating the presence of a modeled fault
In your MATLAB command window, load the transmission data, which is stored in the table
load dfd_Tutorial dataTable
dataTableis an ensemble table that contains 16 members, each of which represents one simulated transmission system. Each row of the table corresponds to one member. Each column of the table corresponds to one variable, such as
FaultCode. All ensemble members have the same variables.
In the app, import
dataTable as a multimember ensemble with the
faultCode. When import is complete, the signal browser
displays the vibration and tacho data, and the dataset browser displays your imported
ensemble. For information on the import process, see Import and Visualize Ensemble Data in Diagnostic Feature Designer.
Compute a time-synchronous average (TSA) signal from your vibration and tacho signals. To do so, in the Feature Designer tab, select Filtering & Averaging > Time-Synchronous Signal Averaging. Set the parameters as shown in the following figure and click OK.
The new signal appears in the app data browser.
For information on TSA signals, see
In the Feature Designer tab, select Time-Domain Features > Signal Features to open the set of available signal features. Select the features for mean, standard deviation, and kurtosis.
View the feature values. In the data browser, select
FeatureTable1. Then, in the plot gallery, click
Feature Table View. These steps open a table containing the feature
values for each member along with the condition variable
Generate a MATLAB function that reproduces the calculations for these features. In the Feature Designer tab, select Export > Generate Function for Features.
Your selection opens a dialog box that allows you to specify the feature table
and the features. The default options
Use All Features are your only choices in this case, as you
have only one feature table and have not performed any ranking.
When you click OK, a function script opens in the MATLAB editor that begins with the following lines.
function [featureTable,outputTable] = diagnosticFeatures(inputData) %DIAGNOSTICFEATURES recreates results in Diagnostic Feature Designer. % % Input: % inputData: A table or a cell array of tables/matrices containing the % data as those imported into the app. % % Output: % featureTable: A table containing all features and condition variables. % outputTable: A table containing the computation results. % % This function computes signals: % Vibration_tsa/Data % % This function computes features: % Vibration_tsa_stats/Data_Kurtosis % Vibration_tsa_stats/Data_Mean % Vibration_tsa_stats/Data_Std % % Organization of the function: % 1. Compute signals/spectra/features % 2. Extract computed features into a table
For more information on the code content, see Anatomy of App-Generated MATLAB Code.
Run the function using
dataTable to create a new feature table
featuretable = diagnosticFeatures(dataTable)
Compare the first eight feature values to the corresponding feature values in the app. At the level of the displayed precision, the values are identical.
16×4 table faultCode Vibration_tsa_stats/Data_Kurtosis Vibration_tsa_stats/Data_Mean Vibration_tsa_stats/Data_Std _________ _________________________________ _____________________________ ____________________________ 0 2.2516 0.022125 0.99955 1 2.2526 -0.027311 0.999 1 2.2571 -0.45475 0.99629 1 2.2526 0.47419 0.999 1 2.2529 0.37326 0.999 1 2.2526 -0.14185 0.999 1 2.2529 0.40644 0.999 1 2.2529 -0.47485 0.99915