スカラーの符号関数 sign(x) は単純だ。正なら +1、負なら −1。信号処理で言えば「理想的なステップ」で、周波数ドメインで完璧なバンドパスの境界を引くことに対応する。こんな単純な関数が、行列へ、制御理論へ、そして量子コンピュータへと繋がっていく。

行列の符号関数とは何か

行列 A に対して sign(A) を定義できる。直感的には、固有値の実部が正なら対応するブロックを +I に、負なら −I に置き換えた行列だ。固有値が虚軸のどちら側にあるかを「スナップショット」で捉える操作、と思えばいい。

これが何に使えるか。連続時間の代数リッカティ方程式（CARE）の解は制御系の安定化フィードバックゲインを与えるが、その解を求めるアルゴリズムの中に sign(A) の計算が本質的に現れる。行列の安定・不安定部分空間を分離するときも同様だ。

Zolotarevの1877年の答え

問題は計算だ。sign(x) は不連続で、多項式ではうまく近似できない（Gibbsの現象がつきまとう）。では有理関数（多項式の比）ならどうか？

1877年、Yegor Zolotarev は「sign 関数に対するミニマックス有理近似」の明示解を与えた。解の零点と極はヤコビの楕円関数で表せる。フィルタ設計に親しんだ人には見覚えがあるはずで、楕円フィルタ（Cauerフィルタ）の極・零点の配置とまったく同じ数学だ。理想的なステップを有理関数で近似する ── という問題が底でつながっていた。

Zolotarev近似には特別な性質がある：高次の近似を低次の近似の合成で作れること。これが後の計算効率に直結する。

Newton法との対応

行列符号関数を求めるNewtonの反復は形が美しい：

A_{k+1} = (A_k + A_k^{-1}) / 2

スカラーで試すと x = 2 から始めて 2 → 1.25 → 1.025 → 1.0003 ... とすごい速さで 1 に向かう。これが二次収束だ。適切なスケーリングを入れたバージョンはZolotarev近似と本質的に等価で、10〜15反復で実用的な精度に達する。「2反復のZolotarev関数で基本行列分解が精度よく計算できる」と示した2014年のSIAM論文がその対応をきれいにまとめている。

量子アルゴリズムへの流れ込み

量子信号処理（QSP）は、量子コンピュータ上でユニタリの固有値に関数を適用する枠組みだ。実装には有理関数や多項式による関数近似が必要で、構造はZolotarevの問題とほぼ同じになる。

今月 arXiv に出た論文「Sign Embedding Quantum Algorithms for Matrix Equations and Matrix Functions」は、行列のSylvester方程式・行列平方根・行列幾何平均・CAREの解を、「符号関数への埋め込み」という統一的な戦略で量子的に解くアルゴリズムをまとめた。ほぼ同時に公開された人間-AI共同発見レポートでは、研究の起点が「jump型の関数には有理近似が特に有効」という人間の直感だったことが記されている。そこからAIが数学的探索を広げて、具体的なアルゴリズムの発見に至ったという経緯だ。

149年分の蓄積

1877年のZolotarevの結果が、楕円フィルタ設計を経由して、行列計算の数値解法になり、量子アルゴリズムの中核に現れる。「sign 関数を有理近似する」というたった一つの発想がこれだけ遠くまで繋がるのは、数学の中でも珍しい部類だと思う。149年分の蓄積が、今年の論文の中に静かに息づいている。

— ランキン

出典

一次情報（arXiv・学術誌）

• Wang et al., “Sign Embedding Quantum Algorithms for Matrix Equations and Matrix Functions,” arXiv:2604.25333 (2026-04)
• Wang et al., “From Meta Idea to Advanced Mathematical Discovery – Human-AI Co-Discovery of Sign-Embedding Quantum Algorithms,” arXiv:2606.24899 (2026-06)
• Nakatsukasa & Freund, “Computing Fundamental Matrix Decompositions Accurately via the Matrix Sign Function in Two Iterations: The Power of Zolotarev’s Functions,” SIAM Review 56(3), 2014

arXiv 掲載のプレプリントは著者の報告値。2604.25333 の査読状況は確認中。