知財求人 - 知財ポータルサイト「IP Force」
特表2020-502906位相アレイシステムのための駆動技術
この文献は図面が300枚以上あるため,図面を表示できません.
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公表特許公報(A)
(11)【公表番号】特表2020-502906(P2020-502906A)
(43)【公表日】2020年1月23日
(54)【発明の名称】位相アレイシステムのための駆動技術
(51)【国際特許分類】
H04R 3/12 20060101AFI20191220BHJP
G06F 3/01 20060101ALI20191220BHJP
【ＦＩ】
H04R3/12 Z
G06F3/01 560
【審査請求】未請求
【予備審査請求】未請求
【全頁数】93
(21)【出願番号】特願2019-530723(P2019-530723)
(86)(22)【出願日】2017年12月13日
(85)【翻訳文提出日】2019年8月6日
(86)【国際出願番号】GB2017053729
(87)【国際公開番号】WO2018109466
(87)【国際公開日】20180621
(31)【優先権主張番号】62/433,544
(32)【優先日】2016年12月13日
(33)【優先権主張国】US
(31)【優先権主張番号】62/436,457
(32)【優先日】2016年12月20日
(33)【優先権主張国】US
(31)【優先権主張番号】62/437,791
(32)【優先日】2016年12月22日
(33)【優先権主張国】US
(31)【優先権主張番号】62/438,512
(32)【優先日】2016年12月23日
(33)【優先権主張国】US
(31)【優先権主張番号】62/594,574
(32)【優先日】2017年12月5日
(33)【優先権主張国】US
(31)【優先権主張番号】62/594,687
(32)【優先日】2017年12月5日
(33)【優先権主張国】US
(31)【優先権主張番号】62/594,836
(32)【優先日】2017年12月5日
(33)【優先権主張国】US
(31)【優先権主張番号】62/594,874
(32)【優先日】2017年12月5日
(33)【優先権主張国】US
(81)【指定国】
AP(BW,GH,GM,KE,LR,LS,MW,MZ,NA,RW,SD,SL,ST,SZ,TZ,UG,ZM,ZW),EA(AM,AZ,BY,KG,KZ,RU,TJ,TM),EP(AL,AT,BE,BG,CH,CY,CZ,DE,DK,EE,ES,FI,FR,GB,GR,HR,HU,IE,IS,IT,LT,LU,LV,MC,MK,MT,NL,NO,PL,PT,RO,RS,SE,SI,SK,SM,TR),OA(BF,BJ,CF,CG,CI,CM,GA,GN,GQ,GW,KM,ML,MR,NE,SN,TD,TG),AE,AG,AL,AM,AO,AT,AU,AZ,BA,BB,BG,BH,BN,BR,BW,BY,BZ,CA,CH,CL,CN,CO,CR,CU,CZ,DE,DJ,DK,DM,DO,DZ,EC,EE,EG,ES,FI,GB,GD,GE,GH,GM,GT,HN,HR,HU,ID,IL,IN,IR,IS,JO,JP,KE,KG,KH,KN,KP,KR,KW,KZ,LA,LC,LK,LR,LS,LU,LY,MA,MD,ME,MG,MK,MN,MW,MX,MY,MZ,NA,NG,NI,NO,NZ,OM,PA,PE,PG,PH,PL,PT,QA,RO,RS,RU,RW,SA,SC,SD,SE,SG,SK,SL,SM,ST,SV,SY,TH,TJ,TM,TN,TR,TT
(71)【出願人】
【識別番号】517084885
【氏名又は名称】ウルトラハプティクス アイピー リミテッド
【氏名又は名称原語表記】ＵＬＴＲＡＨＡＰＴＩＣＳ ＩＰ ＬＴＤ
(74)【代理人】
【識別番号】100109634
【弁理士】
【氏名又は名称】舛谷 威志
(74)【代理人】
【識別番号】100129263
【弁理士】
【氏名又は名称】中尾 洋之
(72)【発明者】
【氏名】ロング，ベンジャミン ジョン オリバー
(72)【発明者】
【氏名】ジョーンズ，フィリップ エドウィンー
【テーマコード（参考）】
5D220
5E555
【Ｆターム（参考）】
5D220AA16
5D220AB06
5E555AA08
5E555AA11
5E555BA02
5E555BB02
5E555BC30
5E555DA24
5E555FA00
(57)【要約】
要約以下のシステムを含む、位相アレイシステムを駆動するための、多様な技術が記載されており、特に、空中触覚、パラメトリックオーディオ、音響浮揚および音響イメージングへの応用を有する音響位相アレイを対象とし：１）一貫した触覚体験を提供するために、空気の変化の影響を軽減することが可能であり；２）空気中でトラップポイントを生産し；３）より一貫した触覚効果の生産のためのベクトルに関して、位相アレイ最適化を定義し；４）コントロールされた音響場を介して、空間内の１つまたは複数のコントロールポイントまたは領域を定義し；５）音響基底関数の構築のために、被約表現法を使用し；６）大量のスループットに対して、複素数値関数の効率的な評価を実行し；７）マトリクスのクリロフ小領域を生成し；並びに８）音響基底関数を作成するために使用するものとは異なる、コントロールポイントおよび／または領域によって、記述される目的を最大にする。
[この文献は図面を表示できません]
【特許請求の範囲】
【請求項1】
音響場においてコントロールポイントを生成すること；
各コントロールポイントに法線ベクトルおよび方向を割り当て、位相アレイソルバーに波形運動量を組み込むこと；および
関心あるポイントで、前記波形運動量の
を含む方法。
各コントロールポイントに法線ベクトルおよび方向を割り当て、位相アレイソルバーに波形運動量を組み込むこと；および
関心あるポイントで、前記波形運動量の
[この文献は図面を表示できません]
成分を解くことを含む方法。
【請求項2】
[この文献は図面を表示できません]
成分が、下記式で定義される[この文献は図面を表示できません]
でそれぞれ特定され、[この文献は図面を表示できません]
上式中、[この文献は図面を表示できません]
は、コントロールポイント[この文献は図面を表示できません]
でのトランスデューサ[この文献は図面を表示できません]
からの進行波の方向を記述するベクトルであり、[この文献は図面を表示できません]
は、前記運動量が測定される位置を横断するコントロールポイント[この文献は図面を表示できません]
の法線ベクトルの方向である請求項１に記載の方法。
【発明の詳細な説明】
【技術分野】
【0001】
関連出願本出願は、以下の８つの米国仮特許出願の利益を主張するものであり、それらの全ては、参照によりその全体が本明細書に組み込まれる：
【0002】
1. 出願番号 ６２／４３３，５４４， 出願日２０１６年１２月１３日
【0003】
2. 出願番号 ６２／４３６，４５７， 出願日２０１６年１２月２０日
【0004】
3. 出願番号 ６２／４３７，７９１， 出願日２０１６年１２月２２日
【0005】
4. 出願番号 ６２／４３８，５１２， 出願日２０１６年１２月２３日
【0006】
5. 出願番号 ６２／５９４，５７４， 出願日２０１７年１２月 ５日
【0007】
6. 出願番号 ６２／５９４，６８７， 出願日２０１７年１２月 ５日
【0008】
7. 出願番号 ６２／５９４，８３６， 出願日２０１７年１２月 ５日
【0009】
8. 出願番号 ６２／５９４，８７４， 出願日２０１７年１２月 ５日
【0010】
開示分野本開示は、概ね、位相アレイシステムを駆動するために改良した技術に関し、主に空中ハプティックス、パラメトリックオーディオ、音響浮揚および音響イメージングを生成するために設計された音響場を生成することに関する。
【背景技術】
【0011】
背景我々が「音響場」と呼ぶ音エネルギーの連続分布は、空中での触覚フィードバックを含む、様々な用途に使用することができる。
【0012】
本音響場では、１つまたは複数のコントロールポイントを明確にすることができる。これらのコントロールポイントは、信号で振幅変調することができ、その結果、空中で振動触覚フィードバックを生成することができる。フィードバックを作成する別の方法は、振幅が変調していないコントロールポイントを作成し、それらを空間的に動かして感じられる時空間変調を生成することである。
【0013】
装置によって作成されるコントロールポイントの最小サイズおよび最大パワーは、生成される１つまたは複数のコントロールポイントの配列レイアウト、形状および相対位置に依存する。これにより、ある場所のコントロールポイントの動作を、別の場所の動作と一致させることが不可能な場合がある。コントロールポイントが空間を通って動的に動かされるとき、さらなる複雑さが起こり得る。 ユーザが経験した点の触覚機能の結果として生じる変化は望ましくなく、しばしば直観的ではなく、結果として混乱を招き、経験的な没入感を喪失する。空気の変化の影響を軽減することができ、一貫した触覚体験を提供するためにそれを逆転させることができるシステムは、商業的に価値があるだろう。
【0014】
さらに、一般的な非線形最適化システムを使用して、空気中にトラップポイントを生成できることが示された。しかし、この目的のために使用されるメカニズムは、位相アレイシステム上の浮揚に適した潜在的な場を評価しないこと、不必要に高価で時間がかかる最適化プロセスを使用すること、位相操作のみに限定される欠点を含む。
【0015】
既存のシステムによって作成したトラップは、それらが可能であるほど強力ではなく、これらの欠点のために正確に拡張することができないので、これらの問題に対処することは重要である。これは、音響的な浮揚を利用して、マルチモーダルアプローチを促進するための、さらなるビジュアルを提供する、システムおよび触覚システムの代替の使用を困難にする。
【0016】
さらに、空間で、１つまたは複数のコントロールポイントを明確にすることによって、音響場をコントロールすることができる。各ポイントには、コントロールポイントにおける所望の振幅に等しい値を割り当てることができる。次に、トランスデューサの物理的セットをコントロールして、コントロールポイントにおいて所望の振幅を示す音場を作り出すことができる。
【0017】
既存のシステムの条件によって、音響位相アレイの自然な指向性成分は、ほとんど注目されていなかった。例えば、目的がエネルギーを狭い領域に集中させることである高周波治療医療システムは、エネルギーを注入するための最適解を決定するためにスカラー圧力（関与する空気分子の数の変化）によって十分に近似され、したがってベクトル解を必要としない。位相アレイベクトルの挙動（空気分子の運動方向を表す）に焦点を合わせることは、ほとんどの場合無視できるので、これを効率的に達成することはまだ優先事項ではない。
【0018】
しかし、空中で触覚フィードバックを作成するには、人間の肌と振動を交換する必要がある。これは完全に波面との接触面と波によって運ばれる運動量に依存し、これは波のベクトルの振る舞いとの関係の結果である。これだけでなく、システムは、トランスデューサ素子間で入射角がますます異なるように焦点が合わせられているので、運動量と圧力との間の仮定された線形関係は崩壊する。このシステムではベクトル成分を無視することは不可能であり、スカラー圧力は空中触覚効果を生み出すための、適切な解を一貫して見つけるには不十分である。この理由から、この運動量の影響をデバイスに認識させることは、商業的に重要である。
【0019】
さらに、空間で１つまたは複数のコントロールポイントを明確にすることによって、音響場をコントロールすることができる。各ポイントには、コントロールポイントにおける所望の振幅に等しい値を割り当てることができる。次に、トランスデューサの物理的セットをコントロールして、コントロールポイントにおいて、所望の振幅を示す音響場を作り出すことができる。
【0020】
固有ベクトルの計算は、行列代数におけるいくつかの非常に具体的な純粋に二次の制約付き問題の解に対応することが示される。空中触覚における特定の問題を考慮して、固有ベクトルがトランスデューサのコントロール係数を解く行列を形成することができれば、その解決策は、既存の技術に対して単純化を示す装置実装の基礎を形成することができる。
【0021】
さらに、基底関数を記述するために簡約表現法を使用するためには、まず音響場の詳細に関して有用な関数の簡単な（低エントロピー）記述を可能にする簡略表現を開発しなければならない。
【0022】
本開示は、理論上、ユーザの視点からアクセス可能なものによるだろう。
【0023】
そのようなアクセス可能な低エントロピー表現を発見したならば、また、それをその等価なトランスデューサ基底（後の解決策において潜在的に多数の基底関数のうちの１つのトランスデューサ音響場の基底を記述するトランスデューサ基底セット）に変換する方法となる。これにより、見込みのあるユーザは音響場内で、多様な異なる構造を設計することが可能になり、焦点、ベッセルビーム、ボトルビーム、ベンディングビームなど、一度に多くの異なる効果の形成を導く現場でのエネルギーの流れを最大限に操作するための特注システムを構築することが可能になる。これらは、コントロールポイントの概念を、より豊かな音響場コントロール方法をカプセル化するコントロール領域に拡張するだろう。
【0024】
触覚における運動量移動を考慮するときに、粒子速度を波の方向の代用として、使用するために対処しなければならない、いくつかの問題があることもまた発見された。
【0025】
さらに、同じ２番目のオペランドを使用して効率的な並列複素数値乗算を実行することは、要件の資料に概説されているものと一致するパフォーマンスメトリクスを持つデバイスを実現するために必要な大量のスループットにとって必要だった。特に、必要な最適化行列の行の計算につながる合計の１つの要素を並列に計算することは、要件を満たすのに十分だが、サブシステムが十分な回数複製された場合に限られる。効率性に対するこの差し迫った必要性は、単純な乗算の上に複雑な超越関数を評価するための、新しいアプローチを生み出した、古いアルゴリズムに対して異なる見方をもたらした。
【0026】
このより広い世界では、この結果得られるアルゴリズムは、処理ハードウェアの要素が係数を無視できない複素数値のデータと効率的に、またはリアルタイムでやり取りしなければならない場合はどこでも適用可能である。これには、ＧＰＳ、ワイヤレス、ラジオおよび衛星送受信システムを含む通信、オーディオ処理、空中触覚を含む位相アレイシステム、飛行時間イメージングおよびデジタル干渉法、散乱効果の計算が含まれるが、これらに限定されない。これは、本発明にとってより広い範囲の機会につながる可能性がある。
【0027】
ＣＯＲＤＩＣとしても知られているＶｏｌｄｅｒのアルゴリズム（Ｖｏｌｄｅｒ、Ｊ．−ＴｈｅＣＯＲＤＩＣ計算手法−１９５９）は、単位複素数値を乗算することができるが、基本ＣＯＲＤＩＣアルゴリズムの計算を超えないと、計算の第２オペランドに係数を含めることができない。
【0028】
この理由で数年前に作者が実装した、１９９４年に発明された別のアルゴリズムは、作者Ｂａｊａｒｄ、ＫｌａおよびＭｕｌｌｅｒにちなんでＢＫＭとして知られており、さらに完全に複素数−複素数乗算の目標を達成することができる（Ｂａｊａｒｄ、Ｊ、Ｋｌａ、S．およびＭｕｌｌｅｒ、Ｊ.−ＢＫＭ：複雑な初等関数のための新しいハードウェアアルゴリズム−１９９４）。しかし、このアルゴリズムは複雑な実装を有するため、そうでない場合ほど多くの注目を集めていない。このアルゴリズムでは、複素数値のオペランドを複素数の対数に減らし、それらを加算してからべき乗化することで、完全な複素数−複素数乗算を実行できる。
【0029】
ＢＫＭと、場合によっては新しいアルゴリズムを使用して、シリアルとパラレルの両方で、この複素数値の乗算と複素数値の除算を実現するためのより効率的な方法がある。
【0030】
１９９９年のＢＫＭアルゴリズムの補遺は、複雑な対数の過度に複雑な実装を扱い、実装をいくらか単純化した。（Ｂａｊａｒｄ、Ｊ.およびＩｍｂｅｒｔ、Ｌ.−複雑な初期機能の評価：ＢＫＭの新バージョン−１９９９年）。１９９４年および１９９９年の補遺のこのＢＫＭ方式は、この文書では単にＢＫＭと呼ばれるアルゴリズムを形成する。
【0031】
実際のハードウェアに組み込むためにアルゴリズムを最適化しようとしている間に、説明したように数学関数を評価することよりも、我々の、またはより一般的には既存のシステムへの統合の容易さを優先する新しい技法のアイデアが明らかになった。この新しい技術は、それが概念実証コードに開発されたときに、元のアルゴリズムよりもハードウェア実装の一般的な採用にとってより有望なものとなる多くの新たな特性を持つことが判明した。加えて、元の技術に対して要求された正確な数値が要求される場合、これらの結果を検索するためにさらなるスケーリングが適用されてもよく、アルゴリズムは全体として競争力を維持したままである。どちらの場合でも、結果は、最先端技術を超える実質的な利点を備えた技術である。
【0032】
単位ノルムを有する二次問題は、その固有問題が二次問題と等価であるマトリクスを構築することによって解決され得ることが知られている。ランク１より大きい二次問題マトリクスは、多くのランク１問題マトリクスを合計すること、および／または既存の問題マトリクスのヌル空間を拡張することによって生成され得る。このマトリクスのクリロフ部分空間を生成することによって、暗示基底関数の最適化された部分集合を生成することができる。クリロフ部分空間を線形システムへの基底セットとして使用することによって、多数のテスト関数を使用して基底セットの正しいまたは近似評価を達成することができる。
【0033】
いくつの固有ベクトルが使用されるか、そしてそれらが選択される場所から、それに続く各固有ベクトルのパワーが減少するため、多くの基底関数の複雑なセットのグレースフルデグラデーションを達成することができる。
【0034】
さらに、音響場形状の結果を定義するための複素数値線形システムは、[この文献は図面を表示できません]
個の基底関数を取り、これらをキャリア周波数における音圧の所望の構成を記述する結果の場に再結合する。
【0035】
この時点までは、コントロールポイントを作動するために、コントロールポイント基底関数を使用したことがある。本明細書に記載しているのは、初期コントロールポイントまたは領域基底関数に異なる一組のコントロールテスト関数を適用するために一般化する方法である。
【0036】
これは、音響場における最初の「パレット」の振る舞いと、線形システムが使用する目的関数との間の関連性を断ち切るため、重要である。例えば、これを使用してベッセル関数の基底セットを作成し、それを使用してコントロールポイントを作動することができる。
【図面の簡単な説明】
【0037】
図面の簡単な説明添付の図面は、以下の詳細な説明と共に、別々の図を通して同一または機能的に類似の要素を指すが、明細書に組み込まれてその一部を形成し、特許請求の範囲に記載の発明を含む概念の実施形態をさらに例示するのに役立つ。これらの実施形態の様々な原理および利点を説明する。
【0038】
【0039】
【0040】
【0041】
【0042】
【0043】
【0044】
【0045】
【0046】
【0047】
【0048】
【0049】
【図31】３１は、所望波とは反対の方向から波に寄与するトランスデューサを除去する効果を示す。
【0050】
【0051】
【0052】
当業者であれば、図中の要素は簡単かつ明確にするために示されており、必ずしも一定の縮尺で描かれていないことを理解するだろう。例えば、図中の要素のいくつかの寸法は、本発明の実施形態の理解を向上させるのを助けるために、他の要素に対して誇張されている場合がある。
【0053】
装置および方法の構成要素は、必要に応じて図面中の従来の記号によって表しており、本発明の実施形態を理解することに関連する特定の詳細のみを示し、開示を曖昧にしないようにする。本明細書の記載の恩恵を受ける当業者は、当業者に自明である。
【発明を実施するための形態】
【0054】
詳細な説明前述の技術、装置および方法は、組み合わせるか、または別々に空中触覚システムで使用できる。
【0055】
１．空中ハプティックデバイスのための動的範囲縮小
【0056】
Ｉ.心理物理学的モデリング
【0057】
システムの振る舞いと触覚出力は、人間が接触を知覚する方法に合わせて、調整しなければならない。例えば、コントロールポイントのサイズが大きくなると、コントロールポイントで作成する物理的な振幅は変化しないが、触覚の知覚も増加する。音声などによるクロスモーダル効果は、それらが相関する場合、触覚フィードバックを向上させるのに役立つ。これは、物理的に定義されたコントロールポイントの間の食い違いを招き、同等の認識を達成するために心理物理学的モデリングステップの使用を必要とする。
【0058】
ＩＩ．心理ハプティックモデルを使った再マッピング
【0059】
再生装置の物理的制限のために、正確な波形は、指定したように再現可能ではないかもしれない。波形は、心理的に定義された知覚空間内に構築され得る。あるいは、それは、実質的に絶対的な物理的スケールでの強度および力の測定値から構成されてもよく、それは触覚テクスチャ、またはさらなるメタデータと共に潜在的に表現される任意の構築波形の記録であり得る。この波形は、「完璧な」空中の触覚デバイスが与えられたときの、コントロールポイントの動作を表す。次に、これは、コントロールポイントの再現の知覚的側面を記述および定量化できる、心理物理学的モデルを使用して知覚的に評価される。この触覚効果を再現しながら、コントロールポイントが動いている可能性がある。運動は、触覚効果の知覚的評価の一部として含まれる。次に、これは、ユーザ位置、相対的幾何学的形状、および利用可能なパワー制限によって制約される空中触覚システムのインスタンスに利用可能な出力の空間に再解釈することができ、システムは実質的に等価なコントロールポイントを作成するように構成する。
【0060】
心理触覚モデルは、完全な装置と実質的に同じ触覚応答を引き起こすために、コントロールポイントの変調波形をどのように修正できるかを決定する能力を、デバイスに与える。この再マッピングを達成するプロセスはまた、結果として生じるハプティック波形が雑音を多くすることを防ぐため、または知覚ハプティック空間に関しては近いがノイズが少ない、より可聴的に許容可能な変調を選択するための心理音響モデリングステップを含むことができる。可聴ノイズの量を増やすことによって、多くの波形をより触覚的に強力にすることができることが知られている。いつかの状況では、他ではこれが望ましいかもしれないが、これは耐えられないほど十分に騒々しいかもしれない。次に、音響心理学的モデルを併用して、音声をハプティックトレードオフ設定にして、出力をバイアスしてより静かまたはより強くすることができる。
【0061】
ＩＩＩ．精神的触覚ユニット
【0062】
振幅が同等であると感じるように意図する心理空間で波形を定義する場合、空中で生成されたコントロールポイントの心理効果を考慮に入れるユニットを設計しなければならない。これらの中で最も注目に値するのは、手に触覚が統合されていることである。各コントロールポイントに含まれる焦点領域は、物理的装置の特性に応じてサイズおよび形状を再び変化させ、そしてこれは、精神触覚ユニットの構成において相殺しなければならない。そのような精神触覚ユニットは、焦点領域のサイズおよび形状に応じて評価する物理的コントロールポイント振幅を変更し、広い面積にわたるタッチ感覚の統合の影響を相殺する。このアプローチは、同じ数の精神触覚ユニットが同等の感覚をコード化することを可能にし、デバイス間の触覚的同等性および物理的インスタンス化を維持する。
【0063】
ＩＶ．動的な波形のマッピング
【0064】
ターゲットデバイスにベストなマッピングを到達させるために、 完全な波形とその軌跡が考慮されるかもしれない。しかし、対話性、したがって短い待ち時間が必要条件であるため、波形の一部または個々のサンプルだけが利用可能である。波形全体をターゲットデバイスの機能に再マッピングすることができるが、これを行うには、波形を事前に知っておく必要があるため、周知の波形の前処理ステップと見なすことができる。代替方法は、振幅波形サンプルの物理的特性を振幅で再マッピングすることである。なぜなら、振幅サンプルの再マッピングは対話型波形を用いてリアルタイムで達成することができ、リアルタイムで触覚フィードバックを提供するのにより適しているからである。
【0065】
多様な関数が、コントロールポイント振幅サンプルのための再マッピングを提供するために、良い候補であると考えられる。これらは、対数線形曲線と共にシグモイド関数を含み、それは他のモダリティにおいて低いダイナミックレンジを有する装置上で、知覚的な感覚入力を納得のいくように再現することができる。これらの曲線は一般にマッピング関数における上限の指定を必要としないが、図１および図２に示すように、利用可能な有限量のパワーがあることを考えると、ダイナミックレンジはより効果的に使用することができる。図１は、任意のダイナミックレンジ１００のハプティック波形を示しており、これは非線形再マッピング機能１１０を介して処理され、コントロールしたダイナミックレンジ１２０の新しい波形を生成することができる。これは任意のダイナミックレンジのハプティック波形を取得し、それをより小さなレンジに変換するために使用できる。波形の振幅は無限大である可能性があるため、これは出力のダイナミックレンジ全体を埋めるわけではなく、未使用の出力範囲になる。図２では、図１と同様に、非線形関数２１０は、振幅波形２００をターゲット装置２２０に利用可能な振幅に再マッピングする。波の広がりを知ることで、再マッピング機能をデバイスに固有の出力範囲に、正確に合わせることができる。これは、出力デバイスの全てのレンジの開発を可能にする。
【0066】
説明のために、以下の例では、空中触覚デバイスにおけるマッピングがいかに望ましくない挙動を防ぐのかを説明するために、ミニマリスト線形スケーリングを使用する。
【0067】
図３に示すグラフ３００は、波形の再現可能な知覚への再マッピングを達成するための第１のステップが、波形およびその限界３００を記述するソルバー関数にデータを供給することである。図３では、実線は[この文献は図面を表示できません]
であり、点線は、[この文献は図面を表示できません]
である。それは、関数を再マッピングする単純な波形の入力である。ここでｙ軸は、デバイスの内部単位である。
【0068】
これは、クリッピングなしで、波形を放射するために再生する必要がある最大振幅を定義することによって達成され、これは、各サンプル時間で、[この文献は図面を表示できません]
と表される。これは、入力波形の振幅[この文献は図面を表示できません]
と共に、入力としてソルバーに入力される。
【0069】
図４に示されているのは、入力データにいかなる変更も適用せずにこの波形を再生しようと試みた結果を示すグラフ４００である。その結果、コントロールポイントで必要な振幅を再現するために利用できるパワーが不足しているために波形がクリップされる。図４では、直線は[この文献は図面を表示できません]
であり、それは振幅[この文献は図面を表示できません]
の結果であり、これは、デバイスのダイナミックレンジに関係なく、図１のユーザ指定振幅を直接出力しようとした結果である。
【0070】
システム内のアクティブトランスデューサの物理的構成を、[この文献は図面を表示できません]
位置[この文献は図面を表示できません]
における所望のコントロールポイント振動を、[この文献は図面を表示できません]
として定義すると、空中触覚システムに対するソルバー関数の補助機能を使用して、関数振動[この文献は図面を表示できません]
を定義することができ、与えられた装置が与えられた時間に、これを物理的に提供することができる場合には、コントロールポイントにおける所望の搬送波振幅、あるいはコントロールポイントの位置において再生できる最も近い振幅のいずれかをもたらす。
【0071】
[この文献は図面を表示できません]
を評価して、この時点で入力振幅をマッピングできる範囲の最大値を求める。コントロールポイントはサンプル間であまり移動しないので、これは正確である必要はなく、コントロールポイントの前の反復からの[この文献は図面を表示できません]
の値が使用できる。丸め誤差および他の雑音源を除去するために、この値をフィルタリングすることもできる。この件では、再マッピング関数は波形の単一サンプルを取り、これを区間[この文献は図面を表示できません]
から区間[この文献は図面を表示できません]
に再マッピングする。図５に示すグラフ５００は、関数５００の再マッピングを通して、[この文献は図面を表示できません]
から[この文献は図面を表示できません]
へ直接マップを産出する。振幅を[この文献は図面を表示できません]
から[この文献は図面を表示できません]
から、値[この文献は図面を表示できません]
で得られる範囲に再マッピングすることで、波形のクリッピングを防止すると、[この文献は図面を表示できません]
となる。
【0072】
Ｖ．複数の同時コントロールポイントの再マッピング
【0073】
図６の２つの正弦波の場合の例のように、ソルバーが波形を正しく表すことができるという保証を提供することができるように、一度に複数のコントロールポイントに対してシステムを使用するには、範囲を同時に評価しなければならない。グラフ６００が示すように、２つのコントロールポイントは、同時に異なるポイントで、異なる振幅[この文献は図面を表示できません]
および[この文献は図面を表示できません]
で生成される。一方はトランスデューサアレイに対して接近して移動し、他方はトランスデューサアレイに対して横方向に移動する。それらは、異なる周波数の２つの正弦波で、振幅が振動するように設定されている。
【0074】
デバイスのパワー制限に違反すると、パワーが不足してクリップされた波形になり、デバイスの制限を超える複数のコントロールポイントがポイント間のクロストークを促進し、それらの触覚効果を混乱させる単一のコントロールポイントとは異なる。図７に示すグラフ７００は、２つの動くコントロールポイントの例を示す。図７は、トランスデューサアレイを使用して２つのコントロールポイントを再作成しようとする試みを示しており、２つのポイント間で異なる共振が発生して置き換えられると、複雑で予測不能な動作になる。
【0075】
図８に示すのは、両方の点に存在するパワー制限を決定するために補助ソルバー評価を使用するグラフ８００であり、示された波形を導き、それは次のように記述することができる。[この文献は図面を表示できません]
ここで、一方のコントロールポイントがより強くなる可能性を持ってデバイスに向かって移動し、他方のコントロールポイントはデバイスの表面を横切って弱くなる可能性を持って離れていく。図８からわかるように、ソルバーソリューションによって生成するダイナミックレンジの制限により、それらの波形は意図したとおりにクロストークなしで表される。したがって、図８は、振幅エンベロープが、エラーなしで、またはコントロールされた周知のエラーレベルで再現されるようにソルバーによって保証されるレベルまで低減された２つのポイントを示している。これは理想的な振幅レベルで振幅誤差の２乗を最小にする構成であるため、それらは同じパワーである。
【0076】
ＶＩ．複数の同時コントロールポイントにわたるパワー再分配
【0077】
振幅範囲は、異なる値に設定でき、パワーレベルを管理するために時間を通して動的に設定できる。この方法では、範囲を操作すると、特定のコントロールポイントに集中している配列出力の比率が決定できる。これは、図９および図１０では、あるコントロールポイントから別のコントロールポイントへのパワーの単純な操作で示す。図１１および図１２で示されるような、コントロールポイント間の動作ダイナミクスを可能にするより複雑な位相外のずれ関係もある。
【0078】
図９に示しているのは、２つの入力コントロールポイントに対する入力波形および範囲を有するグラフ９００である。これらは波形範囲が異なり、各ポイントがデバイスから利用可能なパワーをより多く、またはより少なくできるようにするためのヒントとしても使用される。図１０に示すのは、図９の入力から得られた出力を有するグラフ１０００である。この構成では、アレイによって生成したパワーは、それらがそれぞれアレイから遠ざかっている間およびアレイに向かって動いている間に、一方のコントロールポイントから他方のコントロールポイントに移動する。これにより、図６、図７、図８で予想したように、コントロールポイントの振幅が等しいと予想される場合よりも多くのパワーを使用できる。
【0079】
図１１に示すのは、交互の関数を有するグラフ１１００であり、１つのコントロールポイントまたはコントロールポイント群は、範囲が正弦波に一致することによって効果的にモデル化できる。この結果、各波形は変調システムの各サイドのパワー要件に従って再生し、システムの各サイドはフルアレイパワーにアクセスできることが保証される。その結果、任意の出力クリッピングが発生する可能性なしに、タイムスライスが行われる。図１２に示すのは、図１１の入力から得られた出力を有するが、説明されているグループ分けなしのグラフ１２００である。これにより、パワーバジェットがすべてのコントロールポイント間で不必要に共有しているため、出力ダイナミックレンジが減少する。
【0080】
２．音響浮上のためのトラップポイントの効率的な生成
【0081】
Ｉ．浮遊球状粒子のための最適な目的
【0082】
粒子を浮揚するために必要とする最適化問題は、ゴルコフエネルギーポテンシャルの使用を中心とする。これは、音響場内の粒子の可能な各位置のエネルギーを評価する、スカラー関数である。空間位置で最小または低エネルギー状態を見つけることは、潜在的に安定なトラップポイントを見つけることに対応する。しかし、ゴルコフポテンシャルが考慮に入れていない音響効果、例えば音響波が流体中でバルク運動を引き起こす音響流がある。この理由から、エネルギーポテンシャルの実行可能な最小値であると決定されることが多い変動圧力の領域で、トラップが発生しないことを確実にすることも好ましく、可能性の定義に含まれる。
【0083】
重力および音響流効果のような外力の方向および大きさは、ゴルコフポテンシャルのみからは導き出すことができないので、ゴルコフスカラーポテンシャルのラプラシアンを計算することがしばしば好ましい。これにより、内側に押す力のより安定した配置を生成するために、どの軸を他の軸より優先するかを、さらにコントロールできる。
【0084】
ＩＩ．ゴルコフスポテンシャルの評価
【0085】
極小エネルギーの存在により、球状粒子が浮遊することを可能にするエネルギー場であるゴルコフポテンシャル分野が、（Ｂｒｕｕｓ ２０１２年）として記載される：[この文献は図面を表示できません]
ここで、[この文献は図面を表示できません]
は、音響媒体を通した音のスピードであり、[この文献は図面を表示できません]
は、球状粒子を通した音速であり、[この文献は図面を表示できません]
は、音響媒体の密度であり、[この文献は図面を表示できません]
は、球状粒子の密度であって、時間的平均圧力[この文献は図面を表示できません]
を含む左辺は潜在的なエネルギーの表示であり、粒子速度ベクトル [この文献は図面を表示できません]
を含む右辺は潜在的なエネルギーに保守的に関係する。ゴルコフポテンシャルは、活動的または遠距離場に存在するという仮定をするため、これらの圧力と速度の項は同相の時間調和項と交換することができる。これは、時間調和ゴルコフポテンシャル分野が、次のように書き換えられることを示唆する：[この文献は図面を表示できません]
【0086】
これを多くの情報源に拡張すると、次のようになる：[この文献は図面を表示できません]
【0087】
これに関する問題は、交差項が事実上、これが次のようになることを意味する：[この文献は図面を表示できません]
その結果、二次であり、二乗したトランスデューサの数が増えることとなる。
【0088】
ＩＩＩ．固有問題によるポテンシャルの最適化
【0089】
時間高調波効果を無視し、各トランスデューサに入力する圧力および粒子速度を考慮した潜在的な定義は、次のように書くことができる：[この文献は図面を表示できません]
[この文献は図面を表示できません]
[この文献は図面を表示できません]
ここで、[この文献は図面を表示できません]
は所望のトラップのインデックスであり、余分な[この文献は図面を表示できません]
は、トランデューサ[この文献は図面を表示できません]
および[この文献は図面を表示できません]
からトラップポイント[この文献は図面を表示できません]
へ移動する際に、位相および振幅が変化する位相および振幅である。トランスデューサが非ゼロレベルにパワーを供給することを確実にするために、制約[この文献は図面を表示できません]
を使用して、ポテンシャルウェルを達成するための最適複素トランスデューサ活性化係数は、最大固有値を有する[この文献は図面を表示できません]
個の固有ベクトルによって与えられる。
【0090】
[この文献は図面を表示できません]
をこのように選択すると、目的関数[この文献は図面を表示できません]
が最大化され、これは、[この文献は図面を表示できません]
の定義の負の符号のため、ゴルゴフポテンシャル分野で、最小値が得られる。なぜなら、[この文献は図面を表示できません]
が固有ベクトルの定義であり、[この文献は図面を表示できません]
であることが定義されており、定義により、[この文献は図面を表示できません]
がマトリクスの固有値であり、[この文献は図面を表示できません]
なので、最大化され、これは、存在するのであれば、これは最大値をもたらさなければならない。問題の標準的な記述は、以下のようになる：[この文献は図面を表示できません]
【0091】
しかし、以前に文献で議論したように、ゴルコフ分野のエネルギーを最小にすることは焦点を生成し、実際にはトラップポイントではない。トラップポイントを生成することは、トラップでの圧力変動を最小化することである。圧力変動の定義は複素数値であり、意味のあるように最小化することはできないが、代わりにゼロにする変動量を表すため、これは自明ではない。
【0092】
ＩＶ． 固有値問題内に、ヌル空間を構築して圧力ゼロの生成
【0093】
ゼロ圧力条件は、音響場内の所与の点で、ゼロ圧力に評価される一連の制約として構成することができる。高次元位相空間内に、周知の最適圧力集束解に平行な解ベクトルの一部が存在する場合、圧力は所与の位置でゼロ以外の値に評価することしかできない。トラップポイント[この文献は図面を表示できません]
では、ベクトル[この文献は図面を表示できません]
があり、それは以下のようになる：[この文献は図面を表示できません]
【0094】
それゆえに、トラップポイント[この文献は図面を表示できません]
におけるゼロ圧力に対する制約は、制約式として記述できる：[この文献は図面を表示できません]
【0095】
振幅が、これに対するものであるかどうかは重要ではなく、制約は単位長さである必要があるので、それゆえ以下のようになる：[この文献は図面を表示できません]
【0096】
問題の標準的な記載は、以下のとおりである：[この文献は図面を表示できません]
【0097】
固有値系をＧｏｌｕｂ（１９７３）から導き出し、それを複素数値について言い換えると、以下の議論をもたらす：
【0098】
[この文献は図面を表示できません]
に伴う上記のシステムを挙げると、システムの特性方程式は次のように書くことができる:[この文献は図面を表示できません]
ここで、[この文献は図面を表示できません]
および[この文献は図面を表示できません]
は、ラグランジュ乗数である。これは、[この文献は図面を表示できません]
ベクトルの共役によって微分することができ、システムが得られる：[この文献は図面を表示できません]
【0099】
[この文献は図面を表示できません]
を乗ずると、以下が得られる：[この文献は図面を表示できません]
【0100】
[この文献は図面を表示できません]
が単位長であると仮定すると、[この文献は図面を表示できません]
となり、先の問題定義では[この文献は図面を表示できません]
となり、両方を代入すると以下のようになる：[この文献は図面を表示できません]
【0101】
これを以前の導関数に代入して置換すると以下が得られる：[この文献は図面を表示できません]
【0102】
これは以下のとおりである:[この文献は図面を表示できません]
【0103】
これは固有システムの記述であり、未だ[この文献は図面を表示できません]
は対称ではないため、固有ベクトルの計算は困難である。しかしながら、[この文献は図面を表示できません]
は、べき等[この文献は図面を表示できません]
に示され、固有値は、順序付けによって変わらないので以下が得られる:[この文献は図面を表示できません]
【0104】
これを、より多くの制約ベクトルに拡張しても、左辺をさらに乗算しても問題はない：[この文献は図面を表示できません]
多くの乗算の後、最終的に以下が得られる：[この文献は図面を表示できません]
ここで[この文献は図面を表示できません]
は、ムーア−ペンローズ擬似逆行列を表す。
【0105】
要約すると、初期問題の解は、[この文献は図面を表示できません]
の最大固有値を持つ固有ベクトルを求めることによって見つけることができる。[この文献は図面を表示できません]
の乗算は、解決マトリクスのヌル空間に、基底ベクトルを追加するランク１の変更を追加すると見なすことができる。
【0106】
最終的に、簡単な方法で、ゴルコフポテンシャルのラプラシアンを計算するための、サポートを追加することは有用である。ここで見いだされたゴルコフポテンシャル解はトラップポイントを生み出すが、どの方向から力が加えられるのかをコントロールすることはできず、力がない方向でさえあるかもしれない。これを達成するためには、潜在分野の勾配における最も発散する点に関して最適化し、次いでこの場の勾配を再重み付けする必要がある。これは、一般的に言って、ゴルコフ場のラプラシアンである。
【0107】
Ｖ. ゴルコフ場のラプラシアンを生成するための修正
【0108】
ゴルコフ場のラプラシアンを最小化するための修正は、[この文献は図面を表示できません]
マトリクスのみに適応する必要がある。各トランスデューサからの波面曲率の影響を破棄し、それらを平面波のみと見なすと、定式化を複雑にする特異点が取り除かれ、粒子速度は一定になる。各トランスデューサからの平面波の結果は、以下の式に書くことができる：[この文献は図面を表示できません]
【0109】
マトリクスは、ラプラシアンケース以外では、同等のシステムとして表わされる：[この文献は図面を表示できません]
[この文献は図面を表示できません]
の全ての制約は、[この文献は図面を表示できません]
に移行される。この定式化のラプラシアンを取ると、修正マトリクス定式化が得られる：
[この文献は図面を表示できません]
【0110】
[この文献は図面を表示できません]
マトリクスの対応する要素の代わりにこの値を使用すると、ゴルコフポテンシャルのラプラシアンを最適化するという点で、解決策が得られる。式[この文献は図面を表示できません]
は、以下のように記述することができる：[この文献は図面を表示できません]
【0111】
これは、また、[この文献は図面を表示できません]
方向からの寄与に分割することもできる：[この文献は図面を表示できません]
【0112】
この方法では、重み付けを適用して拘束軸を生成し、各方向から生成された力にバイアスをかけることができる。
【0113】
ＶＩ．デバイスでの実装
【0114】
本明細書におけるシステムは、単純な固有システムソルバーに基づいて、１つまたは複数の浮揚トラップを計算するために使用することができる装置を表す。これは、アイテムを拘束し、トラップし、輸送するために音響浮揚を使用することができる実施形態を作成するために、デバイスファームウェアで潜在的に実施することができる。
【0115】
これが、図１３〜１６に示される。
【0116】
図１３は、同方向の重み付けを伴う、ゴルコフのラプラシアンを示す。左側面１３１０は、トランスデューサアレイに平行なｘ−ｙ平面である。右側面１３２０は、平面トランスデューサアレイに垂直であり、且つそれを通ってスライスするｘ−ｚ平面である。
【0117】
図１４は、重く重み付けされたｚ方向を有するゴルコフのラプラシアンを示す。右側面１４２０は、平面トランスデューサアレイに垂直であり、且つそれを通ってスライスするｘ−ｚ平面である。
【0118】
図１５は、ゼロ圧力状態のゴルコフ最大化を示す。左側面１５１０は、トランスデューサアレイに平行なｘ−ｙ平面であり、それは、一様な強度の輪の周りの不均一な陰影によって示される、らせん状のトラップを記載する。右側面１５２０は、その下の平面トランスデューサアレイに垂直であり、且つそれを通ってスライスするｘ-ｚ平面である。
【0119】
図１６において、さらなる固有モード１６１０、１６２０、１６３０は、やはりコア問題を解決するトラップ解に対応する。これらは、目的関数を最大化しないため、あまり効果的ではない。これは、目的関数の値である修正ゴルコフポテンシャルが、主固有ベクトルに直交する、これらのモードでは小さいからである。達成した目的関数の値は、各固有ベクトルに関連する固有値に反映される。小さい固有値に関連する固有ベクトルは、これらの代替的で効果の低い解を表す。
【0120】
３．空中触覚に適用したベクトル効果に焦点を当てた位相アレイ装置
【0121】
Ｉ．活性音響強度
【0122】
活性音響強度[この文献は図面を表示できません]
は、音響媒体の圧力[この文献は図面を表示できません]
および粒子速度ベクトル[この文献は図面を表示できません]
を用いて定義することができる。音響波は、以下のように記述できる：[この文献は図面を表示できません]
ここで、活性強度は、純粋に実数のベクトルである。波源の近距離場では、[この文献は図面を表示できません]
ではなく、反応場が支配的である。しかし、遠距離場では、１または２つの波長[この文献は図面を表示できません]
を超える。この理由から、[この文献は図面を表示できません]
のみが考慮され、現実感が保証される。
【0123】
下記の単調平面波源について考察する：[この文献は図面を表示できません]
その結果下記が得られる：[この文献は図面を表示できません]
【0124】
ここで、[この文献は図面を表示できません]
は圧力振幅であり、[この文献は図面を表示できません]
は、各方向の粒子速度振幅である。
【0125】
所与の軸に沿って波によって伝達される音の強度、したがってエネルギー束は重要であるので、スカラー寄与を測定するための内積は、以下のようになる：[この文献は図面を表示できません]
【0126】
音波の方程式の形に起因して以下のようになる：[この文献は図面を表示できません]
重ね合わせの原理は、圧力ベクトルと粒子速度ベクトルの両方に当てはまる。[この文献は図面を表示できません]
個の音源を想定し、それぞれを対象の点における平面波として近似すると、次式が得られる：[この文献は図面を表示できません]
効率のためには、線形システムを用いて適切な変数について解くことが望ましい。これは、２つの独立した変数で構成しているため、与えられた方向の音響強度について解くことができないことを意味する。既存の技術は複素数値のスカラー圧力[この文献は図面を表示できません]
解き、それから個々の波のそれぞれの方向は重要ではないか、または、ほぼ平行であると仮定する。代わりに、これは運動量（および運動エネルギー）に直接影響するので、粒子速度ベクトル[この文献は図面を表示できません]
を解き、これは、空中での触覚フィードバックの用途にとって重要である。
【0127】
このことは、波の指向性とその結果生じるエネルギー伝達を反映する解を生成し、複素数値の粒子速度ベクトル[この文献は図面を表示できません]
を利用するために現在の、解のメカニズムを拡張する２つの潜在的な技術につながる。
【0128】
ＩＩ．複素数値に対する解法
【0129】
複素数値に対する解法は、値に対する音響モデルを使用することを含む。このモデルは、コントロールポイントの位置と並んでトランスデューサの位置が与えられて、フィールド内のその点で得られたフェーザが決定されると、照会することができる。これは複素数値であり、[この文献は図面を表示できません]
ここで[この文献は図面を表示できません]
は、音響場内のサンプルポイントのインデックス、[この文献は図面を表示できません]
は、作動トランスデューサのインデックス番号に対応する:[この文献は図面を表示できません]
ここで、[この文献は図面を表示できません]
は、[この文献は図面を表示できません]
の場所で空中に意図したフェーザであり、[この文献は図面を表示できません]
は、トランスデューサ[この文献は図面を表示できません]
からのトランスデューサ活性化係数、空気中での挙動を生成するのに必要な初期位相および振幅が[この文献は図面を表示できません]
で特定される。
【0130】
しかしながら、[この文献は図面を表示できません]
は複素数値であり、位相と振幅仕様の両方を含むため、位相オフセットは自由に修正できる。これを利用して、振幅の最大の可能な部分空間を開く位相値を見つけることができる。これは、システムの主固有ベクトルｙの発見によって達成される。[この文献は図面を表示できません]
【0131】
[この文献は図面を表示できません]
は、音響場[この文献は図面を表示できません]
内の各コントロールポイントについて、音響場の各ポイントで所望される振幅を表し、[この文献は図面を表示できません]
は、内積をマトリックスの各対角内の１にキャンセルするために適用される列当たりの規格化値である。この行列の固有ベクトル[この文献は図面を表示できません]
の正規化された成分は、初期振幅[この文献は図面を表示できません]
と乗算されたときに、場内の圧力点の各々の間で強化を促進するフェーザベクトル[この文献は図面を表示できません]
に対応する。
【0132】
未決定問題の最小ノルム解は、次の形を用いて解くことができる：[この文献は図面を表示できません]
【0133】
この解決に重要なのは、マトリクス[この文献は図面を表示できません]
の逆行列であり、これはコレスキー分解によってシステムを解く最初のステップであり、次のように記述できる：[この文献は図面を表示できません]
これが固有システムのマトリクス[この文献は図面を表示できません]
の導出にも役立つことを意味する。しかし、次の形式でシステムを操作することもできる：[この文献は図面を表示できません]
これは非常に小さな線形システムとそれに続く解のサイズをもたらし、これは任意の手法（ＱＲ分解、コレスキー分解など）で解くことができ、放出されるまでこの小さな状態に保たれる。このＣマトリクスおよびこの[この文献は図面を表示できません]
ベクトルはどちらも、計算されるまで、トランスデューサにインデックスを持たない：[この文献は図面を表示できません]
【0134】
トランスデューサは多数ある可能性があるため、最後の可能性がある時点で[この文献は図面を表示できません]
を乗算して各トランスデューサから放射される波形を乗算し、[この文献は図面を表示できません]
トランスデューサ アクティベーションに変換し直すことがしばしば価値がある。
【0135】
ＩＩＩ．１つまたは複数の異方性コントロールポイントの解決
【0136】
それらの各々に沿って指定された法線ベクトルまたは方向を有するように音響場内にコントロールポイントを作成することは、運動量要件を位相アレイソルバーに組み込むための１つの方法である。各コントロールポイントに法線ベクトルを与えることは、各コントロールポイントにおける圧力について解くのではなく、対象となる点におけるｘ，ｙおよびｚ粒子速度、したがって波の運動量の成分について効果的に解くことを意味する。それを達成するために、線形システム方程式を前から展開して、各粒子速度について３つの成分を、前に定義した圧力については１つの代わりに含める。それらの３つ成分は、[この文献は図面を表示できません]
として記載され、以下のように定義される：[この文献は図面を表示できません]
各ケースで、[この文献は図面を表示できません]
は、コントロールポイント[この文献は図面を表示できません]
におけるトランデューサ[この文献は図面を表示できません]
からの伝搬方向を示すベクトルであり、その位置を横切って運動量が測定される[この文献は図面を表示できません]
は、コントロールポイント[この文献は図面を表示できません]
の法線ベクトル方向である。解決のための核となる線形システムは、以下のようになる：[この文献は図面を表示できません]
これは、各コントロールポイントの音響媒体の粒子の速度の３つの成分を示す。しかし、システムの解法に簡約表現アプローチを使用する場合、つまり必要な作業量を減らすために、Ｃマトリクスを指定する場合、このマトリクスのエントリ数は９倍に増加し、計算するのにより強大な問題となる。先にマトリクスＣの各要素[この文献は図面を表示できません]
をサブマトリクスとする：[この文献は図面を表示できません]
【0137】
これは、固有値マトリクス[この文献は図面を表示できません]
の既存のアプローチを使用した場合にもうまく動作し、同じアプローチを使用して構成することができるが、ユーザ入力の振幅は[この文献は図面を表示できません]
および[この文献は図面を表示できません]
となる。固有システムマトリクスのクロスタームは、別の軸に別の軸の波形を追加した場合の影響を、さらに別のコントロールポイントに考慮したときに、１つのコントロールポイントの位置で、１つの軸で測定を実行していると看做すことができる。[この文献は図面を表示できません]
はここで、成分を３倍多く含むが、[この文献は図面を表示できません]
による乗算は、装置上の空気中で所望の運動量ベクトルを再現するのに必要なトランスデューサ活性化係数にシステムを戻すことができる。
【0138】
この解決策は、空中触覚フィードバックの適用にとって特に興味深い。なぜなら、異なる角度から手に近づく波を作り出す可能性が開かれるからである。これにより、音響波が境界上の皮膚にせん断波を誘導することが可能になり、これは、余分な必要な計算を犠牲にして、空中触覚技術を使用して作り出すことができる全く新しい一連の感覚を開くものである。
【0139】
ＩＶ．音響媒体内の１つまたは複数のコントロールポイントでの粒子速度の解法
【0140】
圧力値を目標とするソルバーから運動量を目標とすることができるものへの変換を達成するための別の技術は、（元の圧力アプローチと同様に）音響粒子速度と完全に連動する線形システムソルバーを生成することであり、また、単調な調和時間のための複素数値である。さらに、これはスカラーである。これは、コントロールポイントにおける最大粒子速度の方向、したがって結果として生じる波の方向[この文献は図面を表示できません]
を計算し、次に乗算して各マトリクスエントリを単一の複素数値のスカラーにすることによって達成することができる。ドット積は、[この文献は図面を表示できません]
にわたって等しくスケールファクタによってすべてのサンプリング方向に影響を与えるので、これは初期マトリクス構成に対する後処理として適用することができる。しかし、結果として得られる波の方向は、ソルバーが機能する前に知る必要があるため、ソルバーの影響を受けない。空中触覚デバイスの構築に適用する場合、このアプローチは、空中であるが比較的近距離のフィードバックが最良であり、現在の解決メカニズムに対する拡張である小型デバイスに最も適しているであろう。
【0141】
これを達成するために、各トランスデューサに対する音響モデルは、音響媒体中の粒子速度を決定するための問題として、再キャストされる。このように、アルファベクトルは、与えられたコントロールポイントについて解く必要があるシステムとして、最大の運動方向に沿って測定される複素数値の粒子速度[この文献は図面を表示できません]
を生じるように再定義し、所与のコントロールポイントの解決を必要とするシステムは、式[この文献は図面を表示できません]
となる。しかし、不幸なことに、[この文献は図面を表示できません]
は、アプリオリに知ることができないが、線形システムマトリクスを構成するのに必要である：[この文献は図面を表示できません]
【0142】
したがって、音響媒体は、値を生成するために修正できる：[この文献は図面を表示できません]
ここで、下記式のとおりである。[この文献は図面を表示できません]
【0143】
これらの値は、ここで、[この文献は図面を表示できません]
マトリクスのドット積を構成するために、使用されなければならない：[この文献は図面を表示できません]
【0144】
しかしながら、それらは、これを達成するための２つの可能性のあるルートがある。その必要なドット積は、以下で定義される ：[この文献は図面を表示できません]
【0145】
これは、以下のように記述することができる：[この文献は図面を表示できません]
ここで、[この文献は図面を表示できません]
はダイアディックプロダクトである。これは、後で法線ベクトルを計算するために、マトリクス要素ごとに９つの別々の総和（対角線の場合は３つ）を維持する必要がある。これは、異方性コントロールポイント点アプローチのそれと多くの類似点を有する[この文献は図面を表示できません]
のための初期マトリクス構成をもたらすという点で有益であり、したがって、両方の技術をサポートする共通の実施インフラストラクチャを可能にする。しかし、多くの実装では、代わりに２段階のモデリングシステムを使用するという方法が、これら９つの別々の合計変数を保持するよりも望ましい場合がある。このようなシステムで、最初の段階では、速度振幅と波の方向のみを計算する単純化されたモデルを実行し、コントロールポイントの合計[この文献は図面を表示できません]
ごとに３つの速度方向を求める。位相計算から分離されているので、これは単純化されたモデルであるこれらの総和は、それぞれ最終の最強振幅の成分を与えるので、第２段階で総和を行う方向を計算するために使用することができる。
【0146】
要約すると以下のとおりである：[この文献は図面を表示できません]
ここで２段階では、モデルは以下のように計算できる：[この文献は図面を表示できません]
プロセスの他の段階を始める前に、トランスデューサーモデル内である。残りの段階は以前とほぼ同じように進行するが、速度が作用する方向が異なるだけである。
【0147】
このアプローチは、圧力ではなく、手への運動量伝達にさらに直接的に対処するための空中触覚フィードバックの応用を可能にする。これにより、感覚の再現性、ひいては知覚行動の測定が向上し、これは、実用化されたテクノロジとしての空中ハプティックスの拡大に不可欠である。
【0148】
Ｖ．添付の図面
【0149】
図１７は、同じ音圧レベル（ＳＰＬ）を有する同じ波が、入射角に応じて表面に異なる効果をもたらすことを示す。左側１７１０では、表面（右）と平行に進行する波（左）がいくらかのエネルギーＥで衝突し、運動量交換の結果としてわずかに少ないエネルギーで反射し離れる。右側１７２０では、同じ表面に対して垂直に進む波はエネルギーを反射および交換せず、したがって同じ量のエネルギーを残す。その結果、触覚効果を与えることは不可能である。
【0150】
図１８では、上の２行１８１０、１８２０に示されているのは、飛行時間法によって作られた２つの異なる焦点の、それぞれｄＢ ＳＰＬ、ｄＢ ＳＩＬおよびｄＢ ＳＶＬでの位相アレイによって発生された音場の測定のシミュレーションであり、ライトグレーは＋１０ｄＢ、白は０ｄＢ、ダークグレーは−１０ｄｂである。0同じカラースキームを有するそれぞれ下列１８３０（ｄＢ ＳＰＬ−ｄＢ ＳＩＬ）、（ｄＢ ＳＰＬ -ｄＢ ＳＶＬ）および（ｄＢ ＳＩＬ−ｄＢ ＳＶＬ）は、各測定方法間の差異を明らかにする。ＳＶＬ場測定シミュレーションは、媒体の粒子の速度に基づいているので、追加のベクトル成分の相殺により、波をさらに相殺することができる。その結果ｄＢ ＳＶＬは常にｄＢ ＳＩＬより低く、ｄＢ ＳＩＬは、ｄＢ ＳＰＬより低いが、それらは、それらのそれぞれの成分の乗算およびＳＩＬ測定値の有効幾何平均のために、対数ｄＢ空間において同様の距離に離れている。
【0151】
４．空中ハプティックシステムにおける二次問題のための独自システムソリューション
【0152】
Ｉ．二次最適化のためのレイリーリッツ定理
【0153】
マトリクス[この文献は図面を表示できません]
が正方の正の半定値エルミート[この文献は図面を表示できません]
マトリクスにスーパーポーズする。マトリクスは、正方のエルミートマトリクスであるため、各固有値は実数でなければならない。これは正の半定値であるので、これらの固有値は[この文献は図面を表示できません]
として正でなければならない。
【0154】
次に、固有値を以下のように順序付ける：[この文献は図面を表示できません]
【0155】
マトリクスの[この文献は図面を表示できません]
は、最適化問題の解であるものとして示される：[この文献は図面を表示できません]
および[この文献は図面を表示できません]
【0156】
これは、レイリーリッツの定理である。したがって、これらのそれぞれからの[この文献は図面を表示できません]
ベクトルは、検出された固有値に対応する固有ベクトルとして見出すことができる。
【0157】
ＩＩ．位相アレイシステムのための二次最適化
【0158】
位相アレイシステムを考慮すると、以下の[この文献は図面を表示できません]
として定義することは有用である：[この文献は図面を表示できません]
[この文献は図面を表示できません]
個の単調音源を用いて、このベクトル解を用いて音源トランスデューサのアレイを、駆動することができる。[この文献は図面を表示できません]
が単位ノルム[この文献は図面を表示できません]
の固有ベクトルになるように、マトリクス[この文献は図面を表示できません]
を[この文献は図面を表示できません]
のエルミートマトリクスとして定義する:[この文献は図面を表示できません]
【0159】
これは、最適化の形成を示す。やらなければいけないことは、マトリクス[この文献は図面を表示できません]
に必要な特定の純粋な二次（マトリクス[この文献は図面を表示できません]
は対称またはエルミートマトリクスである必要があるため）目的関数を入力することである。
【0160】
エルミートマトリクス[この文献は図面を表示できません]
を生成する目的関数には、潜在的に有用な純粋に二次形式の５つの形式があるため、既存の固有ベクトルソルバーを使用して最大化し、空中超音波ハプティックを駆動するために使用できる解を得る：
【0161】
１．最大圧力最適化（全ＳＰＬ）：ａ．圧力の２乗に比例する、[この文献は図面を表示できません]
を最大化する；ｂ.これは、位相アレイシステムを駆動するために一般的に使用される目的関数である。
【0162】
２．最大粒子高速最適化（全ＳＶＬ）は：ａ．粒子高速の２乗に比例する[この文献は図面を表示できません]
を最大化する；ｂ.これは、一貫した触覚知覚を得るために必要であるため、波の方向性を考慮に入れる。
【0163】
３．しかし最大粒子高速最適化は、与えられた方向に沿う（使用可能なＳＶＬ）：ａ. 与えられた方向に沿う粒子高速の２乗に比例する[この文献は図面を表示できません]
を最大化する；ｂ．これは、波の方向性と（皮膚）表面[この文献は図面を表示できません]
に対する法線の影響の両方を考慮に入れる。
【0164】
４．最大音響強度（エネルギー束密度）最適化（全ＳＩＬ）：４．ａ.音響強度に比例する[この文献は図面を表示できません]
を最大化する：ｂ．これは、一貫した触覚知覚を得るために必要であるため、波の方向性を考慮に入れる。
【0165】
５．与えられた方向に沿った最大音響強度（エネルギー束密度）の最適化（全ＳＩＬ）：ａ．所与の方向に沿った音響強度に比例する[この文献は図面を表示できません]
を最大化する；ｂ．これは波の方向性と（皮膚）表面[この文献は図面を表示できません]
に対する法線の影響の両方を考慮に入れる。
【0166】
ＩＩＩ．圧力を最大化するための最適化
【0167】
各コントロールポイント点で圧力を最大化することができる二次形式最適化装置は、量をモデル化することになる:[この文献は図面を表示できません]
各コントロールポイントjでは、目的関数の総和を求める：[この文献は図面を表示できません]
これは、支配的な固有ベクトルが解かれる[この文献は図面を表示できません]
マトリクス要素を以下のようにする：[この文献は図面を表示できません]
ここで[この文献は図面を表示できません]
は、波の圧力の正および実振幅であり、初期放出からのコントロールポイント[この文献は図面を表示できません]
で、トランスデューサ[この文献は図面を表示できません]
から進行する波上の複素位相オフセット[この文献は図面を表示できません]
である。重み付け[この文献は図面を表示できません]
は、コントロールポイント間の相対振幅コントロールを提供するために、二次[この文献は図面を表示できません]
項を再重み付けするために使用される。グローバル振幅スケーリングは、入力パワー[この文献は図面を表示できません]
の単位ノルム条件によって提供される。これを二次重み付け係数およびコントロールループと組み合わせて、個々の圧力振幅に正確に対処することができる。
【0168】
ＩＶ．粒子速度を最大化するための最適化
【0169】
粒子速度を最大化することができる二次形式最適化装置は、量をモデル化することになる:[この文献は図面を表示できません]
各コントロールポイントｊについて行われる。これは、下記の目的関数の総和を求めることとなる：[この文献は図面を表示できません]
これは、支配的な固有ベクトルが解かれる [この文献は図面を表示できません]
マトリクス各要素を以下のようにする：[この文献は図面を表示できません]
ここで[この文献は図面を表示できません]
は、初期放出からコントロールポイント[この文献は図面を表示できません]
で、トランスデューサ[この文献は図面を表示できません]
から進行する波上での誘導粒子速度振幅のｘ方向部分（これは実数であり正数である）および[この文献は図面を表示できません]
複素位相オフセットである。重み付け[この文献は図面を表示できません]
は、コントロールポイント間の相対的な粒子速度振幅制御を提供するために、二次[この文献は図面を表示できません]
項に重みを付けるために使用される。グローバル振幅スケーリングは、入力パワー[この文献は図面を表示できません]
の単位ノルム条件によって提供される。これを二次重み付け係数およびコントロールループと組み合わせると、個々の速度振幅を正確に使用することができる。
【0170】
Ｖ．所与の方向のセットで粒子速度を最大化するための最適化
【0171】
所与の方向のセットにおいて、粒子速度を最大化することができる二次形式最適化装置は、量をモデル化することになる:[この文献は図面を表示できません]
各コントロールポイントｊについて行われる。これは、下記の目的関数の総和を求めることになる：[この文献は図面を表示できません]
これは、支配的な固有ベクトルが解かれるｎ×ｎマトリクス各要素を次のようにする：[この文献は図面を表示できません]
ここで、[この文献は図面を表示できません]
は、誘導粒子速度振幅の[この文献は図面を表示できません]
方向部分（これは実数であり正数である）であり、[この文献は図面を表示できません]
は、初期放出からコントロールポイント[この文献は図面を表示できません]
でトランスデューサ[この文献は図面を表示できません]
から進行する波上での複素位相オフセットである。重み付け[この文献は図面を表示できません]
は、コントロールポイント間の相対振幅コントロールを提供するために、二次項[この文献は図面を表示できません]
を再重み付けするために使用される。グローバル振幅スケーリングは、入力パワー[この文献は図面を表示できません]
に対する単位ノルム条件によって提供される。これを二次重み付け係数およびコントロールループと組み合わせて、個々の振幅を正確にアドレス指定するために使用できる。
【0172】
ＶＩ．活性音響強化
【0173】
活性音響強度[この文献は図面を表示できません]
は、音響媒体の圧力[この文献は図面を表示できません]
および粒子速度ベクトル[この文献は図面を表示できません]
を用いて定義することができる。音響波は、下記のように記述できる：[この文献は図面を表示できません]
ここで、活性強度は、純粋に実数のベクトル量である。波源の近位場では、活性場は、この代わりに[この文献は図面を表示できません]
が支配的である。しかしながら、1つまたは２つの波長を超えた遠位場では、[この文献は図面を表示できません]
となる。この理由により、[この文献は図面を表示できません]
のみが保証現実であると考慮し仮定される。ここで、単調な平面波源を考慮する：[この文献は図面を表示できません]
我々は、下記式を得る:[この文献は図面を表示できません]
ここで、[この文献は図面を表示できません]
は圧振幅であり、[この文献は図面を表示できません]
は各方向で粒子速度振幅である。
【0174】
所与の軸に沿って波によって伝達する音の強度、したがってエネルギー流束もまた重要であるので、スカラー寄与を測定するための内積が生成される:[この文献は図面を表示できません]
【0175】
音響波方程式の形成により：[この文献は図面を表示できません]
重ね合わせの原理は、圧力ベクトルと粒子速度ベクトルの両方に適用される関心のある点で、[この文献は図面を表示できません]
音響源を想定し各平面波として概算すると、方程式が以下のように得られる：[この文献は図面を表示できません]
【0176】
問題の領域では、[この文献は図面を表示できません]
であり、[この文献は図面を表示できません]
が同相であることが周知であり、その場合、共役の配置は問題ではなく、粒子速度の代わりに、共役圧力への共役の順序を逆にすることが可能である。（これは上記の方程式を拡張し、可換性を使用して、各項の複素べき乗の順序を逆にすることで確認できる）推論から、これはエルミートマトリクス[この文献は図面を表示できません]
を生成する方法で記述することができる；[この文献は図面を表示できません]
【0177】
ＶＩＩ．音響強化を最大化するための最適化
【0178】
スカラー値を生成する最初の数量は以下である：[この文献は図面を表示できません]
しかしながら、これは遠場であると想定され、これは以下のようにも記載できる：[この文献は図面を表示できません]
これは圧力和の４乗だが、複素べき乗の構造上、これらのうち２つは同じで、２乗根を使って削除することができる。
【0179】
したがって、音響強度を最大化することができる二次形式オプティマイザは、量をモデル化することになる：[この文献は図面を表示できません]
各コントロールポイント[この文献は図面を表示できません]
についてである。これは、下記の目的関数の相和を得ることになるであろう：[この文献は図面を表示できません]
これは、支配的な固有ベクトルが解かれる[この文献は図面を表示できません]
マトリクス各要素を次のようにする：[この文献は図面を表示できません]
ここで、[この文献は図面を表示できません]
は波前の普通のベクトルの部分方向[この文献は図面を表示できません]
であり、[この文献は図面を表示できません]
は、初期放射からコントロールポイントｊにおいてトランスデューサｑから進行する波上の波の圧力の正および実振幅、および[この文献は図面を表示できません]
は、複素位相オフセットである。重み付け[この文献は図面を表示できません]
は、コントロールポイント間の相対的な粒子速度振幅コントロールを提供するために、圧力の４乗を含む関数の平方根[この文献は図面を表示できません]
を実際に維持する二次項を再重み付けするために使用される。グローバル振幅スケーリングは、入力パワー[この文献は図面を表示できません]
の単位ノルム条件によって提供される。これを二次重み付け係数およびコントロールループと組み合わせると、個々の速度振幅を正確に扱うことができる。
【0180】
ＶＩＩＩ． 所与の方向のセットでの音響強化の最大化のための最適化
【0181】
所与の方向に沿って音響強度を最大化することができる、二次形式最適化装置は、量をモデル化することになる：[この文献は図面を表示できません]
しかし、これを遠位場で想定すると、これは以下で記載できる：[この文献は図面を表示できません]
これは、各コントロールポイント[この文献は図面を表示できません]
についてである。これにより、目的関数の合計は次のようになる：[この文献は図面を表示できません]
これの支配的な[この文献は図面を表示できません]
マトリクスの各要素を生成する固有ベクトルが解かれるべきである：[この文献は図面を表示できません]
ここで、[この文献は図面を表示できません]
は誘導された粒子速度振幅の[この文献は図面を表示できません]
方向の部分であり（これは実数で正数である）、[この文献は図面を表示できません]
は波圧の正、かつ実の振幅であり、[この文献は図面を表示できません]
は、初期放射からコントロールポイント[この文献は図面を表示できません]
でトランスデューサ[この文献は図面を表示できません]
から進行する波上での複素数位相オフセットである。重み付け[この文献は図面を表示できません]
は、コントロールポイント間の相対的な粒子速度振幅コントロールを提供するために、二次の[この文献は図面を表示できません]
項を再重み付けするために使用される。グローバル振幅スケーリングは、入力パワー[この文献は図面を表示できません]
に対する単位ノルム条件によって提供される。これを二次重み付け係数とコントロールループと組み合わせて使用することで、個々の速度振幅を正確に扱うことができる。
【0182】
ＩＸ．重み付けファクターの発見
【0183】
レイリーリッツの定理は、支配的な固有ベクトルを見つけ、トランスデューサの物理的限界までパワーを増加させることによって、上記のマトリクスのいずれも最大化できることを示す。これは、パワー相互を通して到達することができる。コントロールポイントの合計のそれぞれを別々に定義しながら、マトリクス[この文献は図面を表示できません]
を構成することによって、ベクトル[この文献は図面を表示できません]
を乗じる場合、目的関数の評価を追加の副次的な作用として実行することができる。これは、重み付けされたときにシステムの解に適切である振幅を最大化が見つけることを可能にするように、重み付け[この文献は図面を表示できません]
に対する適切な値を見つけるために使用することができる。
【0184】
パワー反復処理は、下記のように記述される：[この文献は図面を表示できません]
ここで：[この文献は図面を表示できません]
【0185】
我々は、各コントロールポイントで目的関数を決定するために、パワー反復処理を使用することができる。
【0186】
定義づけを行う：[この文献は図面を表示できません]
かつ：[この文献は図面を表示できません]
各反復処理で、計算するのではなく以下とすることができる：[この文献は図面を表示できません]
【0187】
つまり、各コントロールポイントでの目的関数は、以下のように評価することができる：[この文献は図面を表示できません]
【0188】
最終的に、以前のように[この文献は図面を表示できません]
の計算をすることができる：[この文献は図面を表示できません]
ここで、最終的な固有ベクトルの解は、上記の通りである。しかし、このように計算すると、[この文献は図面を表示できません]
が、コントロールポイント[この文献は図面を表示できません]
の中間結果を出すことができる。これが所望のコントロールポイント出力を表すレベルを固定するために、一組[この文献は図面を表示できません]
の値が、ユーザによって与えられる目的値と直接比較されることを可能にするために、平均目的関数の値が取られる。このケースでは、重み付け係数は、例えば、Ｎｅｌｄｅｒ−Ｍｅａｄアルゴリズムなどの高価な関数の最適化を目的とした導関数のない方法を使用して、次の反復処理で更新することができる。
【0189】
べき乗法と結合された導関数のない最適化を何度か繰り返した後、中間結果[この文献は図面を表示できません]
を使用してトランスデューサを駆動するべきである。したがって、これにより、システムは、ある程度の誤差を許容しながら、時間の経過と共に、目的関数の所望の値を追跡することが可能になる。このシステムの利点は、それがトランスデューサアレイを駆動するのに、非常に簡単な方法であるということである。最大の不利な点は、多くの計算が必要とされることだが、計算は実質的に複雑な乗算にすぎず、十分なハードウェアで簡単に実装される。
【0190】
Ｘ．添付の図面
【0191】
図１９に示されるものは：左上１９１０：非直交ベクトルを使用して二次粒子速度問題を解くことによって作成されたコントロールポイントであり、[この文献は図面を表示できません]
平面を通るカットとして表示される。右上１９２０：[この文献は図面を表示できません]
平面を通してカットとして表示した、同コントロールポイント。左下１９３０：比較すると、二次圧力問題を解くことによって作成されたコントロールポイントであり、[この文献は図面を表示できません]
平面を通るカットとして表示される。右下１９４０：[この文献は図面を表示できません]
平面を横切る圧力問題の解を通る断面。
【0192】
５．二次問題を通して指定された位相アレイのためのコントロール領域
【0193】
本開示は、ユーザ供給設計パラメータを使用して、基底関数を自動的に生成するための方法を説明する。これらは、オンラインで設定された所与のトランスデューサ基底の空間における、二次関数にわたる効率的な最適化に入力される。これにより、ユーザは、様々な異なる効果を生み出す可能性があるオンライン最適化を通じて、アクセス可能なデザインスペースを手に入れることができる。二次問題のためのコントロール領域であるので、「二次コントロール領域」と呼ばれる線形演算子を使用してトランスデューサ効果にわたって二項積を評価することによって、マトリクスを構築することができる。このマトリクスの主固有ベクトルは、このマトリクスの固有ベクトルである、正規直交解の集合から最良の解として選択されてもよい。最適化は、出力としてトランスデューサ活性化係数−個々のトランスデューサ素子を駆動して音場または所望の特性を有する現象を再現するのに必要な、生の複素数値を生成する。
【0194】
次いで、これらの生の複素数値を使用して、後に使用する線形システムの基底関数をさらに含むことができ、その解によって、ユーザは、それぞれに利用可能な出力パワーを再重み付けしながら、すべての所望の効果を同時に生成することができる。
【0195】
Ｉ．線形コントロール領域
【0196】
線形コントロール領域は、複素数値重み付け関数[この文献は図面を表示できません]
を用いて、圧力などの線形音響量との積分を実行する演算子として定義される：[この文献は図面を表示できません]
これは、音響場を評価するために使用される。[この文献は図面を表示できません]
が音響量における類似の位相分布と対にされる場合、複素共役が、関数評価を最大にするために使用されることに留意されるべきである。この積分に対する任意の近似値（サンプリングなど）を使用することができ、あるいは実際にはこれを使用して、音圧場にわたって任意の線形演算子を近似することができる。それがまた線形であるならば、ベクトル粒子速度を使用してコントロール領域を最適化するようにさらなる積分を構築することができる：[この文献は図面を表示できません]
ここで、[この文献は図面を表示できません]
は、媒体の粒子速度、[この文献は図面を表示できません]
は、各次元が複素数である三次元ベクトル場であり、これは、上記と同じ目的で重み付け関数として機能する。この方法では、音響量のさまざまなパターンが線形的に評価され、位相が[この文献は図面を表示できません]
のスカラー関数またはベクトル関数によって作成された「テンプレート」と一致するときに、結果の特定の位相が最大になる。
【0197】
また、[この文献は図面を表示できません]
は、複数の外見上ばらばらな領域の単一の評価を作成するために、ほぼゼロ（またはゼロベクトルで構成される）である場合があることにも留意すべきである。この定式化は、線形音響量のフィールドにわたって、任意の線形演算子に拡張することができるため、任意の整数近似（サンプリングなど）を使用することができる。
【0198】
位相アレイは、各トランスデューサをコントロールする一連の複素数値の係数によって、記述することができる。マトリクス[この文献は図面を表示できません]
に符号化された音響量を最大にするために、そのような入力[この文献は図面を表示できません]
ベクトルを用いて、固定周波数で音響変換器のアレイを駆動することは、線形代数問題として表すことができ、その意図は以下の通りである：[この文献は図面を表示できません]
【0199】
問題の記述は固有ベクトル、ここでは[この文献は図面を表示できません]
の定義でもあるため、これはマトリクス[この文献は図面を表示できません]
の主要な固有ベクトルをとることによって解くことができる。マトリクス[この文献は図面を表示できません]
を構築することは、まず、問題の音響量の大域空間における基底関数として各トランスデューサによって放出された線形音響量を考慮することによって達成される。これは、[この文献は図面を表示できません]
の絶対位置で表される、１点での圧力やその他の音響量の最大化を達成するための方法である。しかしながら、記述したコントロール領域積分を使用することにより、各個々の変換素子の場に放出される音響量を支配する、[この文献は図面を表示できません]
は、[この文献は図面を表示できません]
と交換され、列ベクトルを生じる:[この文献は図面を表示できません]
【0200】
線形性を呼び出すことによって、各トランスデューサを駆動する複素活性化値を乗算して、ベクトルを導き出すことができる：[この文献は図面を表示できません]
ここで、[この文献は図面を表示できません]
は、トランデューサの数である。この[この文献は図面を表示できません]
が線形場に作用する線形演算子の結果とすると、それはエルミート転置され、それ自身によって乗算されて二次形式を生成することができる、最終的には以下のようになる。[この文献は図面を表示できません]
【0201】
重み付け関数[この文献は図面を表示できません]
を修正し、マトリクス[この文献は図面を表示できません]
の固有ベクトルを解くことによって、様々な異なるビーム、表面および他の音響現象が生成され最大化され、音響内に最適に作動される領域を生成する。
【0202】
先に論じたヌル空間生成技術（「空中ハプティックシステムにおける二次問題のための固有システム解」）もまた、音響浮揚を記述するこの方法を拡張するために使用され、ならびにヌルポイントを記載する構造のポイントのフィールドへ導入することによって、ヌルステアリングおよび空中テクスチャ形成を可能にするために使用される。このプロセスはマトリクスの階数を増加させるため、これを使用したとしても単純化は不可能であることに注意すべきである。
【0203】
等価的に、多数のコントロール領域について、[この文献は図面を表示できません]
個のマトリクスを一緒に加えることによって複数の最適化を実行することができる。これは、コントロールポイントを表すマトリクスを含むこともできる。このことは、ランクが上がったマトリクスを生成し、いくつかの簡略化された方法は使用される本ケースでは不可能である。
【0204】
結果を考慮し、反復的に各点を再重み付けすることによって、音響量の必要な幾何学的構成に収束するように、システムを修正することができる。これは、振幅が出力で構成されない場合に、一般的に必要である。
【0205】
これは、生成した音響量を合計し、次の反復において結果にさらに近づくようにシステムを再重み付けすることによって達成する。例えば、サブスプリットボリュームの積分点または単一の圧力点に対する反復再重み付けは、次のようになる：[この文献は図面を表示できません]
ここで、[この文献は図面を表示できません]
は、本技術における次の粗な反復、[この文献は図面を表示できません]
は、現在の反復、[この文献は図面を表示できません]
は、所望の圧力（これは複素数でもよい）、そして分母は現在の反復によって得られた全圧力を指す。得られた領域の振幅にのみ影響を与え、位相には影響を与えないために、分数は実数のみとしてレンダリングされてもよい。
【0206】
ベッセルビームは、ビームの進行方向に重み付け関数[この文献は図面を表示できません]
を表現し、ビームが前方に移動するにつれて、前方に位相角を修正しながら、単位圧力を生成することによって生成される。後続の各ポイントが通過する角度がビーム軸に沿って位相をサンプリングするのに十分小さい限り、これはポイントサンプルでサンプリングされる。
【0207】
位相のサンプリングは、それが移動する角度において、[この文献は図面を表示できません]
ごとに２つより多い位相サンプルを有することが重要である。これは、経路に沿った増分が[この文献は図面を表示できません]
の場合、波の進行方向が決定できないためである。これはまた、多くの場合、位相がカバーする単色周波数距離の各波長に対して、[この文献は図面を表示できません]
ラジアンの回転を完了するように、位相が前方に増加しなければならないことを意味する。これは、コントロール領域の定義で説明される形式のパス積分への近似である。
【0208】
このように書かれた積分を最大化する結果として、トランスデューサは、アレイから遠ざかる直線で累進的に定義された単純なビームの場合、その理論的能力により最適なベッセル状ビームを生成し、この構成で最適な圧力を生成するだろう。図２０に示されるのは、平面直線トランスデューサアレイに平行な近距離場領域の端部でのフィールドの断面であり、複素数値のトランスデューサ活性化係数の構造を明らかにしている。単純なプログレッシブ位相ビームではベッセル関数近似が得られ（左２０００）、アンダーサンプリングされたヘリカルビームでは双極子の振る舞いを持つベッセル関数が得られ（右２０１０）、双極子は２つの構造を作り出し、それらは互いに位相で反対になるように活性化する。トランスデューサの解像度の限界、有限アレイサイズ、および最適化が回避する必要があるその他の物理的考慮の影響により、ビームは正確なベッセルビームではない可能性があるが、実際のトランスデューサアレイの物理的限界を考慮した最適化により、ベッセルビーム関数の単純な評価よりも効果的に機能する可能性がある。大きなアレイでは、これはベッセル関数に密接に従っており、図２１は、図２０のアレイの中心から半径方向に考慮したときに圧力と共にプロットされたベッセル関数のグラフ２１００を示している。
【0209】
このシステムのさらなる特性は、任意の生成された単色場関数を有する任意の位置、および向きのトランスデューサが使用されることである。これは、定義により平面アレイの場合にのみ機能する、サンプリング方法に対するベッセルビーム生成における重要な利点である。
【0210】
複数の位相外経路を作成することによって、位相シフトおよび位相特異点は、経路の周りにらせん波を巻くことによって暗黙的に作成される。固有経路アプローチに組み込まれている単純な二次最適化を使用しながら、この経路はまた湾曲しており、高度の構成可能性をもたらす。これは、ビーム幅より近いが互いに位相が異なる多くの平行ビームラインにおける別々の位相巻線の定義を通して、ビームの回転を誘導することによって達成する。次に、図２２に示すように、近くの位相をビームの中心線で互いに相殺するように設定することで、高次のベッセルビームを形成できる。図２２では、円の周りの漸進的位相の２つ以上のサンプル（したがって、らせん配置の位相）は、左側２２１０に見られるように角運動量を有するベッセル関数の形成を生じさせる。らせんの周りの位相を特定するために、[この文献は図面を表示できません]
ラジアンを超える値を使用することによって、右側（２２２０）の[この文献は図面を表示できません]
ビームのように高次のベッセルビームを生成することができる。
【0211】
これらの位相線は、互いに素な領域が単一の領域最適化としてどのように定義され得るかの例であるので、点または他の任意の積分近似を用いてサンプリングされ、最適化するボリューム内でのみゼロ以外の値に機能する。高次マシュー型楕円ヘリカルビームは、位相が漸増し楕円の周りに間隔を置いて増加する、位相ビームを導入することによって生成することができる。図２３に示すように、マシュー型ビームなどの高次ベッセルビームと同様の特性を有する楕円形ヘリカルビームは、位相の楕円の円周をパラメータ化することによって作成することができる。[この文献は図面を表示できません]
ラジアンの割合を決定し、楕円弦の円周に比例する位相角を割り当てることは、ビームの最適化のための適切な構成を生成する。左の画像２３１０は、マシュー型ビームの近くで、楕円形に再形成された高次ベッセルを示す。右側の画像２３２０は、中央領域に形成され始めている架橋された、らせん構造を有するおおよそのマシュー型ビームを示す。どちらも[この文献は図面を表示できません]
ラジアンの位相円に基づいており、らせん状の波面を形成する。
【0212】
このようにして生成されたベッセルファミリーのビームは、文献においてより制限的な形態で存在するよりエキゾチックなビームを生成するように修正されてもよい。波の横断に適した方法で前進する物理的にもっともらしい位相関数によって横断される場に任意の経路を作成することによって、さまざまな結果が再現されるかもしれないが、より一般的な形で固有システムを伴う最適化のため、ビームは積分によってのみ表現される必要があるので、初期の帯域幅コストは以前よりも低い。
【0213】
湾曲した経路をとり、位相を前方に動かすと、４５°の円セグメントによって生成された障害物の周りにビームを操縦するために、一般的に使用される湾曲ビームが生成される。図２４に示されているのは、ビームに沿った音圧を等化するために粗い反復法を使用して、ラインに沿った連続に接続された円の４５度セグメントをたどる位相関数から形成された曲げビームの画像２４００である。トランスデューサアレイは、画像の下端に沿って配置される。
【0214】
複数の位相経路を含むさらなる手法は、「ボトルビーム」、すなわち中心から切り取られた部分を有するビームを作り出すことができ、それは大きな物体の浮揚を容易にすることができる。これは、図２５において生成されており、図２５は、それに沿って位相がたどられるようにされてきた、２つのガウス曲線を双晶することの画像２５００を示す。位相関数を作成することによって生成されたボトルビームは、ビームに沿って音圧を均一化するために、粗い反復法を使用してガウス曲線に従う。トランスデューサアレイは、画像の下端に沿って配置される。下部領域の点線の性質は、最適化に入力される位相速度が遅すぎること、および初期区間により速い位相速度を使用することによって、ビームが改善される可能性があることを示唆している。物体が十分に硬い場合には、空間の一部のみを高音圧領域にする必要があるので、図２６のように曲線の一部を使用することもできる。図２６は、不均一な方法で高音圧の領域を最適化することによってさらに大きな物体を浮揚させるために使用することができる、累進相で画定された円の一部から形成された不完全ボトルビームの画像２６００を示す。浮揚すべき物体または物体が低音圧領域を通って流れない、押し出される、または漏れることがないように十分な凝集力を有することを仮定する場合、大きな物体に及ぼされる力を均一に分散させる必要はない。
【0215】
ベッセル状ビームは、経路を位相速度より速く前方に移動させることによって生成することができる。この方法では、図２７に示すように、ビームは通常よりも速い位相速度で最適化され、波長が長くなり、位相速度が音速よりも速く移動するように見えるビームを作成する。図２７に示されるのは、音速よりも速く進む波面の画像２７１０であり、これは、逆説的に同じ周波数でより長い波長を生成するように、角を通って移動する、位相経路を実施することによって生成され得る。暗黙の最適化は、メインビームが交差するビームで構成されていると推論し、位相速度の指定された増加を引き起こす。右側の画像２７２０より見ることができ、波長は、ビームの最初の部分の同じ周波数で約４０％長くなるように操作される。蛇行ビームは、正弦曲線に沿って漸進的位相をサンプリングすることによって構築することができる。図２８に示されているのは、例２８１０であり、そこでは正弦曲線の短い部分が定義され、続いて位相配列上の正弦曲線に続く、漸進的な圧力相から構成される２８２０のより長い２次元蛇行ビームが続く。これらは、漸進的な位相のさらなる部分を単に追加することによる蛇行正弦波ビームの拡大が必要とされる限り（およびトランスデューサの物理的配置によって可能にされる）自然にビームを拡大させることを示す。らせん状に蛇行するビームを生成するために３次元全てに蛇行するビームを作ることによって、アルキメデススクリューシステムと同様の方法で、潜在的には波の放射、波長よりもはるかに大きいフィールド内の固体物体に対して、垂直または反対方向に力を加えるために使用できるビームを生成することができる。これは、それによってビームが物体を浮揚させるために使用される、異なるメカニズムを与える。図２９は、[この文献は図面を表示できません]
平面内の断面を詳細に示す左画像２９１０と、[この文献は図面を表示できません]
平面内のビームの断面を示す右画像２９２０とを示す。[この文献は図面を表示できません]
に曲がっているが、[この文献は図面を表示できません]
にも進行する円を描く経路をたどるこの技術を使用して、真のらせん形状を作り出すことができる。これは、トラクタービームのような仕方で、波長よりはるかに大きい物体に力を加えて浮揚させるために使用することができる。
【0216】
固有システムマトリクスが単一のコントロール領域で定義されている場合、そのマトリクスはランク１であり、複素数値ベクトル[この文献は図面を表示できません]
とそれ自身のみの２項積で形成し、マトリクスを[この文献は図面を表示できません]
と定義する。この特定のケースでは、システムの単一の固有ベクトルはすぐに分かり[この文献は図面を表示できません]
であり、[この文献は図面を表示できません]
を意味する。このケースで目的関数を評価することによって、複素数値のトランスデューサ起動の係数のそれ自身を特定する。これは各トランスデューサのコントロール領域機能を単純に評価することによって定義し、活性化の係数である。
【0217】
従前のとおり、二次形式のために、各反復についての再重み付けを通して、有効性を評価することもまた追求される。これは、各ポイントにおける圧力または各小領域にわたる圧力の積分を評価し、所望の圧力で除算し、新しい関数値を生成して、[この文献は図面を表示できません]
小領域[この文献は図面を表示できません]
の重みを更新することによって達成する。それから[この文献は図面を表示できません]
でもある単一の固有ベクトルはそれから再評価されて、前から固有システム上で繰り返される反復に対する同等のアプローチを生成するが、非常に減少した計算コストで再評価される。
【0218】
媒体の粒子速度でこれを達成するために、粒子速度の試験を使用することができ、これは最大に達するまで方向を変えることになる：[この文献は図面を表示できません]
圧力は、ハイブリッド方法を生成する旧技術を使用して監視することもできる。このハイブリッド技術では、重み関数を指定して以前と同じように圧力を更新するが、効果的なコントロール領域を使用する：[この文献は図面を表示できません]
スカラーだが可変の複素数値[この文献は図面を表示できません]
を得、不変の重み付け法線ベクトル[この文献は図面を表示できません]
と粒子速度[この文献は図面を表示できません]
と対とする。このアプローチは、求められる波の方向を変えることなく、満足のいく圧力を見つけるために最適化することができる。これは、次のセクションで説明されている特定の修正と共に図２８に示されている。
【0219】
ＩＩ．指向性コントロールポイントおよび領域に対する波面フィルタリング
【0220】
指向性出力の作成に必要なさらなる技術は、トランスデューサ出力波面[この文献は図面を表示できません]
と法線ベクトルまたは粒子速度ベクトル[この文献は図面を表示できません]
とが反対である場合に、トランスデューサデータを破棄し、それをゼロに設定することである：[この文献は図面を表示できません]
考えている領域のポイントについてである。
【0221】
図３０は、トランスデューサアレイからさらに定義された２つの速度ベクトルサンプルの左画像３０１０と、解除され定義されたコントロール方向と反対の方向に波面があるトランスデューサを有した、アレイに近接して定義された２つの速度ベクトルサンプルの右画像３０２０とを示す。これら２つの互いに素な領域で指定された速度重み付け関数を使用することは、それぞれの場合において単一の領域を構成し、その領域は、２つの別々のポイントにおける、媒体の粒子速度を表す方向ベクトルに対して最適化される。
【0222】
[この文献は図面を表示できません]
方向の最大量の速度を見つけようとすることは、結果として生じる粒子速度方向が完全に[この文献は図面を表示できません]
方向にあることを、必ずしも意味しないことに留意されたい。図３０は、所望の方向とは反対の波面を提供する、相反するトランスデューサをフィルタリングした結果を示し、図３１は、トランスデューサの波面方向が、波面基準を使用するフィルタリングされてなかった結果である、左画像３１１０と右画像３１２０を示す。相反するトランスデューサは、図３１の場合、自由に干渉することができ、その病状は、定在波構造の形で水平に明らかである。これは、波面方向が粒子速度から欠落している方向情報を伝達するため、粒子速度は波がどの方向から生じたかに関する情報をもたらさず、望ましくない干渉を可能にするからである。
【0223】
これはまた、コントロール領域だけではなく、通常のコントロールポイントの単純な場合にも重要である。コントロールポイントの場合は、事実上単純化されたコントロール領域であり、粒子速度の単一点サンプルと重み係数なしで近似される。どちらの場合でも、「役に立たない」トランスデューサを除外しないことは、「定在波」現象の影響の要素に関して同様の結果をもたらします。これは、複素成分を含む粒子速度ベクトルでは、異なる方向に進む波の違いを特定できないためである：このモデルでは、位相０で左から右へ進む波は、位相[この文献は図面を表示できません]
で右から左へ進む波と数学的に同じである。結果として、これがコントロールポイントまたは領域ベースの生成で正しく処理されない場合、最適化は両方を使用しようとし、意味のあるように寄与しない定在波成分を含む解を生成し、不要な、余分なフィールドのないヌルの領域を作成する。この最適化ステップを用いずに基底関数として線形最適化、固有システム、さらには単純な式に入力した、コントロールポイントがこれらのアーチファクトを示すため、これは重要な点である。
【0224】
記載のこれらの技術は、本明細書で定義されるような単一の固有システムを反復することによって、または単純な場合には単一の点または領域を作成することによって、任意に構成された領域およびビームを作成することができる、音響位相アレイシステムを作成するために装置に統合される。単一のポイントまたは領域の場合には、トランスデューサ素子を駆動するのに必要な振幅および位相が、ポイントまたは領域を評価することによって表される、振幅および位相であることを示す簡単な式を生成するために最適化が適用される。大域的方法の反復および経時的な反復適用は、様々な触覚効果を生み出すために、時間依存的に変更および修正される。これは、それらが機能的に関連しているので、所望の場および導出されたトランスデューサ係数の両方に適用される。
【0225】
ＩＩＩ．可能な請求項１．単調な音響場についての方法であって、
前記音圧場関数と重みづけ関数との積を積分することによって音圧の分布を単一値に見積もることであって、両方の関数が複素数値であり、物理的なトランスデューサの特性から導出されることと、
前記マトリクスのヌル位置に追加されるゼロまたはヌルポイントを適用することと、
これを最適化において最大化するべき目的関数として使用することであって、前記最適化の結果は、各トランスデューサについて前記目的関数の交差乗算項から形成されるマトリクスの固有ベクトルの形式をとることと、
そのようなマトリクスの1つまたはそれ以上を合算して1つまたはそれ以上の固有ベクトルが前記特定の場を生成するような固有システムを生成することを含む方法。
２．単調な音響場についての方法であって、
実際の３次元重み付けベクトル場と粒子速度場と重み付け関数とのドット積を積分することにより単一のベクトル値へと音響粒子速度の分布を見積もることであって、両方の関数が複素数値であり、物理的なトランスデューサの特性から導出されることと、
前記マトリクスのヌル位置に追加されるゼロまたはヌルポイントを適用することと、
これを最適化において最大化するべき目的関数として使用することであって、前記最適化の結果は、各トランスデューサについて前記目的関数の交差乗算項から形成されるマトリクスの固有ベクトルの形式をとることと、
そのようなマトリクスの1つまたはそれ以上を合算して1つまたはそれ以上の固有ベクトルが前記特定の場を生成するような固有システムを生成することを含む方法。
３．単調な音響場についての方法であって、
実際の３次元重み付けベクトル場と粒子速度場と重み付け関数とのドット積を積分することにより単一のベクトル値へと音響粒子速度の分布を見積もることであって、両方の関数が複素数値であり、物理的なトランスデューサの特性から導出されることと、
それぞれ見積もられた積分の副領域で特定される垂直に対して反対方向を有するポイントでの波頭のトランスデューサの寄与をゼロに設定することと、
前記マトリクスのヌル位置に追加されるゼロまたはヌルポイントを適用することと、
これを最適化において最大化するべき目的関数として使用することであって、前記最適化の結果は、各トランスデューサについて前記目的関数の交差乗算項から形成されるマトリクスの固有ベクトルの形式をとることと、
そのようなマトリクスの1つまたはそれ以上を合算して1つまたはそれ以上の固有ベクトルが前記特定の場を生成するような固有システムを生成することを含む方法。
４．単調な音響場についての方法であって、
前記音圧場関数と重み付け関数との積を積分することによって音圧の分布を単一値に見積もることであって、両方の関数が複素数値であり、物理的なトランスデューサの特性から導出されることと、
前記マトリクスのヌル位置に追加されるゼロまたはヌルポイントを適用することと、
これを最適化において最大化するべき目的関数として使用することであって、前記最適化の結果は、各トランスデューサについて前記目的関数の交差乗算項から形成されるマトリクスの固有ベクトルとなることと、
前記トランスデューサの入力信号のための最適な位相および振幅を記述するために各トランスデューサの寄与を使用することを含む方法。
５．実際の３次元重み付けベクトル場と粒子速度場と重み付け関数とのドット積を積分することにより単一のベクトル値へと音響粒子速度の分布を見積もることであって、両方の関数が複素数値であり、物理的なトランスデューサの特性から導出されることと、
前記マトリクスのヌル位置に追加されるゼロまたはヌルポイントを適用することと、
これを最適化において最大化するべき目的関数として使用することであって、前記最適化の結果は、各トランスデューサについて前記目的関数の交差乗算項から形成されるマトリクスの固有ベクトルの形式をとることと、
前記トランスデューサの入力信号のための最適な相および振幅を記述するために各トランスデューサの寄与を使用することを含む方法。
６． 実際の３次元重み付けベクトル場と粒子速度場と重み付け関数とのドット積を積分することにより単一のベクトル値へと音響粒子速度の分布を見積もることであって、両方の関数が複素数値であり、物理的なトランスデューサの特性から導出されることと、
それぞれ見積もられた積分の副領域で特定される垂直に対して反対方向を有するポイントでの波頭のトランスデューサの寄与をゼロに設定することと、
前記マトリクスのヌル位置に追加されるゼロまたはヌルポイントを適用することと、
これを最適化において最大化するべき目的関数として使用することであって、前記最適化の結果は、各トランスデューサについて前記目的関数の交差乗算項から形成されるマトリクスの固有ベクトルの形式をとることと、
前記トランスデューサの入力信号のための最適な相および振幅を記述するために各トランスデューサの寄与を使用することを含む方法。
７．さらに、デバイスが、方法と、音響トランスデューサを駆動するための結果の使用を実装する、請求項１〜６に記載の方法。
８．さらにデバイスが、方法と音響トランスデューサを駆動するために使用する線形基底セットの一部としての結果の使用を実装する、請求項１〜６に記載の方法。
９．音響ビームを生成する、請求項１〜６に記載の方法。
１０．音響量を生成する、請求項１〜６に記載の方法。
１１．音響波面を生成する、請求項１〜６に記載の方法。
１２．音響場についての方法であって、
物理的なトンランスデューサの特性から得られる場の線形結合を使用して1つのポイントにおける単一の値として音響粒子速度を見積もることと、
前記コントロールポイントにおいて、それぞれ見積もられた積分の副領域で特定される垂直に対して反対方向を有するポイントでの波頭のトランスデューサの寄与をゼロに設定することと、
これを前記トランスデューサの入力信号のための相および振幅を見出すための最適化問題における基底関数として使用することと、を含む方法。
１３．音響場についての方法であって、
物理的なトンランスデューサの特性から得られる場の線形結合を使用してポイントにおける単一の値として音響粒子速度を見積もることと、
前記コントロールポイントにおいて、それぞれ見積もられた積分の副領域で特定される垂直に対して反対方向を有するポイントでの波頭のトランスデューサの寄与をゼロに設定することと、
これを最適化において最大化するべき目的関数として使用することであって、前記最適化の結果は、各トランスデューサについて前記目的関数の交差乗算項から形成されるマトリクスの固有ベクトルとなることであって、その後トンランスデューサの入力信号としての位相および振幅を見出すことと、を含む方法。
１４．音響場についての方法であって、
物理的なトンランスデューサの特性から得られる場の線形結合を使用してポイントにおける単一の値として音響粒子速度を見積もることと、
前記コントロールポイントにおいて、それぞれ見積もられた積分の副領域で特定される垂直に対して反対方向を有するポイントでの波頭のトランスデューサの寄与をゼロに設定することと、
トランスデューサの入力信号のため最適な相および振幅を記述するために各トランスデューサの寄与を使用することを含む方法。
１５． 前記トランスデューサは、得られた構造が入力信号として適用され、変調された場合にハプティック効果を生成する、請求項１〜１４に記載の方法。
１６． 前記トランスデューサは、前記得られた結果が入力信号として適用され、時間にわたり変調された場合にハプティック効果を生成する、請求項１〜１４に記載の方法。
【0226】
６．ハードウェアにおける複素数値の超越関数、乗算および分割の評価
【0227】
[この文献は図面を表示できません]
の代わりに[この文献は図面を表示できません]
の間で変換することを目的として、ＢＫＭアルゴリズムは物理的なハードウェア実装においてより有用である。
【0228】
Ｉ．イントロダクション
【0229】
ハードウェアにおける、三角関数および他の特別な関数の計算のための多くの技術が存在する。そのような方法の１つはＢＫＭ法であり、これは複素対数を複素べき乗に交換し、その逆の変換を行い、これらの特別な関数をある表現から別の表現に移動することの副作用として生成する。特に興味深いのは、そのような関数が必要であることがわかっている状況にしばしば遭遇するため、複素数値関数を扱うものである。
【0230】
対数形式とべき乗形式との間を移動する効率的な方法は、単純な加算および減算として、対数空間で乗算および除算を実行することを可能にする。プロセスの性質上、これはべき乗関数に関連する超越関数を容易に生成するためにも、利用される可能性がある。これは単一のアルゴリズムから生成することができるので、１つの要素だけが実装を必要とし、１つの方法だけが実装され最適化される必要があり、開発時間が短縮されるので実装サイズが小さくなる。
【0231】
既存のアルゴリズムであるＢＫＭはこのアプローチを実行するが、提案された方法では、ほとんどまたは完全に除去されている３つの不利な要件を含む：
【0232】
とりわけ、相でなければデータ入力、出力上のステップの前処理と後処理を実装することを困難とするハードウェアの分割。
【0233】
効率的な実施のために、丸めの問題を引き起こす非合理的な価値および要因。
【0234】
浮動小数点型と統合するために、必要とされる追加の演算。
【0235】
次のセクションは、提案した方法がこれらの欠点がないように、どのように設計されてきたかを説明する。
【0236】
ＩＩ．動機と比較範囲の縮小
【0237】
プロセスは、オイラーの式[この文献は図面を表示できません]
に集中しており、[この文献は図面を表示できません]
は、式 [この文献は図面を表示できません]
が複素一般化されているため、さらなる係数である。ＢＫＭアルゴリズムが具体化する複雑な機能はまさにこれであり、実行されるプロセスを[この文献は図面を表示できません]
とする。このアルゴリズムは、収束が保証され得る範囲に、複素数値[この文献は図面を表示できません]
を当てはめるために範囲縮小技術を頻繁に使用する。これは、[この文献は図面を表示できません]
を達成するために、提案された方法を使用して、無理数スケーリングファクターで前処理または後処理することによって評価することがおそらく容易である限りにおいて、適切な実装を有する関数ではない。多くの場合、この無理数スケーリングファクターは、入力と出力が物理的または計算上意味のある有理数の倍数になるように、意図するワークフローを提供するために削除されることがある。
【0238】
複素べき乗を処理するとき、べき乗は実数部と虚数部に分割され、それらは同じ入力データを使用するが、異なる結果を生成する。これは、ある意味ではデカルト座標から極座標への変換に似ている。
【0239】
べき乗化されるべき乗対数の実部だけを考えると、これはデカルト−極変換の[この文献は図面を表示できません]
成分と同様に、すべての係数情報を含む。これは伝統的に実数部の力に引き上げられた無理数[この文献は図面を表示できません]
である。それをアルゴリズムに供給するために、これは最初に範囲縮小を受けなければならず、それは剰余を見つけなければならず、それが収束領域内に収まるように、無理数[この文献は図面を表示できません]
を法で除算ことを意味する。これはハードウェア分割を含み（１９９４年の論文は、ＳＲＴ分割技術を推奨する）、それ故にこれはアルゴリズムと共に分割を実行するのに必要とされる論理の点で、非常に高価な技術を加える。
【0240】
提案されたアルゴリズムの前提は、範囲縮小を簡単にするために、２のべき乗を使用すべきであるということである。この動機で、実部は実部である所与の力に引き上げられた２、つまり底２のべき乗として解釈されるべきである。実数対数入力の整数は、出力上のビットシフトを参照することができ、相対２進対数入力の小数部分だけが収束領域に収まる必要がある複素対数の実部の整数部分は、このようにして四捨五入されて別々に累乗される。
【0241】
対数の虚数部だけをべき乗化すると考えると、これには、直交座標から極座標への変換の[この文献は図面を表示できません]
成分と同様に、すべての引数または角度情報が含まれる。これも伝統的に[この文献は図面を表示できません]
のべき乗とされるが、ここでは[この文献は図面を表示できません]
は問題ではない。虚数対数部分は角度として解釈することができるため、問題は今回[この文献は図面を表示できません]
であるが、収束の範囲を制限するために、回転対称性が利用することができるように[この文献は図面を表示できません]
の整数分数から残りを見つける必要がある。
【0242】
対数の底[この文献は図面を表示できません]
の代わりに、底２を使用することができるだろう。これは、[この文献は図面を表示できません]
という拡張を与える。これはさらに望ましくない無理数の値であるが、底[この文献は図面を表示できません]
を特に使用する場合は、回転数を表すために[この文献は図面を表示できません]
を使用できる。残念なことに、この基底の使用は、収束しない方法を導く。技術と収束範囲の間の最良のトレードオフは、底[この文献は図面を表示できません]
を得る。これは、各整数値が象限を表すことを意味する。４各倍数は、全体の回転である。次いで、ビットごとの演算を使用して、回転数を抽出することができる。象限の数は、実数と虚数および符号の切り替えとして出力に適用できる。この場合もやはり、複素対数の虚数成分の整数部分が四捨五入によって取られ、このようにして別々に処理されるという技法が得られる。この場合もやはり、複素対数の虚数成分の整数部分が四捨五入によって取られ、このようにして別々に処理されるという技術が得られる。
【0243】
これは、図３２の右パネル３２２０に示されるように、対数の小数部分が依然として[この文献は図面を表示できません]
の範囲内でのべき乗化が残り、これは、主アルゴリズムの収束範囲に収まる。図３２では、領域[この文献は図面を表示できません]
で対数−べき乗過程へのインプットのためにプロットされたべき乗での誤差、ここでは、右側は正の実数で、画像の下部は正の虚数である。５０の小数ビットと同様の反復回数では、黒い領域は、２−４５を超える精度に収束した値を表す。左パネル３２１０は、灰色が２−４５よりも大きい誤差に、マッピングする収束領域を示す。右パネル３２２０は、白で、最後に必要な収束領域[この文献は図面を表示できません]
を示す。
【0244】
対数入力の、小数部の主アルゴリズムについては後述する。
【0245】
逆の処理も達成され、そのべき乗化された実数部および虚数部を用いてその値をとり、それをその元の対数に減らす。はじめに、複素数値が属する象限を、４５°オフセットして見つける。これは、どの成分が最大の絶対値；実数または虚数およびその記号を持つかを、決定することによって達成する。これは、これらの象限のどれが、値が属する角度でオフセットしているかを特定する。次に、これを使用して、図３３の中央パネル３３２０に示すように、その面積が実線の正の部分を中心とする、第１象限に値を移動するための回転を計算するために使用する。図３３は、領域内[この文献は図面を表示できません]
のべき乗対数プロセスへの入力に対してプロットされたべき乗の誤差を示し、ここでは、右側は正の実数で、画像の下部は正の虚数である。５０の小数ビットと同様の反復回数では、黒い領域は２−４５を超える精度に収束した値を表す。左パネル３３１０は収束領域を示し、ここでグレーは２−４５より大きい誤差にマッピングされる。真ん中のパネル３３２０は、反転グレースケールを用いて実線を中心とする象限を示す。右パネル３３３０は、白で、最後に必要な収束領域[この文献は図面を表示できません]
を示す。
【0246】
対数出力の虚数部には、出力複素対数の等価虚数整数部があらかじめロードされている。４５°オフセットによって、複雑なブランチカットが実数ラインの負の部分に入るには、さらにテストが必要だが、２の補数システムラップアラウンドが使用されている場合、これは無視でき、必要なビットが与えられると、象限で測定された対数の虚数部は、[この文献は図面を表示できません]
の区間内にあるように、簡単に制限できる。
【0247】
この時点で、値が第１象限にあることは確かであるので、実数部は絶対値で虚数部と同等またはそれ以上であることが周知であり、そして実数部は正であることが周知である。対数の実部の整数部を見つけるために、固定小数点表現の先頭のゼロを数えることができる。固定小数点表現のビット数からこのカウントを減算した値は、事実上、[この文献は図面を表示できません]
であり、ここでｋはビットの表現に依存し、バイナリ対数の整数部分の表現でもあり、これは、対数レジスタの実数部分にすぐに加算できる。図３３の右側のパネルに示すように、実データと虚データの両方の既存データは、実部を区間[この文献は図面を表示できません]
に収めるために必要なビット位置の数だけシフトアップでき、これは、初期カウントの先行ゼロと、ビットの表現に依存するｋによって決まる。これは、出力実数対数の最も近い整数表現ではないかもしれないが、整数であり、プロセスの残りの部分を正しいフォーマットで収束領域に入れることに注意されたい。
【0248】
このプロセスの終わりに、出力複素数対数の整数部分が決定され、出力の小数部分と入力範囲[この文献は図面を表示できません]
は未決定のままであり、これは、主アルゴリズムの収束範囲内に収まる。対数出力の、実質的に小数部分のための主なアルゴリズムの議論は、以下で説明される。
【0249】
この結果は、関数[この文献は図面を表示できません]
が、図３４から分かるように、モデル化されたプロセスとなる。図３４では、べき乗部分３４１０および対数プロセス３４３０の領域のべき乗関数の範囲内で、べき乗部分３４２０および対数３４４０のプロセスの最終的な合成値を構成する対応部分を範囲とともに、プロットされた誤差を十分に小数部分についてのみプロットした。ノイズの多い灰色は、正しい収束が小さな数値誤差によって誤差測定値を支配するような領域を示しており、べき乗関数３４２０と対数３４４０について収束が至る所で達成されていることを示している。
【0250】
結果の底を変更しても、対数とべき乗のすべての特性があることを意味し、つまり、対数の実数部を足したり引いたりして、べき乗の乗算と除算を行うことができる。同様に、角度は同じように振る舞うが、それらは象限の整数倍、つまり回転に対して４、ハードウェアロジックで簡単に操作される値として表されるということだけである。対数形式では、虚数値の加減算は、複素数のべき乗化空間内でゼロを中心に、時計回りおよび反時計回りに回転させるのと同じである。実数部と虚数部が、異なる底で構成されているという事実に関係なく、実数部と虚数部を組み合わせた複素数値の場合、複素数の乗算と除算の実装に使用できる。
【0251】
これらの結果は、必要とされる前処理が、仮数およびべき乗形式によってすでにほとんど達成されるため、浮動小数点表現に自明に拡張することができる。底[この文献は図面を表示できません]
の値が絶対的に必要とされ、ほとんどの場合それらが必要とされない場合、出力対数の実部は[この文献は図面を表示できません]
で虚数は[この文献は図面を表示できません]
で乗算される。同様に、対数データを底[この文献は図面を表示できません]
とするフォーマットで取り込むためには、値を逆数[この文献は図面を表示できません]
とし、変換するために成分ごとに乗算してから、この方法でべき乗化するための初期範囲縮小を行うことができる。
【0252】
ＩＩＩ．一般的な手順
【0253】
新しいアルゴリズムは、対数の小数部分を見つけるか使用するために、ＢＫＭ技術と同様の方法で進行する。[この文献は図面を表示できません]
の値の集合があり、そこから各アルゴリズムステップ[この文献は図面を表示できません]
において[この文献は図面を表示できません]
が選択される。どの[この文献は図面を表示できません]
を使用するかの選択は、切り捨てられた数を、入力レジスタからのビットのサブセットであり得る、セットビット数と比較することによって行われる。[この文献は図面を表示できません]
が選択されると、[この文献は図面を表示できません]
（ビットシフト（潜在的に最上位の未使用ビットに丸めが適用される）と加算／減算）のバランスをとるために、事前計算対数のルックアップテーブルが必要であり、対数値を保持しているレジスタから減算する等価の対数を使用して、べき乗値を保持しているレジスタに乗算を適用する。プロセスが、対数を見つけることであるかべき乗を計算することであるかに応じて、システムを収束構成に配置するために、追加のシフト、追加の計算、またはその両方があり得る。いくつかのＮ回の繰り返しの終わりに、べき乗または対数のいずれかの所望の値は、繰り返しの回数Ｎおよび表現に使用されるビット数に依存する精度の範囲内であることが見出される。
【0254】
ＩＶ．対数のルックアップテーブル
【0255】
提案された方法は、これらの対数がここで、[この文献は図面を表示できません]
べき乗を反映しなければならないという点で既存の技術とは異なる。これにより、対数表の実部には、[この文献は図面を表示できません]
と、虚部[この文献は図面を表示できません]
または[この文献は図面を表示できません]
とをプリロードする必要があり、[この文献は図面を表示できません]
について所望の固定小数点フォーマットである。
【0256】
対数の底として、２と、[この文献は図面を表示できません]
とを選択したが、底の選択を変えると、対数システムのアルゴリズムの最初の反復が異なることはないと思われる。これらの値は、同じハードウェアを使用して、両方のプロセスを達成できる設計を表すために選択された。対数モードにべき乗モードとは異なる異なるハードウェアを使用した別の実装を使用する場合は、これが基本の最良の選択ではない可能性がある。
【0257】
Ｖ．対数の小数部の方法
【0258】
新しいアルゴリズムのための中心的な方法は、いくつかの新しい重要な相違点を除いて、ＢＫＭプロセスによって使用されるものに大部分従う。対数とべき乗の新しい底と、新しい手法に必要な収束領域の変化に対応するために、比較しきい値が新しく導き出された。必要な領域内での収束を達成するために、対数プロセスの最初の反復のために、わずかに修正された比較しきい値もある。
【0259】
ＶＩ．べき乗モード
【0260】
このセクションに記載されているべき乗モードはまた、べき乗値との複素乗算を提供するように修正される。これが達成することになっている場合、出力が正しいことであるならば、これは範囲縮小ステップの前に予めロードされなければならない。これは出力のべき乗値に置き換わるため、この値が必要な場合は使用してはならない。対数の整数部分（値ｚinteger、output）からプロセスの最後まで解を保存して適用するのを待つことも好ましい。これは、処理が行われる中間レジスタに必要な記憶容量を減らすことができるが、これは解の整数部分に必要な追加の記憶容量要件と比較して慮られることである。
【0261】
あるいは、追加のレジスタを作成し、これらの追加のレジスタで同等の演算が行われるようにすることによって、最終のべき乗値との乗算を並行して達成することができる。この方法では、べき乗プロセスは、出力べき乗値をほぼ任意の多くの他の複素数値と、複素数乗算することができる。
【0262】
入力対数の小数部分を入力と仮定すると、収束領域のアルゴリズム[この文献は図面を表示できません]
は次のようになる：
【0263】
１．[この文献は図面を表示できません]
とラベル付けられた、４つの基本レジスタがあると仮定する。並んで、補助複素数乗算に使用されたときの方法の動作を示すために、２つの追加のスレーブ乗算レジスタ[この文献は図面を表示できません]
と[この文献は図面を表示できません]
とがある。これらのレジスタの初期値は以下の通りである：[この文献は図面を表示できません]
ここで、事前乗算の要件がない場合は、[この文献は図面を表示できません]
である。スレーブ乗算レジスタも同様に以下のもので構築することができる：[この文献は図面を表示できません]
【0264】
２．インデックス[この文献は図面を表示できません]
として値[この文献は図面を表示できません]
を繰り返す：【0265】
３.[この文献は図面を表示できません]
を切り捨てて、３ビットになるように[この文献は図面を表示できません]
を作成する；１の符号ビットと２の補数の２つの整数ビットは、範囲が[この文献は図面を表示できません]
で最小の変化は１である。
【0266】
４．[この文献は図面を表示できません]
を切り捨てて、[この文献は図面を表示できません]
を形成するが、これは、１符号ビット、１整数ビット、２の補数で１小数ビットの３ビットを有し、範囲は[この文献は図面を表示できません]
で、最小の変化は０．５である。
【0267】
５.[この文献は図面を表示できません]
を決定するために、３ビット値をテストする：[この文献は図面を表示できません]
【0268】
６．シフトと算術演算を使用して、べき乗レジスタに乗算を適用する:[この文献は図面を表示できません]
[この文献は図面を表示できません]
および：[この文献は図面を表示できません]
[この文献は図面を表示できません]
乗算処理をこれらにも適用するには、[この文献は図面を表示できません]
や[この文献は図面を表示できません]
などの補助レジスタにも同じ操作を行う。
【0269】
７．対数表の対応するエントリをレジスタから減算する：[この文献は図面を表示できません]
[この文献は図面を表示できません]
【0270】
８．対数レジスタを１桁左にシフトして、その内容の値を２倍にする。[この文献は図面を表示できません]
【0271】
９．Ｎに達するまで、次の反復のためにステップ２に戻り、Ｎに達すると、その時点でレジスタはそれらの最終値を含む[この文献は図面を表示できません]
【0272】
および[この文献は図面を表示できません]
【0273】
プロセスの形式を理解したので、[この文献は図面を表示できません]
またはその両方を形成することによって、時には必要な領域で収束する他のテスト手順を見つけることは簡単であるが、我々は、異なるビット数または異なる比較値を使用するもので必要な値テストを最も簡単な周知の形式で指定することで、複雑さを軽減するように努めている。
【0274】
ＶＩＩ．対数モード
【0275】
このセクションに記載されている対数モードはまた、入力値との複素除算を提供するように修正されてもよい。これが達成することになっている場合、追加のレジスタを作成し、これらの追加のレジスタで同等の操作が行われるようにすることによってのみ、並行して発生する可能性がある。これは、元のレジスタが過負荷になる可能性がある補助乗算とは対照的であり、このモードではべき乗レジスタを変更するとアルゴリズムの収束が妨げられるため、ここでは達成できない。しかしながら、補助レジスタを使用することでこれを回避できるため、必要に応じて、このプロセスへの入力値を分母として、対数プロセスで他のほぼ任意の多数の複素数値を複素数除算することができる。
【0276】
１．出力対数の小数部を出力とすると、収束領域[この文献は図面を表示できません]
のアルゴリズムである：
【0277】
２．[この文献は図面を表示できません]
および[この文献は図面を表示できません]
とラベル付けされた、４つの基本レジスタがあると仮定する。並んで、補助複素除算に使用されたときの方法の動作を示すために、２つの追加のスレーブ除算レジスタ[この文献は図面を表示できません]
が存在する。これらのレジスタの初期値は以下のとおりである：[この文献は図面を表示できません]
スレーブ除算レジスタも、同様に構成される：[この文献は図面を表示できません]
レジスタ[この文献は図面を表示できません]
のいずれにも２倍算も−１．０も適用されないことに留意されたい。
【0278】
２．インデックス[この文献は図面を表示できません]
としての値[この文献は図面を表示できません]
を反復する：
【0279】
３．[この文献は図面を表示できません]
を切り捨てて[この文献は図面を表示できません]
を５ビット、１符号ビット、１整数ビット、２の補数の３小数ビットとして形成し、範囲は[この文献は図面を表示できません]
で、最小の変化は０．１２５である。
【0280】
４．[この文献は図面を表示できません]
を切り捨てて、[この文献は図面を表示できません]
１の符号ビット、２つの整数ビット、２の補数の２つの小数ビットの５ビットになるように形成し；範囲は[この文献は図面を表示できません]
で、最小の変化は０．２５である。
【0281】
５．２つの５ビット値をテストし（最初の反復の比較を修正して）、[この文献は図面を表示できません]
を決定する：[この文献は図面を表示できません]
【0282】
６．シフトと算術演算を使用して、べき乗レジスタに乗算を適用する：[この文献は図面を表示できません]
[この文献は図面を表示できません]
および：[この文献は図面を表示できません]
[この文献は図面を表示できません]
【0283】
[この文献は図面を表示できません]
などの補助レジスタにも同じようにして、これらに除算プロセスを適用する。
【0284】
７.[この文献は図面を表示できません]
を、[この文献は図面を表示できません]
に加算する：[この文献は図面を表示できません]
これは、[この文献は図面を表示できません]
に追加しない。
【0285】
８．べき乗レジスタを１桁左にシフトして、その内容の値を２倍にする。[この文献は図面を表示できません]
これは、[この文献は図面を表示できません]
の値を２倍にしない。
【0286】
９．対数表の対応するエントリをレジスタから減算する：[この文献は図面を表示できません]
【0287】
次の反復のため、[この文献は図面を表示できません]
に達するまでステップ２に戻り、このポイントで、でレジスタはそれらの最終値を含む：[この文献は図面を表示できません]
【0288】
および：[この文献は図面を表示できません]
【0289】
プロセスの形式を理解されるように、ある場合には、[この文献は図面を表示できません]
またはその両方を形成することによって、または異なるビット数または異なる比較値を使用して必要な領域で収束する他のテスト手順を見つけることができるが、我々は、必要な値テストを最も簡単な既知の形式で指定することで、複雑さを軽減するように努めている。
【0290】
ＶＩＩＩ．冗長数システム
【0291】
ＢＫＭ法のようなアルゴリズムは、データのローカルフローを遅くする技術をサポートするのに理想的であり、ハードウェアで実施される高速ロジックへの統合を可能にする。各反復での比較のために提示されたビット数の削減と相まって、算術ステップにキャリーセーブのような冗長数システム方式を使用することは、ここに提示されるようなアルゴリズムが高速低面積実装を有することを可能にし、これは、集積回路設計への組み込みに最適である。このような実装の利点には、より高速のクロック速度を達成できること、多数のビットを持つオペランドにも対応できることなどがあり、この設計は複雑な値の操作を含む問題に対する低コストの最新ソリューションの作成に理想的である。
【0292】
直観性を容易にするために、丸めはアルゴリズムの前述の説明では無視されており、本発明では正しい丸めが考慮されていないという意味として解釈されるべきではない。各ステップで適用されるときに丸めの適用を通して、エラーを減少させた結果を達成するために、いかなる算術演算でも直接使用されない最上位ビットを利用することは、ここで説明される結果の精度を高める可能性がある。
【0293】
ＩＸ．高次底のアプローチ
【0294】
ビット単位で評価する他のアルゴリズムと同様の方法で、複数のビットを一度に計算する方法を作成するために、ステップを併合することができる。出力は、２とは異なる底を持つものとして扱われるため、これは高次底アルゴリズムと呼ばれ、つまり、新しい高次底からの１桁がステップごとに計算される。この欠点は、ルックアップテーブルとシフトの選択もマージしなければならず、反復ごとにアルゴリズムが機能するための選択の、より大きなセットを作成することである。
【0295】
Ｘ．実数アルゴリズム
【0296】
多くのシステムは、本明細書に記載のアルゴリズムの実部のみを使用することから恩恵を受けることができる。例えば、ニューラルネットワークアルゴリズムによって、しばしば使用されるシグモイドロジスティック関数は、最初に実数２ｘを決定することによって、またはべき乗化処理を使用して、予め乗算ｅｘされた場合に生成することができる。次に、１を加算してプリセット定数をこの値で除算することが可能であるか、または、元の論文で説明されているＥモードに続いてＬモードを使用し続けると、ロジスティック関数は次のようになる：[この文献は図面を表示できません]
【0297】
ただし、分子ｋとその他の修正子ａおよびｂは、設定可能なままである。
【0298】
ＸＩ．収束性の提示
【0299】
図３５では、べき乗法および対数法の収束は、グラフ３５００によって示される。図３５は、必然的に収束するドメインにおいて可能な各値をテストすることが、ここでグラフ化された最悪の場合の[この文献は図面を表示できません]
エラーをもたらすことを示す。このグラフは、１兆を超える複雑なテスト値（再マッピングされた分数ドメイン内の各実数および虚数の値が、表示されている反復回数についてテストされている）の最悪の場合の動作をまとめたものである。グラフを比較可能にするために、対数の結果は、エラーをテストする前に、べき乗形式に変換した。これは、この方法がこの多数の有効ビットに対して収束的であることを保証する。この誤差は、乗算に丸めを追加することによってさらに減らすことができる。グラフは、精度と反復が増加するにつれて、誤差も同様に減少し、収束を示していることを示している。
【0300】
このグラフ３５００は、表内のビット数およびプロセス内の反復回数と同様の割合で縮小する収束領域内のすべての可能性のあるケースからの最悪ケースの絶対誤差を示す。これは、この固定小数点データのビット数に対する全ての可能な正しい入力に対して、アルゴリズムが正確なビット数で、正しい解に線形に首尾よく収束することを明らかにしている。
【0301】
ＸＩＩ．可能な請求項１．入力された「対数」複素数を得られる「べき乗」複素数に同時的に変換する方法であって、
ａ 前記得られる複素数の係数が前記実入力された複素数部分の実部べき乗に対する二乗からなること；および
ｂ 得られる複素数の引数が、前記虚数入力された複素数部分の実値で乗算された整数べき乗ロテーションの二乗を構成すること；
を含む、方法。
２． 入力された「べき乗」複素数を得られる「対数」複素数に同時的に変換する方法であって、
ａ 入力された複素数の実係数の実バイナリ対数が、得られる複素数の実部を形成すること；および
ｂ 前記入力された複素数の引数において表現された整数べき乗ロテーションの二乗の実数が、得られる複素数の虚部を形成すること；
を含む、方法。
３． 前記変換は、前記複素数のバイナリ表現についてのビットシフト、比較および計算を通して実質的に達成される、請求項１または２に記載の方法。
４． 入力および出力範囲の削減を、実質的に前記入力複素数の成分を別の部分とは別にして丸めることに対応する値の処理により実行する、請求項１に記載の方法。
５． 前記入力および出力範囲の削減を、実質的に前記出力複素数の成分を別の部分とは別にして丸めることに対応する値の処理により実行する、請求項２に記載の方法。
６． 前記出力複素値を、１に等しくない値に初期化して、前記反復の適用が非１の値を修正して本来出力することを意図する値の複素乗算を計算する、請求項１に記載の方法。
７． さらに、複素値を生成し、前記出力複素値に適用される操作がまた、これら更なる値に適用されて、これらのさらなる複素値に対する出力することを本来意図した前記複素値での乗算を与える、請求項１に記載の方法。
８． さらなる複素値を生成し、前記入力複素値に適用される操作がまた、これら更なる値に適用されて、これらのさらなる複素値の前記複素値入力による除算を与える、請求項２に記載の方法。
９． 冗長数システムを使用して、前記残された複素値に対する計算を反復する処理を実装する、請求項１に記載の方法。
１０． 冗長数システムを使用して、前記残された複素値に対する計算を反復する処理を実装する、請求項２に記載の方法。
１１． 請求項１〜１０のいずれかの方法を実行するために設計されたデバイス。
１２． 請求項１〜１０のいずれかの方法を実行するために構成されたデータプロセッサ。
【0302】
７．二次問題のＫＲＹＬＯＶ副空間から構築された音響場
【0303】
本開示は、固有システムとして表現されたときに二次問題を操作して、線形システムへの包含のための、ランク低減基底集合を生成する方法を説明する。縮小ランク基底を使用することにより、（潜在的に低減された）基底関数の最適情報保存集合を、複素値線形システムに含めるために自動的に導き出して、そのような基底関数の最適複素値再重み付けを見つけて、いくつかのテスト関数で説明されているように、目的を達成する。
【0304】
基底関数の縮小されたセットの導出、およびそれらが音場線形最適化問題の実装の適切な段階に、どのように含まれるかについて説明する。
【0305】
Ｉ．線形コントロール領域
【0306】
線形制御領域、これは複素加重重み付け関数、[この文献は図面を表示できません]
を用いて、圧力などの線形音響量との積分を実行する演算子を定義する:[この文献は図面を表示できません]
これは、音響場を評価するために使用する。[この文献は図面を表示できません]
が音響量の類似の位相分布と対になっている場合、複素共役は関数評価を最大化するために使用されることに留意されたい。この積分に対する任意の近似値（サンプリングなど）を使用することができ、あるいは実際には、これを使用して音圧場にわたって、任意の線形演算子を近似することができる。さらに線形であれば、ベクトル粒子速度（各ベクトル成分は場の調和運動により、複素数値である）を使用して、コントロール領域を最適化するようにさらなる積分を構成することができ、ここで：[この文献は図面を表示できません]
である。ここで、[この文献は図面を表示できません]
は媒体の粒子速度、[この文献は図面を表示できません]
は、各次元が複素数である３次元ベクトル場であり、これは、上記と同じ目的で重み付け関数として機能する。この方法では、異なる音響量パターンが線形的に評価されてもよく、ここで位相が[この文献は図面を表示できません]
におけるスカラーまたはベクトル関数によって生成された「テンプレート」と整列するとき、結果の特定の位相は、目的関数を最大化する。
【0307】
また、[この文献は図面を表示できません]
は、複数の外見上ばらばらな領域の単一の評価を作成するために、大部分がゼロ（またはゼロベクトルからなる）であり得ることにも留意されたい。この定式化は、線形音響量のフィールドにわたって、任意の線形演算子に拡張することができるので、任意の整数近似（サンプリングなど）を使用することができる。
【0308】
ＩＩ．二次問題と固有値問題
【0309】
位相アレイは、各トランスデューサをコントロールする、一連の複素数値の係数によって記述することができるマトリクス[この文献は図面を表示できません]
内に符号化された音響量を最大にするために、そのような入力[この文献は図面を表示できません]
ベクトルを用いて固定周波数で音響変換器のアレイを駆動することは、以下の意図をもった線形代数問題として表すことができる。[この文献は図面を表示できません]
【0310】
これは、固有ベクトル、ここでは[この文献は図面を表示できません]
の定義でもあるので、これはマトリクス[この文献は図面を表示できません]
の主要な固有ベクトルをとることによって、解くことができる。マトリクス[この文献は図面を表示できません]
を構築することは、まず、問題の音響量の大域空間における基底関数として、各トランスデューサによって放出された、線形音響量を考慮することによって達成することができる。これは、絶対位置[この文献は図面を表示できません]
によって記述する、単一ポイントで、圧力および他の音響量の最大化を達成するための方法である。しかし、記載されたコントロール領域積分の使用により、各個々の変換素子の場に放出される音響量を支配する[この文献は図面を表示できません]
は、[この文献は図面を表示できません]
と交換され、列ベクトルをもたらす：[この文献は図面を表示できません]
【0311】
線形性を呼び出すことによって、各トランスデューサを駆動する複素活性化値を乗算して、ベクトルを導き出すことができる：[この文献は図面を表示できません]
ここで、ｎはトランデューサの数である。このｍが線形場に作用する線形演算子の結果であると仮定すると、それはエルミート転置され、それ自身によって乗算されて二次形式を生成することができる。最後に下記式が得られる。[この文献は図面を表示できません]
【0312】
重み付け関数[この文献は図面を表示できません]
を修正し、マトリクス[この文献は図面を表示できません]
の固有ベクトルを解くことによって、様々な異なるビーム、表面、および他の音響現象が生成され最大化され、音場内に最適に作動される領域を作り出すことができる。
【0313】
前述のヌル空間生成技術（「空中ハプティックシステムにおける二次問題のための固有システム解」）もまた、音響浮上を記述するためにこの方法を拡張するために、ならびにヌルステアリングおよびミッド形成を可能にするために使用され、構造を記述するための位置で、フィールドにヌルポイントを導入することによって、ヌルステアリングおよび空中テクスチャの形成を可能にする。このプロセスはマトリクスのランクを増加するため、簡素化は使用されることが不可能であることに留意されたい。
【0314】
同様に、複数の最適化は、いくつかのコントロール領域について、[この文献は図面を表示できません]
個のマトリクスを互いに加算することによって、実行することができる。これは、コントロールポイントを表すマトリクスも含み得る。これは、増加したランクを持つマトリクスを生成することにより、簡素化の方法は使用するケースでは、不可能である。
【0315】
ＩＩＩ．高ランクマトリクスに対するランク低減近似
【0316】
キーポイントは、ランクの低い近似値が生成されるように、ランクの高いマトリクスを意図的に作成することができることである。例えば、独立したポイントを何百も持つシステムは、独立した基底関数で記述されるものが非常に少なく、忠実度を高め、複雑さを減らし、基底関数を必要に応じて自動的に作成することを可能にするよりも、高次元であるように見える基底セットの実現を作成するように、自動的に再編成される。
【0317】
既存のアプローチは、各コントロールポイント、またはコントロール領域のために、基底関数を生成する。各コントロールポイントまたは領域を最適化するための固有システムを作成し、それらを合計するか、またはヌルを追加することによって、固有システムはランクを獲得する。基底関数マトリクスおよび／またはヌルを追加することを繰り返されるプロセスを通して、高度なランクが最終マトリクスによって得られることが可能である。線形システムは、リソースに対する制約のために、最大で[この文献は図面を表示できません]
個のテスト関数を使用して、最大で[この文献は図面を表示できません]
個の基底関数を処理することしかできないかもしれない。この最適化マトリクスの高いランクはこの場合[この文献は図面を表示できません]
であり、ここで、個々の基底関数を生成するための[この文献は図面を表示できません]
個の最適化マトリクスは、追加で併合された。このアプローチを使用して、[この文献は図面を表示できません]
および／または[この文献は図面を表示できません]
であるという他の不可能な状況は、有限処理能力を有するこのシステムで解決可能な[この文献は図面を表示できません]
個の基底関数を有する、縮小基底集合で表される問題に帰着される。
【0318】
この縮小基底集合は、基底関数に対する複数の最適化からなる最終マトリクスの[この文献は図面を表示できません]
個の最も主要な固有ベクトルと、実質的に同じマトリクス部分空間を占める、クリロフ部分空間を計算することによって、見出される。クリロフ部分空間は、主要な固有ベクトルが見つけられ、次いで残りのマトリクスから取り除かれるアーノルディ反復法のような、標準的な技術を用いて見つけられ、それから[この文献は図面を表示できません]
個のそのような固有ベクトルが見つかるまでこれが繰り返される。
【0319】
ＩＶ．混合アプローチ
【0320】
基底関数が上位マトリクスから抽出されるこのアプローチは、既存の特定の基底関数アプローチと混合して適合させることができる。例えば、２つのコントロール領域またはポイントを必須項目として指定し、線形システムへの基本セット入力に追加することができる。それから、それらが二次最適化、したがってそれ自身の最適化マトリクスによって特定されるということを考えると、残りの基底関数は、単一のマトリクスに加えることができ、そのマトリクスは、線形システムに入力された基底関数のセットの残りの数を構成する、最も支配的な[この文献は図面を表示できません]
の固有ベクトルを見つけるために、このクリロフ部分空間法を用いて分解される。この方法で、最初に指定された２つの基底関数は、線形システムに直接インポートされるため保証されているが、セットの残りの関数は近似する。このシステムはさらに用いられ、ここで、底関数の集合数は、上位ランクの最適化マトリクスによって、それぞれ記述される基底関数の互いに素な集合から特定され、そして基底関数の集合数が取り出される。縮小基底関数を導出するために使用できる、さらなる技術は、パワーしきい値を超えるまで、マトリクスから基底関数を取り続けることである。このしきい値は、結合された基底関数から構成されるマトリクスごとに異なる可能性があるため、異なるパワーの固有ベクトルを生成し、この場合、固有ベクトルを選択することで固有ベクトルを生成でき、各固有ベクトルのパワーは、実際の優先順位の値と対応する固有値の大きさの積に基づいて、各マトリクスから計算される。
【0321】
８．テスト関数を使用した線形システムソリューション
【0322】
基底関数セットが基底関数セットとは別に定義される、特定の目標を達成するかどうかを「テストする」方法は、価値があるだろう。これらの目標が達成されたかどうかの定量的評価は、最適化に提供され、これは、これらのテスト関数によって定義される音場内の所望の特性を生成するために使用される。
【0323】
テスト関数および複素数値線形システム最適化における、それらの役割について説明する。
【0324】
Ｉ．線形コントロール領域
【0325】
線形コントロール領域、これは、複素数値重み付け関数[この文献は図面を表示できません]
を用いて、圧力などの線形音響量との積分を実行する演算子として定義する：[この文献は図面を表示できません]
これは、音響場で評価するために使われる。[この文献は図面を表示できません]
が、音響量の類似の位相分布と対になっている場合、複素共役は、関数評価を最大化するために使用されることに留意されたい。この積分に対する任意の近似値（サンプリングなど）を使用することができ、または実際には、これを使用して音圧場にわたって任意の線形演算子を近似することができる。ベクトル粒子速度（ここで各ベクトル成分は、場の調和運動のための複素数値である）を用いて、コントロール領域を最適化するように、さらに積分されてもよい、ここで：[この文献は図面を表示できません]
[この文献は図面を表示できません]
は、媒体の粒子速度、[この文献は図面を表示できません]
は、それぞれの次元が複素数値の３次元ベクトル場であり上記と同じ目的で重み付け関数として機能する。この方法では、異なる音響量パターンが線形的に評価されてもよく、ここで位相が[この文献は図面を表示できません]
におけるスカラーまたはベクトル関数によって生成された「テンプレート」と整列するとき、結果の特定の位相は目的関数を最大化する。
【0326】
また、[この文献は図面を表示できません]
は、複数の外見上ばらばらな領域の単一の評価を作成するために、大部分がゼロ（またはゼロベクトルからなる）であり得ることにも留意されたい。この定式化は、線形音響量のフィールドにわたって、任意の線形演算子に拡張することができるので、任意の整数近似（サンプリングなど）を使用することができる。
【0327】
ＩＩ．基底関数セット
【0328】
入力として[この文献は図面を表示できません]
個の基底関数のセットが与えられると、[この文献は図面を表示できません]
個の基底関数の最適な線形結合を見つけるために、線形システムを通して[この文献は図面を表示できません]
個のテスト関数が使用され、[この文献は図面を表示できません]
個のテスト関数の所望の客観的評価を達成する。これら[この文献は図面を表示できません]
個のテスト関数は、線形コントロール領域評価の形をとることができる。上記の目的[この文献は図面を表示できません]
の形式を使用すると、次のようになる：[この文献は図面を表示できません]
または、単一のコントロールポイント定義を通して、次のように評価することもできる：[この文献は図面を表示できません]
ここで、上線は複素共役を表し、[この文献は図面を表示できません]
は、三次元位置を表し、[この文献は図面を表示できません]
は全てのトランスデューサ素子をパラメータ化し、[この文献は図面を表示できません]
はトランスデューサの圧力場の共役を乗じた、トランスデューサ活性化である。法線方向が定義されている場合、粒子速度によって定義されるコントロールポイントの等価形式も、評価子として使用できる：[この文献は図面を表示できません]
ただし、間違った方向を向いているトランスデューサ[この文献は図面を表示できません]
も、[この文献は図面を表示できません]
の場合、その寄与のセットをゼロに設定し、省略することができる）。
【0329】
両方の単一コントロールポイント定義において、既存のアプローチは、コントロール領域評価の単位重み付け関数の場合を伴う、特別な単一ポイント積分評価と見なすことができる。ここでの目的関数は、余分な副作用のない単純な音響量である。
【0330】
ＩＩＩ．主要な固有ベクトルの発見
【0331】
[この文献は図面を表示できません]
個のテスト関数の各々を用いて、[この文献は図面を表示できません]
個の基底関数の各々を評価することにより、複素数値である矩形マトリクスが得られる。前述のように、主要な固有ベクトルは、[この文献は図面を表示できません]
個のテスト関数によって定義される空間を最も良く張る[この文献は図面を表示できません]
個の基底関数の組み合わせを表す。しかし、これは必ずしも正方マトリクスではないため、特異値分解がより適切であるが、唯一必要な固有ベクトルは支配的なものである。その支配的な固有ベクトルが、システム基底関数の支配的な組み合わせの要素を表す、正方エルミート行列を、構築することができる。固有システムは、ここで、以下のように記載される：[この文献は図面を表示できません]
ここで、[この文献は図面を表示できません]
は、位相が[この文献は図面を表示できません]
個のテスト関数のパラメータ空間のカバーを生成する係数を記載する。したがって、優勢固有ベクトルは、単一の組み合わせベクトルで達成可能な空間の最良のカバーを表す成分[この文献は図面を表示できません]
を有し、その位相はこのカバーを最適化する。
【0332】
テスト関数の数が基底関数の数と等しいため、[この文献は図面を表示できません]
である場合でも、既存の固有システムアプローチは必ずしも適用可能ではない場合があることに、留意されたい。これは、マトリクスがエルミーションではなくなる可能性があるため、上記の手順を完了せずに、レイリー−リッツ（Ｒａｙｌｅｉｇｈ−Ｒｉｔｚ）法をこのマトリクスに直接適用すると、方法が主要の固有ベクトルに収束できなくなる可能性があるためである。
【0333】
位相データが指定されていないか、または不完全である場合であっても、主要な固有ベクトルを位相情報のために使用すべきである。固有ベクトルは、各反復の終わりに、振幅データを除去するように修正され、位相がゼロに平均化され、前の解から別の時間ステップで位相進行が制限され、または、既存する定常の位相の関係を有する、いくつかの基底関数に対して１組の位相解を有し、このセットの補数は、それら自体と既存の位相の定常セットとの間の位相で最適化されている。この方法では、この矩形マトリクスシステムは、基底集合を生成する関数が、基底集合の目的を評価するために使用される関数と同じ関数である条件について、特化した既存の固有システムアプローチと同様の効果を達成するために使用される。
【0334】
ＩＶ．線形システムソリューション
【0335】
複素数値線形システムは、修正された役割、すなわちテスト関数の所与の値を最もよく再現する、基底関数の集合を見つけるという役割を持つ。ここで、線形システムは、下記の矩形マトリックスの形態を取る：[この文献は図面を表示できません]
【0336】
マトリクスが矩形であるため、テスト関数の数が基底関数の数よりも多く、その結果、複素数値の最小二乗システムになるか、またはテスト関数の数が基底関数の数よりも少なくなる可能性があり、複素数値の最小ノルムシステムが得られる。いずれの場合も、標準的な分解を使用して解ベクトルを見つけることができ、その解は、解ベクトルの各成分が、試験関数の所与の値を表現する場合に音響場を生じさせる基底関数のそれぞれに適用される係数である。
【0337】
これらのテスト関数は、あるポイントにおける音圧または方向性粒子速度、あるいは空間内のある領域にわたる圧力、および方向性粒子速度の混合を表すことができる。これらの関数のどのような組み合わせも[この文献は図面を表示できません]
個のテスト関数を構成することができ、それは線形システムの最終ステップで解かれる。
【0338】
基底関数の係数が分かれば、トランスデューサを作動させるための最終値は、次のように拡張される：[この文献は図面を表示できません]
【0339】
単一の単調搬送波周波数（状況の変化により多かれ少なかれ有効であり得る）の仮定を与える入力信号は、この[この文献は図面を表示できません]
の成分を各変換器に適用することによって構成される。
【0340】
９．結論
【0341】
前述の実施形態の様々な特徴は、改良した触角のシステムの多数のバリエーションを生産するために選択され、そして組み合わされてもよい。
【0342】
前述の明細書において、特定の実施形態が説明された。しかしながら、当業者であれば、特許請求の範囲に記載されている本発明の範囲から逸脱することなく、様々な修正および変更が可能であることを理解するであろう。従って、明細書および図面は限定的な意味ではなく、例示的な意味で見なされるべきであり、そのような修正はすべて本教示の範囲内に含まれることが意図される。
【0343】
利益、利点、問題に対する解決策、および何らかの利益、利点、または解決策を生じさせるか、またはより顕著にする可能性があるあらゆる要素は、いずれかもしくは全ての請求項の重要な、必要の、本質的な特徴、要素として解釈されるべきではない。本発明は、本願の係属中に行われたいかなる修正、および発行された特許請求の範囲のすべての均等物を含む、添付の特許請求の範囲によってのみ定義される。
【0344】
さらに、本明細書では、第１および第２、上下などの関係用語は、必ずしもそのような実態、または動作間の実際のそのような関係または順序を必要とせず、暗示することなく、ある実態または動作を、別の実体または動作から区別するためにのみ、使用される。用語「構成する」、「構成している」、「有する」、「有している」、「含む」、「含んでいる」、「包含する」、「包含している」またはそれらの任意の他のバリエーションは、非排他的包含に及ぶものであることを意図しており、要素のリストを構成する、有する、含む、包含するプロセス、方法、物品、または装置は、それらの要素だけを含むのではなく、そのようなプロセス、方法、物品、または装置に明示的に挙げられないか、または明示的にリストされた固有の他の要素を含む。「構成する」、「有する」、「含む」、「包含する」などの要素は、プロセス、方法、記事、または要素を構成する、有する、含む、包含する装置の中に、追加の同一の要素が存在することを排除するものではない。「ａ」および「ａｎ」という用語は、本明細書で、他に明示的に述べられていない限り、１つまたはそれ以上として定義される。用語「実質的に」、「本質的に」、「およそ」、「約」またはそれらの任意の他の類型は、当業者によって理解されるものに近いとものとして定義される。本明細書で使用する「結合」という用語は、必ずしも直接ではなく、必ずしも機械的にではないが、接続されていると定義される。特定の方法で「構成」されているデバイスまたは構造は、少なくともそのように構成されているが、記載されていない方法で構成することもできる。
【0345】
開示の要約は、読者が技術的開示の性質を迅速に確認することを可能にするために提供され、これは、特許請求の範囲または意味を、解釈または限定するために使用するものではないとの理解のもとに提出される。加えて、前述の詳細な説明では、開示を簡素化する目的で、多様な特徴が、様々な実施形態において、一緒にグループとしてまとめられていることが理解できる。本開示の方法は、請求された実施形態が、各請求項に明示的に記載されているよりも、多くの特徴を必要とするという意図を反映していると解釈されるべきではない。そうではなく、添付の特許請求の範囲が反映するように、発明の主題は、単一の開示された実施形態の全ての特徴より少ない特徴にある。したがって、添付の特許請求の範囲は、ここで詳細な説明に組み込まれ、各特許請求の範囲は、別々に特許請求される主題としてそれ自体で存在する。【図1】
[この文献は図面を表示できません]
【図2】
[この文献は図面を表示できません]
【図3】
[この文献は図面を表示できません]
【図4】
[この文献は図面を表示できません]
【図5】
[この文献は図面を表示できません]
【図6】
[この文献は図面を表示できません]
【図7】
[この文献は図面を表示できません]
【図8】
[この文献は図面を表示できません]
【図9】
[この文献は図面を表示できません]
【図10】
[この文献は図面を表示できません]
【図11】
[この文献は図面を表示できません]
【図12】
[この文献は図面を表示できません]
【図13】
[この文献は図面を表示できません]
【図14】
[この文献は図面を表示できません]
【図15】
[この文献は図面を表示できません]
【図16】
[この文献は図面を表示できません]
【図17】
[この文献は図面を表示できません]
【図18】
[この文献は図面を表示できません]
【図19】
[この文献は図面を表示できません]
【図20】
[この文献は図面を表示できません]
【図21】
[この文献は図面を表示できません]
【図22】
[この文献は図面を表示できません]
【図23】
[この文献は図面を表示できません]
【図24】
[この文献は図面を表示できません]
【図25】
[この文献は図面を表示できません]
【図26】
[この文献は図面を表示できません]
【図27】
[この文献は図面を表示できません]
【図28】
[この文献は図面を表示できません]
【図29】
[この文献は図面を表示できません]
【図30】
[この文献は図面を表示できません]
【図31】
[この文献は図面を表示できません]
【図32】
[この文献は図面を表示できません]
【図33】
[この文献は図面を表示できません]
【図34】
[この文献は図面を表示できません]
【図35】
[この文献は図面を表示できません]
【手続補正書】
【提出日】2019年10月31日
【手続補正1】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項1】
ｉ）既知の相対的な位置および配向を有する位相トランスデューサアレイから音響場を生成すること；
ｉｉ）複数のコントロールポイントであって、それぞれ前記位相トランスデューサアレイに相対して既知の空間的関係を有する前記複数のコントロールポイントを定義すること；
ｉｉｉ）前記位相トランスデューサアレイにより生成される量の物理場において既知の方向に所望する場の値の実質的な生成を最適化することであって、前記量がベクトル量であること；
を含む方法。
ｉｉ）複数のコントロールポイントであって、それぞれ前記位相トランスデューサアレイに相対して既知の空間的関係を有する前記複数のコントロールポイントを定義すること；
ｉｉｉ）前記位相トランスデューサアレイにより生成される量の物理場において既知の方向に所望する場の値の実質的な生成を最適化することであって、前記量がベクトル量であること；
を含む方法。
【請求項2】
さらに、固有値システムの合計を使用して最適化マトリクスを生成することをさらに含む、請求項１に記載の方法。
【請求項3】
さらに、少なくとも１つの場において、少なくとも１つのヌルポイントを使用して最適化マトリクスをさらに生成することを含む、請求項１に記載の方法。
【請求項4】
前記所望する場の値の少なくとも１つが、ホスト媒体内での既知のポイントの粒子の速度を実質的に記述する、請求項１に記載の方法。
【請求項5】
前記所望する場の値の１つが、波の運動量を実質的に記述する、請求項１に記載の方法。
【請求項6】
さらに、複数の基底関数を計算するための第１の値を使用して最適化を生成する、請求項１に記載の方法。
【請求項7】
複数の前記基底関数は、前記物理場の第２の値を特定する少なくとも１つのコントロールポイントを実現するために使用され、前記第２の値は、前記第１の値とは異なる、請求項６に記載の方法。
【請求項8】
ｉ）既知の相対的な位置および配向を有する位相トランスデューサアレイから音響場を生成すること；
ｉｉ）複数のコントロールポイントであって、それぞれ前記位相トランスデューサアレイに相対して既知の空間的関係を有する前記複数のコントロールポイントを定義すること；
ｉｉｉ）前記位相トランスデューサアレイにより生成された量の物理場における既知の方向に沿って所望する場の値の絶対的な大きさを実質的に最大化するように最適化することであって、前記量がベクトル量であること；
を含む方法。
ｉｉ）複数のコントロールポイントであって、それぞれ前記位相トランスデューサアレイに相対して既知の空間的関係を有する前記複数のコントロールポイントを定義すること；
ｉｉｉ）前記位相トランスデューサアレイにより生成された量の物理場における既知の方向に沿って所望する場の値の絶対的な大きさを実質的に最大化するように最適化することであって、前記量がベクトル量であること；
を含む方法。
【請求項9】
さらに、固有値システムの合計を使用して最適化マトリクスを生成することをさらに含む、請求項１に記載の方法。
【請求項10】
さらに、少なくとも１つの場において、少なくとも１つのヌルポイントを使用して最適化マトリクスをさらに生成することを含む、請求項８に記載の方法。
【請求項11】
前記所望する場の値の少なくとも１つが、音響強度を記述する、請求項８に記載の方法。
【請求項12】
前記所望する場の値の１つが、エネルギーフラックスベクトルを実質的に記述する、請求項８に記載の方法。
【請求項13】
前記所望する場の値の少なくとも１つが、ホスト媒体内での既知のポイントの粒子の速度を実質的に記述する、請求項８に記載の方法。
【請求項14】
前記所望する場の値の１つが、波の運動量を実質的に記述する、請求項８に記載の方法。
【請求項15】
さらに、複数の基底関数を計算するための第１の値を使用して最適化を生成する、請求項８に記載の方法。
【請求項16】
複数の前記基底関数は、前記物理場の第２の値を特定する少なくとも１つのコントロールポイントを実現するために使用され、前記第２の値は、前記第１の値とは異なる、請求項１５に記載の方法。
【請求項17】
ｉ）既知の相対的な位置および配向を有する位相トランスデューサアレイから音響場を生成すること；
ｉｉ）複数のコントロールポイントであって、それぞれ前記位相トランスデューサアレイに相対して既知の空間的関係を有する前記複数のコントロールポイントを定義すること；および
ｉｉｉ）前記位相トランスデューサアレイにより生成された少なくとも１つの物理場における場の値の所望する領域での重み付け積分を実質的に最大化するように最適化すること；
を含む方法。
ｉｉ）複数のコントロールポイントであって、それぞれ前記位相トランスデューサアレイに相対して既知の空間的関係を有する前記複数のコントロールポイントを定義すること；および
ｉｉｉ）前記位相トランスデューサアレイにより生成された少なくとも１つの物理場における場の値の所望する領域での重み付け積分を実質的に最大化するように最適化すること；
を含む方法。
【請求項18】
さらに、固有値システムの合計を使用して最適化マトリクスを生成することをさらに含む、請求項１７に記載の方法。
【請求項19】
さらに、少なくとも１つの場において、少なくとも１つのヌルポイントを使用して最適化マトリクスをさらに生成することを含む、請求項１７に記載の方法。
【請求項20】
前記所望する場の値の少なくとも１つが、音響強度を記述する、請求項１７に記載の方法。
【請求項21】
前記所望する場の値の１つが、エネルギーフラックスベクトルを実質的に記述する、請求項１７に記載の方法。
【請求項22】
前記所望する場の値の少なくとも１つが、音響圧を実質的に記述する、請求項１７に記載の方法。
【請求項23】
前記所望する場の値の少なくとも１つが、ホスト媒体内での既知のポイントの粒子の速度を実質的に記述する、請求項１７に記載の方法。
【請求項24】
前記所望する場の値の１つが、波の運動量を実質的に記述する、請求項１７に記載の方法。
【請求項25】
さらに、複数の基底関数を計算するための第１の値を使用して最適化を生成する、請求項１７に記載の方法。
【請求項26】
複数の前記基底関数は、前記物理場の第２の値を特定する少なくとも１つのコントロールポイントを実現するために使用され、前記第２の値は、前記第１の値とは異なる、請求項２５に記載の方法。
【請求項27】
複数の複素値入力を有するハードウェアを操作すること；
前記複数の複素値入力のそれぞれを、べき乗関連の処理において複素値出力に変換すること；
前記複素値入力の前記実部の２を底とする（バイナリ）べき乗の大きさを有する係数と；
完全な回転の整数乗の二乗として定義される角度ユニットで乗算した前記複素値入力の前記虚部である角度引数と
を含む方法。
前記複数の複素値入力のそれぞれを、べき乗関連の処理において複素値出力に変換すること；
前記複素値入力の前記実部の２を底とする（バイナリ）べき乗の大きさを有する係数と；
完全な回転の整数乗の二乗として定義される角度ユニットで乗算した前記複素値入力の前記虚部である角度引数と
を含む方法。
【請求項28】
前記変換は、前記複素数のバイナリ表現に対するビットシフト、比較および計算を通して達成される、請求項２７に記載の方法。
【請求項29】
さらに、前記複数の複素値入力のそれぞれの成分の実質的な丸めに対応する値を、さらなる部分と別別に処理して、入力および出力範囲の削減を達成することを含む請求項２７に記載の方法。
【請求項30】
さらに、１ではない値に出力複素値を初期化して、反復の適用が、非１の値を変更して出力することを意図する元々の値による複素乗算を計算することを含む、請求項２７に記載の方法。
【請求項31】
さらに、複素値を生成して、前記複素値出力に適用される操作をまたさらなる複素値にも適用して、出力することを意図する元々の複素値による前記さらなる複素値の乗算を与える、請求項２７に記載の方法。
【請求項32】
さらに、冗長数システムを使用して残りの複素値に適用する反復計算の処理を実装することを含む、請求項２７に記載の方法。
【請求項33】
複数の複素値入力を有するハードウェアを操作すること；
前記複数の複素値入力のそれぞれを、対数を求めることに関連する処理において複複素値出力に変換することであって、前記出力複素値出力が、
前記複素値入力の係数の２を底とする（バイナリ）対数に等しい実部と；
前記複素値入力の角度引数において出現する角度ユニットの数に等しい虚部であって、前記角度ユニットが前記完全な回転の整数乗の二乗である、
方法。
前記複数の複素値入力のそれぞれを、対数を求めることに関連する処理において複複素値出力に変換することであって、前記出力複素値出力が、
前記複素値入力の係数の２を底とする（バイナリ）対数に等しい実部と；
前記複素値入力の角度引数において出現する角度ユニットの数に等しい虚部であって、前記角度ユニットが前記完全な回転の整数乗の二乗である、
方法。
【請求項34】
前記変換は、、前記複素数のバイナリ表現に対するビットシフト、比較および計算を通して達成される、請求項３３に記載の方法。
【請求項35】
さらに、実質的に象限の回転および前記複数の複素値入力のそれぞれの成分をさらなる部分とは別々にビットシフトすることに対応する値の処理を行い、入力および出力の削減を達成する、請求項３３に記載の方法。
【請求項36】
さらに、非ゼロの値を有する出力複素値を生成して、前記複素値出力に適用される操作をまたさらなる複素値にも適用して、出力することを意図する元々の複素値による複素加算を計算する、請求項３３に記載の方法。
【請求項37】
さらに、複素値を生成して、前記複素値入力に適用される操作をまたさらなる複素値にも適用して、前記さらなる複素値の前記複素値入力による除算を与える、請求項３３に記載の方法。
【請求項38】
さらに、冗長数システムを使用して残りの複素値に適用する反復計算の処理を実装することを含む、請求項３３に記載の方法。
【国際調査報告】
[この文献は図面を表示できません]
[この文献は図面を表示できません]
[この文献は図面を表示できません]
[この文献は図面を表示できません]