ドキュメンテーション

最新のリリースでは、このページがまだ翻訳されていません。 このページの最新版は英語でご覧になれます。

MATLAB Function ブロックの使用によるクラスの予測

この例では、ラベル予測用に MATLAB Function ブロックを Simulink® に追加する方法を示します。この MATLAB Function ブロックは、ストリーミング データを受け入れて、学習済みのサポート ベクター マシン (SVM) 分類モデルを使用することによりラベルと分類スコアを予測します。MATLAB Function ブロックの使用法の詳細については、MATLAB Function ブロックを使用するモデルの作成 (Simulink)を参照してください。

Statistics and Machine Learning Toolbox™ に含まれている ionosphere データセットには、レーダー反射の品質 (Y) と予測子データ (X) が含まれています。レーダー反射の品質は良好 ('g') または不良 ('b') のいずれかです。

ionosphere データセットを読み込みます。標本サイズを調べます。

load ionosphere
n = numel(Y)
n =

   351

MATLAB Function ブロックで cell 配列を返すことはできません。このため、レーダー反射が良好な場合は要素が 1、それ以外の場合は 0 である logical ベクトルに応答変数を変換します。

Y = strcmp(Y,'g');

レーダー反射は連続的に検出されるものとし、また、はじめの 300 個の観測値を入手しており、残りの 51 個はまだ入手していないとします。現在の標本と将来の標本にデータを分割します。

prsntX = X(1:300,:);
prsntY = Y(1:300);
ftrX = X(301:end,:);
ftrY = Y(301:end);

現在利用できるすべてのデータを使用して SVM モデルに学習をさせます。予測子データを標準化するよう指定します。

Mdl = fitcsvm(prsntX,prsntY,'Standardize',true);

MdlClassificationSVMモデルです。コマンド ラインでは、Mdl を使用して新しい観測値について予測を行うことができます。しかし、コード生成用の関数で入力引数として Mdl を使用することはできません。

saveCompactModelを使用して、関数内に読み込めるように Mdl を準備します。

saveCompactModel(Mdl,'SVMIonosphere');

saveCompactModelMdl をコンパクトにしてから MAT ファイル SVMIonosphere.mat に保存します。

レーダー反射が良好な品質であるかどうかを予測する svmIonospherePredict.m という名前の関数を定義します。この関数は以下を行う必要があります。

  • コード生成命令 %#codegen を関数内に含める。

  • レーダー反射の予測子データを受け入れる。データは、行数以外は X と同等でなければなりません。

  • loadCompactModelを使用して SVMIonosphere.mat を読み込む。

  • 予測されたラベルと、レーダー反射の品質を良好であると予測する分類スコア (つまり、陽性クラス スコア) を返す。

function [label,score] = svmIonospherePredict(X) %#codegen
%svmIonospherePredict Predict radar-return quality using SVM model
%   svmIonospherePredict predicts labels and estimates classification
%   scores of the radar returns in the numeric matrix of preditor data X
%   using the compact SVM model in the file SVMIonosphere.mat.  Rows of X
%   correspond to observations and columns to predictor variables.  label
%   is the predicted label and score is the confidence measure for
%   classifying the radar-return quality as good.
%
% Copyright 2016 The MathWorks Inc.
Mdl = loadCompactModel('SVMIonosphere');
[label,bothscores] = predict(Mdl,X);
score = bothscores(:,2);
end

svmIonospherePredict.mmlr/help/toolbox/stats/examples フォルダーにあります。mlrmatlabroot の値です。

Simulink® モデル slexSVMIonospherePredictExample.slx を読み込みます。

SimMdlName = 'slexSVMIonospherePredictExample';
open_system(SimMdlName);

Simulink® モデルが Figure に表示されます。入力ノードでレーダー反射が検出されると、svmIonospherePredict.m にディスパッチする MATLAB Function ブロックにその観測値が送られます。ラベルとスコアを予測した後で、これらの値がワークスペースに返され、モデル内の値が一度に 1 つずつ表示されます。slexSVMIonospherePredictExample.slx を読み込む際、MATLAB® は radarReturnInput という名前の必要なデータセットも読み込みます。しかし、この例では必要なデータセットを構築する方法を示します。

このモデルでは、次のフィールドが含まれている radarReturnInput という名前の構造体配列を入力データとして受け入れることを想定しています。

  • time - 観測値がモデルに入力された時点。この例では、0 ~ 50 の整数を期間に含めます。方向は予測子データ内の観測値に対応しなければなりません。したがって、この例の場合は time が列ベクトルでなければなりません。

  • signals - values フィールドと dimensions フィールドが含まれている、入力データを説明する 1 行 1 列の構造体配列。values は予測子データの行列、dimensions は予測子変数の個数です。

将来のレーダー反射用に適切な構造体配列を作成します。

radarReturnInput.time = (0:50)';
radarReturnInput.signals(1).values = ftrX;
radarReturnInput.signals(1).dimensions = size(ftrX,2);

ここで名前を radarReturnInput から変更して、新しい名前をモデル内で指定できます。ただし、Simulink® では、前述したフィールド名が構造体配列に含まれていると想定しています。

学習から取り分けておいたデータ、つまり radarReturnInput のデータを使用してモデルをシミュレートします。

sim(SimMdlName);

radarReturnInput 内のすべての観測値が一度に 1 つずつ処理された後で、モデルが Figure に表示されます。X(351,:) について予測されたラベルは 1、陽性クラス スコアは 1.431 です。変数 toutyout および svmlogsout がワークスペースに表示されます。youtsvmlogsout は、予測したラベルおよびスコアが含まれている SimulinkData.Dataset オブジェクトです。詳細は、記録されたシミュレーション データのデータ形式 (Simulink)を参照してください。

シミュレーション データをシミュレーション ログから抽出します。

labelsSL = svmlogsout.getElement(1).Values.Data;
scoresSL = svmlogsout.getElement(2).Values.Data;

labelsSL は、予測したラベルが格納されている 51 行 1 列の数値ベクトルです。labelsSL(j) は、将来の標本におけるレーダー反射 j の品質が良好であると SVM モデルが予測した場合は 1、それ以外の場合は 0 になります。scoresSL は、陽性クラス スコア、つまり判定境界からの符号付き距離が格納されている 51 行 1 列の数値ベクトルです。陽性スコアは予測したラベルの値 1 に、陰性スコアは 0 に対応します。

predictを使用してコマンド ラインでラベルと陽性クラス スコアを予測します。

[labelCMD,scoresCMD] = predict(Mdl,ftrX);
scoresCMD = scoresCMD(:,2);

labelCMDscoresCMDlabelsSL および scoresSL に相応します。

slexSVMIonospherePredictExample が返した将来の標本の陽性クラス スコアを、コマンド ラインで predict を呼び出すことにより得られたスコアと比較します。

err = sum((scoresCMD - scoresSL).^2);
err < eps
ans =

  logical

   1

スコアの集合間における偏差の二乗和は無視できる値です。

Simulink® Coder™ のライセンスもある場合、Simulink® の slexSVMIonospherePredictExample.slx またはコマンド ラインから rtwbuild を使用して C コードを生成できます。詳細は、モデル用 C コードの生成 (Simulink Coder)を参照してください。

参考

| | | |

関連するトピック