7720 歯のピニオンギヤで駆動する 2 MW 風力タービン高速シャフトから収集されたデータを使用します[ 1] 。まず、[ https://github.com/mathworks/WindTurbineHighSpeedBearingPrognosis-Data ] ( https://github.com/mathworks/WindTurbineHighSpeedBearingPrognosis-Data ) からリポジトリ全体を zip ファイルとしてダウンロードして、本スクリプトと同じディレクトリに保存してください。以下のコマンドを使用してファイルを解凍します。こちらのデータでは計測間隔は 1 日間隔です。
88
99``` matlab:Code
10- if exist('WindTurbineHighSpeedBearingPrognosis-Data-master .zip', 'file')
11- unzip('WindTurbineHighSpeedBearingPrognosis-Data-master .zip')
10+ if exist('WindTurbineHighSpeedBearingPrognosis-Data-main .zip', 'file')
11+ unzip('WindTurbineHighSpeedBearingPrognosis-Data-main .zip')
1212end
1313```
1414
1717まず風力タービンデータに対して fileDatastore を作成します。使用するデータには振動とタコメータの信号が含まれています。各 mat ファイルからcutomreader 関数で読み取りますが、1 秒分の振動データに加えてファイル名から日付情報も合わせて取得します。
1818
1919``` matlab:Code
20- dir = 'WindTurbineHighSpeedBearingPrognosis-Data-master ';
20+ dir = 'WindTurbineHighSpeedBearingPrognosis-Data-main ';
2121ds = fileDatastore(fullfile('.', dir,'*.mat'),'ReadFcn',@customreader,'UniformRead',true);
2222data = readall(ds)
2323```
@@ -34,6 +34,10 @@ data = readall(ds)
3434| 8| 2013/03/14| 97656x1 timetable|
3535| 9| 2013/03/15| 97656x1 timetable|
3636| 10| 2013/03/16| 97656x1 timetable|
37+ | 11| 2013/03/17| 97656x1 timetable|
38+ | 12| 2013/03/17| 97656x1 timetable|
39+ | 13| 2013/03/18| 97656x1 timetable|
40+ | 14| 2013/03/20| 97656x1 timetable|
3741
3842# Data Reduction
3943
@@ -80,11 +84,11 @@ vib = mfile.vibration;
8084
8185% Extract the first 1s of signal and convert it to timetable
8286fs = 97656; % Sampling rate
83- tVibration = timetable(vib(1:fs),'SamplingRate ',fs);
87+ tVibration = timetable(vib(1:fs),'SampleRate ',fs);
8488data.vibration = {tVibration};
8589
8690end
8791
8892```
8993
90- * Copyright 2018 The MathWorks, Inc.*
94+ * Copyright 2018-2022 The MathWorks, Inc.*
0 commit comments