浮動小数点から固定小数点への変換

この例では、MATLAB® で浮動小数点設計を開始し、MATLAB で効率的な固定小数点設計に収束するまで反復処理してから、生成された固定小数点設計の数値精度を検証する方法を示します。

再構成可能なプラットフォーム用の信号処理アプリケーションには、一般に浮動小数点演算を使用して指定されるアルゴリズムが必要です。しかし、電力、コスト、パフォーマンスの理由から、DSP コアについてはソフトウェアで、FPGA では特殊用途のハードウェアとして、通常は固定小数点演算で実装されます。固定小数点の変換は非常に難しく、時間がかかることがあり、一般に設計と実装にかかる合計時間の 25 ～ 50 パーセントを占めます。自動化ツールを使用すると、変換プロセスが簡単になり、時間を短縮できます。

ソフトウェア実装の場合、与えられた計算精度の制約の中でコードサイズと実行時間を最小化する、最適化された固定小数点の仕様を定義することが目標になります。この最適化は、ターゲットプロセッサでサポートされるそれぞれのデータ型に応じて、2 進小数点の位置の変更 (スケーリング用) とデータの語長の選択を行うことで達成されます。

ハードウェア実装の場合、アーキテクチャ全体を最適化できます。実装が効率的であれば、使用面積と電力消費の両方が最小化されます。そのため、変換プロセスの目標は、通常は演算子の語長を最小化することが中心になります。

浮動小数点から固定小数点へのワークフローは、MATLAB から HDL へのワークフロー入門で説明されているように、現在は HDL ワークフローアドバイザーに統合されています。

はじめに

HDL Coder™ の浮動小数点から固定小数点への変換ワークフローには次の手順が含まれます。

浮動小数点設計がコード生成と互換性があることを確認する。
テストベンチのシミュレーションに基づいて固定小数点型を計算する。
推奨された型を適用して、読み取りおよび追跡が可能な固定小数点 MATLAB コードを生成する。
生成された固定小数点設計を検証する。
生成された固定小数点コードの数値精度を元の浮動小数点コードと比較する。

MATLAB 設計

この例で使用する MATLAB コードはシンプルな 2 次直接型 2 転置フィルターです。この例には、このフィルターを使用する MATLAB テストベンチも含まれています。

design_name = 'mlhdlc_df2t_filter';
testbench_name = 'mlhdlc_df2t_filter_tb';

MATLAB 設計を調べます。

type(design_name);

%#codegen
function y = mlhdlc_df2t_filter(x)

%   Copyright 2011-2015 The MathWorks, Inc.

persistent z;
if isempty(z)
    % Filter states as a column vector
    z = zeros(2,1);
end

% Filter coefficients as constants
b = [0.29290771484375   0.585784912109375  0.292907714843750];
a = [1.0                0.0                0.171600341796875];

y    =  b(1)*x + z(1);
z(1) = (b(2)*x + z(2)) - a(2) * y;
z(2) =  b(3)*x - a(3) * y;

end

浮動小数点から固定小数点へのワークフローでは、完全なテストベンチを用意することが推奨されます。推奨された固定小数点データ型の品質は、テストベンチが設計のダイナミックレンジを適切な精度でどれだけカバーできるかによって異なります。

浮動小数点設計とテストベンチの要件の詳細については、生成された固定小数点ファイルの操作の「浮動小数点設計の構造」の節を参照してください。

type(testbench_name);


%

%   Copyright 2011-2015 The MathWorks, Inc.


Fs = 256;              % Sampling frequency
Ts = 1/Fs;             % Sample time
t = 0:Ts:1-Ts;         % Time vector from 0 to 1 second
f1 = Fs/2;             % Target frequency of chirp set to Nyquist
in = sin(pi*f1*t.^2);  % Linear chirp from 0 to Fs/2 Hz in 1 second
out = zeros(size(in)); % Output the same size as the input

for ii=1:length(in)
    out(ii) = mlhdlc_df2t_filter(in(ii));
end

