ブレインマシンインタフェースの閉ループ性能を研究する実験プラットフォーム

1。準備を飛ぶ実験のセットアップの最初のステップはない不随意運動が壊れるように神経細胞の記録の安定性をフライを準備し、ハエの頭は、視覚刺激装置と向きが正しいことを伴います。フライの準備を開始するには、氷上で冷却してから、翼を押したままになまってカクテルスティックを使うと、顕微鏡スライド上に両面テープの部分にフライの背面を固定する。 次に、スライドに翼を添付すると、フライトモータの作用をブロックするように蜜ろうを適用するために電気焼灼針を使用してください。このステップでは、ハエが手続き中にウォームアップしないように、迅速かつ正確な処理が必要です。 今、顕微鏡下で、鉗子で各脚を保持し、身体に最も近い関節にそれらを遮断するための小さなハサミを使う。口吻のためにこれを繰り返します。乾燥からハエを防ぐために、穴はワックスでシールする必要があります。 次に、翼のいずれかを切断し、その側にフライを入れてください。平均棍をカバーする帽子を残したまま、翼の残りの部分を削除し、ワックスで穴を塞ぐ。他の翼についてもこの手順を繰り返します。 定義された方法で、ターゲットのニューロンを刺激するために、ハエの頭には、適切にコンピュータのモニタで整列する必要があります。これを行うには、フライの首が配置されるノッチカットと一方の端にあるフライのボディと付属物のための広いスペースを持っているカスタマイズされたホルダーが必要になります。 接着場所の腹部をしながらそれを押し、ノッチでその首にホルダーの上にフライを置きます。あなたが顕微鏡を通してハエの頭の正面を見ることができるように、今すぐスタンドにフライホルダーを配置。 赤い光でハエの表示、光学現象は擬似瞳孔がそれぞれの目で見ることができると呼ばれる。擬似瞳孔は刺激（フランチェスキーニ1975）とハエの頭を揃えるために使用できる参照フレームを提供します。擬似瞳孔は、以下の画像の挿入によって示される特定の形状を、想定している場合は、ハエの頭の向きが完全に定義されています。 ホルダーに接着剤、それをするためにワックスを使用して正しく配向させるハエの頭をマイクロマニピュレーターを使用して、と。 次に、フラットダウン胸部を押すとホルダーにワックス。これは、電極が、ハエの脳に挿入できるように後頭部のカプセルを開くことができます。 慎重に右側頭部カプセルのキューティクルにウィンドウをカットするマイクロメスまたは微細な注射針を使用してください。キューティクルの下の神経組織の右側をカットしないように注意してください。キューティクルの一部が削除されると、リンガー溶液数滴を追加します。 lobulaプレートをカバーする可能性のある浮遊毛、脂肪沈着、または筋肉組織を削除するには、ピンセットを使用してください。 lobulaプレートは、その後部の表面を覆う銀色の気管の特徴的な分岐パターンによって識別することができます。 参照電極を位置決めするための左後頭部のカプセルのキューティクルに小さな穴を開けます。フライは準備して、記録電極を配置する方法を確認することができます。 2。記録電極の配置フライの準備で、H1のニューロンからの信号を配置し、記録に進んでみましょう。記録電極は、H1のニューロンに近接して配置する必要があります。 H1ニューロンは主に、その受容野（クラップら、2001）に提示水平バックツーフロントの動きに応答します。 記録電極を配置するには、視覚的なランドマークとして気管を使用してください。最初に、一番上の気管の間に電極を配置。 それは、音響信号に記録された電位に変換するオーディオアンプを使用するのに役立ちます。それぞれの個々のスパイクが特徴のクリック音に変換されます。近い電極は、個々のニューロンに到達する、明確なクリック音がなります。 その動きの好みによってH1のニューロンを識別するには、水平方向の動きでそれを刺激する。代わりに、記録電極と、視覚刺激と記録に移動することができます。 3。視覚刺激と記録閉ループ実験では、ターンテーブルの動きのためのロボットの補償でH1ニューロンの結果のように刺激するようにセットアップされます。開始するには、二つのCRTコンピュータモニターの前でハエを置きます。ハエの視覚システムが人間の10倍の速さなので、モニターは、毎秒200フレームを表示する必要があります。モニターや他の電気機器は、電磁的測定神経信号で外部ノイズを最小限に抑えるためにシールドする必要があります。 ハエの方向に対して45度 – + /でモニターの中心を合わせます。フライアイ赤道から分かるように、各モニタは、水平​​で+ / -25度の角度をsubtends、および+垂直面で/ -19度。 コンピュータのモニタに同期した入力が実験のために変更されている小さな、二輪ASUROロボットに搭載された2台のビデオカメラによって提供されます。 その壁垂直方向、黒と白のストライプの模様が並んでいます円筒状の領域内でターンテーブルの上にロボットを配置します。水平面内のターンテーブルを回転させることによって、ロボットの動きは自由の一程度に制限されています。 最初にターンテーブルとロボットの両方が安静時です。ターンテーブルが移動して起動すると、その回転が同じ方向に、ロボットを搬送し、ビデオカメラは、ロボットと競技場の縞模様の間の相対運動を記録する。 45度 – ロボットのバッテリ駆動のビデオカメラは、+ /の向きにマウントされています。彼らは、ハエの目の前にコンピュータのモニタのフレームレートに合わせて、毎秒200の画像をキャプチャ。 640 × 480（グレースケール）の解像度で毎秒200フレームでのコンピュータのモニタに表示される画像を記録します。 フライは、縞模様の動きを見ながら、録画バンドを渡される（例えば、300〜2kHzの）少なくとも10 kHzのサンプリングレートを使用してデジタル画像取り込みボードとの電気信号をフィルタリング。 しきい値が適用される帯域通過するためにバックグラウンドアクティビティからスパイクを分離する電気信号をフィルタリング。因果、ハーフガウシアンフィルタは、H1のセルの円滑なスパイク活動の推定値を得るためにスパイクと畳み込まれている。 ブレインマシンインターフェースのループを閉じるには、制御アルゴリズムは、ロボットの車輪を駆動する2個のDCモータを制御するために、Bluetoothインターフェースを介してフィードバックされるロボットの速度にH1細胞のスパイクレートを変換するために使用されます。 純粋な正弦波は、ターンテーブルの速度プロファイルとして選択されています。正弦波はターンテーブルのみをその好ましい方向に沿ってH1のニューロンを刺激する方向に回転するように- DC相殺されている。ターンテーブルの動きのためのロボットの補償でH1のニューロンの結果の刺激。 図1：閉ループセットアップ 。我々のセットアップでは、左側のH1細胞のスパイク活動は、ターンテーブルに搭載されたロボットの動きを制御するために使用されます。ロボットとターンテーブルとの間の相対運動の結果として生成された視覚的なイメージの動きは、高速度カメラを介してキャプチャされ、フライの前に2つのCRTモニターに表示されます。左半球からH1スパイク活動は、ロボットのための補償の速度を計算するために、制御則を使用してリアルタイムのスパイクレートを推定するために使用されます。ロボットの反回転は、閉ループ制御の間にその場で観測された視覚的な画像の動きを安定させます。 4。代表的な成果と結果正しくセットアップすると、視覚的な安定化は、ロボットの反回転がコンピュータのモニタ上にほとんど、あるいはまったくパターンの動きで、その結果、ターンテーブルの回転と一致したときに達成されます。システムの全体的なパフォーマンスは、ループを閉じるために使用されている制御アルゴリズムに依存します。 我々はテストの最初のアルゴリズムは、更新されたロボットの速度は、ロボット、ωR、そしてターンテーブル、ωP間角速度の差に比例する比例制御（図2）です。ターンテーブルの信号、ωP、のための静的ゲイン、Kpを、と入力周波数の異なる値は、コントローラの性能をテストするために選択されています。 図2：比例制御 。 （a）の好ましい方向にロボットとターンテーブルとの間の相対運動は、スパイクの速度を与えるために、H1のセルを刺激し、F.このスパイクの速度は、速度誤差、Eに変換され、比例制御が推定するために使用されます。更新されたロボットの速度、VR（t +1）を。スパイクの速度、Fは、それがより小さいか自発スパイク率、Fspontより大きいかどうかに基づいて、速度誤差、Eに変換されます。エラーの変換を高速化するためにスパイクの速度は、速度のしきい値の非直線性スパイクを考慮するためにコサイン（区間[π、0]）上のFを投影することによって行われます。定数70と150は、ターンテーブルの最小値と最大角速度にロボットの8ビット入力の速度を一致させるために使用されています。 （b）に比例コントローラを使用して閉ループ系を示すブロック図。システムへの入力は、ターンテーブルの角速度、ωP（T）、および対応するロボットの応答ωRのsinosuidal変調（T +1）が記録されている。 ωpとωrのサンプルトレースは、K、P = 1、ωPは0.6 Hzの入力周波数（図3参照）のためにここに示されています。ロボットは（緑の）遅れでターンテーブルを（青で）続くとピーク振幅より小さい。 H1細胞を刺激パターンの動きの水平成分は、（赤で）を以下に示します。 図3：閉ループ応答 。 （）入力周波数での角度のターンテーブルの速度（青）、ωP、そしてロボット（緑）、ωR、= 0.6 Hzの。 （b）の水平方向の光の流れが記録された画像（赤）から計算として示されています。ピラミッド型のルーカスKanade法（3ピラミッドのレベル）が連続する画像フレーム間の光の流れ場を計算するために使用されます。水平方向の角速度は、水平方向の単位ベクトルiの上に、流れ場内の個々のベクトルの投影を合計して計算されます。 H1セルを阻害する、水平、前面から背面への動きはND（ヌル方向）と呼ばれている間H1細胞を興奮させる水平方向のバックツーフロントの動きは、PD（好ましい方向）と呼ばれています。 ターンテーブルの信号、ωP、のための入力周波数は、0.03から3 Hzと対応するロボットの信号、ωRとの間で選ばれて記録されます。両方の信号を高速フーリエ変換によって周波数領域に変換される（図4参照）と、振幅と位相の値は入力周波数で計算されます。 図4：ターンテーブル、ωP、およびロボット、ωR、 の周波数応答角速度信号が入力周波数での振幅と位相成分を計算するために高速フーリエ変換（FFT）法を用いて周波数領域に変換されています。 FFTの位相成分は、図に示されていません。 K pと比例コントローラの大きさのプロットは= 1に示すテストの入力周波数におけるシステムの応答（図5参照）ボード線図。コントローラの性能は、一般的に増加する周波数で減少する。 1 Hzでわずかに増加の利得は1つだけH1セルが動的に（出力）範囲は主に水平方向のバックツーフロントの動きをカバーを使用のため、ロボットの信号の振動の結果としてです。 図5：比例制御器の性能 。比例制御のために振幅と位相のボード線図（8ハエで平均）、静的ゲインKP = 1.0（）ボード線図の振幅のプロットは、ほぼローパスフィルタ特性に従います。 1 Hzでわずかに増加の利得は1つだけH1セルが動的に（出力）範囲は主にフロントの動きに戻って水平方向のカバーを使用のため、ロボットの信号、ωR、の振動の結果としてです。ロボットの信号、ωRの振動数は、これらの周波数でわずかに増加した利得につながる増加する入力周波数で減少する。 （b）のボード線図、位相プロットは、3 Hzで≤1 Hzとアプローチ不安定入力周波数に対して180 °未満である。特定の入力周波数を超えて、コントローラはロボットの運動学の不安定が原因となります。この不安定性は、ハエの視覚システム（Warzecha ら、1999）の既知の最適な応答の範囲を超えて発生します。 位相のボード線図（図5 – bを参照）入力周波数<0.6 Hzのためのπより小さいコントローラの位相遅れを示しています。コントローラは、周波数に対して安定であることを示して<0.6 Hzと≥1 Hzの入力周波数が不安定。 静的KP（青）と比例コントローラの性能はK pの値が時間間隔で計算したピークスパイク率、F 最高 、に基づいて、すべての50ミリ秒に更新される適応コントローラ（赤で）、と比較した[T – 500ミリ – T]（図6参照）。大規模な統合の時間窓の結果として、比例コントローラは、テストパラメータの範囲のための適応コントローラ（図7参照）よりもパフォーマンスが良くなります。 500ミリ秒の初期の時間積分のウィンドウは、私たちが使用しているロボットのプラットフォームに関連する技術的な理由で選ばれました。適応コントローラは、比例制御（図7 – bを参照）と同様の位相特性を持っていた。 図6：適応型コントローラのゲイン 。適応コントローラは、閉ループ制御の間に継続的にゲインを推定するのに対し、比例制御は、静的ゲイン、Kpを使用しています。ダイナミックゲイン、Kpは、間隔における最大スパイク率、Fmaxと、T – 500ミリ≤τ≤tに反比例する図は、Fmaxのが時間の経過と推定される3つのインスタンスを示しています。 Fmaxのの推定値に基づいて、Kpは時間ウィンドウ3（オレンジ）の間の時間ウィンドウ2（緑）と最低の時に最高です。 図7：比例対適応コントローラ 。 （KP = 1）比例と適応制御のために振幅と位相のボード線図（a）はKpの値ごとに50msの適応コントローラの更新は最後の500ミリ秒以上のピークスパイクレート推定Fmaxとに基づいて。比例制御のゲインのプロット（青）は、すべての入力周波数でより良い実行していることを示す適応コントローラ（赤）よりも高いです。 （b）の両方のコントローラのボード線位相プロットは、f = 0.3 Hzの有意差と似ています。両方のコントローラは、3 Hzで不安定に近づく。大幅に異なるゲインと位相の値は、アスタリスク（ウィルコクソンの順位和法、P = 0.05）で示されます。 ターンテーブル周りの格子パターンを除去し、ラボ環境は、フライH1 -セルのための自然の視覚入力の近似値として使用されました。平均では、自然主義的視覚的なイメージの空間周波数の広い範囲のため、おそらく格子の視覚入力（赤で）でのものよりわずかに高いゲインを示した自然な視覚入力（青色の）のために大きさのプロットを（図8参照）ボード線図活用されている。回折格子対自然な視覚入力のためのボード線図の位相プロットの特性は（図8 – bを参照）と同様であった。 図8：パターン対のラボ環境を剥奪 。閉ループ下で取り除いたパターン（赤）VSラボ環境（青）映像が表示されたら比例制御のための振幅と位相のボード線図（）のラボ環境の画像が使用されているときの大きさのプロット時に最低限のパターンよりもわずかに高くなってボード線図このような刺激の下で良いパフォーマンスを示す（F = 0.1 Hzのを除いて）使用されています。 （b）両方の視覚的な条件下での位相のプロットは3 Hzで近づいて不安定の両方で、同じパターンに従ってボード線図。大幅に異なるゲインと位相の値は、アスタリスク（ウィルコクソンの順位和法、P = 0.05）で示されます。