(57)【要約】要約以下のシステムを含む、位相アレイシステムを駆動するための、多様な技術が記載されており、特に、空中触覚、パラメトリックオーディオ、音響浮揚および音響イメージングへの応用を有する音響位相アレイを対象とし：１）一貫した触覚体験を提供するために、空気の変化の影響を軽減することが可能であり；２）空気中でトラップポイントを生産し；３）より一貫した触覚効果の生産のためのベクトルに関して、位相アレイ最適化を定義し；４）コントロールされた音響場を介して、空間内の１つまたは複数のコントロールポイントまたは領域を定義し；５）音響基底関数の構築のために、被約表現法を使用し；６）大量のスループットに対して、複素数値関数の効率的な評価を実行し；７）マトリクスのクリロフ小領域を生成し；並びに８）音響基底関数を作成するために使用するものとは異なる、コントロールポイントおよび／または領域によって、記述される目的を最大にする。
2020502906.tif 000002
【特許請求の範囲】
【請求項１】
音響場においてコントロールポイントを生成すること；
各コントロールポイントに法線ベクトルおよび方向を割り当て、位相アレイソルバーに波形運動量を組み込むこと；および
関心あるポイントで、前記波形運動量の
2020502906.tif 000537
成分を解くこと
を含む方法。
【請求項２】
2020502906.tif 000538
成分が、下記式で定義される
2020502906.tif 000539
でそれぞれ特定され、
2020502906.tif 000540
上式中、
2020502906.tif 000541
は、コントロールポイント
2020502906.tif 000542
でのトランスデューサ
2020502906.tif 000543
からの進行波の方向を記述するベクトルであり、
2020502906.tif 000544
は、前記運動量が測定される位置を横断するコントロールポイント
2020502906.tif 000545
の法線ベクトルの方向である
請求項１に記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
関連出願
本出願は、以下の８つの米国仮特許出願の利益を主張するものであり、それらの全ては、参照によりその全体が本明細書に組み込まれる：
【０００２】
1.  出願番号  ６２／４３３，５４４，  出願日２０１６年１２月１３日
【０００３】
2.  出願番号  ６２／４３６，４５７，  出願日２０１６年１２月２０日
【０００４】
3.  出願番号  ６２／４３７，７９１，  出願日２０１６年１２月２２日
【０００５】
4.  出願番号  ６２／４３８，５１２，  出願日２０１６年１２月２３日
【０００６】
5.  出願番号  ６２／５９４，５７４，  出願日２０１７年１２月  ５日
【０００７】
6.  出願番号  ６２／５９４，６８７，  出願日２０１７年１２月  ５日
【０００８】
7.  出願番号  ６２／５９４，８３６，  出願日２０１７年１２月  ５日
【０００９】
8.  出願番号  ６２／５９４，８７４，  出願日２０１７年１２月  ５日
【００１０】
開示分野
本開示は、概ね、位相アレイシステムを駆動するために改良した技術に関し、主に空中ハプティックス、パラメトリックオーディオ、音響浮揚および音響イメージングを生成するために設計された音響場を生成することに関する。
【背景技術】
【００１１】
背景
我々が「音響場」と呼ぶ音エネルギーの連続分布は、空中での触覚フィードバックを含む、様々な用途に使用することができる。
【００１２】
本音響場では、１つまたは複数のコントロールポイントを明確にすることができる。これらのコントロールポイントは、信号で振幅変調することができ、その結果、空中で振動触覚フィードバックを生成することができる。フィードバックを作成する別の方法は、振幅が変調していないコントロールポイントを作成し、それらを空間的に動かして感じられる時空間変調を生成することである。
【００１３】
装置によって作成されるコントロールポイントの最小サイズおよび最大パワーは、生成される１つまたは複数のコントロールポイントの配列レイアウト、形状および相対位置に依存する。これにより、ある場所のコントロールポイントの動作を、別の場所の動作と一致させることが不可能な場合がある。コントロールポイントが空間を通って動的に動かされるとき、さらなる複雑さが起こり得る。 ユーザが経験した点の触覚機能の結果として生じる変化は望ましくなく、しばしば直観的ではなく、結果として混乱を招き、経験的な没入感を喪失する。空気の変化の影響を軽減することができ、一貫した触覚体験を提供するためにそれを逆転させることができるシステムは、商業的に価値があるだろう。
【００１４】
さらに、一般的な非線形最適化システムを使用して、空気中にトラップポイントを生成できることが示された。しかし、この目的のために使用されるメカニズムは、位相アレイシステム上の浮揚に適した潜在的な場を評価しないこと、不必要に高価で時間がかかる最適化プロセスを使用すること、位相操作のみに限定される欠点を含む。
【００１５】
既存のシステムによって作成したトラップは、それらが可能であるほど強力ではなく、これらの欠点のために正確に拡張することができないので、これらの問題に対処することは重要である。これは、音響的な浮揚を利用して、マルチモーダルアプローチを促進するための、さらなるビジュアルを提供する、システムおよび触覚システムの代替の使用を困難にする。
【００１６】
さらに、空間で、１つまたは複数のコントロールポイントを明確にすることによって、音響場をコントロールすることができる。各ポイントには、コントロールポイントにおける所望の振幅に等しい値を割り当てることができる。次に、トランスデューサの物理的セットをコントロールして、コントロールポイントにおいて所望の振幅を示す音場を作り出すことができる。
【００１７】
既存のシステムの条件によって、音響位相アレイの自然な指向性成分は、ほとんど注目されていなかった。例えば、目的がエネルギーを狭い領域に集中させることである高周波治療医療システムは、エネルギーを注入するための最適解を決定するためにスカラー圧力（関与する空気分子の数の変化）によって十分に近似され、したがってベクトル解を必要としない。位相アレイベクトルの挙動（空気分子の運動方向を表す）に焦点を合わせることは、ほとんどの場合無視できるので、これを効率的に達成することはまだ優先事項ではない。
【００１８】
しかし、空中で触覚フィードバックを作成するには、人間の肌と振動を交換する必要がある。これは完全に波面との接触面と波によって運ばれる運動量に依存し、これは波のベクトルの振る舞いとの関係の結果である。これだけでなく、システムは、トランスデューサ素子間で入射角がますます異なるように焦点が合わせられているので、運動量と圧力との間の仮定された線形関係は崩壊する。このシステムではベクトル成分を無視することは不可能であり、スカラー圧力は空中触覚効果を生み出すための、適切な解を一貫して見つけるには不十分である。この理由から、この運動量の影響をデバイスに認識させることは、商業的に重要である。
【００１９】
さらに、空間で１つまたは複数のコントロールポイントを明確にすることによって、音響場をコントロールすることができる。各ポイントには、コントロールポイントにおける所望の振幅に等しい値を割り当てることができる。次に、トランスデューサの物理的セットをコントロールして、コントロールポイントにおいて、所望の振幅を示す音響場を作り出すことができる。
【００２０】
固有ベクトルの計算は、行列代数におけるいくつかの非常に具体的な純粋に二次の制約付き問題の解に対応することが示される。空中触覚における特定の問題を考慮して、固有ベクトルがトランスデューサのコントロール係数を解く行列を形成することができれば、その解決策は、既存の技術に対して単純化を示す装置実装の基礎を形成することができる。
【００２１】
さらに、基底関数を記述するために簡約表現法を使用するためには、まず音響場の詳細に関して有用な関数の簡単な（低エントロピー）記述を可能にする簡略表現を開発しなければならない。
【００２２】
本開示は、理論上、ユーザの視点からアクセス可能なものによるだろう。
【００２３】
そのようなアクセス可能な低エントロピー表現を発見したならば、また、それをその等価なトランスデューサ基底（後の解決策において潜在的に多数の基底関数のうちの１つのトランスデューサ音響場の基底を記述するトランスデューサ基底セット）に変換する方法となる。これにより、見込みのあるユーザは音響場内で、多様な異なる構造を設計することが可能になり、焦点、ベッセルビーム、ボトルビーム、ベンディングビームなど、一度に多くの異なる効果の形成を導く現場でのエネルギーの流れを最大限に操作するための特注システムを構築することが可能になる。これらは、コントロールポイントの概念を、より豊かな音響場コントロール方法をカプセル化するコントロール領域に拡張するだろう。
【００２４】
触覚における運動量移動を考慮するときに、粒子速度を波の方向の代用として、使用するために対処しなければならない、いくつかの問題があることもまた発見された。
【００２５】
さらに、同じ２番目のオペランドを使用して効率的な並列複素数値乗算を実行することは、要件の資料に概説されているものと一致するパフォーマンスメトリクスを持つデバイスを実現するために必要な大量のスループットにとって必要だった。特に、必要な最適化行列の行の計算につながる合計の１つの要素を並列に計算することは、要件を満たすのに十分だが、サブシステムが十分な回数複製された場合に限られる。効率性に対するこの差し迫った必要性は、単純な乗算の上に複雑な超越関数を評価するための、新しいアプローチを生み出した、古いアルゴリズムに対して異なる見方をもたらした。
【００２６】
このより広い世界では、この結果得られるアルゴリズムは、処理ハードウェアの要素が係数を無視できない複素数値のデータと効率的に、またはリアルタイムでやり取りしなければならない場合はどこでも適用可能である。これには、ＧＰＳ、ワイヤレス、ラジオおよび衛星送受信システムを含む通信、オーディオ処理、空中触覚を含む位相アレイシステム、飛行時間イメージングおよびデジタル干渉法、散乱効果の計算が含まれるが、これらに限定されない。これは、本発明にとってより広い範囲の機会につながる可能性がある。
【００２７】
ＣＯＲＤＩＣとしても知られているＶｏｌｄｅｒのアルゴリズム（Ｖｏｌｄｅｒ、Ｊ．−ＴｈｅＣＯＲＤＩＣ計算手法−１９５９）は、単位複素数値を乗算することができるが、基本ＣＯＲＤＩＣアルゴリズムの計算を超えないと、計算の第２オペランドに係数を含めることができない。
【００２８】
この理由で数年前に作者が実装した、１９９４年に発明された別のアルゴリズムは、作者Ｂａｊａｒｄ、ＫｌａおよびＭｕｌｌｅｒにちなんでＢＫＭとして知られており、さらに完全に複素数−複素数乗算の目標を達成することができる（Ｂａｊａｒｄ、Ｊ、Ｋｌａ、S．およびＭｕｌｌｅｒ、Ｊ.−ＢＫＭ：複雑な初等関数のための新しいハードウェアアルゴリズム−１９９４）。しかし、このアルゴリズムは複雑な実装を有するため、そうでない場合ほど多くの注目を集めていない。このアルゴリズムでは、複素数値のオペランドを複素数の対数に減らし、それらを加算してからべき乗化することで、完全な複素数−複素数乗算を実行できる。
【００２９】
ＢＫＭと、場合によっては新しいアルゴリズムを使用して、シリアルとパラレルの両方で、この複素数値の乗算と複素数値の除算を実現するためのより効率的な方法がある。
【００３０】
１９９９年のＢＫＭアルゴリズムの補遺は、複雑な対数の過度に複雑な実装を扱い、実装をいくらか単純化した。（Ｂａｊａｒｄ、Ｊ.およびＩｍｂｅｒｔ、Ｌ.−複雑な初期機能の評価：ＢＫＭの新バージョン−１９９９年）。１９９４年および１９９９年の補遺のこのＢＫＭ方式は、この文書では単にＢＫＭと呼ばれるアルゴリズムを形成する。
【００３１】
実際のハードウェアに組み込むためにアルゴリズムを最適化しようとしている間に、説明したように数学関数を評価することよりも、我々の、またはより一般的には既存のシステムへの統合の容易さを優先する新しい技法のアイデアが明らかになった。この新しい技術は、それが概念実証コードに開発されたときに、元のアルゴリズムよりもハードウェア実装の一般的な採用にとってより有望なものとなる多くの新たな特性を持つことが判明した。加えて、元の技術に対して要求された正確な数値が要求される場合、これらの結果を検索するためにさらなるスケーリングが適用されてもよく、アルゴリズムは全体として競争力を維持したままである。どちらの場合でも、結果は、最先端技術を超える実質的な利点を備えた技術である。
【００３２】
単位ノルムを有する二次問題は、その固有問題が二次問題と等価であるマトリクスを構築することによって解決され得ることが知られている。ランク１より大きい二次問題マトリクスは、多くのランク１問題マトリクスを合計すること、および／または既存の問題マトリクスのヌル空間を拡張することによって生成され得る。このマトリクスのクリロフ部分空間を生成することによって、暗示基底関数の最適化された部分集合を生成することができる。クリロフ部分空間を線形システムへの基底セットとして使用することによって、多数のテスト関数を使用して基底セットの正しいまたは近似評価を達成することができる。
【００３３】
いくつの固有ベクトルが使用されるか、そしてそれらが選択される場所から、それに続く各固有ベクトルのパワーが減少するため、多くの基底関数の複雑なセットのグレースフルデグラデーションを達成することができる。
【００３４】
さらに、音響場形状の結果を定義するための複素数値線形システムは、
2020502906.tif 000003
個の基底関数を取り、これらをキャリア周波数における音圧の所望の構成を記述する結果の場に再結合する。
【００３５】
この時点までは、コントロールポイントを作動するために、コントロールポイント基底関数を使用したことがある。本明細書に記載しているのは、初期コントロールポイントまたは領域基底関数に異なる一組のコントロールテスト関数を適用するために一般化する方法である。
【００３６】
これは、音響場における最初の「パレット」の振る舞いと、線形システムが使用する目的関数との間の関連性を断ち切るため、重要である。例えば、これを使用してベッセル関数の基底セットを作成し、それを使用してコントロールポイントを作動することができる。
【図面の簡単な説明】
【００３７】
図面の簡単な説明
添付の図面は、以下の詳細な説明と共に、別々の図を通して同一または機能的に類似の要素を指すが、明細書に組み込まれてその一部を形成し、特許請求の範囲に記載の発明を含む概念の実施形態をさらに例示するのに役立つ。これらの実施形態の様々な原理および利点を説明する。
【００３８】
【図１】図１は、コントロールしたダイナミックレンジの新しい波形を生成するために非線形再マッピング機能を通して、処理することができる任意のダイナミックレンジの触覚波形である。
【図２】図２は、コントロールしたダイナミックレンジの新しい波形を生成するために非線形再マッピング機能を通して、処理することができる任意のダイナミックレンジの触覚波形である。
【００３９】
【図３】図３は、異なる範囲のコントロール技術を適用した異なる構成の、コントロールポイントから生じる有効振幅変調波形である。
【図４】図４は、異なる範囲のコントロール技術を適用した異なる構成の、コントロールポイントから生じる有効振幅変調波形である。
【図５】図５は、異なる範囲のコントロール技術を適用した異なる構成の、コントロールポイントから生じる有効振幅変調波形である。
【図６】図６は、異なる範囲のコントロール技術を適用した異なる構成の、コントロールポイントから生じる有効振幅変調波形である。
【図７】図７は、異なる範囲のコントロール技術を適用した異なる構成の、コントロールポイントから生じる有効振幅変調波形である。
【図８】図８は、異なる範囲のコントロール技術を適用した異なる構成の、コントロールポイントから生じる有効振幅変調波形である。
【図９】図９は、異なる範囲のコントロール技術を適用した異なる構成の、コントロールポイントから生じる有効振幅変調波形である。
【図１０】図１０は、異なる範囲のコントロール技術を適用した異なる構成の、コントロールポイントから生じる有効振幅変調波形である。
【図１１】図１１は、異なる範囲のコントロール技術を適用した異なる構成の、コントロールポイントから生じる有効振幅変調波形である。
【図１２】図１２は、異なる範囲のコントロール技術を適用した異なる構成の、コントロールポイントから生じる有効振幅変調波形である。
【００４０】
【図１３】図１３は、ゴルコフ関連関数の最適化から生じる音圧の投影図である。
【図１４】図１４は、ゴルコフ関連関数の最適化から生じる音圧の投影図である。
【図１５】図１５は、ゴルコフ関連関数の最適化から生じる音圧の投影図である。
【図１６】図１６は、ゴルコフ関連関数の最適化から生じる音圧の投影図である。
【００４１】
【図１７】図１７は、同じ音圧レベルを有する同じ波のグラフであるが、ベクトル成分の効果が最大の触覚効果を引き起こし得るが、触覚効果を生じないシナリオである。
【００４２】
【図１８】図１８は、ｄＢ ＳＰＬ、ｄＢ ＳＩＬおよびｄＢ ＳＶＬの位相アレイによって生成された、音響場の測定値のシミュレーションである。
【００４３】
【図１９】図１９は、二次問題を解くことによって作成された、コントロールポイントプロットを示す。
【００４４】
【図２０】図２０は、平坦な直線形トランスデューサアレイに平行な近距離場領域の端部における音場を通るカットを示す。
【００４５】
【図２１】図２１は、ベッセル関数のグラフを示す。
【００４６】
【図２２】図２２は、高次ベッセルビームおよびマシュービームへの近似の断面図を示す。
【図２３】図２３は、高次ベッセルビームおよびマシュービームへの近似の断面図を示す。
【００４７】
【図２４】図２４は、曲げビーム、ボトルビーム、およびボトルビームの部分断面を示す。
【図２５】図２５は、曲げビーム、ボトルビーム、およびボトルビームの部分断面を示す。
【図２６】図２６は、曲げビーム、ボトルビーム、およびボトルビームの部分断面を示す。
【００４８】
【図２７】図２７は、波長が長くなったビーム、スネークビーム、およびヘリカルビームを示す。
【図２８】図２８は、波長が長くなったビーム、スネークビーム、およびヘリカルビームを示す。
【図２９】図２９は、波長が長くなったビーム、スネークビーム、およびヘリカルビームを示す。
【００４９】
【図３０】図３０は、所望波とは反対の方向から波に寄与するトランスデューサを除去する効果を示す。
【図３１】３１は、所望波とは反対の方向から波に寄与するトランスデューサを除去する効果を示す。
【００５０】
【図３２】図３２は、再マッピングなしの対数法とべき乗法の収束領域を示す。
【図３３】図３３は、再マッピングなしの対数法とべき乗法の収束領域を示す。
【図３４】図３４は、再マッピングなしの対数法とべき乗法の収束領域を示す。
【００５１】
【図３５】図３５は、対数法とべき乗法の収束グラフを示す。
【００５２】
当業者であれば、図中の要素は簡単かつ明確にするために示されており、必ずしも一定の縮尺で描かれていないことを理解するだろう。例えば、図中の要素のいくつかの寸法は、本発明の実施形態の理解を向上させるのを助けるために、他の要素に対して誇張されている場合がある。 
【００５３】
装置および方法の構成要素は、必要に応じて図面中の従来の記号によって表しており、本発明の実施形態を理解することに関連する特定の詳細のみを示し、開示を曖昧にしないようにする。本明細書の記載の恩恵を受ける当業者は、当業者に自明である。
【発明を実施するための形態】
【００５４】
詳細な説明
前述の技術、装置および方法は、組み合わせるか、または別々に空中触覚システムで使用できる。
【００５５】
１．空中ハプティックデバイスのための動的範囲縮小
【００５６】
Ｉ.心理物理学的モデリング
【００５７】
システムの振る舞いと触覚出力は、人間が接触を知覚する方法に合わせて、調整しなければならない。例えば、コントロールポイントのサイズが大きくなると、コントロールポイントで作成する物理的な振幅は変化しないが、触覚の知覚も増加する。音声などによるクロスモーダル効果は、それらが相関する場合、触覚フィードバックを向上させるのに役立つ。これは、物理的に定義されたコントロールポイントの間の食い違いを招き、同等の認識を達成するために心理物理学的モデリングステップの使用を必要とする。
【００５８】
ＩＩ．心理ハプティックモデルを使った再マッピング
【００５９】
再生装置の物理的制限のために、正確な波形は、指定したように再現可能ではないかもしれない。波形は、心理的に定義された知覚空間内に構築され得る。あるいは、それは、実質的に絶対的な物理的スケールでの強度および力の測定値から構成されてもよく、それは触覚テクスチャ、またはさらなるメタデータと共に潜在的に表現される任意の構築波形の記録であり得る。この波形は、「完璧な」空中の触覚デバイスが与えられたときの、コントロールポイントの動作を表す。次に、これは、コントロールポイントの再現の知覚的側面を記述および定量化できる、心理物理学的モデルを使用して知覚的に評価される。この触覚効果を再現しながら、コントロールポイントが動いている可能性がある。運動は、触覚効果の知覚的評価の一部として含まれる。次に、これは、ユーザ位置、相対的幾何学的形状、および利用可能なパワー制限によって制約される空中触覚システムのインスタンスに利用可能な出力の空間に再解釈することができ、システムは実質的に等価なコントロールポイントを作成するように構成する。
【００６０】
心理触覚モデルは、完全な装置と実質的に同じ触覚応答を引き起こすために、コントロールポイントの変調波形をどのように修正できるかを決定する能力を、デバイスに与える。この再マッピングを達成するプロセスはまた、結果として生じるハプティック波形が雑音を多くすることを防ぐため、または知覚ハプティック空間に関しては近いがノイズが少ない、より可聴的に許容可能な変調を選択するための心理音響モデリングステップを含むことができる。可聴ノイズの量を増やすことによって、多くの波形をより触覚的に強力にすることができることが知られている。いつかの状況では、他ではこれが望ましいかもしれないが、これは耐えられないほど十分に騒々しいかもしれない。次に、音響心理学的モデルを併用して、音声をハプティックトレードオフ設定にして、出力をバイアスしてより静かまたはより強くすることができる。
【００６１】
ＩＩＩ．精神的触覚ユニット
【００６２】
振幅が同等であると感じるように意図する心理空間で波形を定義する場合、空中で生成されたコントロールポイントの心理効果を考慮に入れるユニットを設計しなければならない。これらの中で最も注目に値するのは、手に触覚が統合されていることである。各コントロールポイントに含まれる焦点領域は、物理的装置の特性に応じてサイズおよび形状を再び変化させ、そしてこれは、精神触覚ユニットの構成において相殺しなければならない。そのような精神触覚ユニットは、焦点領域のサイズおよび形状に応じて評価する物理的コントロールポイント振幅を変更し、広い面積にわたるタッチ感覚の統合の影響を相殺する。このアプローチは、同じ数の精神触覚ユニットが同等の感覚をコード化することを可能にし、デバイス間の触覚的同等性および物理的インスタンス化を維持する。
【００６３】
ＩＶ．動的な波形のマッピング
【００６４】
ターゲットデバイスにベストなマッピングを到達させるために、 完全な波形とその軌跡が考慮されるかもしれない。しかし、対話性、したがって短い待ち時間が必要条件であるため、波形の一部または個々のサンプルだけが利用可能である。波形全体をターゲットデバイスの機能に再マッピングすることができるが、これを行うには、波形を事前に知っておく必要があるため、周知の波形の前処理ステップと見なすことができる。代替方法は、振幅波形サンプルの物理的特性を振幅で再マッピングすることである。なぜなら、振幅サンプルの再マッピングは対話型波形を用いてリアルタイムで達成することができ、リアルタイムで触覚フィードバックを提供するのにより適しているからである。
【００６５】
多様な関数が、コントロールポイント振幅サンプルのための再マッピングを提供するために、良い候補であると考えられる。これらは、対数線形曲線と共にシグモイド関数を含み、それは他のモダリティにおいて低いダイナミックレンジを有する装置上で、知覚的な感覚入力を納得のいくように再現することができる。これらの曲線は一般にマッピング関数における上限の指定を必要としないが、図１および図２に示すように、利用可能な有限量のパワーがあることを考えると、ダイナミックレンジはより効果的に使用することができる。図１は、任意のダイナミックレンジ１００のハプティック波形を示しており、これは非線形再マッピング機能１１０を介して処理され、コントロールしたダイナミックレンジ１２０の新しい波形を生成することができる。これは任意のダイナミックレンジのハプティック波形を取得し、それをより小さなレンジに変換するために使用できる。波形の振幅は無限大である可能性があるため、これは出力のダイナミックレンジ全体を埋めるわけではなく、未使用の出力範囲になる。図２では、図１と同様に、非線形関数２１０は、振幅波形２００をターゲット装置２２０に利用可能な振幅に再マッピングする。波の広がりを知ることで、再マッピング機能をデバイスに固有の出力範囲に、正確に合わせることができる。これは、出力デバイスの全てのレンジの開発を可能にする。
【００６６】
説明のために、以下の例では、空中触覚デバイスにおけるマッピングがいかに望ましくない挙動を防ぐのかを説明するために、ミニマリスト線形スケーリングを使用する。
【００６７】
図３に示すグラフ３００は、波形の再現可能な知覚への再マッピングを達成するための第１のステップが、波形およびその限界３００を記述するソルバー関数にデータを供給することである。図３では、実線は
2020502906.tif 000004
であり、点線は、
2020502906.tif 000005
である。それは、関数を再マッピングする単純な波形の入力である。ここでｙ軸は、デバイスの内部単位である。
【００６８】
これは、クリッピングなしで、波形を放射するために再生する必要がある最大振幅を定義することによって達成され、これは、各サンプル時間で、
2020502906.tif 000006
と表される。これは、入力波形の振幅
2020502906.tif 000007
と共に、入力としてソルバーに入力される。
【００６９】
図４に示されているのは、入力データにいかなる変更も適用せずにこの波形を再生しようと試みた結果を示すグラフ４００である。その結果、コントロールポイントで必要な振幅を再現するために利用できるパワーが不足しているために波形がクリップされる。図４では、直線は
2020502906.tif 000008
であり、それは振幅
2020502906.tif 000009
の結果であり、これは、デバイスのダイナミックレンジに関係なく、図１のユーザ指定振幅を直接出力しようとした結果である。
【００７０】
システム内のアクティブトランスデューサの物理的構成を、
2020502906.tif 000010
位置
2020502906.tif 000011
における所望のコントロールポイント振動を、
2020502906.tif 000012
として定義すると、空中触覚システムに対するソルバー関数の補助機能を使用して、関数振動
2020502906.tif 000013
を定義することができ、与えられた装置が与えられた時間に、これを物理的に提供することができる場合には、コントロールポイントにおける所望の搬送波振幅、あるいはコントロールポイントの位置において再生できる最も近い振幅のいずれかをもたらす。
【００７１】
2020502906.tif 000014
を評価して、この時点で入力振幅をマッピングできる範囲の最大値を求める。コントロールポイントはサンプル間であまり移動しないので、これは正確である必要はなく、コントロールポイントの前の反復からの
2020502906.tif 000015
の値が使用できる。丸め誤差および他の雑音源を除去するために、この値をフィルタリングすることもできる。この件では、再マッピング関数は波形の単一サンプルを取り、これを区間
2020502906.tif 000016
から区間
2020502906.tif 000017
に再マッピングする。図５に示すグラフ５００は、関数５００の再マッピングを通して、
2020502906.tif 000018
から
2020502906.tif 000019
へ直接マップを産出する。振幅を
2020502906.tif 000020
から
2020502906.tif 000021
から、値
2020502906.tif 000022
で得られる範囲に再マッピングすることで、波形のクリッピングを防止すると、
2020502906.tif 000023
となる。
【００７２】
Ｖ．複数の同時コントロールポイントの再マッピング
【００７３】
図６の２つの正弦波の場合の例のように、ソルバーが波形を正しく表すことができるという保証を提供することができるように、一度に複数のコントロールポイントに対してシステムを使用するには、範囲を同時に評価しなければならない。グラフ６００が示すように、２つのコントロールポイントは、同時に異なるポイントで、異なる振幅
2020502906.tif 000024
および
2020502906.tif 000025
で生成される。一方はトランスデューサアレイに対して接近して移動し、他方はトランスデューサアレイに対して横方向に移動する。それらは、異なる周波数の２つの正弦波で、振幅が振動するように設定されている。
【００７４】
デバイスのパワー制限に違反すると、パワーが不足してクリップされた波形になり、デバイスの制限を超える複数のコントロールポイントがポイント間のクロストークを促進し、それらの触覚効果を混乱させる単一のコントロールポイントとは異なる。図７に示すグラフ７００は、２つの動くコントロールポイントの例を示す。図７は、トランスデューサアレイを使用して２つのコントロールポイントを再作成しようとする試みを示しており、２つのポイント間で異なる共振が発生して置き換えられると、複雑で予測不能な動作になる。
【００７５】
図８に示すのは、両方の点に存在するパワー制限を決定するために補助ソルバー評価を使用するグラフ８００であり、示された波形を導き、それは次のように記述することができる。
2020502906.tif 000026
ここで、一方のコントロールポイントがより強くなる可能性を持ってデバイスに向かって移動し、他方のコントロールポイントはデバイスの表面を横切って弱くなる可能性を持って離れていく。図８からわかるように、ソルバーソリューションによって生成するダイナミックレンジの制限により、それらの波形は意図したとおりにクロストークなしで表される。したがって、図８は、振幅エンベロープが、エラーなしで、またはコントロールされた周知のエラーレベルで再現されるようにソルバーによって保証されるレベルまで低減された２つのポイントを示している。これは理想的な振幅レベルで振幅誤差の２乗を最小にする構成であるため、それらは同じパワーである。
【００７６】
ＶＩ．複数の同時コントロールポイントにわたるパワー再分配
【００７７】
振幅範囲は、異なる値に設定でき、パワーレベルを管理するために時間を通して動的に設定できる。この方法では、範囲を操作すると、特定のコントロールポイントに集中している配列出力の比率が決定できる。これは、図９および図１０では、あるコントロールポイントから別のコントロールポイントへのパワーの単純な操作で示す。図１１および図１２で示されるような、コントロールポイント間の動作ダイナミクスを可能にするより複雑な位相外のずれ関係もある。
【００７８】
図９に示しているのは、２つの入力コントロールポイントに対する入力波形および範囲を有するグラフ９００である。これらは波形範囲が異なり、各ポイントがデバイスから利用可能なパワーをより多く、またはより少なくできるようにするためのヒントとしても使用される。図１０に示すのは、図９の入力から得られた出力を有するグラフ１０００である。この構成では、アレイによって生成したパワーは、それらがそれぞれアレイから遠ざかっている間およびアレイに向かって動いている間に、一方のコントロールポイントから他方のコントロールポイントに移動する。これにより、図６、図７、図８で予想したように、コントロールポイントの振幅が等しいと予想される場合よりも多くのパワーを使用できる。
【００７９】
図１１に示すのは、交互の関数を有するグラフ１１００であり、１つのコントロールポイントまたはコントロールポイント群は、範囲が正弦波に一致することによって効果的にモデル化できる。この結果、各波形は変調システムの各サイドのパワー要件に従って再生し、システムの各サイドはフルアレイパワーにアクセスできることが保証される。その結果、任意の出力クリッピングが発生する可能性なしに、タイムスライスが行われる。図１２に示すのは、図１１の入力から得られた出力を有するが、説明されているグループ分けなしのグラフ１２００である。これにより、パワーバジェットがすべてのコントロールポイント間で不必要に共有しているため、出力ダイナミックレンジが減少する。
【００８０】
２．音響浮上のためのトラップポイントの効率的な生成
【００８１】
Ｉ．浮遊球状粒子のための最適な目的
【００８２】
粒子を浮揚するために必要とする最適化問題は、ゴルコフエネルギーポテンシャルの使用を中心とする。これは、音響場内の粒子の可能な各位置のエネルギーを評価する、スカラー関数である。空間位置で最小または低エネルギー状態を見つけることは、潜在的に安定なトラップポイントを見つけることに対応する。しかし、ゴルコフポテンシャルが考慮に入れていない音響効果、例えば音響波が流体中でバルク運動を引き起こす音響流がある。この理由から、エネルギーポテンシャルの実行可能な最小値であると決定されることが多い変動圧力の領域で、トラップが発生しないことを確実にすることも好ましく、可能性の定義に含まれる。
【００８３】
重力および音響流効果のような外力の方向および大きさは、ゴルコフポテンシャルのみからは導き出すことができないので、ゴルコフスカラーポテンシャルのラプラシアンを計算することがしばしば好ましい。これにより、内側に押す力のより安定した配置を生成するために、どの軸を他の軸より優先するかを、さらにコントロールできる。
【００８４】
ＩＩ．ゴルコフスポテンシャルの評価
【００８５】
極小エネルギーの存在により、球状粒子が浮遊することを可能にするエネルギー場であるゴルコフポテンシャル分野が、（Ｂｒｕｕｓ ２０１２年）として記載される：
2020502906.tif 000027
ここで、
2020502906.tif 000028
は、音響媒体を通した音のスピードであり、
2020502906.tif 000029
は、球状粒子を通した音速であり、
2020502906.tif 000030
は、音響媒体の密度であり、
2020502906.tif 000031
は、球状粒子の密度であって、時間的平均圧力
2020502906.tif 000032
を含む左辺は潜在的なエネルギーの表示であり、粒子速度ベクトル 
2020502906.tif 000033
を含む右辺は潜在的なエネルギーに保守的に関係する。ゴルコフポテンシャルは、活動的または遠距離場に存在するという仮定をするため、これらの圧力と速度の項は同相の時間調和項と交換することができる。これは、時間調和ゴルコフポテンシャル分野が、次のように書き換えられることを示唆する：
2020502906.tif 000034
【００８６】
これを多くの情報源に拡張すると、次のようになる：
2020502906.tif 000035
【００８７】
これに関する問題は、交差項が事実上、これが次のようになることを意味する：
2020502906.tif 000036
その結果、二次であり、二乗したトランスデューサの数が増えることとなる。
【００８８】
ＩＩＩ．固有問題によるポテンシャルの最適化
【００８９】
時間高調波効果を無視し、各トランスデューサに入力する圧力および粒子速度を考慮した潜在的な定義は、次のように書くことができる：
2020502906.tif 000037
2020502906.tif 000038
2020502906.tif 000039
ここで、
2020502906.tif 000040
は所望のトラップのインデックスであり、余分な
2020502906.tif 000041
は、トランデューサ
2020502906.tif 000042
および
2020502906.tif 000043
からトラップポイント
2020502906.tif 000044
へ移動する際に、位相および振幅が変化する位相および振幅である。トランスデューサが非ゼロレベルにパワーを供給することを確実にするために、制約
2020502906.tif 000045
を使用して、ポテンシャルウェルを達成するための最適複素トランスデューサ活性化係数は、最大固有値を有する
2020502906.tif 000046
個の固有ベクトルによって与えられる。
【００９０】
2020502906.tif 000047
をこのように選択すると、目的関数
2020502906.tif 000048
が最大化され、これは、
2020502906.tif 000049
の定義の負の符号のため、ゴルゴフポテンシャル分野で、最小値が得られる。なぜなら、
2020502906.tif 000050
が固有ベクトルの定義であり、
2020502906.tif 000051
であることが定義されており、定義により、
2020502906.tif 000052
がマトリクスの固有値であり、
2020502906.tif 000053
なので、最大化され、これは、存在するのであれば、これは最大値をもたらさなければならない。問題の標準的な記述は、以下のようになる：
2020502906.tif 000054
【００９１】
しかし、以前に文献で議論したように、ゴルコフ分野のエネルギーを最小にすることは焦点を生成し、実際にはトラップポイントではない。トラップポイントを生成することは、トラップでの圧力変動を最小化することである。圧力変動の定義は複素数値であり、意味のあるように最小化することはできないが、代わりにゼロにする変動量を表すため、これは自明ではない。
【００９２】
ＩＶ．  固有値問題内に、ヌル空間を構築して圧力ゼロの生成
【００９３】
ゼロ圧力条件は、音響場内の所与の点で、ゼロ圧力に評価される一連の制約として構成することができる。高次元位相空間内に、周知の最適圧力集束解に平行な解ベクトルの一部が存在する場合、圧力は所与の位置でゼロ以外の値に評価することしかできない。トラップポイント
2020502906.tif 000055
では、ベクトル
2020502906.tif 000056
があり、それは以下のようになる：
2020502906.tif 000057
【００９４】
それゆえに、トラップポイント
2020502906.tif 000058
におけるゼロ圧力に対する制約は、制約式として記述できる：
2020502906.tif 000059
【００９５】
振幅が、これに対するものであるかどうかは重要ではなく、制約は単位長さである必要があるので、それゆえ以下のようになる：
2020502906.tif 000060
【００９６】
問題の標準的な記載は、以下のとおりである：
2020502906.tif 000061
【００９７】
固有値系をＧｏｌｕｂ（１９７３）から導き出し、それを複素数値について言い換えると、以下の議論をもたらす：
【００９８】
2020502906.tif 000062
に伴う上記のシステムを挙げると、システムの特性方程式は次のように書くことができる:
2020502906.tif 000063
ここで、
2020502906.tif 000064
および
2020502906.tif 000065
は、ラグランジュ乗数である。これは、
2020502906.tif 000066
ベクトルの共役によって微分することができ、システムが得られる：
2020502906.tif 000067
【００９９】
2020502906.tif 000068
を乗ずると、以下が得られる：
2020502906.tif 000069
【０１００】
2020502906.tif 000070
が単位長であると仮定すると、
2020502906.tif 000071
となり、先の問題定義では
2020502906.tif 000072
となり、両方を代入すると以下のようになる：
2020502906.tif 000073
【０１０１】
これを以前の導関数に代入して置換すると以下が得られる：
2020502906.tif 000074
【０１０２】
これは以下のとおりである:
2020502906.tif 000075
【０１０３】
これは固有システムの記述であり、未だ
2020502906.tif 000076
は対称ではないため、固有ベクトルの計算は困難である。しかしながら、
2020502906.tif 000077
は、べき等
2020502906.tif 000078
に示され、固有値は、順序付けによって変わらないので以下が得られる:
2020502906.tif 000079
【０１０４】
これを、より多くの制約ベクトルに拡張しても、左辺をさらに乗算しても問題はない：
2020502906.tif 000080
多くの乗算の後、最終的に以下が得られる：
2020502906.tif 000081
ここで
2020502906.tif 000082
は、ムーア−ペンローズ擬似逆行列を表す。
【０１０５】
要約すると、初期問題の解は、
2020502906.tif 000083
の最大固有値を持つ固有ベクトルを求めることによって見つけることができる。
2020502906.tif 000084
の乗算は、解決マトリクスのヌル空間に、基底ベクトルを追加するランク１の変更を追加すると見なすことができる。
【０１０６】
最終的に、簡単な方法で、ゴルコフポテンシャルのラプラシアンを計算するための、サポートを追加することは有用である。ここで見いだされたゴルコフポテンシャル解はトラップポイントを生み出すが、どの方向から力が加えられるのかをコントロールすることはできず、力がない方向でさえあるかもしれない。これを達成するためには、潜在分野の勾配における最も発散する点に関して最適化し、次いでこの場の勾配を再重み付けする必要がある。これは、一般的に言って、ゴルコフ場のラプラシアンである。
【０１０７】
Ｖ.  ゴルコフ場のラプラシアンを生成するための修正
【０１０８】
ゴルコフ場のラプラシアンを最小化するための修正は、
2020502906.tif 000085
マトリクスのみに適応する必要がある。各トランスデューサからの波面曲率の影響を破棄し、それらを平面波のみと見なすと、定式化を複雑にする特異点が取り除かれ、粒子速度は一定になる。各トランスデューサからの平面波の結果は、以下の式に書くことができる：
2020502906.tif 000086
【０１０９】
マトリクスは、ラプラシアンケース以外では、同等のシステムとして表わされる：
2020502906.tif 000087
2020502906.tif 000088
の全ての制約は、
2020502906.tif 000089
に移行される。
この定式化のラプラシアンを取ると、修正マトリクス定式化が得られる：
2020502906.tif 000090
【０１１０】
2020502906.tif 000091
マトリクスの対応する要素の代わりにこの値を使用すると、ゴルコフポテンシャルのラプラシアンを最適化するという点で、解決策が得られる。式
2020502906.tif 000092
は、以下のように記述することができる：
2020502906.tif 000093
【０１１１】
これは、また、
2020502906.tif 000094
方向からの寄与に分割することもできる：
2020502906.tif 000095
【０１１２】
この方法では、重み付けを適用して拘束軸を生成し、各方向から生成された力にバイアスをかけることができる。
【０１１３】
ＶＩ．デバイスでの実装
【０１１４】
本明細書におけるシステムは、単純な固有システムソルバーに基づいて、１つまたは複数の浮揚トラップを計算するために使用することができる装置を表す。これは、アイテムを拘束し、トラップし、輸送するために音響浮揚を使用することができる実施形態を作成するために、デバイスファームウェアで潜在的に実施することができる。
【０１１５】
これが、図１３〜１６に示される。
【０１１６】
図１３は、同方向の重み付けを伴う、ゴルコフのラプラシアンを示す。左側面１３１０は、トランスデューサアレイに平行なｘ−ｙ平面である。右側面１３２０は、平面トランスデューサアレイに垂直であり、且つそれを通ってスライスするｘ−ｚ平面である。
【０１１７】
図１４は、重く重み付けされたｚ方向を有するゴルコフのラプラシアンを示す。右側面１４２０は、平面トランスデューサアレイに垂直であり、且つそれを通ってスライスするｘ−ｚ平面である。
【０１１８】
図１５は、ゼロ圧力状態のゴルコフ最大化を示す。左側面１５１０は、トランスデューサアレイに平行なｘ−ｙ平面であり、それは、一様な強度の輪の周りの不均一な陰影によって示される、らせん状のトラップを記載する。右側面１５２０は、その下の平面トランスデューサアレイに垂直であり、且つそれを通ってスライスするｘ-ｚ平面である。
【０１１９】
図１６において、さらなる固有モード１６１０、１６２０、１６３０は、やはりコア問題を解決するトラップ解に対応する。これらは、目的関数を最大化しないため、あまり効果的ではない。これは、目的関数の値である修正ゴルコフポテンシャルが、主固有ベクトルに直交する、これらのモードでは小さいからである。達成した目的関数の値は、各固有ベクトルに関連する固有値に反映される。小さい固有値に関連する固有ベクトルは、これらの代替的で効果の低い解を表す。
【０１２０】
３．空中触覚に適用したベクトル効果に焦点を当てた位相アレイ装置
【０１２１】
Ｉ．活性音響強度
【０１２２】
活性音響強度
2020502906.tif 000096
は、音響媒体の圧力
2020502906.tif 000097
および粒子速度ベクトル
2020502906.tif 000098
を用いて定義することができる。音響波は、以下のように記述できる：
2020502906.tif 000099
ここで、活性強度は、純粋に実数のベクトルである。波源の近距離場では、
2020502906.tif 000100
ではなく、反応場が支配的である。しかし、遠距離場では、１または２つの波長
2020502906.tif 000101
を超える。この理由から、
2020502906.tif 000102
のみが考慮され、現実感が保証される。
【０１２３】
下記の単調平面波源について考察する：
2020502906.tif 000103
その結果下記が得られる：
2020502906.tif 000104
【０１２４】
ここで、
2020502906.tif 000105
は圧力振幅であり、
2020502906.tif 000106
は、各方向の粒子速度振幅である。
【０１２５】
所与の軸に沿って波によって伝達される音の強度、したがってエネルギー束は重要であるので、スカラー寄与を測定するための内積は、以下のようになる：
2020502906.tif 000107
【０１２６】
音波の方程式の形に起因して以下のようになる：
2020502906.tif 000108
重ね合わせの原理は、圧力ベクトルと粒子速度ベクトルの両方に当てはまる。
2020502906.tif 000109
個の音源を想定し、それぞれを対象の点における平面波として近似すると、次式が得られる：
2020502906.tif 000110
効率のためには、線形システムを用いて適切な変数について解くことが望ましい。これは、２つの独立した変数で構成しているため、与えられた方向の音響強度について解くことができないことを意味する。既存の技術は複素数値のスカラー圧力
2020502906.tif 000111
解き、それから個々の波のそれぞれの方向は重要ではないか、または、ほぼ平行であると仮定する。代わりに、これは運動量（および運動エネルギー）に直接影響するので、粒子速度ベクトル
2020502906.tif 000112
を解き、これは、空中での触覚フィードバックの用途にとって重要である。
【０１２７】
このことは、波の指向性とその結果生じるエネルギー伝達を反映する解を生成し、複素数値の粒子速度ベクトル
2020502906.tif 000113
を利用するために現在の、解のメカニズムを拡張する２つの潜在的な技術につながる。
【０１２８】
ＩＩ．複素数値に対する解法
【０１２９】
複素数値に対する解法は、値に対する音響モデルを使用することを含む。このモデルは、コントロールポイントの位置と並んでトランスデューサの位置が与えられて、フィールド内のその点で得られたフェーザが決定されると、照会することができる。これは複素数値であり、
2020502906.tif 000114
ここで
2020502906.tif 000115
は、音響場内のサンプルポイントのインデックス、
2020502906.tif 000116
は、作動トランスデューサのインデックス番号に対応する:
2020502906.tif 000117
ここで、
2020502906.tif 000118
は、
2020502906.tif 000119
の場所で空中に意図したフェーザであり、
2020502906.tif 000120
は、トランスデューサ
2020502906.tif 000121
からのトランスデューサ活性化係数、空気中での挙動を生成するのに必要な初期位相および振幅が
2020502906.tif 000122
で特定される。
【０１３０】
しかしながら、
2020502906.tif 000123
は複素数値であり、位相と振幅仕様の両方を含むため、位相オフセットは自由に修正できる。これを利用して、振幅の最大の可能な部分空間を開く位相値を見つけることができる。これは、システムの主固有ベクトルｙの発見によって達成される。
2020502906.tif 000124
【０１３１】
2020502906.tif 000125
は、音響場
2020502906.tif 000126
内の各コントロールポイントについて、音響場の各ポイントで所望される振幅を表し、
2020502906.tif 000127
は、内積をマトリックスの各対角内の１にキャンセルするために適用される列当たりの規格化値である。この行列の固有ベクトル
2020502906.tif 000128
の正規化された成分は、初期振幅
2020502906.tif 000129
と乗算されたときに、場内の圧力点の各々の間で強化を促進するフェーザベクトル
2020502906.tif 000130
に対応する。
【０１３２】
未決定問題の最小ノルム解は、次の形を用いて解くことができる：
2020502906.tif 000131
【０１３３】
この解決に重要なのは、マトリクス
2020502906.tif 000132
の逆行列であり、これはコレスキー分解によってシステムを解く最初のステップであり、次のように記述できる：
2020502906.tif 000133
これが固有システムのマトリクス
2020502906.tif 000134
の導出にも役立つことを意味する。しかし、次の形式でシステムを操作することもできる：
2020502906.tif 000135
これは非常に小さな線形システムとそれに続く解のサイズをもたらし、これは任意の手法（ＱＲ分解、コレスキー分解など）で解くことができ、放出されるまでこの小さな状態に保たれる。このＣマトリクスおよびこの
2020502906.tif 000136
ベクトルはどちらも、計算されるまで、トランスデューサにインデックスを持たない：
2020502906.tif 000137
【０１３４】
トランスデューサは多数ある可能性があるため、最後の可能性がある時点で
2020502906.tif 000138
を乗算して各トランスデューサから放射される波形を乗算し、
2020502906.tif 000139
トランスデューサ アクティベーションに変換し直すことがしばしば価値がある。
【０１３５】
ＩＩＩ．１つまたは複数の異方性コントロールポイントの解決
【０１３６】
それらの各々に沿って指定された法線ベクトルまたは方向を有するように音響場内にコントロールポイントを作成することは、運動量要件を位相アレイソルバーに組み込むための１つの方法である。各コントロールポイントに法線ベクトルを与えることは、各コントロールポイントにおける圧力について解くのではなく、対象となる点におけるｘ，ｙおよびｚ粒子速度、したがって波の運動量の成分について効果的に解くことを意味する。それを達成するために、線形システム方程式を前から展開して、各粒子速度について３つの成分を、前に定義した圧力については１つの代わりに含める。それらの３つ成分は、
2020502906.tif 000140
として記載され、以下のように定義される：
2020502906.tif 000141
各ケースで、
2020502906.tif 000142
は、コントロールポイント
2020502906.tif 000143
におけるトランデューサ
2020502906.tif 000144
からの伝搬方向を示すベクトルであり、その位置を横切って運動量が測定される
2020502906.tif 000145
は、コントロールポイント
2020502906.tif 000146
の法線ベクトル方向である。解決のための核となる線形システムは、以下のようになる：
2020502906.tif 000147
これは、各コントロールポイントの音響媒体の粒子の速度の３つの成分を示す。しかし、システムの解法に簡約表現アプローチを使用する場合、つまり必要な作業量を減らすために、Ｃマトリクスを指定する場合、このマトリクスのエントリ数は９倍に増加し、計算するのにより強大な問題となる。先にマトリクスＣの各要素
2020502906.tif 000148
をサブマトリクスとする：
2020502906.tif 000149
【０１３７】
これは、固有値マトリクス
2020502906.tif 000150
の既存のアプローチを使用した場合にもうまく動作し、同じアプローチを使用して構成することができるが、ユーザ入力の振幅は
2020502906.tif 000151
および
2020502906.tif 000152
となる。固有システムマトリクスのクロスタームは、別の軸に別の軸の波形を追加した場合の影響を、さらに別のコントロールポイントに考慮したときに、１つのコントロールポイントの位置で、１つの軸で測定を実行していると看做すことができる。
2020502906.tif 000153
はここで、成分を３倍多く含むが、
2020502906.tif 000154
による乗算は、装置上の空気中で所望の運動量ベクトルを再現するのに必要なトランスデューサ活性化係数にシステムを戻すことができる。
【０１３８】
この解決策は、空中触覚フィードバックの適用にとって特に興味深い。なぜなら、異なる角度から手に近づく波を作り出す可能性が開かれるからである。これにより、音響波が境界上の皮膚にせん断波を誘導することが可能になり、これは、余分な必要な計算を犠牲にして、空中触覚技術を使用して作り出すことができる全く新しい一連の感覚を開くものである。
【０１３９】
ＩＶ．音響媒体内の１つまたは複数のコントロールポイントでの粒子速度の解法
【０１４０】
圧力値を目標とするソルバーから運動量を目標とすることができるものへの変換を達成するための別の技術は、（元の圧力アプローチと同様に）音響粒子速度と完全に連動する線形システムソルバーを生成することであり、また、単調な調和時間のための複素数値である。さらに、これはスカラーである。これは、コントロールポイントにおける最大粒子速度の方向、したがって結果として生じる波の方向
2020502906.tif 000155
を計算し、次に乗算して各マトリクスエントリを単一の複素数値のスカラーにすることによって達成することができる。ドット積は、
2020502906.tif 000156
にわたって等しくスケールファクタによってすべてのサンプリング方向に影響を与えるので、これは初期マトリクス構成に対する後処理として適用することができる。しかし、結果として得られる波の方向は、ソルバーが機能する前に知る必要があるため、ソルバーの影響を受けない。空中触覚デバイスの構築に適用する場合、このアプローチは、空中であるが比較的近距離のフィードバックが最良であり、現在の解決メカニズムに対する拡張である小型デバイスに最も適しているであろう。
【０１４１】
これを達成するために、各トランスデューサに対する音響モデルは、音響媒体中の粒子速度を決定するための問題として、再キャストされる。このように、アルファベクトルは、与えられたコントロールポイントについて解く必要があるシステムとして、最大の運動方向に沿って測定される複素数値の粒子速度
2020502906.tif 000157
を生じるように再定義し、所与のコントロールポイントの解決を必要とするシステムは、式
2020502906.tif 000158
となる。しかし、不幸なことに、
2020502906.tif 000159
は、アプリオリに知ることができないが、線形システムマトリクスを構成するのに必要である：
2020502906.tif 000160
【０１４２】
したがって、音響媒体は、値を生成するために修正できる：
2020502906.tif 000161
ここで、下記式のとおりである。
2020502906.tif 000162
【０１４３】
これらの値は、ここで、
2020502906.tif 000163
マトリクスのドット積を構成するために、使用されなければならない：
2020502906.tif 000164
【０１４４】
しかしながら、それらは、これを達成するための２つの可能性のあるルートがある。その必要なドット積は、以下で定義される ：
2020502906.tif 000165
【０１４５】
これは、以下のように記述することができる：
2020502906.tif 000166
ここで、
2020502906.tif 000167
はダイアディックプロダクトである。これは、後で法線ベクトルを計算するために、マトリクス要素ごとに９つの別々の総和（対角線の場合は３つ）を維持する必要がある。これは、異方性コントロールポイント点アプローチのそれと多くの類似点を有する
2020502906.tif 000168
のための初期マトリクス構成をもたらすという点で有益であり、したがって、両方の技術をサポートする共通の実施インフラストラクチャを可能にする。しかし、多くの実装では、代わりに２段階のモデリングシステムを使用するという方法が、これら９つの別々の合計変数を保持するよりも望ましい場合がある。このようなシステムで、最初の段階では、速度振幅と波の方向のみを計算する単純化されたモデルを実行し、コントロールポイントの合計
2020502906.tif 000169
ごとに３つの速度方向を求める。位相計算から分離されているので、これは単純化されたモデルであるこれらの総和は、それぞれ最終の最強振幅の成分を与えるので、第２段階で総和を行う方向を計算するために使用することができる。
【０１４６】
要約すると以下のとおりである：
2020502906.tif 000170
ここで２段階では、モデルは以下のように計算できる：
2020502906.tif 000171
プロセスの他の段階を始める前に、トランスデューサーモデル内である。残りの段階は以前とほぼ同じように進行するが、速度が作用する方向が異なるだけである。
【０１４７】
このアプローチは、圧力ではなく、手への運動量伝達にさらに直接的に対処するための空中触覚フィードバックの応用を可能にする。これにより、感覚の再現性、ひいては知覚行動の測定が向上し、これは、実用化されたテクノロジとしての空中ハプティックスの拡大に不可欠である。
【０１４８】
Ｖ．添付の図面
【０１４９】
図１７は、同じ音圧レベル（ＳＰＬ）を有する同じ波が、入射角に応じて表面に異なる効果をもたらすことを示す。左側１７１０では、表面（右）と平行に進行する波（左）がいくらかのエネルギーＥで衝突し、運動量交換の結果としてわずかに少ないエネルギーで反射し離れる。右側１７２０では、同じ表面に対して垂直に進む波はエネルギーを反射および交換せず、したがって同じ量のエネルギーを残す。その結果、触覚効果を与えることは不可能である。
【０１５０】
図１８では、上の２行１８１０、１８２０に示されているのは、飛行時間法によって作られた２つの異なる焦点の、それぞれｄＢ ＳＰＬ、ｄＢ ＳＩＬおよびｄＢ ＳＶＬでの位相アレイによって発生された音場の測定のシミュレーションであり、ライトグレーは＋１０ｄＢ、白は０ｄＢ、ダークグレーは−１０ｄｂである。0同じカラースキームを有するそれぞれ下列１８３０（ｄＢ ＳＰＬ−ｄＢ ＳＩＬ）、（ｄＢ ＳＰＬ -ｄＢ ＳＶＬ）および（ｄＢ ＳＩＬ−ｄＢ ＳＶＬ）は、各測定方法間の差異を明らかにする。ＳＶＬ場測定シミュレーションは、媒体の粒子の速度に基づいているので、追加のベクトル成分の相殺により、波をさらに相殺することができる。その結果ｄＢ  ＳＶＬは常にｄＢ  ＳＩＬより低く、ｄＢ  ＳＩＬは、ｄＢ  ＳＰＬより低いが、それらは、それらのそれぞれの成分の乗算およびＳＩＬ測定値の有効幾何平均のために、対数ｄＢ空間において同様の距離に離れている。
【０１５１】
４．空中ハプティックシステムにおける二次問題のための独自システムソリューション
【０１５２】
Ｉ．二次最適化のためのレイリーリッツ定理
【０１５３】
マトリクス
2020502906.tif 000172
が正方の正の半定値エルミート
2020502906.tif 000173
マトリクスにスーパーポーズする。マトリクスは、正方のエルミートマトリクスであるため、各固有値は実数でなければならない。これは正の半定値であるので、これらの固有値は
2020502906.tif 000174
として正でなければならない。
【０１５４】
次に、固有値を以下のように順序付ける：
2020502906.tif 000175
【０１５５】
マトリクスの
2020502906.tif 000176
は、最適化問題の解であるものとして示される：
2020502906.tif 000177
および
2020502906.tif 000178
【０１５６】
これは、レイリーリッツの定理である。したがって、これらのそれぞれからの
2020502906.tif 000179
ベクトルは、検出された固有値に対応する固有ベクトルとして見出すことができる。
【０１５７】
ＩＩ．位相アレイシステムのための二次最適化
【０１５８】
位相アレイシステムを考慮すると、以下の
2020502906.tif 000180
として定義することは有用である：
2020502906.tif 000181
2020502906.tif 000182
個の単調音源を用いて、このベクトル解を用いて音源トランスデューサのアレイを、駆動することができる。
2020502906.tif 000183
が単位ノルム
2020502906.tif 000184
の固有ベクトルになるように、マトリクス
2020502906.tif 000185
を
2020502906.tif 000186
のエルミートマトリクスとして定義する:
2020502906.tif 000187
【０１５９】
これは、最適化の形成を示す。やらなければいけないことは、マトリクス
2020502906.tif 000188
に必要な特定の純粋な二次（マトリクス
2020502906.tif 000189
は対称またはエルミートマトリクスである必要があるため）目的関数を入力することである。
【０１６０】
エルミートマトリクス
2020502906.tif 000190
を生成する目的関数には、潜在的に有用な純粋に二次形式の５つの形式があるため、既存の固有ベクトルソルバーを使用して最大化し、空中超音波ハプティックを駆動するために使用できる解を得る：
【０１６１】
１．最大圧力最適化（全ＳＰＬ）：ａ．圧力の２乗に比例する、
2020502906.tif 000191
を最大化する；
ｂ.これは、位相アレイシステムを駆動するために一般的に使用される目的関数である。
【０１６２】
２．最大粒子高速最適化（全ＳＶＬ）は：ａ．粒子高速の２乗に比例する
2020502906.tif 000192
を最大化する；ｂ.これは、一貫した触覚知覚を得るために必要であるため、波の方向性を考慮に入れる。 
【０１６３】
３．しかし最大粒子高速最適化は、与えられた方向に沿う（使用可能なＳＶＬ）：ａ. 与えられた方向に沿う粒子高速の２乗に比例する
2020502906.tif 000193
を最大化する；ｂ．これは、波の方向性と（皮膚）表面
2020502906.tif 000194
に対する法線の影響の両方を考慮に入れる。
【０１６４】
４．最大音響強度（エネルギー束密度）最適化（全ＳＩＬ）：４．ａ.音響強度に比例する
2020502906.tif 000195
を最大化する：ｂ．これは、一貫した触覚知覚を得るために必要であるため、波の方向性を考慮に入れる。
【０１６５】
５．与えられた方向に沿った最大音響強度（エネルギー束密度）の最適化（全ＳＩＬ）：ａ．所与の方向に沿った音響強度に比例する
2020502906.tif 000196
を最大化する；ｂ．これは波の方向性と（皮膚）表面
2020502906.tif 000197
に対する法線の影響の両方を考慮に入れる。
【０１６６】
ＩＩＩ．圧力を最大化するための最適化
【０１６７】
各コントロールポイント点で圧力を最大化することができる二次形式最適化装置は、量をモデル化することになる:
2020502906.tif 000198
各コントロールポイントjでは、目的関数の総和を求める：
2020502906.tif 000199
これは、支配的な固有ベクトルが解かれる
2020502906.tif 000200
マトリクス要素を以下のようにする：
2020502906.tif 000201
ここで
2020502906.tif 000202
は、波の圧力の正および実振幅であり、初期放出からのコントロールポイント
2020502906.tif 000203
で、トランスデューサ
2020502906.tif 000204
から進行する波上の複素位相オフセット
2020502906.tif 000205
である。重み付け
2020502906.tif 000206
は、コントロールポイント間の相対振幅コントロールを提供するために、二次
2020502906.tif 000207
項を再重み付けするために使用される。グローバル振幅スケーリングは、入力パワー
2020502906.tif 000208
の単位ノルム条件によって提供される。これを二次重み付け係数およびコントロールループと組み合わせて、個々の圧力振幅に正確に対処することができる。
【０１６８】
ＩＶ．粒子速度を最大化するための最適化
【０１６９】
粒子速度を最大化することができる二次形式最適化装置は、量をモデル化することになる:
2020502906.tif 000209
各コントロールポイントｊについて行われる。これは、下記の目的関数の総和を求めることとなる：
2020502906.tif 000210
これは、支配的な固有ベクトルが解かれる 
2020502906.tif 000211
マトリクス各要素を以下のようにする：
2020502906.tif 000212
ここで
2020502906.tif 000213
は、初期放出からコントロールポイント
2020502906.tif 000214
で、トランスデューサ
2020502906.tif 000215
から進行する波上での誘導粒子速度振幅のｘ方向部分（これは実数であり正数である）および
2020502906.tif 000216
複素位相オフセットである。重み付け
2020502906.tif 000217
は、コントロールポイント間の相対的な粒子速度振幅制御を提供するために、二次
2020502906.tif 000218
項に重みを付けるために使用される。グローバル振幅スケーリングは、入力パワー
2020502906.tif 000219
の単位ノルム条件によって提供される。これを二次重み付け係数およびコントロールループと組み合わせると、個々の速度振幅を正確に使用することができる。
【０１７０】
Ｖ．所与の方向のセットで粒子速度を最大化するための最適化
【０１７１】
所与の方向のセットにおいて、粒子速度を最大化することができる二次形式最適化装置は、量をモデル化することになる:
2020502906.tif 000220
各コントロールポイントｊについて行われる。これは、下記の目的関数の総和を求めることになる：
2020502906.tif 000221
これは、支配的な固有ベクトルが解かれるｎ×ｎマトリクス各要素を次のようにする：
2020502906.tif 000222
ここで、
2020502906.tif 000223
は、誘導粒子速度振幅の
2020502906.tif 000224
方向部分（これは実数であり正数である）であり、
2020502906.tif 000225
は、初期放出からコントロールポイント
2020502906.tif 000226
でトランスデューサ
2020502906.tif 000227
から進行する波上での複素位相オフセットである。重み付け
2020502906.tif 000228
は、コントロールポイント間の相対振幅コントロールを提供するために、二次項
2020502906.tif 000229
を再重み付けするために使用される。グローバル振幅スケーリングは、入力パワー
2020502906.tif 000230
に対する単位ノルム条件によって提供される。これを二次重み付け係数およびコントロールループと組み合わせて、個々の振幅を正確にアドレス指定するために使用できる。
【０１７２】
ＶＩ．活性音響強化
【０１７３】
活性音響強度
2020502906.tif 000231
は、音響媒体の圧力
2020502906.tif 000232
および粒子速度ベクトル
2020502906.tif 000233
を用いて定義することができる。音響波は、下記のように記述できる：
2020502906.tif 000234
ここで、活性強度は、純粋に実数のベクトル量である。波源の近位場では、活性場は、この代わりに
2020502906.tif 000235
が支配的である。しかしながら、1つまたは２つの波長を超えた遠位場では、
2020502906.tif 000236
となる。この理由により、
2020502906.tif 000237
のみが保証現実であると考慮し仮定される。ここで、単調な平面波源を考慮する：
2020502906.tif 000238
我々は、下記式を得る:
2020502906.tif 000239
ここで、
2020502906.tif 000240
は圧振幅であり、
2020502906.tif 000241
は各方向で粒子速度振幅である。
【０１７４】
所与の軸に沿って波によって伝達する音の強度、したがってエネルギー流束もまた重要であるので、スカラー寄与を測定するための内積が生成される:
2020502906.tif 000242
【０１７５】
音響波方程式の形成により：
2020502906.tif 000243
重ね合わせの原理は、圧力ベクトルと粒子速度ベクトルの両方に適用される関心のある点で、
2020502906.tif 000244
音響源を想定し各平面波として概算すると、方程式が以下のように得られる：
2020502906.tif 000245
【０１７６】
問題の領域では、
2020502906.tif 000246
であり、
2020502906.tif 000247
が同相であることが周知であり、その場合、共役の配置は問題ではなく、粒子速度の代わりに、共役圧力への共役の順序を逆にすることが可能である。（これは上記の方程式を拡張し、可換性を使用して、各項の複素べき乗の順序を逆にすることで確認できる）推論から、これはエルミートマトリクス
2020502906.tif 000248
を生成する方法で記述することができる；
2020502906.tif 000249
【０１７７】
ＶＩＩ．音響強化を最大化するための最適化
【０１７８】
スカラー値を生成する最初の数量は以下である：
2020502906.tif 000250
しかしながら、これは遠場であると想定され、これは以下のようにも記載できる：
2020502906.tif 000251
これは圧力和の４乗だが、複素べき乗の構造上、これらのうち２つは同じで、２乗根を使って削除することができる。
【０１７９】
したがって、音響強度を最大化することができる二次形式オプティマイザは、量をモデル化することになる：
2020502906.tif 000252
各コントロールポイント
2020502906.tif 000253
についてである。これは、下記の目的関数の相和を得ることになるであろう：
2020502906.tif 000254
これは、支配的な固有ベクトルが解かれる
2020502906.tif 000255
マトリクス各要素を次のようにする：
2020502906.tif 000256
ここで、
2020502906.tif 000257
は波前の普通のベクトルの部分方向
2020502906.tif 000258
であり、
2020502906.tif 000259
は、初期放射からコントロールポイントｊにおいてトランスデューサｑから進行する波上の波の圧力の正および実振幅、および
2020502906.tif 000260
は、複素位相オフセットである。重み付け
2020502906.tif 000261
は、コントロールポイント間の相対的な粒子速度振幅コントロールを提供するために、圧力の４乗を含む関数の平方根
2020502906.tif 000262
を実際に維持する二次項を再重み付けするために使用される。グローバル振幅スケーリングは、入力パワー
2020502906.tif 000263
の単位ノルム条件によって提供される。これを二次重み付け係数およびコントロールループと組み合わせると、個々の速度振幅を正確に扱うことができる。
【０１８０】
ＶＩＩＩ． 所与の方向のセットでの音響強化の最大化のための最適化
【０１８１】
所与の方向に沿って音響強度を最大化することができる、二次形式最適化装置は、量をモデル化することになる：
2020502906.tif 000264
しかし、これを遠位場で想定すると、これは以下で記載できる：
2020502906.tif 000265
これは、各コントロールポイント
2020502906.tif 000266
についてである。これにより、目的関数の合計は次のようになる：
2020502906.tif 000267
これの支配的な
2020502906.tif 000268
マトリクスの各要素を生成する固有ベクトルが解かれるべきである：
2020502906.tif 000269
ここで、
2020502906.tif 000270
は誘導された粒子速度振幅の
2020502906.tif 000271
方向の部分であり（これは実数で正数である）、
2020502906.tif 000272
は波圧の正、かつ実の振幅であり、
2020502906.tif 000273
は、初期放射からコントロールポイント
2020502906.tif 000274
でトランスデューサ
2020502906.tif 000275
から進行する波上での複素数位相オフセットである。重み付け
2020502906.tif 000276
は、コントロールポイント間の相対的な粒子速度振幅コントロールを提供するために、二次の
2020502906.tif 000277
項を再重み付けするために使用される。グローバル振幅スケーリングは、入力パワー
2020502906.tif 000278
に対する単位ノルム条件によって提供される。これを二次重み付け係数とコントロールループと組み合わせて使用することで、個々の速度振幅を正確に扱うことができる。
【０１８２】
ＩＸ．重み付けファクターの発見
【０１８３】
レイリーリッツの定理は、支配的な固有ベクトルを見つけ、トランスデューサの物理的限界までパワーを増加させることによって、上記のマトリクスのいずれも最大化できることを示す。これは、パワー相互を通して到達することができる。コントロールポイントの合計のそれぞれを別々に定義しながら、マトリクス
2020502906.tif 000279
を構成することによって、ベクトル
2020502906.tif 000280
を乗じる場合、目的関数の評価を追加の副次的な作用として実行することができる。これは、重み付けされたときにシステムの解に適切である振幅を最大化が見つけることを可能にするように、重み付け
2020502906.tif 000281
に対する適切な値を見つけるために使用することができる。
【０１８４】
パワー反復処理は、下記のように記述される：
2020502906.tif 000282
ここで：
2020502906.tif 000283
【０１８５】
我々は、各コントロールポイントで目的関数を決定するために、パワー反復処理を使用することができる。
【０１８６】
定義づけを行う：
2020502906.tif 000284
かつ：
2020502906.tif 000285
各反復処理で、計算するのではなく以下とすることができる：
2020502906.tif 000286
【０１８７】
つまり、各コントロールポイントでの目的関数は、以下のように評価することができる：
2020502906.tif 000287
【０１８８】
最終的に、以前のように
2020502906.tif 000288
の計算をすることができる：
2020502906.tif 000289
ここで、最終的な固有ベクトルの解は、上記の通りである。しかし、このように計算すると、
2020502906.tif 000290
が、コントロールポイント
2020502906.tif 000291
の中間結果を出すことができる。これが所望のコントロールポイント出力を表すレベルを固定するために、一組
2020502906.tif 000292
の値が、ユーザによって与えられる目的値と直接比較されることを可能にするために、平均目的関数の値が取られる。このケースでは、重み付け係数は、例えば、Ｎｅｌｄｅｒ−Ｍｅａｄアルゴリズムなどの高価な関数の最適化を目的とした導関数のない方法を使用して、次の反復処理で更新することができる。
【０１８９】
べき乗法と結合された導関数のない最適化を何度か繰り返した後、中間結果
2020502906.tif 000293
を使用してトランスデューサを駆動するべきである。したがって、これにより、システムは、ある程度の誤差を許容しながら、時間の経過と共に、目的関数の所望の値を追跡することが可能になる。このシステムの利点は、それがトランスデューサアレイを駆動するのに、非常に簡単な方法であるということである。最大の不利な点は、多くの計算が必要とされることだが、計算は実質的に複雑な乗算にすぎず、十分なハードウェアで簡単に実装される。
【０１９０】
Ｘ．添付の図面
【０１９１】
図１９に示されるものは：左上１９１０：非直交ベクトルを使用して二次粒子速度問題を解くことによって作成されたコントロールポイントであり、
2020502906.tif 000294
平面を通るカットとして表示される。右上１９２０：
2020502906.tif 000295
平面を通してカットとして表示した、同コントロールポイント。左下１９３０：比較すると、二次圧力問題を解くことによって作成されたコントロールポイントであり、
2020502906.tif 000296
平面を通るカットとして表示される。右下１９４０：
2020502906.tif 000297
平面を横切る圧力問題の解を通る断面。
【０１９２】
５．二次問題を通して指定された位相アレイのためのコントロール領域
【０１９３】
本開示は、ユーザ供給設計パラメータを使用して、基底関数を自動的に生成するための方法を説明する。これらは、オンラインで設定された所与のトランスデューサ基底の空間における、二次関数にわたる効率的な最適化に入力される。これにより、ユーザは、様々な異なる効果を生み出す可能性があるオンライン最適化を通じて、アクセス可能なデザインスペースを手に入れることができる。二次問題のためのコントロール領域であるので、「二次コントロール領域」と呼ばれる線形演算子を使用してトランスデューサ効果にわたって二項積を評価することによって、マトリクスを構築することができる。このマトリクスの主固有ベクトルは、このマトリクスの固有ベクトルである、正規直交解の集合から最良の解として選択されてもよい。最適化は、出力としてトランスデューサ活性化係数−個々のトランスデューサ素子を駆動して音場または所望の特性を有する現象を再現するのに必要な、生の複素数値を生成する。
【０１９４】
次いで、これらの生の複素数値を使用して、後に使用する線形システムの基底関数をさらに含むことができ、その解によって、ユーザは、それぞれに利用可能な出力パワーを再重み付けしながら、すべての所望の効果を同時に生成することができる。
【０１９５】
Ｉ．線形コントロール領域
【０１９６】
線形コントロール領域は、複素数値重み付け関数
2020502906.tif 000298
を用いて、圧力などの線形音響量との積分を実行する演算子として定義される：
2020502906.tif 000299
これは、音響場を評価するために使用される。
2020502906.tif 000300
が音響量における類似の位相分布と対にされる場合、複素共役が、関数評価を最大にするために使用されることに留意されるべきである。この積分に対する任意の近似値（サンプリングなど）を使用することができ、あるいは実際にはこれを使用して、音圧場にわたって任意の線形演算子を近似することができる。それがまた線形であるならば、ベクトル粒子速度を使用してコントロール領域を最適化するようにさらなる積分を構築することができる：
2020502906.tif 000301
ここで、
2020502906.tif 000302
は、媒体の粒子速度、
2020502906.tif 000303
は、各次元が複素数である三次元ベクトル場であり、これは、上記と同じ目的で重み付け関数として機能する。この方法では、音響量のさまざまなパターンが線形的に評価され、位相が
2020502906.tif 000304
のスカラー関数またはベクトル関数によって作成された「テンプレート」と一致するときに、結果の特定の位相が最大になる。
【０１９７】
また、
2020502906.tif 000305
は、複数の外見上ばらばらな領域の単一の評価を作成するために、ほぼゼロ（またはゼロベクトルで構成される）である場合があることにも留意すべきである。この定式化は、線形音響量のフィールドにわたって、任意の線形演算子に拡張することができるため、任意の整数近似（サンプリングなど）を使用することができる。
【０１９８】
位相アレイは、各トランスデューサをコントロールする一連の複素数値の係数によって、記述することができる。マトリクス
2020502906.tif 000306
に符号化された音響量を最大にするために、そのような入力
2020502906.tif 000307
ベクトルを用いて、固定周波数で音響変換器のアレイを駆動することは、線形代数問題として表すことができ、その意図は以下の通りである：
2020502906.tif 000308
【０１９９】
問題の記述は固有ベクトル、ここでは
2020502906.tif 000309
の定義でもあるため、これはマトリクス
2020502906.tif 000310
の主要な固有ベクトルをとることによって解くことができる。マトリクス
2020502906.tif 000311
を構築することは、まず、問題の音響量の大域空間における基底関数として各トランスデューサによって放出された線形音響量を考慮することによって達成される。これは、
2020502906.tif 000312
の絶対位置で表される、１点での圧力やその他の音響量の最大化を達成するための方法である。しかしながら、記述したコントロール領域積分を使用することにより、各個々の変換素子の場に放出される音響量を支配する、
2020502906.tif 000313
は、
2020502906.tif 000314
と交換され、列ベクトルを生じる:
2020502906.tif 000315
【０２００】
線形性を呼び出すことによって、各トランスデューサを駆動する複素活性化値を乗算して、ベクトルを導き出すことができる：
2020502906.tif 000316
ここで、
2020502906.tif 000317
は、トランデューサの数である。この
2020502906.tif 000318
が線形場に作用する線形演算子の結果とすると、それはエルミート転置され、それ自身によって乗算されて二次形式を生成することができる、最終的には以下のようになる。
2020502906.tif 000319
【０２０１】
重み付け関数
2020502906.tif 000320
を修正し、マトリクス
2020502906.tif 000321
の固有ベクトルを解くことによって、様々な異なるビーム、表面および他の音響現象が生成され最大化され、音響内に最適に作動される領域を生成する。
【０２０２】
先に論じたヌル空間生成技術（「空中ハプティックシステムにおける二次問題のための固有システム解」）もまた、音響浮揚を記述するこの方法を拡張するために使用され、ならびにヌルポイントを記載する構造のポイントのフィールドへ導入することによって、ヌルステアリングおよび空中テクスチャ形成を可能にするために使用される。このプロセスはマトリクスの階数を増加させるため、これを使用したとしても単純化は不可能であることに注意すべきである。
【０２０３】
等価的に、多数のコントロール領域について、
2020502906.tif 000322
個のマトリクスを一緒に加えることによって複数の最適化を実行することができる。これは、コントロールポイントを表すマトリクスを含むこともできる。このことは、ランクが上がったマトリクスを生成し、いくつかの簡略化された方法は使用される本ケースでは不可能である。
【０２０４】
結果を考慮し、反復的に各点を再重み付けすることによって、音響量の必要な幾何学的構成に収束するように、システムを修正することができる。これは、振幅が出力で構成されない場合に、一般的に必要である。
【０２０５】
これは、生成した音響量を合計し、次の反復において結果にさらに近づくようにシステムを再重み付けすることによって達成する。例えば、サブスプリットボリュームの積分点または単一の圧力点に対する反復再重み付けは、次のようになる：
2020502906.tif 000323
ここで、
2020502906.tif 000324
は、本技術における次の粗な反復、
2020502906.tif 000325
は、現在の反復、
2020502906.tif 000326
は、所望の圧力（これは複素数でもよい）、そして分母は現在の反復によって得られた全圧力を指す。得られた領域の振幅にのみ影響を与え、位相には影響を与えないために、分数は実数のみとしてレンダリングされてもよい。
【０２０６】
ベッセルビームは、ビームの進行方向に重み付け関数
2020502906.tif 000327
を表現し、ビームが前方に移動するにつれて、前方に位相角を修正しながら、単位圧力を生成することによって生成される。後続の各ポイントが通過する角度がビーム軸に沿って位相をサンプリングするのに十分小さい限り、これはポイントサンプルでサンプリングされる。
【０２０７】
位相のサンプリングは、それが移動する角度において、
2020502906.tif 000328
ごとに２つより多い位相サンプルを有することが重要である。これは、経路に沿った増分が
2020502906.tif 000329
の場合、波の進行方向が決定できないためである。これはまた、多くの場合、位相がカバーする単色周波数距離の各波長に対して、
2020502906.tif 000330
ラジアンの回転を完了するように、位相が前方に増加しなければならないことを意味する。これは、コントロール領域の定義で説明される形式のパス積分への近似である。
【０２０８】
このように書かれた積分を最大化する結果として、トランスデューサは、アレイから遠ざかる直線で累進的に定義された単純なビームの場合、その理論的能力により最適なベッセル状ビームを生成し、この構成で最適な圧力を生成するだろう。図２０に示されるのは、平面直線トランスデューサアレイに平行な近距離場領域の端部でのフィールドの断面であり、複素数値のトランスデューサ活性化係数の構造を明らかにしている。単純なプログレッシブ位相ビームではベッセル関数近似が得られ（左２０００）、アンダーサンプリングされたヘリカルビームでは双極子の振る舞いを持つベッセル関数が得られ（右２０１０）、双極子は２つの構造を作り出し、それらは互いに位相で反対になるように活性化する。トランスデューサの解像度の限界、有限アレイサイズ、および最適化が回避する必要があるその他の物理的考慮の影響により、ビームは正確なベッセルビームではない可能性があるが、実際のトランスデューサアレイの物理的限界を考慮した最適化により、ベッセルビーム関数の単純な評価よりも効果的に機能する可能性がある。大きなアレイでは、これはベッセル関数に密接に従っており、図２１は、図２０のアレイの中心から半径方向に考慮したときに圧力と共にプロットされたベッセル関数のグラフ２１００を示している。
【０２０９】
このシステムのさらなる特性は、任意の生成された単色場関数を有する任意の位置、および向きのトランスデューサが使用されることである。これは、定義により平面アレイの場合にのみ機能する、サンプリング方法に対するベッセルビーム生成における重要な利点である。
【０２１０】
複数の位相外経路を作成することによって、位相シフトおよび位相特異点は、経路の周りにらせん波を巻くことによって暗黙的に作成される。固有経路アプローチに組み込まれている単純な二次最適化を使用しながら、この経路はまた湾曲しており、高度の構成可能性をもたらす。これは、ビーム幅より近いが互いに位相が異なる多くの平行ビームラインにおける別々の位相巻線の定義を通して、ビームの回転を誘導することによって達成する。次に、図２２に示すように、近くの位相をビームの中心線で互いに相殺するように設定することで、高次のベッセルビームを形成できる。図２２では、円の周りの漸進的位相の２つ以上のサンプル（したがって、らせん配置の位相）は、左側２２１０に見られるように角運動量を有するベッセル関数の形成を生じさせる。らせんの周りの位相を特定するために、
2020502906.tif 000331
ラジアンを超える値を使用することによって、右側（２２２０）の
2020502906.tif 000332
ビームのように高次のベッセルビームを生成することができる。
【０２１１】
これらの位相線は、互いに素な領域が単一の領域最適化としてどのように定義され得るかの例であるので、点または他の任意の積分近似を用いてサンプリングされ、最適化するボリューム内でのみゼロ以外の値に機能する。高次マシュー型楕円ヘリカルビームは、位相が漸増し楕円の周りに間隔を置いて増加する、位相ビームを導入することによって生成することができる。図２３に示すように、マシュー型ビームなどの高次ベッセルビームと同様の特性を有する楕円形ヘリカルビームは、位相の楕円の円周をパラメータ化することによって作成することができる。
2020502906.tif 000333
ラジアンの割合を決定し、楕円弦の円周に比例する位相角を割り当てることは、ビームの最適化のための適切な構成を生成する。左の画像２３１０は、マシュー型ビームの近くで、楕円形に再形成された高次ベッセルを示す。右側の画像２３２０は、中央領域に形成され始めている架橋された、らせん構造を有するおおよそのマシュー型ビームを示す。どちらも
2020502906.tif 000334
ラジアンの位相円に基づいており、らせん状の波面を形成する。
【０２１２】
このようにして生成されたベッセルファミリーのビームは、文献においてより制限的な形態で存在するよりエキゾチックなビームを生成するように修正されてもよい。波の横断に適した方法で前進する物理的にもっともらしい位相関数によって横断される場に任意の経路を作成することによって、さまざまな結果が再現されるかもしれないが、より一般的な形で固有システムを伴う最適化のため、ビームは積分によってのみ表現される必要があるので、初期の帯域幅コストは以前よりも低い。
【０２１３】
湾曲した経路をとり、位相を前方に動かすと、４５°の円セグメントによって生成された障害物の周りにビームを操縦するために、一般的に使用される湾曲ビームが生成される。図２４に示されているのは、ビームに沿った音圧を等化するために粗い反復法を使用して、ラインに沿った連続に接続された円の４５度セグメントをたどる位相関数から形成された曲げビームの画像２４００である。トランスデューサアレイは、画像の下端に沿って配置される。
【０２１４】
複数の位相経路を含むさらなる手法は、「ボトルビーム」、すなわち中心から切り取られた部分を有するビームを作り出すことができ、それは大きな物体の浮揚を容易にすることができる。これは、図２５において生成されており、図２５は、それに沿って位相がたどられるようにされてきた、２つのガウス曲線を双晶することの画像２５００を示す。位相関数を作成することによって生成されたボトルビームは、ビームに沿って音圧を均一化するために、粗い反復法を使用してガウス曲線に従う。トランスデューサアレイは、画像の下端に沿って配置される。下部領域の点線の性質は、最適化に入力される位相速度が遅すぎること、および初期区間により速い位相速度を使用することによって、ビームが改善される可能性があることを示唆している。物体が十分に硬い場合には、空間の一部のみを高音圧領域にする必要があるので、図２６のように曲線の一部を使用することもできる。図２６は、不均一な方法で高音圧の領域を最適化することによってさらに大きな物体を浮揚させるために使用することができる、累進相で画定された円の一部から形成された不完全ボトルビームの画像２６００を示す。浮揚すべき物体または物体が低音圧領域を通って流れない、押し出される、または漏れることがないように十分な凝集力を有することを仮定する場合、大きな物体に及ぼされる力を均一に分散させる必要はない。
【０２１５】
ベッセル状ビームは、経路を位相速度より速く前方に移動させることによって生成することができる。この方法では、図２７に示すように、ビームは通常よりも速い位相速度で最適化され、波長が長くなり、位相速度が音速よりも速く移動するように見えるビームを作成する。図２７に示されるのは、音速よりも速く進む波面の画像２７１０であり、これは、逆説的に同じ周波数でより長い波長を生成するように、角を通って移動する、位相経路を実施することによって生成され得る。暗黙の最適化は、メインビームが交差するビームで構成されていると推論し、位相速度の指定された増加を引き起こす。右側の画像２７２０より見ることができ、波長は、ビームの最初の部分の同じ周波数で約４０％長くなるように操作される。蛇行ビームは、正弦曲線に沿って漸進的位相をサンプリングすることによって構築することができる。図２８に示されているのは、例２８１０であり、そこでは正弦曲線の短い部分が定義され、続いて位相配列上の正弦曲線に続く、漸進的な圧力相から構成される２８２０のより長い２次元蛇行ビームが続く。これらは、漸進的な位相のさらなる部分を単に追加することによる蛇行正弦波ビームの拡大が必要とされる限り（およびトランスデューサの物理的配置によって可能にされる）自然にビームを拡大させることを示す。らせん状に蛇行するビームを生成するために３次元全てに蛇行するビームを作ることによって、アルキメデススクリューシステムと同様の方法で、潜在的には波の放射、波長よりもはるかに大きいフィールド内の固体物体に対して、垂直または反対方向に力を加えるために使用できるビームを生成することができる。これは、それによってビームが物体を浮揚させるために使用される、異なるメカニズムを与える。図２９は、
2020502906.tif 000335
平面内の断面を詳細に示す左画像２９１０と、
2020502906.tif 000336
平面内のビームの断面を示す右画像２９２０とを示す。
2020502906.tif 000337
に曲がっているが、
2020502906.tif 000338
にも進行する円を描く経路をたどるこの技術を使用して、真のらせん形状を作り出すことができる。これは、トラクタービームのような仕方で、波長よりはるかに大きい物体に力を加えて浮揚させるために使用することができる。
【０２１６】
固有システムマトリクスが単一のコントロール領域で定義されている場合、そのマトリクスはランク１であり、複素数値ベクトル
2020502906.tif 000339
とそれ自身のみの２項積で形成し、マトリクスを
2020502906.tif 000340
と定義する。この特定のケースでは、システムの単一の固有ベクトルはすぐに分かり
2020502906.tif 000341
であり、
2020502906.tif 000342
を意味する。このケースで目的関数を評価することによって、複素数値のトランスデューサ起動の係数のそれ自身を特定する。これは各トランスデューサのコントロール領域機能を単純に評価することによって定義し、活性化の係数である。
【０２１７】
従前のとおり、二次形式のために、各反復についての再重み付けを通して、有効性を評価することもまた追求される。これは、各ポイントにおける圧力または各小領域にわたる圧力の積分を評価し、所望の圧力で除算し、新しい関数値を生成して、
2020502906.tif 000343
小領域
2020502906.tif 000344
の重みを更新することによって達成する。それから
2020502906.tif 000345
でもある単一の固有ベクトルはそれから再評価されて、前から固有システム上で繰り返される反復に対する同等のアプローチを生成するが、非常に減少した計算コストで再評価される。
【０２１８】
媒体の粒子速度でこれを達成するために、粒子速度の試験を使用することができ、これは最大に達するまで方向を変えることになる：
2020502906.tif 000346
圧力は、ハイブリッド方法を生成する旧技術を使用して監視することもできる。このハイブリッド技術では、重み関数を指定して以前と同じように圧力を更新するが、効果的なコントロール領域を使用する：
2020502906.tif 000347
スカラーだが可変の複素数値
2020502906.tif 000348
を得、不変の重み付け法線ベクトル
2020502906.tif 000349
と粒子速度
2020502906.tif 000350
と対とする。このアプローチは、求められる波の方向を変えることなく、満足のいく圧力を見つけるために最適化することができる。これは、次のセクションで説明されている特定の修正と共に図２８に示されている。
【０２１９】
ＩＩ．指向性コントロールポイントおよび領域に対する波面フィルタリング
【０２２０】
指向性出力の作成に必要なさらなる技術は、トランスデューサ出力波面
2020502906.tif 000351
と法線ベクトルまたは粒子速度ベクトル
2020502906.tif 000352
とが反対である場合に、トランスデューサデータを破棄し、それをゼロに設定することである：
2020502906.tif 000353
考えている領域のポイントについてである。
【０２２１】
図３０は、トランスデューサアレイからさらに定義された２つの速度ベクトルサンプルの左画像３０１０と、解除され定義されたコントロール方向と反対の方向に波面があるトランスデューサを有した、アレイに近接して定義された２つの速度ベクトルサンプルの右画像３０２０とを示す。これら２つの互いに素な領域で指定された速度重み付け関数を使用することは、それぞれの場合において単一の領域を構成し、その領域は、２つの別々のポイントにおける、媒体の粒子速度を表す方向ベクトルに対して最適化される。
【０２２２】
2020502906.tif 000354
方向の最大量の速度を見つけようとすることは、結果として生じる粒子速度方向が完全に
2020502906.tif 000355
方向にあることを、必ずしも意味しないことに留意されたい。図３０は、所望の方向とは反対の波面を提供する、相反するトランスデューサをフィルタリングした結果を示し、図３１は、トランスデューサの波面方向が、波面基準を使用するフィルタリングされてなかった結果である、左画像３１１０と右画像３１２０を示す。相反するトランスデューサは、図３１の場合、自由に干渉することができ、その病状は、定在波構造の形で水平に明らかである。これは、波面方向が粒子速度から欠落している方向情報を伝達するため、粒子速度は波がどの方向から生じたかに関する情報をもたらさず、望ましくない干渉を可能にするからである。
【０２２３】
これはまた、コントロール領域だけではなく、通常のコントロールポイントの単純な場合にも重要である。コントロールポイントの場合は、事実上単純化されたコントロール領域であり、粒子速度の単一点サンプルと重み係数なしで近似される。どちらの場合でも、「役に立たない」トランスデューサを除外しないことは、「定在波」現象の影響の要素に関して同様の結果をもたらします。これは、複素成分を含む粒子速度ベクトルでは、異なる方向に進む波の違いを特定できないためである：このモデルでは、位相０で左から右へ進む波は、位相
2020502906.tif 000356
で右から左へ進む波と数学的に同じである。結果として、これがコントロールポイントまたは領域ベースの生成で正しく処理されない場合、最適化は両方を使用しようとし、意味のあるように寄与しない定在波成分を含む解を生成し、不要な、余分なフィールドのないヌルの領域を作成する。この最適化ステップを用いずに基底関数として線形最適化、固有システム、さらには単純な式に入力した、コントロールポイントがこれらのアーチファクトを示すため、これは重要な点である。
【０２２４】
記載のこれらの技術は、本明細書で定義されるような単一の固有システムを反復することによって、または単純な場合には単一の点または領域を作成することによって、任意に構成された領域およびビームを作成することができる、音響位相アレイシステムを作成するために装置に統合される。単一のポイントまたは領域の場合には、トランスデューサ素子を駆動するのに必要な振幅および位相が、ポイントまたは領域を評価することによって表される、振幅および位相であることを示す簡単な式を生成するために最適化が適用される。大域的方法の反復および経時的な反復適用は、様々な触覚効果を生み出すために、時間依存的に変更および修正される。これは、それらが機能的に関連しているので、所望の場および導出されたトランスデューサ係数の両方に適用される。
【０２２５】
ＩＩＩ．可能な請求項
１．単調な音響場についての方法であって、
前記音圧場関数と重みづけ関数との積を積分することによって音圧の分布を単一値に見積もることであって、両方の関数が複素数値であり、物理的なトランスデューサの特性から導出されることと、
前記マトリクスのヌル位置に追加されるゼロまたはヌルポイントを適用することと、
これを最適化において最大化するべき目的関数として使用することであって、前記最適化の結果は、各トランスデューサについて前記目的関数の交差乗算項から形成されるマトリクスの固有ベクトルの形式をとることと、
そのようなマトリクスの1つまたはそれ以上を合算して1つまたはそれ以上の固有ベクトルが前記特定の場を生成するような固有システムを生成することを含む方法。
２．単調な音響場についての方法であって、
実際の３次元重み付けベクトル場と粒子速度場と重み付け関数とのドット積を積分することにより単一のベクトル値へと音響粒子速度の分布を見積もることであって、両方の関数が複素数値であり、物理的なトランスデューサの特性から導出されることと、
前記マトリクスのヌル位置に追加されるゼロまたはヌルポイントを適用することと、
これを最適化において最大化するべき目的関数として使用することであって、前記最適化の結果は、各トランスデューサについて前記目的関数の交差乗算項から形成されるマトリクスの固有ベクトルの形式をとることと、
そのようなマトリクスの1つまたはそれ以上を合算して1つまたはそれ以上の固有ベクトルが前記特定の場を生成するような固有システムを生成することを含む方法。
３．単調な音響場についての方法であって、
実際の３次元重み付けベクトル場と粒子速度場と重み付け関数とのドット積を積分することにより単一のベクトル値へと音響粒子速度の分布を見積もることであって、両方の関数が複素数値であり、物理的なトランスデューサの特性から導出されることと、
それぞれ見積もられた積分の副領域で特定される垂直に対して反対方向を有するポイントでの波頭のトランスデューサの寄与をゼロに設定することと、
前記マトリクスのヌル位置に追加されるゼロまたはヌルポイントを適用することと、
これを最適化において最大化するべき目的関数として使用することであって、前記最適化の結果は、各トランスデューサについて前記目的関数の交差乗算項から形成されるマトリクスの固有ベクトルの形式をとることと、
そのようなマトリクスの1つまたはそれ以上を合算して1つまたはそれ以上の固有ベクトルが前記特定の場を生成するような固有システムを生成することを含む方法。
４．単調な音響場についての方法であって、
前記音圧場関数と重み付け関数との積を積分することによって音圧の分布を単一値に見積もることであって、両方の関数が複素数値であり、物理的なトランスデューサの特性から導出されることと、
前記マトリクスのヌル位置に追加されるゼロまたはヌルポイントを適用することと、
これを最適化において最大化するべき目的関数として使用することであって、前記最適化の結果は、各トランスデューサについて前記目的関数の交差乗算項から形成されるマトリクスの固有ベクトルとなることと、
前記トランスデューサの入力信号のための最適な位相および振幅を記述するために各トランスデューサの寄与を使用することを含む方法。
５．実際の３次元重み付けベクトル場と粒子速度場と重み付け関数とのドット積を積分することにより単一のベクトル値へと音響粒子速度の分布を見積もることであって、両方の関数が複素数値であり、物理的なトランスデューサの特性から導出されることと、
前記マトリクスのヌル位置に追加されるゼロまたはヌルポイントを適用することと、
これを最適化において最大化するべき目的関数として使用することであって、前記最適化の結果は、各トランスデューサについて前記目的関数の交差乗算項から形成されるマトリクスの固有ベクトルの形式をとることと、
前記トランスデューサの入力信号のための最適な相および振幅を記述するために各トランスデューサの寄与を使用することを含む方法。
６．  実際の３次元重み付けベクトル場と粒子速度場と重み付け関数とのドット積を積分することにより単一のベクトル値へと音響粒子速度の分布を見積もることであって、両方の関数が複素数値であり、物理的なトランスデューサの特性から導出されることと、
それぞれ見積もられた積分の副領域で特定される垂直に対して反対方向を有するポイントでの波頭のトランスデューサの寄与をゼロに設定することと、
前記マトリクスのヌル位置に追加されるゼロまたはヌルポイントを適用することと、
これを最適化において最大化するべき目的関数として使用することであって、前記最適化の結果は、各トランスデューサについて前記目的関数の交差乗算項から形成されるマトリクスの固有ベクトルの形式をとることと、
前記トランスデューサの入力信号のための最適な相および振幅を記述するために各トランスデューサの寄与を使用することを含む方法。
７．さらに、デバイスが、方法と、音響トランスデューサを駆動するための結果の使用を実装する、請求項１〜６に記載の方法。
８．さらにデバイスが、方法と音響トランスデューサを駆動するために使用する線形基底セットの一部としての結果の使用を実装する、請求項１〜６に記載の方法。
９．音響ビームを生成する、請求項１〜６に記載の方法。
１０．音響量を生成する、請求項１〜６に記載の方法。
１１．音響波面を生成する、請求項１〜６に記載の方法。
１２．音響場についての方法であって、
物理的なトンランスデューサの特性から得られる場の線形結合を使用して1つのポイントにおける単一の値として音響粒子速度を見積もることと、
前記コントロールポイントにおいて、それぞれ見積もられた積分の副領域で特定される垂直に対して反対方向を有するポイントでの波頭のトランスデューサの寄与をゼロに設定することと、
これを前記トランスデューサの入力信号のための相および振幅を見出すための最適化問題における基底関数として使用することと、を含む方法。
１３．音響場についての方法であって、
物理的なトンランスデューサの特性から得られる場の線形結合を使用してポイントにおける単一の値として音響粒子速度を見積もることと、
前記コントロールポイントにおいて、それぞれ見積もられた積分の副領域で特定される垂直に対して反対方向を有するポイントでの波頭のトランスデューサの寄与をゼロに設定することと、
これを最適化において最大化するべき目的関数として使用することであって、前記最適化の結果は、各トランスデューサについて前記目的関数の交差乗算項から形成されるマトリクスの固有ベクトルとなることであって、その後トンランスデューサの入力信号としての位相および振幅を見出すことと、を含む方法。
１４．音響場についての方法であって、
物理的なトンランスデューサの特性から得られる場の線形結合を使用してポイントにおける単一の値として音響粒子速度を見積もることと、
前記コントロールポイントにおいて、それぞれ見積もられた積分の副領域で特定される垂直に対して反対方向を有するポイントでの波頭のトランスデューサの寄与をゼロに設定することと、
トランスデューサの入力信号のため最適な相および振幅を記述するために各トランスデューサの寄与を使用することを含む方法。
１５．  前記トランスデューサは、得られた構造が入力信号として適用され、変調された場合にハプティック効果を生成する、請求項１〜１４に記載の方法。
１６．  前記トランスデューサは、前記得られた結果が入力信号として適用され、時間にわたり変調された場合にハプティック効果を生成する、請求項１〜１４に記載の方法。
【０２２６】
６．ハードウェアにおける複素数値の超越関数、乗算および分割の評価
【０２２７】
2020502906.tif 000357
の代わりに
2020502906.tif 000358
の間で変換することを目的として、ＢＫＭアルゴリズムは物理的なハードウェア実装においてより有用である。
【０２２８】
Ｉ．イントロダクション
【０２２９】
ハードウェアにおける、三角関数および他の特別な関数の計算のための多くの技術が存在する。そのような方法の１つはＢＫＭ法であり、これは複素対数を複素べき乗に交換し、その逆の変換を行い、これらの特別な関数をある表現から別の表現に移動することの副作用として生成する。特に興味深いのは、そのような関数が必要であることがわかっている状況にしばしば遭遇するため、複素数値関数を扱うものである。
【０２３０】
対数形式とべき乗形式との間を移動する効率的な方法は、単純な加算および減算として、対数空間で乗算および除算を実行することを可能にする。プロセスの性質上、これはべき乗関数に関連する超越関数を容易に生成するためにも、利用される可能性がある。これは単一のアルゴリズムから生成することができるので、１つの要素だけが実装を必要とし、１つの方法だけが実装され最適化される必要があり、開発時間が短縮されるので実装サイズが小さくなる。
【０２３１】
既存のアルゴリズムであるＢＫＭはこのアプローチを実行するが、提案された方法では、ほとんどまたは完全に除去されている３つの不利な要件を含む：
【０２３２】
とりわけ、相でなければデータ入力、出力上のステップの前処理と後処理を実装することを困難とするハードウェアの分割。
【０２３３】
効率的な実施のために、丸めの問題を引き起こす非合理的な価値および要因。
【０２３４】
浮動小数点型と統合するために、必要とされる追加の演算。
【０２３５】
次のセクションは、提案した方法がこれらの欠点がないように、どのように設計されてきたかを説明する。
【０２３６】
ＩＩ．動機と比較範囲の縮小
【０２３７】
プロセスは、オイラーの式
2020502906.tif 000359
に集中しており、
2020502906.tif 000360
は、式 
2020502906.tif 000361
が複素一般化されているため、さらなる係数である。ＢＫＭアルゴリズムが具体化する複雑な機能はまさにこれであり、実行されるプロセスを
2020502906.tif 000362
とする。このアルゴリズムは、収束が保証され得る範囲に、複素数値
2020502906.tif 000363
を当てはめるために範囲縮小技術を頻繁に使用する。これは、
2020502906.tif 000364
を達成するために、提案された方法を使用して、無理数スケーリングファクターで前処理または後処理することによって評価することがおそらく容易である限りにおいて、適切な実装を有する関数ではない。多くの場合、この無理数スケーリングファクターは、入力と出力が物理的または計算上意味のある有理数の倍数になるように、意図するワークフローを提供するために削除されることがある。
【０２３８】
複素べき乗を処理するとき、べき乗は実数部と虚数部に分割され、それらは同じ入力データを使用するが、異なる結果を生成する。これは、ある意味ではデカルト座標から極座標への変換に似ている。
【０２３９】
べき乗化されるべき乗対数の実部だけを考えると、これはデカルト−極変換の
2020502906.tif 000365
成分と同様に、すべての係数情報を含む。これは伝統的に実数部の力に引き上げられた無理数
2020502906.tif 000366
である。それをアルゴリズムに供給するために、これは最初に範囲縮小を受けなければならず、それは剰余を見つけなければならず、それが収束領域内に収まるように、無理数
2020502906.tif 000367
を法で除算ことを意味する。これはハードウェア分割を含み（１９９４年の論文は、ＳＲＴ分割技術を推奨する）、それ故にこれはアルゴリズムと共に分割を実行するのに必要とされる論理の点で、非常に高価な技術を加える。
【０２４０】
提案されたアルゴリズムの前提は、範囲縮小を簡単にするために、２のべき乗を使用すべきであるということである。この動機で、実部は実部である所与の力に引き上げられた２、つまり底２のべき乗として解釈されるべきである。実数対数入力の整数は、出力上のビットシフトを参照することができ、相対２進対数入力の小数部分だけが収束領域に収まる必要がある複素対数の実部の整数部分は、このようにして四捨五入されて別々に累乗される。
【０２４１】
対数の虚数部だけをべき乗化すると考えると、これには、直交座標から極座標への変換の
2020502906.tif 000368
成分と同様に、すべての引数または角度情報が含まれる。これも伝統的に
2020502906.tif 000369
のべき乗とされるが、ここでは
2020502906.tif 000370
は問題ではない。虚数対数部分は角度として解釈することができるため、問題は今回
2020502906.tif 000371
であるが、収束の範囲を制限するために、回転対称性が利用することができるように
2020502906.tif 000372
の整数分数から残りを見つける必要がある。
【０２４２】
対数の底
2020502906.tif 000373
の代わりに、底２を使用することができるだろう。これは、
2020502906.tif 000374
という拡張を与える。これはさらに望ましくない無理数の値であるが、底
2020502906.tif 000375
を特に使用する場合は、回転数を表すために
2020502906.tif 000376
を使用できる。残念なことに、この基底の使用は、収束しない方法を導く。技術と収束範囲の間の最良のトレードオフは、底
2020502906.tif 000377
を得る。これは、各整数値が象限を表すことを意味する。４各倍数は、全体の回転である。次いで、ビットごとの演算を使用して、回転数を抽出することができる。象限の数は、実数と虚数および符号の切り替えとして出力に適用できる。この場合もやはり、複素対数の虚数成分の整数部分が四捨五入によって取られ、このようにして別々に処理されるという技法が得られる。この場合もやはり、複素対数の虚数成分の整数部分が四捨五入によって取られ、このようにして別々に処理されるという技術が得られる。
【０２４３】
これは、図３２の右パネル３２２０に示されるように、対数の小数部分が依然として
2020502906.tif 000378
の範囲内でのべき乗化が残り、これは、主アルゴリズムの収束範囲に収まる。図３２では、領域
2020502906.tif 000379
で対数−べき乗過程へのインプットのためにプロットされたべき乗での誤差、ここでは、右側は正の実数で、画像の下部は正の虚数である。５０の小数ビットと同様の反復回数では、黒い領域は、２−４５を超える精度に収束した値を表す。左パネル３２１０は、灰色が２−４５よりも大きい誤差に、マッピングする収束領域を示す。右パネル３２２０は、白で、最後に必要な収束領域
2020502906.tif 000380
を示す。
【０２４４】
対数入力の、小数部の主アルゴリズムについては後述する。
【０２４５】
逆の処理も達成され、そのべき乗化された実数部および虚数部を用いてその値をとり、それをその元の対数に減らす。はじめに、複素数値が属する象限を、４５°オフセットして見つける。これは、どの成分が最大の絶対値；実数または虚数およびその記号を持つかを、決定することによって達成する。これは、これらの象限のどれが、値が属する角度でオフセットしているかを特定する。次に、これを使用して、図３３の中央パネル３３２０に示すように、その面積が実線の正の部分を中心とする、第１象限に値を移動するための回転を計算するために使用する。図３３は、領域内
2020502906.tif 000381
のべき乗対数プロセスへの入力に対してプロットされたべき乗の誤差を示し、ここでは、右側は正の実数で、画像の下部は正の虚数である。５０の小数ビットと同様の反復回数では、黒い領域は２−４５を超える精度に収束した値を表す。左パネル３３１０は収束領域を示し、ここでグレーは２−４５より大きい誤差にマッピングされる。真ん中のパネル３３２０は、反転グレースケールを用いて実線を中心とする象限を示す。右パネル３３３０は、白で、最後に必要な収束領域
2020502906.tif 000382
を示す。
【０２４６】
対数出力の虚数部には、出力複素対数の等価虚数整数部があらかじめロードされている。４５°オフセットによって、複雑なブランチカットが実数ラインの負の部分に入るには、さらにテストが必要だが、２の補数システムラップアラウンドが使用されている場合、これは無視でき、必要なビットが与えられると、象限で測定された対数の虚数部は、
2020502906.tif 000383
の区間内にあるように、簡単に制限できる。
【０２４７】
この時点で、値が第１象限にあることは確かであるので、実数部は絶対値で虚数部と同等またはそれ以上であることが周知であり、そして実数部は正であることが周知である。対数の実部の整数部を見つけるために、固定小数点表現の先頭のゼロを数えることができる。固定小数点表現のビット数からこのカウントを減算した値は、事実上、
2020502906.tif 000384
であり、ここでｋはビットの表現に依存し、バイナリ対数の整数部分の表現でもあり、これは、対数レジスタの実数部分にすぐに加算できる。図３３の右側のパネルに示すように、実データと虚データの両方の既存データは、実部を区間
2020502906.tif 000385
に収めるために必要なビット位置の数だけシフトアップでき、これは、初期カウントの先行ゼロと、ビットの表現に依存するｋによって決まる。これは、出力実数対数の最も近い整数表現ではないかもしれないが、整数であり、プロセスの残りの部分を正しいフォーマットで収束領域に入れることに注意されたい。
【０２４８】
このプロセスの終わりに、出力複素数対数の整数部分が決定され、出力の小数部分と入力範囲
2020502906.tif 000386
は未決定のままであり、これは、主アルゴリズムの収束範囲内に収まる。対数出力の、実質的に小数部分のための主なアルゴリズムの議論は、以下で説明される。
【０２４９】
この結果は、関数
2020502906.tif 000387
が、図３４から分かるように、モデル化されたプロセスとなる。図３４では、べき乗部分３４１０および対数プロセス３４３０の領域のべき乗関数の範囲内で、べき乗部分３４２０および対数３４４０のプロセスの最終的な合成値を構成する対応部分を範囲とともに、プロットされた誤差を十分に小数部分についてのみプロットした。ノイズの多い灰色は、正しい収束が小さな数値誤差によって誤差測定値を支配するような領域を示しており、べき乗関数３４２０と対数３４４０について収束が至る所で達成されていることを示している。
【０２５０】
結果の底を変更しても、対数とべき乗のすべての特性があることを意味し、つまり、対数の実数部を足したり引いたりして、べき乗の乗算と除算を行うことができる。同様に、角度は同じように振る舞うが、それらは象限の整数倍、つまり回転に対して４、ハードウェアロジックで簡単に操作される値として表されるということだけである。対数形式では、虚数値の加減算は、複素数のべき乗化空間内でゼロを中心に、時計回りおよび反時計回りに回転させるのと同じである。実数部と虚数部が、異なる底で構成されているという事実に関係なく、実数部と虚数部を組み合わせた複素数値の場合、複素数の乗算と除算の実装に使用できる。
【０２５１】
これらの結果は、必要とされる前処理が、仮数およびべき乗形式によってすでにほとんど達成されるため、浮動小数点表現に自明に拡張することができる。底
2020502906.tif 000388
の値が絶対的に必要とされ、ほとんどの場合それらが必要とされない場合、出力対数の実部は
2020502906.tif 000389
で虚数は
2020502906.tif 000390
で乗算される。同様に、対数データを底
2020502906.tif 000391
とするフォーマットで取り込むためには、値を逆数
2020502906.tif 000392
とし、変換するために成分ごとに乗算してから、この方法でべき乗化するための初期範囲縮小を行うことができる。
【０２５２】
ＩＩＩ．一般的な手順
【０２５３】
新しいアルゴリズムは、対数の小数部分を見つけるか使用するために、ＢＫＭ技術と同様の方法で進行する。
2020502906.tif 000393
の値の集合があり、そこから各アルゴリズムステップ
2020502906.tif 000394
において
2020502906.tif 000395
が選択される。どの
2020502906.tif 000396
を使用するかの選択は、切り捨てられた数を、入力レジスタからのビットのサブセットであり得る、セットビット数と比較することによって行われる。
2020502906.tif 000397
が選択されると、
2020502906.tif 000398
（ビットシフト（潜在的に最上位の未使用ビットに丸めが適用される）と加算／減算）のバランスをとるために、事前計算対数のルックアップテーブルが必要であり、対数値を保持しているレジスタから減算する等価の対数を使用して、べき乗値を保持しているレジスタに乗算を適用する。プロセスが、対数を見つけることであるかべき乗を計算することであるかに応じて、システムを収束構成に配置するために、追加のシフト、追加の計算、またはその両方があり得る。いくつかのＮ回の繰り返しの終わりに、べき乗または対数のいずれかの所望の値は、繰り返しの回数Ｎおよび表現に使用されるビット数に依存する精度の範囲内であることが見出される。
【０２５４】
ＩＶ．対数のルックアップテーブル
【０２５５】
提案された方法は、これらの対数がここで、
2020502906.tif 000399
べき乗を反映しなければならないという点で既存の技術とは異なる。これにより、対数表の実部には、
2020502906.tif 000400
と、虚部
2020502906.tif 000401
または
2020502906.tif 000402
とをプリロードする必要があり、
2020502906.tif 000403
について所望の固定小数点フォーマットである。
【０２５６】
対数の底として、２と、
2020502906.tif 000404
とを選択したが、底の選択を変えると、対数システムのアルゴリズムの最初の反復が異なることはないと思われる。これらの値は、同じハードウェアを使用して、両方のプロセスを達成できる設計を表すために選択された。対数モードにべき乗モードとは異なる異なるハードウェアを使用した別の実装を使用する場合は、これが基本の最良の選択ではない可能性がある。
【０２５７】
Ｖ．対数の小数部の方法
【０２５８】
新しいアルゴリズムのための中心的な方法は、いくつかの新しい重要な相違点を除いて、ＢＫＭプロセスによって使用されるものに大部分従う。対数とべき乗の新しい底と、新しい手法に必要な収束領域の変化に対応するために、比較しきい値が新しく導き出された。必要な領域内での収束を達成するために、対数プロセスの最初の反復のために、わずかに修正された比較しきい値もある。
【０２５９】
ＶＩ．べき乗モード
【０２６０】
このセクションに記載されているべき乗モードはまた、べき乗値との複素乗算を提供するように修正される。これが達成することになっている場合、出力が正しいことであるならば、これは範囲縮小ステップの前に予めロードされなければならない。これは出力のべき乗値に置き換わるため、この値が必要な場合は使用してはならない。対数の整数部分（値ｚinteger、output）からプロセスの最後まで解を保存して適用するのを待つことも好ましい。これは、処理が行われる中間レジスタに必要な記憶容量を減らすことができるが、これは解の整数部分に必要な追加の記憶容量要件と比較して慮られることである。
【０２６１】
あるいは、追加のレジスタを作成し、これらの追加のレジスタで同等の演算が行われるようにすることによって、最終のべき乗値との乗算を並行して達成することができる。この方法では、べき乗プロセスは、出力べき乗値をほぼ任意の多くの他の複素数値と、複素数乗算することができる。
【０２６２】
入力対数の小数部分を入力と仮定すると、収束領域のアルゴリズム
2020502906.tif 000405
は次のようになる：
【０２６３】
１．
2020502906.tif 000406
とラベル付けられた、４つの基本レジスタがあると仮定する。並んで、補助複素数乗算に使用されたときの方法の動作を示すために、２つの追加のスレーブ乗算レジスタ
2020502906.tif 000407
と
2020502906.tif 000408
とがある。これらのレジスタの初期値は以下の通りである：
2020502906.tif 000409
ここで、事前乗算の要件がない場合は、
2020502906.tif 000410
である。スレーブ乗算レジスタも同様に以下のもので構築することができる：
2020502906.tif 000411
【０２６４】
２．インデックス
2020502906.tif 000412
として値
2020502906.tif 000413
を繰り返す：
【０２６５】
３.
2020502906.tif 000414
を切り捨てて、３ビットになるように
2020502906.tif 000415
を作成する；１の符号ビットと２の補数の２つの整数ビットは、範囲が
2020502906.tif 000416
で最小の変化は１である。
【０２６６】
４．
2020502906.tif 000417
を切り捨てて、
2020502906.tif 000418
を形成するが、これは、１符号ビット、１整数ビット、２の補数で１小数ビットの３ビットを有し、範囲は
2020502906.tif 000419
で、最小の変化は０．５である。 
【０２６７】
５.
2020502906.tif 000420
を決定するために、３ビット値をテストする：
2020502906.tif 000421
【０２６８】
６．シフトと算術演算を使用して、べき乗レジスタに乗算を適用する:
2020502906.tif 000422
2020502906.tif 000423
および：
2020502906.tif 000424
2020502906.tif 000425
乗算処理をこれらにも適用するには、
2020502906.tif 000426
や
2020502906.tif 000427
などの補助レジスタにも同じ操作を行う。
【０２６９】
７．対数表の対応するエントリをレジスタから減算する：
2020502906.tif 000428
2020502906.tif 000429
【０２７０】
８．対数レジスタを１桁左にシフトして、その内容の値を２倍にする。
2020502906.tif 000430
【０２７１】
９．Ｎに達するまで、次の反復のためにステップ２に戻り、Ｎに達すると、その時点でレジスタはそれらの最終値を含む
2020502906.tif 000431
【０２７２】
および
2020502906.tif 000432
【０２７３】
プロセスの形式を理解したので、
2020502906.tif 000433
またはその両方を形成することによって、時には必要な領域で収束する他のテスト手順を見つけることは簡単であるが、我々は、異なるビット数または異なる比較値を使用するもので必要な値テストを最も簡単な周知の形式で指定することで、複雑さを軽減するように努めている。
【０２７４】
ＶＩＩ．対数モード
【０２７５】
このセクションに記載されている対数モードはまた、入力値との複素除算を提供するように修正されてもよい。これが達成することになっている場合、追加のレジスタを作成し、これらの追加のレジスタで同等の操作が行われるようにすることによってのみ、並行して発生する可能性がある。これは、元のレジスタが過負荷になる可能性がある補助乗算とは対照的であり、このモードではべき乗レジスタを変更するとアルゴリズムの収束が妨げられるため、ここでは達成できない。しかしながら、補助レジスタを使用することでこれを回避できるため、必要に応じて、このプロセスへの入力値を分母として、対数プロセスで他のほぼ任意の多数の複素数値を複素数除算することができる。
【０２７６】
１．出力対数の小数部を出力とすると、収束領域
2020502906.tif 000434
のアルゴリズムである：
【０２７７】
２．
2020502906.tif 000435
および
2020502906.tif 000436
とラベル付けされた、４つの基本レジスタがあると仮定する。並んで、補助複素除算に使用されたときの方法の動作を示すために、２つの追加のスレーブ除算レジスタ
2020502906.tif 000437
が存在する。これらのレジスタの初期値は以下のとおりである：
2020502906.tif 000438
スレーブ除算レジスタも、同様に構成される：
2020502906.tif 000439
レジスタ
2020502906.tif 000440
のいずれにも２倍算も−１．０も適用されないことに留意されたい。
【０２７８】
２．インデックス
2020502906.tif 000441
としての値
2020502906.tif 000442
を反復する：
【０２７９】
３．
2020502906.tif 000443
を切り捨てて
2020502906.tif 000444
を５ビット、１符号ビット、１整数ビット、２の補数の３小数ビットとして形成し、範囲は
2020502906.tif 000445
で、最小の変化は０．１２５である。
【０２８０】
４．
2020502906.tif 000446
を切り捨てて、
2020502906.tif 000447
１の符号ビット、２つの整数ビット、２の補数の２つの小数ビットの５ビットになるように形成し；範囲は
2020502906.tif 000448
で、最小の変化は０．２５である。
【０２８１】
５．２つの５ビット値をテストし（最初の反復の比較を修正して）、
2020502906.tif 000449
を決定する：
2020502906.tif 000450
【０２８２】
６．シフトと算術演算を使用して、べき乗レジスタに乗算を適用する：
2020502906.tif 000451
2020502906.tif 000452
および：
2020502906.tif 000453
2020502906.tif 000454
【０２８３】
2020502906.tif 000455
などの補助レジスタにも同じようにして、これらに除算プロセスを適用する。
【０２８４】
７.
2020502906.tif 000456
を、
2020502906.tif 000457
に加算する：
2020502906.tif 000458
これは、
2020502906.tif 000459
に追加しない。
【０２８５】
８．べき乗レジスタを１桁左にシフトして、その内容の値を２倍にする。
2020502906.tif 000460
これは、
2020502906.tif 000461
の値を２倍にしない。
【０２８６】
９．対数表の対応するエントリをレジスタから減算する：
2020502906.tif 000462
【０２８７】
次の反復のため、
2020502906.tif 000463
に達するまでステップ２に戻り、このポイントで、でレジスタはそれらの最終値を含む：
2020502906.tif 000464
【０２８８】
および：
2020502906.tif 000465
【０２８９】
プロセスの形式を理解されるように、ある場合には、
2020502906.tif 000466
またはその両方を形成することによって、または異なるビット数または異なる比較値を使用して必要な領域で収束する他のテスト手順を見つけることができるが、我々は、必要な値テストを最も簡単な既知の形式で指定することで、複雑さを軽減するように努めている。
【０２９０】
ＶＩＩＩ．冗長数システム
【０２９１】
ＢＫＭ法のようなアルゴリズムは、データのローカルフローを遅くする技術をサポートするのに理想的であり、ハードウェアで実施される高速ロジックへの統合を可能にする。各反復での比較のために提示されたビット数の削減と相まって、算術ステップにキャリーセーブのような冗長数システム方式を使用することは、ここに提示されるようなアルゴリズムが高速低面積実装を有することを可能にし、これは、集積回路設計への組み込みに最適である。このような実装の利点には、より高速のクロック速度を達成できること、多数のビットを持つオペランドにも対応できることなどがあり、この設計は複雑な値の操作を含む問題に対する低コストの最新ソリューションの作成に理想的である。
【０２９２】
直観性を容易にするために、丸めはアルゴリズムの前述の説明では無視されており、本発明では正しい丸めが考慮されていないという意味として解釈されるべきではない。各ステップで適用されるときに丸めの適用を通して、エラーを減少させた結果を達成するために、いかなる算術演算でも直接使用されない最上位ビットを利用することは、ここで説明される結果の精度を高める可能性がある。
【０２９３】
ＩＸ．高次底のアプローチ
【０２９４】
ビット単位で評価する他のアルゴリズムと同様の方法で、複数のビットを一度に計算する方法を作成するために、ステップを併合することができる。出力は、２とは異なる底を持つものとして扱われるため、これは高次底アルゴリズムと呼ばれ、つまり、新しい高次底からの１桁がステップごとに計算される。この欠点は、ルックアップテーブルとシフトの選択もマージしなければならず、反復ごとにアルゴリズムが機能するための選択の、より大きなセットを作成することである。
【０２９５】
Ｘ．実数アルゴリズム
【０２９６】
多くのシステムは、本明細書に記載のアルゴリズムの実部のみを使用することから恩恵を受けることができる。例えば、ニューラルネットワークアルゴリズムによって、しばしば使用されるシグモイドロジスティック関数は、最初に実数２ｘを決定することによって、またはべき乗化処理を使用して、予め乗算ｅｘされた場合に生成することができる。次に、１を加算してプリセット定数をこの値で除算することが可能であるか、または、元の論文で説明されているＥモードに続いてＬモードを使用し続けると、ロジスティック関数は次のようになる：
2020502906.tif 000467
【０２９７】
ただし、分子ｋとその他の修正子ａおよびｂは、設定可能なままである。
【０２９８】
ＸＩ．収束性の提示
【０２９９】
図３５では、べき乗法および対数法の収束は、グラフ３５００によって示される。図３５は、必然的に収束するドメインにおいて可能な各値をテストすることが、ここでグラフ化された最悪の場合の
2020502906.tif 000468
エラーをもたらすことを示す。このグラフは、１兆を超える複雑なテスト値（再マッピングされた分数ドメイン内の各実数および虚数の値が、表示されている反復回数についてテストされている）の最悪の場合の動作をまとめたものである。グラフを比較可能にするために、対数の結果は、エラーをテストする前に、べき乗形式に変換した。これは、この方法がこの多数の有効ビットに対して収束的であることを保証する。この誤差は、乗算に丸めを追加することによってさらに減らすことができる。グラフは、精度と反復が増加するにつれて、誤差も同様に減少し、収束を示していることを示している。
【０３００】
このグラフ３５００は、表内のビット数およびプロセス内の反復回数と同様の割合で縮小する収束領域内のすべての可能性のあるケースからの最悪ケースの絶対誤差を示す。これは、この固定小数点データのビット数に対する全ての可能な正しい入力に対して、アルゴリズムが正確なビット数で、正しい解に線形に首尾よく収束することを明らかにしている。
【０３０１】
ＸＩＩ．可能な請求項
１．入力された「対数」複素数を得られる「べき乗」複素数に同時的に変換する方法であって、
ａ  前記得られる複素数の係数が前記実入力された複素数部分の実部べき乗に対する二乗からなること；および
ｂ  得られる複素数の引数が、前記虚数入力された複素数部分の実値で乗算された整数べき乗ロテーションの二乗を構成すること；
を含む、方法。
２．  入力された「べき乗」複素数を得られる「対数」複素数に同時的に変換する方法であって、
ａ  入力された複素数の実係数の実バイナリ対数が、得られる複素数の実部を形成すること；および
ｂ  前記入力された複素数の引数において表現された整数べき乗ロテーションの二乗の実数が、得られる複素数の虚部を形成すること；
を含む、方法。
３．  前記変換は、前記複素数のバイナリ表現についてのビットシフト、比較および計算を通して実質的に達成される、請求項１または２に記載の方法。
４．  入力および出力範囲の削減を、実質的に前記入力複素数の成分を別の部分とは別にして丸めることに対応する値の処理により実行する、請求項１に記載の方法。
５．  前記入力および出力範囲の削減を、実質的に前記出力複素数の成分を別の部分とは別にして丸めることに対応する値の処理により実行する、請求項２に記載の方法。
６．  前記出力複素値を、１に等しくない値に初期化して、前記反復の適用が非１の値を修正して本来出力することを意図する値の複素乗算を計算する、請求項１に記載の方法。
７．  さらに、複素値を生成し、前記出力複素値に適用される操作がまた、これら更なる値に適用されて、これらのさらなる複素値に対する出力することを本来意図した前記複素値での乗算を与える、請求項１に記載の方法。
８．  さらなる複素値を生成し、前記入力複素値に適用される操作がまた、これら更なる値に適用されて、これらのさらなる複素値の前記複素値入力による除算を与える、請求項２に記載の方法。
９．  冗長数システムを使用して、前記残された複素値に対する計算を反復する処理を実装する、請求項１に記載の方法。
１０．  冗長数システムを使用して、前記残された複素値に対する計算を反復する処理を実装する、請求項２に記載の方法。
１１．  請求項１〜１０のいずれかの方法を実行するために設計されたデバイス。
１２．  請求項１〜１０のいずれかの方法を実行するために構成されたデータプロセッサ。
【０３０２】
７．二次問題のＫＲＹＬＯＶ副空間から構築された音響場
【０３０３】
本開示は、固有システムとして表現されたときに二次問題を操作して、線形システムへの包含のための、ランク低減基底集合を生成する方法を説明する。縮小ランク基底を使用することにより、（潜在的に低減された）基底関数の最適情報保存集合を、複素値線形システムに含めるために自動的に導き出して、そのような基底関数の最適複素値再重み付けを見つけて、いくつかのテスト関数で説明されているように、目的を達成する。
【０３０４】
基底関数の縮小されたセットの導出、およびそれらが音場線形最適化問題の実装の適切な段階に、どのように含まれるかについて説明する。
【０３０５】
Ｉ．線形コントロール領域
【０３０６】
線形制御領域、これは複素加重重み付け関数、
2020502906.tif 000469
を用いて、圧力などの線形音響量との積分を実行する演算子を定義する:
2020502906.tif 000470
これは、音響場を評価するために使用する。
2020502906.tif 000471
が音響量の類似の位相分布と対になっている場合、複素共役は関数評価を最大化するために使用されることに留意されたい。この積分に対する任意の近似値（サンプリングなど）を使用することができ、あるいは実際には、これを使用して音圧場にわたって、任意の線形演算子を近似することができる。さらに線形であれば、ベクトル粒子速度（各ベクトル成分は場の調和運動により、複素数値である）を使用して、コントロール領域を最適化するようにさらなる積分を構成することができ、ここで：
2020502906.tif 000472
である。ここで、
2020502906.tif 000473
は媒体の粒子速度、
2020502906.tif 000474
は、各次元が複素数である３次元ベクトル場であり、これは、上記と同じ目的で重み付け関数として機能する。この方法では、異なる音響量パターンが線形的に評価されてもよく、ここで位相が
2020502906.tif 000475
におけるスカラーまたはベクトル関数によって生成された「テンプレート」と整列するとき、結果の特定の位相は、目的関数を最大化する。
【０３０７】
また、
2020502906.tif 000476
は、複数の外見上ばらばらな領域の単一の評価を作成するために、大部分がゼロ（またはゼロベクトルからなる）であり得ることにも留意されたい。この定式化は、線形音響量のフィールドにわたって、任意の線形演算子に拡張することができるので、任意の整数近似（サンプリングなど）を使用することができる。
【０３０８】
ＩＩ．二次問題と固有値問題
【０３０９】
位相アレイは、各トランスデューサをコントロールする、一連の複素数値の係数によって記述することができるマトリクス
2020502906.tif 000477
内に符号化された音響量を最大にするために、そのような入力
2020502906.tif 000478
ベクトルを用いて固定周波数で音響変換器のアレイを駆動することは、以下の意図をもった線形代数問題として表すことができる。
2020502906.tif 000479
【０３１０】
これは、固有ベクトル、ここでは
2020502906.tif 000480
の定義でもあるので、これはマトリクス
2020502906.tif 000481
の主要な固有ベクトルをとることによって、解くことができる。マトリクス
2020502906.tif 000482
を構築することは、まず、問題の音響量の大域空間における基底関数として、各トランスデューサによって放出された、線形音響量を考慮することによって達成することができる。これは、絶対位置
2020502906.tif 000483
によって記述する、単一ポイントで、圧力および他の音響量の最大化を達成するための方法である。しかし、記載されたコントロール領域積分の使用により、各個々の変換素子の場に放出される音響量を支配する
2020502906.tif 000484
は、
2020502906.tif 000485
と交換され、列ベクトルをもたらす：
2020502906.tif 000486
【０３１１】
線形性を呼び出すことによって、各トランスデューサを駆動する複素活性化値を乗算して、ベクトルを導き出すことができる：
2020502906.tif 000487
ここで、ｎはトランデューサの数である。このｍが線形場に作用する線形演算子の結果であると仮定すると、それはエルミート転置され、それ自身によって乗算されて二次形式を生成することができる。最後に下記式が得られる。
2020502906.tif 000488
【０３１２】
重み付け関数
2020502906.tif 000489
を修正し、マトリクス
2020502906.tif 000490
の固有ベクトルを解くことによって、様々な異なるビーム、表面、および他の音響現象が生成され最大化され、音場内に最適に作動される領域を作り出すことができる。
【０３１３】
前述のヌル空間生成技術（「空中ハプティックシステムにおける二次問題のための固有システム解」）もまた、音響浮上を記述するためにこの方法を拡張するために、ならびにヌルステアリングおよびミッド形成を可能にするために使用され、構造を記述するための位置で、フィールドにヌルポイントを導入することによって、ヌルステアリングおよび空中テクスチャの形成を可能にする。このプロセスはマトリクスのランクを増加するため、簡素化は使用されることが不可能であることに留意されたい。
【０３１４】
同様に、複数の最適化は、いくつかのコントロール領域について、
2020502906.tif 000491
個のマトリクスを互いに加算することによって、実行することができる。これは、コントロールポイントを表すマトリクスも含み得る。これは、増加したランクを持つマトリクスを生成することにより、簡素化の方法は使用するケースでは、不可能である。
【０３１５】
ＩＩＩ．高ランクマトリクスに対するランク低減近似
【０３１６】
キーポイントは、ランクの低い近似値が生成されるように、ランクの高いマトリクスを意図的に作成することができることである。例えば、独立したポイントを何百も持つシステムは、独立した基底関数で記述されるものが非常に少なく、忠実度を高め、複雑さを減らし、基底関数を必要に応じて自動的に作成することを可能にするよりも、高次元であるように見える基底セットの実現を作成するように、自動的に再編成される。
【０３１７】
既存のアプローチは、各コントロールポイント、またはコントロール領域のために、基底関数を生成する。各コントロールポイントまたは領域を最適化するための固有システムを作成し、それらを合計するか、またはヌルを追加することによって、固有システムはランクを獲得する。基底関数マトリクスおよび／またはヌルを追加することを繰り返されるプロセスを通して、高度なランクが最終マトリクスによって得られることが可能である。線形システムは、リソースに対する制約のために、最大で
2020502906.tif 000492
個のテスト関数を使用して、最大で
2020502906.tif 000493
個の基底関数を処理することしかできないかもしれない。この最適化マトリクスの高いランクはこの場合
2020502906.tif 000494
であり、ここで、個々の基底関数を生成するための
2020502906.tif 000495
個の最適化マトリクスは、追加で併合された。このアプローチを使用して、
2020502906.tif 000496
および／または
2020502906.tif 000497
であるという他の不可能な状況は、有限処理能力を有するこのシステムで解決可能な
2020502906.tif 000498
個の基底関数を有する、縮小基底集合で表される問題に帰着される。
【０３１８】
この縮小基底集合は、基底関数に対する複数の最適化からなる最終マトリクスの
2020502906.tif 000499
個の最も主要な固有ベクトルと、実質的に同じマトリクス部分空間を占める、クリロフ部分空間を計算することによって、見出される。クリロフ部分空間は、主要な固有ベクトルが見つけられ、次いで残りのマトリクスから取り除かれるアーノルディ反復法のような、標準的な技術を用いて見つけられ、それから
2020502906.tif 000500
個のそのような固有ベクトルが見つかるまでこれが繰り返される。
【０３１９】
ＩＶ．混合アプローチ
【０３２０】
基底関数が上位マトリクスから抽出されるこのアプローチは、既存の特定の基底関数アプローチと混合して適合させることができる。例えば、２つのコントロール領域またはポイントを必須項目として指定し、線形システムへの基本セット入力に追加することができる。それから、それらが二次最適化、したがってそれ自身の最適化マトリクスによって特定されるということを考えると、残りの基底関数は、単一のマトリクスに加えることができ、そのマトリクスは、線形システムに入力された基底関数のセットの残りの数を構成する、最も支配的な
2020502906.tif 000501
の固有ベクトルを見つけるために、このクリロフ部分空間法を用いて分解される。この方法で、最初に指定された２つの基底関数は、線形システムに直接インポートされるため保証されているが、セットの残りの関数は近似する。このシステムはさらに用いられ、ここで、底関数の集合数は、上位ランクの最適化マトリクスによって、それぞれ記述される基底関数の互いに素な集合から特定され、そして基底関数の集合数が取り出される。縮小基底関数を導出するために使用できる、さらなる技術は、パワーしきい値を超えるまで、マトリクスから基底関数を取り続けることである。このしきい値は、結合された基底関数から構成されるマトリクスごとに異なる可能性があるため、異なるパワーの固有ベクトルを生成し、この場合、固有ベクトルを選択することで固有ベクトルを生成でき、各固有ベクトルのパワーは、実際の優先順位の値と対応する固有値の大きさの積に基づいて、各マトリクスから計算される。
【０３２１】
８．テスト関数を使用した線形システムソリューション
【０３２２】
基底関数セットが基底関数セットとは別に定義される、特定の目標を達成するかどうかを「テストする」方法は、価値があるだろう。これらの目標が達成されたかどうかの定量的評価は、最適化に提供され、これは、これらのテスト関数によって定義される音場内の所望の特性を生成するために使用される。
【０３２３】
テスト関数および複素数値線形システム最適化における、それらの役割について説明する。
【０３２４】
Ｉ．線形コントロール領域 
【０３２５】
線形コントロール領域、これは、複素数値重み付け関数
2020502906.tif 000502
を用いて、圧力などの線形音響量との積分を実行する演算子として定義する：
2020502906.tif 000503
これは、音響場で評価するために使われる。
2020502906.tif 000504
が、音響量の類似の位相分布と対になっている場合、複素共役は、関数評価を最大化するために使用されることに留意されたい。この積分に対する任意の近似値（サンプリングなど）を使用することができ、または実際には、これを使用して音圧場にわたって任意の線形演算子を近似することができる。ベクトル粒子速度（ここで各ベクトル成分は、場の調和運動のための複素数値である）を用いて、コントロール領域を最適化するように、さらに積分されてもよい、ここで：
2020502906.tif 000505
2020502906.tif 000506
は、媒体の粒子速度、
2020502906.tif 000507
は、それぞれの次元が複素数値の３次元ベクトル場であり上記と同じ目的で重み付け関数として機能する。この方法では、異なる音響量パターンが線形的に評価されてもよく、ここで位相が
2020502906.tif 000508
におけるスカラーまたはベクトル関数によって生成された「テンプレート」と整列するとき、結果の特定の位相は目的関数を最大化する。
【０３２６】
また、
2020502906.tif 000509
は、複数の外見上ばらばらな領域の単一の評価を作成するために、大部分がゼロ（またはゼロベクトルからなる）であり得ることにも留意されたい。この定式化は、線形音響量のフィールドにわたって、任意の線形演算子に拡張することができるので、任意の整数近似（サンプリングなど）を使用することができる。
【０３２７】
ＩＩ．基底関数セット
【０３２８】
入力として
2020502906.tif 000510
個の基底関数のセットが与えられると、
2020502906.tif 000511
個の基底関数の最適な線形結合を見つけるために、線形システムを通して
2020502906.tif 000512
個のテスト関数が使用され、
2020502906.tif 000513
個のテスト関数の所望の客観的評価を達成する。これら
2020502906.tif 000514
個のテスト関数は、線形コントロール領域評価の形をとることができる。上記の目的
2020502906.tif 000515
の形式を使用すると、次のようになる：
2020502906.tif 000516
または、単一のコントロールポイント定義を通して、次のように評価することもできる：
2020502906.tif 000517
ここで、上線は複素共役を表し、
2020502906.tif 000518
は、三次元位置を表し、
2020502906.tif 000519
は全てのトランスデューサ素子をパラメータ化し、
2020502906.tif 000520
はトランスデューサの圧力場の共役を乗じた、トランスデューサ活性化である。法線方向が定義されている場合、粒子速度によって定義されるコントロールポイントの等価形式も、評価子として使用できる：
2020502906.tif 000521
ただし、間違った方向を向いているトランスデューサ
2020502906.tif 000522
も、
2020502906.tif 000523
の場合、その寄与のセットをゼロに設定し、省略することができる）。
【０３２９】
両方の単一コントロールポイント定義において、既存のアプローチは、コントロール領域評価の単位重み付け関数の場合を伴う、特別な単一ポイント積分評価と見なすことができる。ここでの目的関数は、余分な副作用のない単純な音響量である。
【０３３０】
ＩＩＩ．主要な固有ベクトルの発見
【０３３１】
2020502906.tif 000524
個のテスト関数の各々を用いて、
2020502906.tif 000525
個の基底関数の各々を評価することにより、複素数値である矩形マトリクスが得られる。前述のように、主要な固有ベクトルは、
2020502906.tif 000526
個のテスト関数によって定義される空間を最も良く張る
2020502906.tif 000527
個の基底関数の組み合わせを表す。しかし、これは必ずしも正方マトリクスではないため、特異値分解がより適切であるが、唯一必要な固有ベクトルは支配的なものである。その支配的な固有ベクトルが、システム基底関数の支配的な組み合わせの要素を表す、正方エルミート行列を、構築することができる。固有システムは、ここで、以下のように記載される：
2020502906.tif 000528
ここで、
2020502906.tif 000529
は、位相が
2020502906.tif 000530
個のテスト関数のパラメータ空間のカバーを生成する係数を記載する。したがって、優勢固有ベクトルは、単一の組み合わせベクトルで達成可能な空間の最良のカバーを表す成分
2020502906.tif 000531
を有し、その位相はこのカバーを最適化する。
【０３３２】
テスト関数の数が基底関数の数と等しいため、
2020502906.tif 000532
である場合でも、既存の固有システムアプローチは必ずしも適用可能ではない場合があることに、留意されたい。これは、マトリクスがエルミーションではなくなる可能性があるため、上記の手順を完了せずに、レイリー−リッツ（Ｒａｙｌｅｉｇｈ−Ｒｉｔｚ）法をこのマトリクスに直接適用すると、方法が主要の固有ベクトルに収束できなくなる可能性があるためである。
【０３３３】
位相データが指定されていないか、または不完全である場合であっても、主要な固有ベクトルを位相情報のために使用すべきである。固有ベクトルは、各反復の終わりに、振幅データを除去するように修正され、位相がゼロに平均化され、前の解から別の時間ステップで位相進行が制限され、または、既存する定常の位相の関係を有する、いくつかの基底関数に対して１組の位相解を有し、このセットの補数は、それら自体と既存の位相の定常セットとの間の位相で最適化されている。この方法では、この矩形マトリクスシステムは、基底集合を生成する関数が、基底集合の目的を評価するために使用される関数と同じ関数である条件について、特化した既存の固有システムアプローチと同様の効果を達成するために使用される。
【０３３４】
ＩＶ．線形システムソリューション
【０３３５】
複素数値線形システムは、修正された役割、すなわちテスト関数の所与の値を最もよく再現する、基底関数の集合を見つけるという役割を持つ。ここで、線形システムは、下記の矩形マトリックスの形態を取る：
2020502906.tif 000533
【０３３６】
マトリクスが矩形であるため、テスト関数の数が基底関数の数よりも多く、その結果、複素数値の最小二乗システムになるか、またはテスト関数の数が基底関数の数よりも少なくなる可能性があり、複素数値の最小ノルムシステムが得られる。いずれの場合も、標準的な分解を使用して解ベクトルを見つけることができ、その解は、解ベクトルの各成分が、試験関数の所与の値を表現する場合に音響場を生じさせる基底関数のそれぞれに適用される係数である。
【０３３７】
これらのテスト関数は、あるポイントにおける音圧または方向性粒子速度、あるいは空間内のある領域にわたる圧力、および方向性粒子速度の混合を表すことができる。これらの関数のどのような組み合わせも
2020502906.tif 000534
個のテスト関数を構成することができ、それは線形システムの最終ステップで解かれる。
【０３３８】
基底関数の係数が分かれば、トランスデューサを作動させるための最終値は、次のように拡張される：
2020502906.tif 000535
【０３３９】
単一の単調搬送波周波数（状況の変化により多かれ少なかれ有効であり得る）の仮定を与える入力信号は、この
2020502906.tif 000536
の成分を各変換器に適用することによって構成される。
【０３４０】
９．結論
【０３４１】
前述の実施形態の様々な特徴は、改良した触角のシステムの多数のバリエーションを生産するために選択され、そして組み合わされてもよい。
【０３４２】
前述の明細書において、特定の実施形態が説明された。しかしながら、当業者であれば、特許請求の範囲に記載されている本発明の範囲から逸脱することなく、様々な修正および変更が可能であることを理解するであろう。従って、明細書および図面は限定的な意味ではなく、例示的な意味で見なされるべきであり、そのような修正はすべて本教示の範囲内に含まれることが意図される。
【０３４３】
利益、利点、問題に対する解決策、および何らかの利益、利点、または解決策を生じさせるか、またはより顕著にする可能性があるあらゆる要素は、いずれかもしくは全ての請求項の重要な、必要の、本質的な特徴、要素として解釈されるべきではない。本発明は、本願の係属中に行われたいかなる修正、および発行された特許請求の範囲のすべての均等物を含む、添付の特許請求の範囲によってのみ定義される。
【０３４４】
さらに、本明細書では、第１および第２、上下などの関係用語は、必ずしもそのような実態、または動作間の実際のそのような関係または順序を必要とせず、暗示することなく、ある実態または動作を、別の実体または動作から区別するためにのみ、使用される。用語「構成する」、「構成している」、「有する」、「有している」、「含む」、「含んでいる」、「包含する」、「包含している」またはそれらの任意の他のバリエーションは、非排他的包含に及ぶものであることを意図しており、要素のリストを構成する、有する、含む、包含するプロセス、方法、物品、または装置は、それらの要素だけを含むのではなく、そのようなプロセス、方法、物品、または装置に明示的に挙げられないか、または明示的にリストされた固有の他の要素を含む。「構成する」、「有する」、「含む」、「包含する」などの要素は、プロセス、方法、記事、または要素を構成する、有する、含む、包含する装置の中に、追加の同一の要素が存在することを排除するものではない。「ａ」および「ａｎ」という用語は、本明細書で、他に明示的に述べられていない限り、１つまたはそれ以上として定義される。用語「実質的に」、「本質的に」、「およそ」、「約」またはそれらの任意の他の類型は、当業者によって理解されるものに近いとものとして定義される。本明細書で使用する「結合」という用語は、必ずしも直接ではなく、必ずしも機械的にではないが、接続されていると定義される。特定の方法で「構成」されているデバイスまたは構造は、少なくともそのように構成されているが、記載されていない方法で構成することもできる。
【０３４５】
開示の要約は、読者が技術的開示の性質を迅速に確認することを可能にするために提供され、これは、特許請求の範囲または意味を、解釈または限定するために使用するものではないとの理解のもとに提出される。加えて、前述の詳細な説明では、開示を簡素化する目的で、多様な特徴が、様々な実施形態において、一緒にグループとしてまとめられていることが理解できる。本開示の方法は、請求された実施形態が、各請求項に明示的に記載されているよりも、多くの特徴を必要とするという意図を反映していると解釈されるべきではない。そうではなく、添付の特許請求の範囲が反映するように、発明の主題は、単一の開示された実施形態の全ての特徴より少ない特徴にある。したがって、添付の特許請求の範囲は、ここで詳細な説明に組み込まれ、各特許請求の範囲は、別々に特許請求される主題としてそれ自体で存在する。
【図１】
2020502906.tif 000546
【図２】
2020502906.tif 000547
【図３】
2020502906.tif 000548
【図４】
2020502906.tif 000549
【図５】
2020502906.tif 000550
【図６】
2020502906.tif 000551
【図７】
2020502906.tif 000552
【図８】
2020502906.tif 000553
【図９】
2020502906.tif 000554
【図１０】
2020502906.tif 000555
【図１１】
2020502906.tif 000556
【図１２】
2020502906.tif 000557
【図１３】
2020502906.tif 000558
【図１４】
2020502906.tif 000559
【図１５】
2020502906.tif 000560
【図１６】
2020502906.tif 000561
【図１７】
2020502906.tif 000562
【図１８】
2020502906.tif 000563
【図１９】
2020502906.tif 000564
【図２０】
2020502906.tif 000565
【図２１】
2020502906.tif 000566
【図２２】
2020502906.tif 000567
【図２３】
2020502906.tif 000568
【図２４】
2020502906.tif 000569
【図２５】
2020502906.tif 000570
【図２６】
2020502906.tif 000571
【図２７】
2020502906.tif 000572
【図２８】
2020502906.tif 000573
【図２９】
2020502906.tif 000574
【図３０】
2020502906.tif 000575
【図３１】
2020502906.tif 000576
【図３２】
2020502906.tif 000577
【図３３】
2020502906.tif 000578
【図３４】
2020502906.tif 000579
【図３５】
2020502906.tif 000580
【手続補正書】
【提出日】令和1年10月31日(2019.10.31)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
ｉ）既知の相対的な位置および配向を有する位相トランスデューサアレイから音響場を生成すること；
ｉｉ）複数のコントロールポイントであって、それぞれ前記位相トランスデューサアレイに相対して既知の空間的関係を有する前記複数のコントロールポイントを定義すること；
ｉｉｉ）前記位相トランスデューサアレイにより生成される量の物理場において既知の方向に所望する場の値の実質的な生成を最適化することであって、前記量がベクトル量であること；
を含む方法。
【請求項２】
さらに、固有値システムの合計を使用して最適化マトリクスを生成することをさらに含む、請求項１に記載の方法。
【請求項３】
さらに、少なくとも１つの場において、少なくとも１つのヌルポイントを使用して最適化マトリクスをさらに生成することを含む、請求項１に記載の方法。
【請求項４】
前記所望する場の値の少なくとも１つが、ホスト媒体内での既知のポイントの粒子の速度を実質的に記述する、請求項１に記載の方法。
【請求項５】
前記所望する場の値の１つが、波の運動量を実質的に記述する、請求項１に記載の方法。
【請求項６】
さらに、複数の基底関数を計算するための第１の値を使用して最適化を生成する、請求項１に記載の方法。
【請求項７】
複数の前記基底関数は、前記物理場の第２の値を特定する少なくとも１つのコントロールポイントを実現するために使用され、前記第２の値は、前記第１の値とは異なる、請求項６に記載の方法。
【請求項８】
ｉ）既知の相対的な位置および配向を有する位相トランスデューサアレイから音響場を生成すること；
ｉｉ）複数のコントロールポイントであって、それぞれ前記位相トランスデューサアレイに相対して既知の空間的関係を有する前記複数のコントロールポイントを定義すること；
ｉｉｉ）前記位相トランスデューサアレイにより生成された量の物理場における既知の方向に沿って所望する場の値の絶対的な大きさを実質的に最大化するように最適化することであって、前記量がベクトル量であること；
を含む方法。
【請求項９】
さらに、固有値システムの合計を使用して最適化マトリクスを生成することをさらに含む、請求項１に記載の方法。
【請求項１０】
さらに、少なくとも１つの場において、少なくとも１つのヌルポイントを使用して最適化マトリクスをさらに生成することを含む、請求項８に記載の方法。
【請求項１１】
前記所望する場の値の少なくとも１つが、音響強度を記述する、請求項８に記載の方法。
【請求項１２】
前記所望する場の値の１つが、エネルギーフラックスベクトルを実質的に記述する、請求項８に記載の方法。
【請求項１３】
前記所望する場の値の少なくとも１つが、ホスト媒体内での既知のポイントの粒子の速度を実質的に記述する、請求項８に記載の方法。
【請求項１４】
前記所望する場の値の１つが、波の運動量を実質的に記述する、請求項８に記載の方法。
【請求項１５】
さらに、複数の基底関数を計算するための第１の値を使用して最適化を生成する、請求項８に記載の方法。
【請求項１６】
複数の前記基底関数は、前記物理場の第２の値を特定する少なくとも１つのコントロールポイントを実現するために使用され、前記第２の値は、前記第１の値とは異なる、請求項１５に記載の方法。
【請求項１７】
ｉ）既知の相対的な位置および配向を有する位相トランスデューサアレイから音響場を生成すること；
ｉｉ）複数のコントロールポイントであって、それぞれ前記位相トランスデューサアレイに相対して既知の空間的関係を有する前記複数のコントロールポイントを定義すること；および
ｉｉｉ）前記位相トランスデューサアレイにより生成された少なくとも１つの物理場における場の値の所望する領域での重み付け積分を実質的に最大化するように最適化すること；
を含む方法。
【請求項１８】
さらに、固有値システムの合計を使用して最適化マトリクスを生成することをさらに含む、請求項１７に記載の方法。
【請求項１９】
さらに、少なくとも１つの場において、少なくとも１つのヌルポイントを使用して最適化マトリクスをさらに生成することを含む、請求項１７に記載の方法。
【請求項２０】
前記所望する場の値の少なくとも１つが、音響強度を記述する、請求項１７に記載の方法。
【請求項２１】
前記所望する場の値の１つが、エネルギーフラックスベクトルを実質的に記述する、請求項１７に記載の方法。
【請求項２２】
前記所望する場の値の少なくとも１つが、音響圧を実質的に記述する、請求項１７に記載の方法。
【請求項２３】
前記所望する場の値の少なくとも１つが、ホスト媒体内での既知のポイントの粒子の速度を実質的に記述する、請求項１７に記載の方法。
【請求項２４】
前記所望する場の値の１つが、波の運動量を実質的に記述する、請求項１７に記載の方法。
【請求項２５】
さらに、複数の基底関数を計算するための第１の値を使用して最適化を生成する、請求項１７に記載の方法。
【請求項２６】
複数の前記基底関数は、前記物理場の第２の値を特定する少なくとも１つのコントロールポイントを実現するために使用され、前記第２の値は、前記第１の値とは異なる、請求項２５に記載の方法。
【請求項２７】
複数の複素値入力を有するハードウェアを操作すること；
前記複数の複素値入力のそれぞれを、べき乗関連の処理において複素値出力に変換すること；
前記複素値入力の前記実部の２を底とする（バイナリ）べき乗の大きさを有する係数と；
完全な回転の整数乗の二乗として定義される角度ユニットで乗算した前記複素値入力の前記虚部である角度引数と
を含む方法。
【請求項２８】
前記変換は、前記複素数のバイナリ表現に対するビットシフト、比較および計算を通して達成される、請求項２７に記載の方法。
【請求項２９】
さらに、前記複数の複素値入力のそれぞれの成分の実質的な丸めに対応する値を、さらなる部分と別別に処理して、入力および出力範囲の削減を達成することを含む請求項２７に記載の方法。
【請求項３０】
さらに、１ではない値に出力複素値を初期化して、反復の適用が、非１の値を変更して出力することを意図する元々の値による複素乗算を計算することを含む、請求項２７に記載の方法。
【請求項３１】
さらに、複素値を生成して、前記複素値出力に適用される操作をまたさらなる複素値にも適用して、出力することを意図する元々の複素値による前記さらなる複素値の乗算を与える、請求項２７に記載の方法。
【請求項３２】
さらに、冗長数システムを使用して残りの複素値に適用する反復計算の処理を実装することを含む、請求項２７に記載の方法。
【請求項３３】
複数の複素値入力を有するハードウェアを操作すること；
前記複数の複素値入力のそれぞれを、対数を求めることに関連する処理において複複素値出力に変換することであって、前記出力複素値出力が、
前記複素値入力の係数の２を底とする（バイナリ）対数に等しい実部と；
前記複素値入力の角度引数において出現する角度ユニットの数に等しい虚部であって、前記角度ユニットが前記完全な回転の整数乗の二乗である、
方法。
【請求項３４】
前記変換は、、前記複素数のバイナリ表現に対するビットシフト、比較および計算を通して達成される、請求項３３に記載の方法。
【請求項３５】
さらに、実質的に象限の回転および前記複数の複素値入力のそれぞれの成分をさらなる部分とは別々にビットシフトすることに対応する値の処理を行い、入力および出力の削減を達成する、請求項３３に記載の方法。
【請求項３６】
さらに、非ゼロの値を有する出力複素値を生成して、前記複素値出力に適用される操作をまたさらなる複素値にも適用して、出力することを意図する元々の複素値による複素加算を計算する、請求項３３に記載の方法。
【請求項３７】
さらに、複素値を生成して、前記複素値入力に適用される操作をまたさらなる複素値にも適用して、前記さらなる複素値の前記複素値入力による除算を与える、請求項３３に記載の方法。
【請求項３８】
さらに、冗長数システムを使用して残りの複素値に適用する反復計算の処理を実装することを含む、請求項３３に記載の方法。
【国際調査報告】
2020502906.tif 000581
2020502906.tif 000582
2020502906.tif 000583
2020502906.tif 000584