% Plot
figure('Name', [mfilename, '_plot']);
subplot(2,1,1);
plot(in);
xlabel('Time')
ylabel('Amplitude')
title('Input Signal (with Noise)')

subplot(2,1,2);
plot(out);
xlabel('Time')
ylabel('Amplitude')
title('Output Signal (filtered)')

設計のシミュレーション

コードの生成前にテストベンチを使用して設計をシミュレートし、実行時エラーがないことを確認します。

mlhdlc_df2t_filter_tb

HDL Coder プロジェクトの新規作成

新しいプロジェクトを作成するには、次のコマンドを入力します。

coder -hdlcoder -new flt2fix_project

次に、'mlhdlc_filter.m' ファイルを [MATLAB 関数] としてプロジェクトに追加し、'mlhdlc_filter_tb.m' を [MATLAB テストベンチ] として追加します。

MATLAB HDL Coder プロジェクトの作成と入力に関する詳細なチュートリアルについては、MATLAB から HDL へのワークフロー入門を参照してください。

固定小数点コード生成ワークフロー

浮動小数点から固定小数点への変換ワークフローでは以下を実行できます。

浮動小数点設計がコード生成に適合していることを確認する
シミュレーションデータと語長の設定に基づいて固定小数点型を推奨する
推奨された固定小数点型をユーザーが手動で調整できるように許可する
推奨された固定小数点型を検証する
生成された固定小数点 MATLAB コードが目的の数値精度になっていることを確認する

手順 1: ワークフローアドバイザーを起動します。

[ワークフローアドバイザー] ボタンをクリックして HDL ワークフローアドバイザーを起動します。
[固定小数点の変換] オプションで [ビルド時に固定小数点に変換] を選択します。

手順 2: 入力タイプの定義

この手順では、入力の型を手動で定義するか、テストベンチを指定して実行することで定義できます。

[実行] をクリックして、この手順を実行します。

シミュレーション後に確認すると、入力変数 x が double 型のスカラー double(1x1) として定義されていることがわかります。

手順 3: シミュレーションの実行

[固定小数点の変換] の手順をクリックします。

前の手順で定義した入力の型で設計がコンパイルされ、コンパイルが正常に完了すると、設計内のすべての関数に対する推定された型が変数テーブルに表示されます。

この手順では、シミュレーション中に設計内のすべての変数の最小値と最大値が収集されるように、元の設計がインストルメント化されます。

[解析] ボタンをクリックします。

テーブルの [シミュレーションの最小値] と [シミュレーションの最大値] にシミュレーション範囲が取り込まれていることがわかります。固定小数点型は、既定の語長の設定に基づいて推奨されます。

この段階で、すべての変数に対する計算されたシミュレーション範囲に基づいて以下を計算できます。

語長を固定の設定としたときの対応する小数部の長さ。
小数部の長さを固定の設定としたときの対応する語長。

型のテーブルには、浮動小数点 MATLAB 設計に存在するそれぞれの変数について、以下の情報が関数別にまとめて表示されます。

シミュレーションの最小値: シミュレーション中に変数に代入される最小値。
シミュレーションの最大値: シミュレーション中に変数に代入される最大値。
0 または正の整数: シミュレーション中に代入される値がすべて整数かどうか。

型の推奨手順では、上記の情報をユーザーが指定した語長の設定と組み合わせて使用して、各変数の固定小数点型を推奨します。

[解析] ボタンのドロップダウンメニューで [ヒストグラムデータを記録] オプションを有効にして、ヒストグラムデータの記録を有効にすることもできます。

ヒストグラムの表示では、変数に対するシミュレーションデータのダイナミックレンジに関する情報が簡潔に示されます。x 軸はビットの重みに対応し、y 軸は出現回数を表します。このグラフに推奨された数値型の情報が重ねて表示され、ここで情報を編集できます。白い境界ボックスを左右に動かすと、2 進小数点の位置が変わります。左右の端を動かすと、それに応じて小数部の長さまたは語長が変わります。推奨された型に対する変更は、プロジェクトにすべて保存されます。