嗜好判別装置、嗜好判別方法、嗜好判別プログラム、及び脳波解析方法
【課題】 被験者の負担にならない簡易脳波計を用いて脳波を測定することで、被験者の嗜好および／または生理的状態を簡単に判別することができる嗜好判別装置、嗜好判別方法、嗜好判別プログラム、及び脳波解析方法を提供する。
【解決手段】 刺激発生部７は、被験者が好む匂い刺激を発生する。簡易脳波計２は、装着した被験者の脳波を測定する。処理部５内のデータ行列作成部５１は、測定した脳波データについて高速フーリエ変換による周波数解析を行う。データ行列作成部５１は、高速フーリエ変換を行った脳波データから、行が時間、列が周波数成分のデータ行列を作成する。特徴値抽出部５２は、作成したデータ行列から、因子分析により特徴値の抽出を行う。制御部５３は、抽出した特徴値について、データベース６にあらかじめ格納されている、分類テーブルを参照して、好む匂いか嫌いな匂いかの嗜好を判定する。
【解決手段】 刺激発生部７は、被験者が好む匂い刺激を発生する。簡易脳波計２は、装着した被験者の脳波を測定する。処理部５内のデータ行列作成部５１は、測定した脳波データについて高速フーリエ変換による周波数解析を行う。データ行列作成部５１は、高速フーリエ変換を行った脳波データから、行が時間、列が周波数成分のデータ行列を作成する。特徴値抽出部５２は、作成したデータ行列から、因子分析により特徴値の抽出を行う。制御部５３は、抽出した特徴値について、データベース６にあらかじめ格納されている、分類テーブルを参照して、好む匂いか嫌いな匂いかの嗜好を判定する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、脳波から人間の生理的状態および／または刺激に対する人間の嗜好を判別する嗜好判別装置、嗜好判別方法、嗜好判別プログラム、及び脳波解析方法に関する。
【背景技術】
【０００２】
近年、個人の嗜好を考慮したサービスが提供されている。それに伴い、個人の嗜好、感性の取得に関する研究が盛んに行われている。嗜好を取得するために、アンケートなどの主観的指標を用いる方法や、生体信号等の生理的指標を用いる方法がある。特に生理的指標を用いる方法の場合、物理量を数値化することができるため、客観的で工学的なアプローチが可能である。
【０００３】
生体信号のひとつである脳波を指標に用い、人間の生理状態を推定する方法はよく知られている。しかし、脳波は微弱信号であるため、瞬き等でノイズがのりやすく、解析が難しい。更に、脳波は個人差が大きいため、従来の解析法では、例えば好きなにおいでの脳波の傾向、分布が一様に定まらなかった。
また、現在主に使われている脳波計は、装置自体が大型であり、装着する被験者にかかる負担が大きく、脳波計を装着しない場合の脳波と異なる測定値が出る可能性がある。
本発明者は、嗜好の判定のひとつとして、匂い刺激を与えたときの被験者の脳波を解析して特徴値を算出し、その被験者が好む匂いと好まない匂いとで特徴値が異なり、分類出来ることを見出した（非特許文献１参照）。
【０００４】
【非特許文献１】横松恵理子、伊藤伸一、満倉靖恵、曹建庭、福見稔、「因子分析を用いた嗜好の取得」、平成19年電気学会産業応用部門大会、社団法人電気学会、p.II-407-410
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかし、非特許文献１の技術において、被験者の脳波の特徴値を抽出しているが、特徴値の傾向が各人により異なり、例えば同じ好きなにおいを嗅いでいるときでも、各人によって特徴値の数値が異なる。よって、ある嗜好における一般的な傾向が一様に定まっていないために、未知の嗜好に対する被験者の脳波から、嗜好を判定することが困難であった。
本発明はこのような点に鑑みてなされたもので、被験者の負担にならない簡易脳波計を用いて脳波を測定することで、被験者の嗜好および／または生理的状態を簡単に判別することができる嗜好判別装置、嗜好判別方法、嗜好判別プログラム、及び脳波解析方法を提供することを目的としている。
【課題を解決するための手段】
【０００６】
上記した目的を達成するために、請求項１記載の嗜好判別装置は、脳波から、人間の生理的状態および／または刺激に対する人間の嗜好を判別する嗜好判別装置において、脳波を検出する検出手段と、前記検出手段により検出された脳波を周波数解析し、前記脳波のデータ行列を作成するデータ行列作成手段と、前記データ行列作成手段で作成されたデータ行列を多変量解析し、前記脳波の特徴値を抽出する特徴値抽出手段と、当該特徴値のデータ集合を複数個格納するデータベースとを有し、前記データベースに格納されている前記特徴値のデータ集合に基づいて、前記特徴値抽出手段により抽出された前記特徴値から、当該脳波における人間の生理的状態および／または刺激に対する人間の嗜好を判断することを特徴とする。
【０００７】
請求項２記載の発明は、請求項１記載の嗜好判別装置において、前記データ行列作成手段が、周波数解析した脳波から、行を時間、列を周波数とするデータ行列を作成することを特徴とする。
【０００８】
請求項３記載の発明は、請求項１記載の嗜好判別装置において、前記多変量解析が、因子分析であることを特徴とする。
【０００９】
請求項４記載の発明は、請求項１記載の嗜好判別装置において、前記特徴値抽出手段により抽出される特徴値が、前記因子分析により算出された因子得点であることを特徴とする。
【００１０】
請求項５記載の発明は、請求項１記載の嗜好判別装置において、前記データベースが、人間のある生理的状態および／またはある刺激に対する人間の嗜好に対応した各特徴値のデータ集合を格納することを特徴とする。
【００１１】
請求項６記載の発明は、請求項１記載の嗜好判別装置において、前記検出手段が、国際10-20法で定めるFp1の部位で脳波を測定し検出することを特徴とする。
【００１２】
請求項７記載の嗜好判別方法は、脳波から、人間の生理的状態および／または刺激に対する人間の嗜好を判別する嗜好判別方法において、脳波を検出するステップと、検出された前記脳波を周波数解析し、前記脳波に対するデータ行列を作成するステップと、データ行列を多変量解析し、前記脳波の特徴値を抽出するステップと、特徴値のデータ集合に基づいて、抽出された前記特徴値から、当該脳波における嗜好を判断するステップとを有することを特徴とする。
【００１３】
請求項８記載の嗜好判別プログラムは、脳波から、人間の生理的状態および／または刺激に対する人間の嗜好を判別するプログラムであって、脳波を検出するステップと、検出された前記脳波を周波数解析し、前記脳波に対するデータ行列を作成するステップと、データ行列を多変量解析し、前記脳波の特徴値を抽出するステップと、複数の特定の嗜好に対応した特徴値のデータ範囲に基づいて、抽出された前記特徴値から、当該脳波における嗜好および／または生理的状態を判断するステップとを実行させる。
【００１４】
請求項９記載の脳波解析装置は、脳波から、人間の生理的状態および／または刺激に対する人間の嗜好を解析する脳波解析装置において、脳波を検出する検出手段と、前記検出手段により検出された脳波を周波数解析し、前記脳波のデータ行列を作成するデータ行列作成手段と、前記データ行列作成手段で作成されたデータ行列を多変量解析し、前記脳波の特徴値を抽出する特徴値抽出手段とを有することを特徴とする。
【００１５】
請求項１０記載の脳波解析方法は、脳波から、人間の生理的状態および／または刺激に対する人間の嗜好を解析する脳波解析方法において、脳波を検出するステップと、検出された前記脳波を周波数解析し、前記脳波に対するデータ行列を作成するステップと、データ行列を多変量解析し、前記脳波の特徴値を抽出するステップとを有することを特徴とする。
【発明の効果】
【００１６】
本発明の嗜好判別装置によれば、多変量解析を行うことにより、データベースに格納されている特徴値のデータ集合が、各種の嗜好、生理的状態によりそれぞれ傾向を持っているため、検出した未知の脳波から抽出した特徴値の嗜好あるいは生理的状態を簡単に判別することができる。
【００１７】
また、請求項２記載の本発明によれば、行を時間、列を周波数とするデータ行列を作成することで、被験者に因らず一般的な嗜好の分布を測定することが出来る。
【００１８】
また、請求項３記載の本発明によれば、脳波のデータ行列を因子分析を用いて解析することにより、より嗜好の判別がしやすい脳波の特徴値を抽出することができる。
【００１９】
また、請求項４記載の本発明によれば、特徴値が因子分析により算出された因子得点であるため、人間の嗜好および生理的状態が簡単に判別することができる。
【００２０】
また、請求項５記載の本発明によれば、データベースに特徴値のデータ集合を格納しているので、未知の脳波の特徴値に対して、その脳波の示す嗜好を判別することが出来る。
【００２１】
また、請求項６記載の本発明によれば、Fp1の部位で脳波を測定し検出することで、被験者が簡易脳波計を装着すればよく、測定時の脳波計装着ストレスによる脳波のノイズを軽減することが出来る。
【００２２】
本発明の嗜好判別方法によれば、多変量解析を行うことにより、データベースに格納されている特徴値のデータ集合が、各種の嗜好、生理的状態によりそれぞれ傾向を持っているため、検出した未知の脳波から抽出した特徴値の嗜好あるいは生理的状態を簡単に判別することができる。
【００２３】
本発明の嗜好判別プログラムによれば、多変量解析を行うことにより、データベースに格納されている特徴値のデータ集合が、各種の嗜好、生理的状態によりそれぞれ傾向を持っているため、検出した未知の脳波から抽出した特徴値の嗜好あるいは生理的状態を簡単に判別することができる。
【００２４】
本発明の脳波解析装置によれば、多変量解析を行うことにより、データベースに格納されている特徴値のデータ集合が、各種の嗜好、生理的状態によりそれぞれ傾向を持つことができる。
【００２５】
本発明の脳波解析方法によれば、多変量解析を行うことにより、データベースに格納されている特徴値のデータ集合が、各種の嗜好、生理的状態によりそれぞれ傾向を持つことができる。
【発明を実施するための最良の形態】
【００２６】
以下、図面を参照して本発明の実施の形態を説明する。
図１は、本発明の一実施形態に係る嗜好判別装置１００の概略構成を示すブロック図である。
簡易脳波計２（検出手段）は被験者の脳波を測定する。簡易脳波計２は、被験者の頭部、国際10-20法におけるFp1に当たる箇所（左前部前頭葉）に、センサバンドで固定された電極が接する構造となっている。Fp1に当たる箇所は、本装置における嗜好判別に十分な脳波が発生する箇所である。Fp１は一般に感情や情動に関する反応を計測することが出来ることが知られている。
簡易脳波計２での測定は、電気的に不活性な耳朶（本発明は左耳朶）に電極を装着し、電気的に活性な頭皮上にもう一方の電極を装着して二極間の電位差の変動を導出する単極誘導法によって計測する。単極誘導法は、脳の各部位の電気的活動を観察するのに適している。
【００２７】
メモリ３は、ワーク領域となるもので、その機能はＲＡＭなどにより実現できる。
情報記憶媒体４（コンピュータにより読み取り可能な媒体）は、プログラムやデータなどを格納するものであり、その機能は、光ディスク（ＣＤ、ＤＶＤ）、光磁気ディスク（ＭＯ）、磁気ディスク、ハードディスク、磁気テープ、或いはメモリ（ＲＯＭ）などにより実現できる。処理部５は、情報記憶媒体４に格納されるプログラム（データ）に基づいて本実施形態の種々の処理を行う。即ち情報記憶媒体４には、本実施形態の各部としてコンピュータを機能させるためのプログラム（各部の処理をコンピュータに実行させるためのプログラム）が記憶される。
【００２８】
処理部５はデータ行列作成部５１、特徴値抽出部５２、制御部５３を含む。処理部５（プロセッサ）は、情報記憶媒体４に記憶されるプログラム等に基づいて、所与の処理を行う。また、処理部５はメモリ３をワーク領域として各種処理を行う。処理部５の機能は各種プロセッサ（ＣＰＵ、ＤＳＰ等）、ＡＳＩＣ（ゲートアレイ等）などのハードウェアや、プログラムにより実現できる。
データベース６は、処理部５で計算した特徴値のデータ集合を格納、蓄積する。特徴値の数値範囲により、「好む匂い」「嫌いな匂い」「薬を飲んでいる」「飲酒している」等といった分類テーブルを有している。
【００２９】
刺激発生部７は被験者の脳波に影響する何らかの刺激を発生させる。本実施例では匂い刺激を発生させる。なお、刺激発生部７は、制御部５３により制御される構成でもよい。
【実施例１】
【００３０】
脳波の特徴値の抽出、データベースへの格納における嗜好判別装置１００の動作の一例として、被験者が好む匂い刺激を与えた場合について図２を参照して説明する。
まず、刺激発生部７は、被験者が好む匂い刺激を発生する（ステップＳ１０）。簡易脳波計２は、装着した被験者の脳波を測定する（ステップＳ１１）。この時の刺激発生は、無臭状態が15秒、次に刺激発生部７で匂い刺激を発生する時間が30秒、その後匂い刺激を取り除いた無臭状態が15秒の計１分間である。
次に、処理部５内のデータ行列作成部５１は、測定した脳波データについて高速フーリエ変換による周波数解析を行う（ステップＳ１３）。高速フーリエ変換は、サンプリング周波数128Hzで1秒間毎に行う。出力される周波数範囲は1Hz間隔で1-24Hzであるが、本装置は4-22Hzのバンドパスフィルタを用いるため、有効となる周波数範囲は4-22Hzである。また、脳波データは、計測の状態が安定したと考えられる、計測開始5秒後から計測終了5秒前の50秒間を用いる。
次にデータ行列作成部５１は、高速フーリエ変換を行った脳波データから、行が時間、列が周波数成分のｎ×ｍデータ行列を作成する(ステップＳ１５)。
【００３１】
次に、特徴値抽出部５２は、ステップＳ１５でデータ行列作成部５１が作成したデータ行列から、因子分析により特徴値の抽出を行う(ステップＳ１７)。因子分析については後述する。
処理部５は、抽出した特徴値をデータベース６に「好む匂い」として格納する（ステップＳ１９）。
【００３２】
図３に因子分析のモデルを示す。因子分析は、多変量間の相関関係を、少数の潜在的な因子によって説明するための統計的手法である。つまり、あるいくつかの変量の中にみられる変動の共通部分について明らかにしようとする手法である。変量とは実際に測定された変動、共通因子とはn個1組の変量を表すのに共通に用いられる仮説的変動を表すものであり、独自因子とはn個の変量のそれぞれに対応して考えられている固有の変動を表すものである。
本発明では因子分析を用いることで、測定した脳波データにおけるノイズを効果的に取り除き、嗜好の分布を一般化することができるので、嗜好の判別を容易にすることができる。
【００３３】
ステップＳ１７の特徴値抽出方法である因子分析での解析手法について、図３を参照して説明する。まず、ステップＳ１５で作成されたデータ行列xij(変量jについての個体iの測定値)について、変量ｊの平均を用い、（１）式に従い正規化する（ステップＳ２０１）。なお、本実施例での変量は時間、個体は周波数である。
【００３４】
【数１】
【００３５】
次に、正規化したデータ行列Zから（２）式より相関行列Rを求める。相関行列Rの各成分は（３）式より求められる（ステップＳ２０３）。
【００３６】
【数２】
【００３７】
次に、重相関係数の平方（ＳＭＣ）を用いて（４）式より共通性の推定値を求め、ステップＳ２０３で求めた相関係数の対角成分に代入する（ステップＳ２０５）。ここで、r(j)は相関行列Rの第j列成分からなるベクトルのj番目の成分1を除いたベクトル、R(j)はもとの相関行列Rから，そのj行とj列にあたる成分をすべて取り除いた行列である。
【００３８】
【数３】
【００３９】
重相関係数は、2つ以上の変量の１次結合としてできる合成変量と、他のもう1つの変量との相関を表し、特に本発明では、重相関係数が最大となるように1次結合の重みが選ばれている。
次に、共通性の推定値が対角成分に代入された相関行列を固有値分解し、主因子法を用いて（５）式により因子を抽出する（ステップＳ２０７）。
【００４０】
【数４】
【００４１】
求めた因子のうち、寄与率が６０％となる第５因子までを用いて、（６）式を最大にするバリマックス回転を行う（ステップＳ２０９）。
【００４２】
【数５】
【００４３】
そして、最小二乗法により因子得点Fを求めることで（ステップＳ２１１）、特徴値を抽出する。
【００４４】
【数６】
【００４５】
本発明では、求めた因子のうち、第１因子または第２因子が匂い刺激の影響を表す因子であると仮定した。
図５Ａ、Ｂおよび図６Ａ、Ｂは、４人の被験者についての、周波数に対する因子得点の定常状態からの変化量を示す。因子得点は第１因子（実線）および第２因子（点線）について示している。図５は好む匂い刺激を与えたときの結果、図６は嫌いな匂い刺激を与えたときの結果である。
図５および図６より、好む匂い刺激を与えたときは第２因子の因子得点が6〜8Hzの間で急激に低下しており、下に凸となっている。それに対し、嫌いな匂い刺激を与えたときは、第２因子の因子得点が同じ6〜8Hzの間で上に凸となっている。また、被験者毎に比較すると、好む匂い刺激を与えた場合と嫌いな匂い刺激を与えた場合とで、周波数軸に対してほぼ対称となっている。
【００４６】
好む匂い刺激における傾向や分布を明確にするために、図５で表示した第１因子及び第２因子の値をニューラルネットワークに入力した場合の、中間ユニットの値２つをＸ軸とＹ軸にして、４人の被験者について図7に表示する。図７で、neutralは定常状態、fragranceは好む匂い刺激を与えたときの分布である。
図７より、好む匂い刺激を与えたときは、各被験者ともほぼ同様の位置に分布しており、好む匂い刺激を与えられたときの脳波に対して本発明の解析手法を用いれば、被験者に因らない一般化した傾向を判別することができる。
【実施例２】
【００４７】
次に、脳波から被験者の嗜好を判定する嗜好判別装置１００の動作の一例として、被験者の嗜好が未知の、匂い刺激を与えた場合について図９を参照して説明する。
刺激発生部７は匂い刺激を発生する（ステップＳ３０）。簡易脳波計２は、装着した被験者の脳波を測定する（ステップＳ３１）。この時の刺激発生時間は実施例１と同様である。
処理部５は、測定した脳波データを解析して、因子分析により特徴値を抽出する（ステップＳ３３）。解析動作の詳細は、図２のステップＳ１３〜Ｓ１７と同様であるので説明を省略する。
処理部５の制御部５３は、ステップＳ３３で特徴値抽出部５２が抽出した特徴値について、データベース６にあらかじめ格納されている、特徴値の数値範囲と「好む匂い」または「嫌いな匂い」が対応した分類テーブルを参照して、嗜好を判定する（ステップＳ３５）。
【００４８】
図１０に、被験者１〜７について、未知の嗜好の匂い刺激を測定して嗜好を判定した際の、嗜好判定装置１００の認識率を示す。図１０より、特に嫌いな匂いに対する認識率が高いことが判る。好む匂いの認識率は、２人の被験者で認識率がよくなかったが、これは「好む匂い」のサンプルが１種類しかなかったため、個人の匂いに対する好みの違いが表れてしまったと推察される。
【産業上の利用可能性】
【００４９】
以上説明したように、本実施形態によれば、匂い刺激を与えた被験者の脳波を高速フーリエ変換し、行が時間、列が周波数成分のデータ行列を作成して、因子分析により特徴値である因子得点を求める。これにより、被験者に因らず一般的な嗜好の傾向、分布を測定することが出来る。また、未知の嗜好に対して被験者に匂い刺激を与えた場合に、データベースにあらかじめ格納している嗜好の傾向を参照して、被験者の嗜好を判定することが出来る。
【００５０】
以上、本発明の実施形態について図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、本発明の要旨を逸脱しない範囲の設計変更等も含まれる。
【００５１】
例えば、上記実施形態では、匂いの嗜好の判定について説明したが、本発明はその他に、食欲のわく匂い、薬を飲んでいる、飲酒している、メガネの矯正視力が合っていない、食欲がある、喫煙している、排泄をしたい、といった嗜好、生理的状態についても、算出した特徴値のデータ集合がそれぞれ傾向を持っている。そのため、未知の脳波について特徴値を抽出すれば、それぞれの嗜好、生理的状態を判定することができる。
【００５２】
例えば、車の運転をしているドライバーの脳波について本発明を適用すれば、飲酒の有無が判定できるので、ドライバーがごまかしても飲酒を発覚でき、飲酒運転の撲滅に役立つ。また、メガネの矯正時に装用者の脳波を測定して本発明を適用すれば、装用者にとって正しい矯正になっているのかを判定することが出来る。
【図面の簡単な説明】
【００５３】
【図１】本発明の一実施形態に係る嗜好判別装置１００の概略構成を示すブロック図である。
【図２】同上の実施形態における脳波の特徴値抽出の動作を示すフローチャートである。
【図３】同上の実施形態における因子分析のモデルである。
【図４】同上の実施形態における因子分析の演算を示すフローチャートである。
【図５Ａ】実施例１における、好む匂い刺激を与えた時の、周波数と因子得点の変化量を示すグラフである。
【図５Ｂ】実施例１における、好む匂い刺激を与えた時の、周波数と因子得点の変化量を示すグラフである。
【図６Ａ】同上の実施例１における、嫌いな匂い刺激を与えた時の、周波数と因子得点の変化量を示すグラフである。
【図６Ｂ】同上の実施例１における、嫌いな匂い刺激を与えた時の、周波数と因子得点の変化量を示すグラフである。
【図７】同上の実施例１における、好む匂い刺激を与えた時の分布を示すグラフである。
【符号の説明】
【００５４】
１００ 嗜好判別装置
２ 簡易脳波計
３ メモリ
４ 情報記憶媒体
５ 処理部
５１ データ行列作成部
５２ 特徴値抽出部
５３ 制御部
６ データベース
【技術分野】
【０００１】
本発明は、脳波から人間の生理的状態および／または刺激に対する人間の嗜好を判別する嗜好判別装置、嗜好判別方法、嗜好判別プログラム、及び脳波解析方法に関する。
【背景技術】
【０００２】
近年、個人の嗜好を考慮したサービスが提供されている。それに伴い、個人の嗜好、感性の取得に関する研究が盛んに行われている。嗜好を取得するために、アンケートなどの主観的指標を用いる方法や、生体信号等の生理的指標を用いる方法がある。特に生理的指標を用いる方法の場合、物理量を数値化することができるため、客観的で工学的なアプローチが可能である。
【０００３】
生体信号のひとつである脳波を指標に用い、人間の生理状態を推定する方法はよく知られている。しかし、脳波は微弱信号であるため、瞬き等でノイズがのりやすく、解析が難しい。更に、脳波は個人差が大きいため、従来の解析法では、例えば好きなにおいでの脳波の傾向、分布が一様に定まらなかった。
また、現在主に使われている脳波計は、装置自体が大型であり、装着する被験者にかかる負担が大きく、脳波計を装着しない場合の脳波と異なる測定値が出る可能性がある。
本発明者は、嗜好の判定のひとつとして、匂い刺激を与えたときの被験者の脳波を解析して特徴値を算出し、その被験者が好む匂いと好まない匂いとで特徴値が異なり、分類出来ることを見出した（非特許文献１参照）。
【０００４】
【非特許文献１】横松恵理子、伊藤伸一、満倉靖恵、曹建庭、福見稔、「因子分析を用いた嗜好の取得」、平成19年電気学会産業応用部門大会、社団法人電気学会、p.II-407-410
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかし、非特許文献１の技術において、被験者の脳波の特徴値を抽出しているが、特徴値の傾向が各人により異なり、例えば同じ好きなにおいを嗅いでいるときでも、各人によって特徴値の数値が異なる。よって、ある嗜好における一般的な傾向が一様に定まっていないために、未知の嗜好に対する被験者の脳波から、嗜好を判定することが困難であった。
本発明はこのような点に鑑みてなされたもので、被験者の負担にならない簡易脳波計を用いて脳波を測定することで、被験者の嗜好および／または生理的状態を簡単に判別することができる嗜好判別装置、嗜好判別方法、嗜好判別プログラム、及び脳波解析方法を提供することを目的としている。
【課題を解決するための手段】
【０００６】
上記した目的を達成するために、請求項１記載の嗜好判別装置は、脳波から、人間の生理的状態および／または刺激に対する人間の嗜好を判別する嗜好判別装置において、脳波を検出する検出手段と、前記検出手段により検出された脳波を周波数解析し、前記脳波のデータ行列を作成するデータ行列作成手段と、前記データ行列作成手段で作成されたデータ行列を多変量解析し、前記脳波の特徴値を抽出する特徴値抽出手段と、当該特徴値のデータ集合を複数個格納するデータベースとを有し、前記データベースに格納されている前記特徴値のデータ集合に基づいて、前記特徴値抽出手段により抽出された前記特徴値から、当該脳波における人間の生理的状態および／または刺激に対する人間の嗜好を判断することを特徴とする。
【０００７】
請求項２記載の発明は、請求項１記載の嗜好判別装置において、前記データ行列作成手段が、周波数解析した脳波から、行を時間、列を周波数とするデータ行列を作成することを特徴とする。
【０００８】
請求項３記載の発明は、請求項１記載の嗜好判別装置において、前記多変量解析が、因子分析であることを特徴とする。
【０００９】
請求項４記載の発明は、請求項１記載の嗜好判別装置において、前記特徴値抽出手段により抽出される特徴値が、前記因子分析により算出された因子得点であることを特徴とする。
【００１０】
請求項５記載の発明は、請求項１記載の嗜好判別装置において、前記データベースが、人間のある生理的状態および／またはある刺激に対する人間の嗜好に対応した各特徴値のデータ集合を格納することを特徴とする。
【００１１】
請求項６記載の発明は、請求項１記載の嗜好判別装置において、前記検出手段が、国際10-20法で定めるFp1の部位で脳波を測定し検出することを特徴とする。
【００１２】
請求項７記載の嗜好判別方法は、脳波から、人間の生理的状態および／または刺激に対する人間の嗜好を判別する嗜好判別方法において、脳波を検出するステップと、検出された前記脳波を周波数解析し、前記脳波に対するデータ行列を作成するステップと、データ行列を多変量解析し、前記脳波の特徴値を抽出するステップと、特徴値のデータ集合に基づいて、抽出された前記特徴値から、当該脳波における嗜好を判断するステップとを有することを特徴とする。
【００１３】
請求項８記載の嗜好判別プログラムは、脳波から、人間の生理的状態および／または刺激に対する人間の嗜好を判別するプログラムであって、脳波を検出するステップと、検出された前記脳波を周波数解析し、前記脳波に対するデータ行列を作成するステップと、データ行列を多変量解析し、前記脳波の特徴値を抽出するステップと、複数の特定の嗜好に対応した特徴値のデータ範囲に基づいて、抽出された前記特徴値から、当該脳波における嗜好および／または生理的状態を判断するステップとを実行させる。
【００１４】
請求項９記載の脳波解析装置は、脳波から、人間の生理的状態および／または刺激に対する人間の嗜好を解析する脳波解析装置において、脳波を検出する検出手段と、前記検出手段により検出された脳波を周波数解析し、前記脳波のデータ行列を作成するデータ行列作成手段と、前記データ行列作成手段で作成されたデータ行列を多変量解析し、前記脳波の特徴値を抽出する特徴値抽出手段とを有することを特徴とする。
【００１５】
請求項１０記載の脳波解析方法は、脳波から、人間の生理的状態および／または刺激に対する人間の嗜好を解析する脳波解析方法において、脳波を検出するステップと、検出された前記脳波を周波数解析し、前記脳波に対するデータ行列を作成するステップと、データ行列を多変量解析し、前記脳波の特徴値を抽出するステップとを有することを特徴とする。
【発明の効果】
【００１６】
本発明の嗜好判別装置によれば、多変量解析を行うことにより、データベースに格納されている特徴値のデータ集合が、各種の嗜好、生理的状態によりそれぞれ傾向を持っているため、検出した未知の脳波から抽出した特徴値の嗜好あるいは生理的状態を簡単に判別することができる。
【００１７】
また、請求項２記載の本発明によれば、行を時間、列を周波数とするデータ行列を作成することで、被験者に因らず一般的な嗜好の分布を測定することが出来る。
【００１８】
また、請求項３記載の本発明によれば、脳波のデータ行列を因子分析を用いて解析することにより、より嗜好の判別がしやすい脳波の特徴値を抽出することができる。
【００１９】
また、請求項４記載の本発明によれば、特徴値が因子分析により算出された因子得点であるため、人間の嗜好および生理的状態が簡単に判別することができる。
【００２０】
また、請求項５記載の本発明によれば、データベースに特徴値のデータ集合を格納しているので、未知の脳波の特徴値に対して、その脳波の示す嗜好を判別することが出来る。
【００２１】
また、請求項６記載の本発明によれば、Fp1の部位で脳波を測定し検出することで、被験者が簡易脳波計を装着すればよく、測定時の脳波計装着ストレスによる脳波のノイズを軽減することが出来る。
【００２２】
本発明の嗜好判別方法によれば、多変量解析を行うことにより、データベースに格納されている特徴値のデータ集合が、各種の嗜好、生理的状態によりそれぞれ傾向を持っているため、検出した未知の脳波から抽出した特徴値の嗜好あるいは生理的状態を簡単に判別することができる。
【００２３】
本発明の嗜好判別プログラムによれば、多変量解析を行うことにより、データベースに格納されている特徴値のデータ集合が、各種の嗜好、生理的状態によりそれぞれ傾向を持っているため、検出した未知の脳波から抽出した特徴値の嗜好あるいは生理的状態を簡単に判別することができる。
【００２４】
本発明の脳波解析装置によれば、多変量解析を行うことにより、データベースに格納されている特徴値のデータ集合が、各種の嗜好、生理的状態によりそれぞれ傾向を持つことができる。
【００２５】
本発明の脳波解析方法によれば、多変量解析を行うことにより、データベースに格納されている特徴値のデータ集合が、各種の嗜好、生理的状態によりそれぞれ傾向を持つことができる。
【発明を実施するための最良の形態】
【００２６】
以下、図面を参照して本発明の実施の形態を説明する。
図１は、本発明の一実施形態に係る嗜好判別装置１００の概略構成を示すブロック図である。
簡易脳波計２（検出手段）は被験者の脳波を測定する。簡易脳波計２は、被験者の頭部、国際10-20法におけるFp1に当たる箇所（左前部前頭葉）に、センサバンドで固定された電極が接する構造となっている。Fp1に当たる箇所は、本装置における嗜好判別に十分な脳波が発生する箇所である。Fp１は一般に感情や情動に関する反応を計測することが出来ることが知られている。
簡易脳波計２での測定は、電気的に不活性な耳朶（本発明は左耳朶）に電極を装着し、電気的に活性な頭皮上にもう一方の電極を装着して二極間の電位差の変動を導出する単極誘導法によって計測する。単極誘導法は、脳の各部位の電気的活動を観察するのに適している。
【００２７】
メモリ３は、ワーク領域となるもので、その機能はＲＡＭなどにより実現できる。
情報記憶媒体４（コンピュータにより読み取り可能な媒体）は、プログラムやデータなどを格納するものであり、その機能は、光ディスク（ＣＤ、ＤＶＤ）、光磁気ディスク（ＭＯ）、磁気ディスク、ハードディスク、磁気テープ、或いはメモリ（ＲＯＭ）などにより実現できる。処理部５は、情報記憶媒体４に格納されるプログラム（データ）に基づいて本実施形態の種々の処理を行う。即ち情報記憶媒体４には、本実施形態の各部としてコンピュータを機能させるためのプログラム（各部の処理をコンピュータに実行させるためのプログラム）が記憶される。
【００２８】
処理部５はデータ行列作成部５１、特徴値抽出部５２、制御部５３を含む。処理部５（プロセッサ）は、情報記憶媒体４に記憶されるプログラム等に基づいて、所与の処理を行う。また、処理部５はメモリ３をワーク領域として各種処理を行う。処理部５の機能は各種プロセッサ（ＣＰＵ、ＤＳＰ等）、ＡＳＩＣ（ゲートアレイ等）などのハードウェアや、プログラムにより実現できる。
データベース６は、処理部５で計算した特徴値のデータ集合を格納、蓄積する。特徴値の数値範囲により、「好む匂い」「嫌いな匂い」「薬を飲んでいる」「飲酒している」等といった分類テーブルを有している。
【００２９】
刺激発生部７は被験者の脳波に影響する何らかの刺激を発生させる。本実施例では匂い刺激を発生させる。なお、刺激発生部７は、制御部５３により制御される構成でもよい。
【実施例１】
【００３０】
脳波の特徴値の抽出、データベースへの格納における嗜好判別装置１００の動作の一例として、被験者が好む匂い刺激を与えた場合について図２を参照して説明する。
まず、刺激発生部７は、被験者が好む匂い刺激を発生する（ステップＳ１０）。簡易脳波計２は、装着した被験者の脳波を測定する（ステップＳ１１）。この時の刺激発生は、無臭状態が15秒、次に刺激発生部７で匂い刺激を発生する時間が30秒、その後匂い刺激を取り除いた無臭状態が15秒の計１分間である。
次に、処理部５内のデータ行列作成部５１は、測定した脳波データについて高速フーリエ変換による周波数解析を行う（ステップＳ１３）。高速フーリエ変換は、サンプリング周波数128Hzで1秒間毎に行う。出力される周波数範囲は1Hz間隔で1-24Hzであるが、本装置は4-22Hzのバンドパスフィルタを用いるため、有効となる周波数範囲は4-22Hzである。また、脳波データは、計測の状態が安定したと考えられる、計測開始5秒後から計測終了5秒前の50秒間を用いる。
次にデータ行列作成部５１は、高速フーリエ変換を行った脳波データから、行が時間、列が周波数成分のｎ×ｍデータ行列を作成する(ステップＳ１５)。
【００３１】
次に、特徴値抽出部５２は、ステップＳ１５でデータ行列作成部５１が作成したデータ行列から、因子分析により特徴値の抽出を行う(ステップＳ１７)。因子分析については後述する。
処理部５は、抽出した特徴値をデータベース６に「好む匂い」として格納する（ステップＳ１９）。
【００３２】
図３に因子分析のモデルを示す。因子分析は、多変量間の相関関係を、少数の潜在的な因子によって説明するための統計的手法である。つまり、あるいくつかの変量の中にみられる変動の共通部分について明らかにしようとする手法である。変量とは実際に測定された変動、共通因子とはn個1組の変量を表すのに共通に用いられる仮説的変動を表すものであり、独自因子とはn個の変量のそれぞれに対応して考えられている固有の変動を表すものである。
本発明では因子分析を用いることで、測定した脳波データにおけるノイズを効果的に取り除き、嗜好の分布を一般化することができるので、嗜好の判別を容易にすることができる。
【００３３】
ステップＳ１７の特徴値抽出方法である因子分析での解析手法について、図３を参照して説明する。まず、ステップＳ１５で作成されたデータ行列xij(変量jについての個体iの測定値)について、変量ｊの平均を用い、（１）式に従い正規化する（ステップＳ２０１）。なお、本実施例での変量は時間、個体は周波数である。
【００３４】
【数１】
【００３５】
次に、正規化したデータ行列Zから（２）式より相関行列Rを求める。相関行列Rの各成分は（３）式より求められる（ステップＳ２０３）。
【００３６】
【数２】
【００３７】
次に、重相関係数の平方（ＳＭＣ）を用いて（４）式より共通性の推定値を求め、ステップＳ２０３で求めた相関係数の対角成分に代入する（ステップＳ２０５）。ここで、r(j)は相関行列Rの第j列成分からなるベクトルのj番目の成分1を除いたベクトル、R(j)はもとの相関行列Rから，そのj行とj列にあたる成分をすべて取り除いた行列である。
【００３８】
【数３】
【００３９】
重相関係数は、2つ以上の変量の１次結合としてできる合成変量と、他のもう1つの変量との相関を表し、特に本発明では、重相関係数が最大となるように1次結合の重みが選ばれている。
次に、共通性の推定値が対角成分に代入された相関行列を固有値分解し、主因子法を用いて（５）式により因子を抽出する（ステップＳ２０７）。
【００４０】
【数４】
【００４１】
求めた因子のうち、寄与率が６０％となる第５因子までを用いて、（６）式を最大にするバリマックス回転を行う（ステップＳ２０９）。
【００４２】
【数５】
【００４３】
そして、最小二乗法により因子得点Fを求めることで（ステップＳ２１１）、特徴値を抽出する。
【００４４】
【数６】
【００４５】
本発明では、求めた因子のうち、第１因子または第２因子が匂い刺激の影響を表す因子であると仮定した。
図５Ａ、Ｂおよび図６Ａ、Ｂは、４人の被験者についての、周波数に対する因子得点の定常状態からの変化量を示す。因子得点は第１因子（実線）および第２因子（点線）について示している。図５は好む匂い刺激を与えたときの結果、図６は嫌いな匂い刺激を与えたときの結果である。
図５および図６より、好む匂い刺激を与えたときは第２因子の因子得点が6〜8Hzの間で急激に低下しており、下に凸となっている。それに対し、嫌いな匂い刺激を与えたときは、第２因子の因子得点が同じ6〜8Hzの間で上に凸となっている。また、被験者毎に比較すると、好む匂い刺激を与えた場合と嫌いな匂い刺激を与えた場合とで、周波数軸に対してほぼ対称となっている。
【００４６】
好む匂い刺激における傾向や分布を明確にするために、図５で表示した第１因子及び第２因子の値をニューラルネットワークに入力した場合の、中間ユニットの値２つをＸ軸とＹ軸にして、４人の被験者について図7に表示する。図７で、neutralは定常状態、fragranceは好む匂い刺激を与えたときの分布である。
図７より、好む匂い刺激を与えたときは、各被験者ともほぼ同様の位置に分布しており、好む匂い刺激を与えられたときの脳波に対して本発明の解析手法を用いれば、被験者に因らない一般化した傾向を判別することができる。
【実施例２】
【００４７】
次に、脳波から被験者の嗜好を判定する嗜好判別装置１００の動作の一例として、被験者の嗜好が未知の、匂い刺激を与えた場合について図９を参照して説明する。
刺激発生部７は匂い刺激を発生する（ステップＳ３０）。簡易脳波計２は、装着した被験者の脳波を測定する（ステップＳ３１）。この時の刺激発生時間は実施例１と同様である。
処理部５は、測定した脳波データを解析して、因子分析により特徴値を抽出する（ステップＳ３３）。解析動作の詳細は、図２のステップＳ１３〜Ｓ１７と同様であるので説明を省略する。
処理部５の制御部５３は、ステップＳ３３で特徴値抽出部５２が抽出した特徴値について、データベース６にあらかじめ格納されている、特徴値の数値範囲と「好む匂い」または「嫌いな匂い」が対応した分類テーブルを参照して、嗜好を判定する（ステップＳ３５）。
【００４８】
図１０に、被験者１〜７について、未知の嗜好の匂い刺激を測定して嗜好を判定した際の、嗜好判定装置１００の認識率を示す。図１０より、特に嫌いな匂いに対する認識率が高いことが判る。好む匂いの認識率は、２人の被験者で認識率がよくなかったが、これは「好む匂い」のサンプルが１種類しかなかったため、個人の匂いに対する好みの違いが表れてしまったと推察される。
【産業上の利用可能性】
【００４９】
以上説明したように、本実施形態によれば、匂い刺激を与えた被験者の脳波を高速フーリエ変換し、行が時間、列が周波数成分のデータ行列を作成して、因子分析により特徴値である因子得点を求める。これにより、被験者に因らず一般的な嗜好の傾向、分布を測定することが出来る。また、未知の嗜好に対して被験者に匂い刺激を与えた場合に、データベースにあらかじめ格納している嗜好の傾向を参照して、被験者の嗜好を判定することが出来る。
【００５０】
以上、本発明の実施形態について図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、本発明の要旨を逸脱しない範囲の設計変更等も含まれる。
【００５１】
例えば、上記実施形態では、匂いの嗜好の判定について説明したが、本発明はその他に、食欲のわく匂い、薬を飲んでいる、飲酒している、メガネの矯正視力が合っていない、食欲がある、喫煙している、排泄をしたい、といった嗜好、生理的状態についても、算出した特徴値のデータ集合がそれぞれ傾向を持っている。そのため、未知の脳波について特徴値を抽出すれば、それぞれの嗜好、生理的状態を判定することができる。
【００５２】
例えば、車の運転をしているドライバーの脳波について本発明を適用すれば、飲酒の有無が判定できるので、ドライバーがごまかしても飲酒を発覚でき、飲酒運転の撲滅に役立つ。また、メガネの矯正時に装用者の脳波を測定して本発明を適用すれば、装用者にとって正しい矯正になっているのかを判定することが出来る。
【図面の簡単な説明】
【００５３】
【図１】本発明の一実施形態に係る嗜好判別装置１００の概略構成を示すブロック図である。
【図２】同上の実施形態における脳波の特徴値抽出の動作を示すフローチャートである。
【図３】同上の実施形態における因子分析のモデルである。
【図４】同上の実施形態における因子分析の演算を示すフローチャートである。
【図５Ａ】実施例１における、好む匂い刺激を与えた時の、周波数と因子得点の変化量を示すグラフである。
【図５Ｂ】実施例１における、好む匂い刺激を与えた時の、周波数と因子得点の変化量を示すグラフである。
【図６Ａ】同上の実施例１における、嫌いな匂い刺激を与えた時の、周波数と因子得点の変化量を示すグラフである。
【図６Ｂ】同上の実施例１における、嫌いな匂い刺激を与えた時の、周波数と因子得点の変化量を示すグラフである。
【図７】同上の実施例１における、好む匂い刺激を与えた時の分布を示すグラフである。
【符号の説明】
【００５４】
１００ 嗜好判別装置
２ 簡易脳波計
３ メモリ
４ 情報記憶媒体
５ 処理部
５１ データ行列作成部
５２ 特徴値抽出部
５３ 制御部
６ データベース
【特許請求の範囲】
【請求項１】
脳波から、人間の生理的状態および／または刺激に対する人間の嗜好を判別する嗜好判別装置において、
脳波を検出する検出手段と、
前記検出手段により検出された脳波を周波数解析し、前記脳波のデータ行列を作成するデータ行列作成手段と、
前記データ行列作成手段で作成されたデータ行列を多変量解析し、前記脳波の特徴値を抽出する特徴値抽出手段と、
当該特徴値のデータ集合を複数個格納するデータベースとを有し、
前記データベースに格納されている前記特徴値のデータ集合に基づいて、前記特徴値抽出手段により抽出された前記特徴値から、当該脳波における人間の生理的状態および／または刺激に対する人間の嗜好を判断する
ことを特徴とする嗜好判別装置。
【請求項２】
前記データ行列作成手段が、周波数解析した脳波から、行を時間、列を周波数とするデータ行列を作成する
ことを特徴とする請求項１記載の嗜好判別装置。
【請求項３】
前記多変量解析が、因子分析である
ことを特徴とする請求項１記載の嗜好判別装置。
【請求項４】
前記特徴値抽出手段により抽出される特徴値が、前記因子分析により算出された因子得点である
ことを特徴とする請求項１〜３に記載の嗜好判別装置。
【請求項５】
前記データベースが、人間のある生理的状態および／またはある刺激に対する人間の嗜好に対応した各特徴値のデータ集合を格納する
ことを特徴とする請求項１記載の嗜好判別装置。
【請求項６】
前記検出手段が、国際10-20法で定めるFp1の部位で脳波を測定し検出することを特徴とする請求項1記載の嗜好判別装置。
【請求項７】
脳波から、人間の生理的状態および／または刺激に対する人間の嗜好を判別する嗜好判別方法において、
脳波を検出するステップと、
検出された前記脳波を周波数解析し、前記脳波に対するデータ行列を作成するステップと、
データ行列を多変量解析し、前記脳波の特徴値を抽出するステップと、
特徴値のデータ集合に基づいて、抽出された前記特徴値から、当該脳波における嗜好を判断するステップと
を有することを特徴とする嗜好判別方法。
【請求項８】
脳波から、人間の生理的状態および／または刺激に対する人間の嗜好を判別するプログラムであって、
脳波を検出するステップと、
検出された前記脳波を周波数解析し、前記脳波に対するデータ行列を作成するステップと、
データ行列を多変量解析し、前記脳波の特徴値を抽出するステップと、
複数の特定の嗜好に対応した特徴値のデータ範囲に基づいて、抽出された前記特徴値から、当該脳波における嗜好および／または生理的状態を判断するステップと
を実行させるための嗜好判別プログラム。
【請求項９】
脳波から、人間の生理的状態および／または刺激に対する人間の嗜好を解析する脳波解析装置において、
脳波を検出する検出手段と、
前記検出手段により検出された脳波を周波数解析し、前記脳波のデータ行列を作成するデータ行列作成手段と、
前記データ行列作成手段で作成されたデータ行列を多変量解析し、前記脳波の特徴値を抽出する特徴値抽出手段と
を有することを特徴とする脳波解析装置。
【請求項１０】
脳波から、人間の生理的状態および／または刺激に対する人間の嗜好を解析する脳波解析方法において、
脳波を検出するステップと、
検出された前記脳波を周波数解析し、前記脳波に対するデータ行列を作成するステップと、
データ行列を多変量解析し、前記脳波の特徴値を抽出するステップと
を有することを特徴とする脳波解析方法。
【請求項１】
脳波から、人間の生理的状態および／または刺激に対する人間の嗜好を判別する嗜好判別装置において、
脳波を検出する検出手段と、
前記検出手段により検出された脳波を周波数解析し、前記脳波のデータ行列を作成するデータ行列作成手段と、
前記データ行列作成手段で作成されたデータ行列を多変量解析し、前記脳波の特徴値を抽出する特徴値抽出手段と、
当該特徴値のデータ集合を複数個格納するデータベースとを有し、
前記データベースに格納されている前記特徴値のデータ集合に基づいて、前記特徴値抽出手段により抽出された前記特徴値から、当該脳波における人間の生理的状態および／または刺激に対する人間の嗜好を判断する
ことを特徴とする嗜好判別装置。
【請求項２】
前記データ行列作成手段が、周波数解析した脳波から、行を時間、列を周波数とするデータ行列を作成する
ことを特徴とする請求項１記載の嗜好判別装置。
【請求項３】
前記多変量解析が、因子分析である
ことを特徴とする請求項１記載の嗜好判別装置。
【請求項４】
前記特徴値抽出手段により抽出される特徴値が、前記因子分析により算出された因子得点である
ことを特徴とする請求項１〜３に記載の嗜好判別装置。
【請求項５】
前記データベースが、人間のある生理的状態および／またはある刺激に対する人間の嗜好に対応した各特徴値のデータ集合を格納する
ことを特徴とする請求項１記載の嗜好判別装置。
【請求項６】
前記検出手段が、国際10-20法で定めるFp1の部位で脳波を測定し検出することを特徴とする請求項1記載の嗜好判別装置。
【請求項７】
脳波から、人間の生理的状態および／または刺激に対する人間の嗜好を判別する嗜好判別方法において、
脳波を検出するステップと、
検出された前記脳波を周波数解析し、前記脳波に対するデータ行列を作成するステップと、
データ行列を多変量解析し、前記脳波の特徴値を抽出するステップと、
特徴値のデータ集合に基づいて、抽出された前記特徴値から、当該脳波における嗜好を判断するステップと
を有することを特徴とする嗜好判別方法。
【請求項８】
脳波から、人間の生理的状態および／または刺激に対する人間の嗜好を判別するプログラムであって、
脳波を検出するステップと、
検出された前記脳波を周波数解析し、前記脳波に対するデータ行列を作成するステップと、
データ行列を多変量解析し、前記脳波の特徴値を抽出するステップと、
複数の特定の嗜好に対応した特徴値のデータ範囲に基づいて、抽出された前記特徴値から、当該脳波における嗜好および／または生理的状態を判断するステップと
を実行させるための嗜好判別プログラム。
【請求項９】
脳波から、人間の生理的状態および／または刺激に対する人間の嗜好を解析する脳波解析装置において、
脳波を検出する検出手段と、
前記検出手段により検出された脳波を周波数解析し、前記脳波のデータ行列を作成するデータ行列作成手段と、
前記データ行列作成手段で作成されたデータ行列を多変量解析し、前記脳波の特徴値を抽出する特徴値抽出手段と
を有することを特徴とする脳波解析装置。
【請求項１０】
脳波から、人間の生理的状態および／または刺激に対する人間の嗜好を解析する脳波解析方法において、
脳波を検出するステップと、
検出された前記脳波を周波数解析し、前記脳波に対するデータ行列を作成するステップと、
データ行列を多変量解析し、前記脳波の特徴値を抽出するステップと
を有することを特徴とする脳波解析方法。
【図２】
【図３】
【図４】
【図５Ａ】
【図５Ｂ】
【図６Ａ】
【図６Ｂ】
【図１】
【図７】
【図３】
【図４】
【図５Ａ】
【図５Ｂ】
【図６Ａ】
【図６Ｂ】
【図１】
【図７】
【公開番号】特開２０１０−１３１３２８（Ｐ２０１０−１３１３２８Ａ）
【公開日】平成２２年６月１７日（２０１０．６．１７）
【国際特許分類】
【出願番号】特願２００８−３１２６５１（Ｐ２００８−３１２６５１）
【出願日】平成２０年１２月８日（２００８．１２．８）
【出願人】（５０４１３２８８１）国立大学法人東京農工大学 (595)
【Ｆターム（参考）】
【公開日】平成２２年６月１７日（２０１０．６．１７）
【国際特許分類】
【出願日】平成２０年１２月８日（２００８．１２．８）
【出願人】（５０４１３２８８１）国立大学法人東京農工大学 (595)
【Ｆターム（参考）】
[ Back to top ]