スプレッドシートで第2種変形ベッセル関数 K_n(x) を求めたいけれど、どの関数を使えばいいかわからない。そんな場面はありませんか。
熱伝導や電磁場の減衰計算、拡散現象の解析では、ベッセル関数の値が必要になります。手計算や紙の数表で求めるのは現実的ではないですよね。
GoogleスプレッドシートのBESSELK関数を使えば、セルに数式を入力するだけで第2種変形ベッセル関数の値を一発で計算できますよ。この記事では、BESSELK関数の構文から実務での使いどころ、エラー対処法まで丁寧に解説します。
スプレッドシートのBESSELK関数とは？第2種変形ベッセル関数を計算する関数
BESSELK関数は、第2種変形ベッセル関数 K_n(x) の値を返すスプレッドシートの関数です。読み方は「ベッセル・ケイ」です。
ベッセル関数は、円筒座標系の微分方程式を解くときに登場する特殊関数です。工学や物理学の分野で広く使われています。
GoogleスプレッドシートではExcelと同じ仕様で動作します。Excel 2007以降と互換性があるため、Excel経験者にも違和感なく使えますよ。
BESSELK関数の数学的な意味
BESSELK関数が扱う「第2種変形ベッセル関数」は、変形ベッセル方程式の特異解として定義される関数です。
K_n(x) の特徴は、x が大きくなるにつれて指数関数的に減衰することです。これが変形ベッセル関数 I_n(x) と大きく異なる点で、BESSELI は増加しますが、BESSELK は単調に減衰します。
もう一つ重要な性質として、x → 0 のとき K_n(x) は正の無限大に発散します。この「原点で特異」という性質が、BESSELI（原点で有限）との使い分けの基準です。
どんなときに使う？
BESSELK関数は、主に次のような場面で活躍します。
- 円筒座標系で減衰を伴う熱伝導方程式を解くとき
- 電磁場解析で外部領域（無限遠で減衰する）の電場を計算するとき
- 拡散現象のモデリングで境界条件を設定するとき
- 粒子物理学における散乱問題を扱うとき
理工系の研究やエンジニアリング業務で使うことが多い関数ですね。
BESSELIとBESSELKの関係
BESSELI（第1種変形ベッセル関数 I_n(x)）とBESSELK（第2種変形ベッセル関数 K_n(x)）は、変形ベッセル方程式の2つの独立な解です。この2つがペアで変形ベッセル方程式の一般解を構成します。
- I_n(x): 原点で有限、無限大に増加する解（正則解）
- K_n(x): 原点で特異（発散）、無限遠で減衰する解（特異解）
境界条件に応じてどちらか、あるいは両方を使います。セットで覚えておくと便利ですよ。
BESSELK関数の構文と引数
基本構文
=BESSELK(x, n)引数は2つで、どちらも必須です。省略するとエラーになります。
引数の説明
| 引数 | 必須/省略可 | 説明 |
|---|---|---|
| x | 必須 | 関数を評価する値（正の実数のみ） |
| n | 必須 | ベッセル関数の次数（0以上の整数） |
引数 x には、関数に代入したい数値を指定します。BESSELK関数では x は正の数のみ指定できます。 0以下の値を指定すると #NUM! エラーになるので注意してください。これがBESSELI（x は任意の実数を指定可能）と異なる重要なポイントです。
引数 n には、ベッセル関数の次数（じすう）を指定します。0以上の整数を入力してください。小数を渡すと小数点以下が自動的に切り捨てられます。たとえば n に 2.7 を指定すると、2として計算されますよ。
BESSELK関数の基本的な使い方
実際にスプレッドシートでBESSELK関数を使ってみましょう。代表的なパターンを3つ紹介します。
数値を直接指定する
セルに次のように入力します。
=BESSELK(1.5, 1)この数式は、x=1.5、次数 n=1 のときの第2種変形ベッセル関数 K_1(1.5) の値を返します。結果は約 0.2774 になります。
数式バーに直接数値を打ち込むだけなので、ちょっと値を確認したいときに便利ですよ。
セル参照を使う
A1セルに x の値、B1セルに次数 n を入力しておけば、次のように書けます。
=BESSELK(A1, B1)A1 に 1.5、B1 に 1 を入力すると、先ほどと同じ約 0.2774 という結果になります。パラメータを変えて繰り返し計算するときは、セル参照のほうが圧倒的にラクです。
次数ごとの値を一覧で確認する
x の値を固定して、次数 n を 0 から順に変えていくと、関数の振る舞いがわかりやすくなります。
| x の値 | 次数 n | 数式 | 結果（概算） |
|---|---|---|---|
| 1.5 | 0 | =BESSELK(1.5, 0) | 0.2138 |
| 1.5 | 1 | =BESSELK(1.5, 1) | 0.2774 |
| 1.5 | 2 | =BESSELK(1.5, 2) | 0.5836 |
| 1.5 | 3 | =BESSELK(1.5, 3) | 1.9953 |
次数が大きくなるほど値が大きくなる（x が小さい領域では）特徴がよくわかりますね。この一覧をスプレッドシートで作っておくと、パラメータの傾向をつかみやすいですよ。
BESSELK関数の実務活用パターン
ここからはスプレッドシートならではの便利な使い方を紹介します。表形式での一括計算と、x 値の振る舞い確認の2パターンです。
ARRAYFORMULAで複数の x 値を一気に計算する
スプレッドシートの強みである ARRAYFORMULA関数（範囲全体に数式を一気に展開する関数）と組み合わせると、複数の入力値に対して BESSELK を一発で計算できます。
A2:A6 に x の値（0.5, 1.0, 1.5, 2.0, 2.5）を、B2:B6 に次数（すべて 0）を入れておきます。C2 に次の数式を入れると、C2:C6 まで一気に展開されますよ。
=ARRAYFORMULA(BESSELK(A2:A6, B2:B6))| x | n | K_n(x)（結果） |
|---|---|---|
| 0.5 | 0 | 0.9244 |
| 1.0 | 0 | 0.4210 |
| 1.5 | 0 | 0.2138 |
| 2.0 | 0 | 0.1139 |
| 2.5 | 0 | 0.0623 |
x が大きくなるにつれて K_0(x) が急速に減衰していく様子が、表からも一目でわかります。I_0(x) とは逆の傾向ですね。手で1セルずつ数式を書くより、ARRAYFORMULA でまとめて流すほうが圧倒的に速いですよ。
BESSELIとBESSELKのペアで一般解を確認する
変形ベッセル方程式の一般解は C1・I_n(x) + C2・K_n(x) の形になります。たとえば C1=1、C2=1 の場合に、x=1 〜 3 での一般解の値を表で確認するには次のようにします。
A列に x の値（1.0, 1.5, 2.0, 2.5, 3.0）、B列に =BESSELI(A2, 0)、C列に =BESSELK(A2, 0)、D列に =B2+C2 を入力して下にコピーします。
| x | I_0(x) | K_0(x) | 一般解（I+K） |
|---|---|---|---|
| 1.0 | 1.2661 | 0.4210 | 1.6871 |
| 1.5 | 1.6467 | 0.2138 | 1.8605 |
| 2.0 | 2.2796 | 0.1139 | 2.3935 |
| 2.5 | 3.2898 | 0.0623 | 3.3521 |
| 3.0 | 4.8808 | 0.0347 | 4.9155 |
x が大きくなると K_0(x) の寄与が急速に小さくなっていきますね。このようにBESSELIとBESSELKを組み合わせて使うと、変形ベッセル方程式の解の構造がスプレッドシートで直感的に確認できますよ。
BESSELI/BESSELJ/BESSELK/BESSELY の使い分け
スプレッドシートにはBESSELKのほかに、3つのベッセル関数が用意されています。それぞれの違いを整理しておきましょう。
4関数の違いを比較表で整理
| 関数名 | 正式名称 | 数学記号 | 特性 | 主な用途 |
|---|---|---|---|---|
| BESSELJ | 第1種ベッセル関数 | J_n(x) | 振動的・原点で有限 | 振動・波動の正則解 |
| BESSELY | 第2種ベッセル関数 | Y_n(x) | x→0 で発散・振動的 | 振動・波動の特異解 |
| BESSELI | 第1種変形ベッセル関数 | I_n(x) | 単調増加（指数発散） | 熱伝導・拡散の正則解 |
| BESSELK | 第2種変形ベッセル関数 | K_n(x) | 単調減少（指数減衰） | 熱伝導・拡散の特異解 |
大きく分けると2グループあります。
- BESSELJ / BESSELY: 通常のベッセル関数。振動や波動を扱う問題で使う
- BESSELI / BESSELK: 変形ベッセル関数。熱伝導や拡散など指数的な振る舞いの問題で使う
BESSELI と BESSELK はペアで使うことが多いので、セットで覚えておくと便利ですよ。
場面別の選び方
実務で「どれを使えばいい？」と迷ったときは、解きたい現象の性質で判断しましょう。
- 振動・波動現象（弦の振動、円形膜の振動、波動方程式）→ BESSELJ / BESSELY
- 熱伝導・拡散現象（円柱の温度分布、拡散方程式）→ BESSELI / BESSELK
- 原点で有限の解が欲しい（例: 円柱の中心温度が有限） → BESSELI
- 無限遠で減衰する解が欲しい（例: 外部領域の電場）→ BESSELK
どの関数を使うかは、最終的には解いている微分方程式の境界条件で決まります。迷ったら関連現象から逆引きしてみてくださいね。
Excel版の各関数についても詳しく知りたいときは、BESSELK関数・BESSELI関数・BESSELJ関数 の各記事も参考になりますよ。
BESSELK関数のよくあるエラーと対処法
BESSELK関数で発生しやすいエラーと、その対処法をまとめました。
#NUM! エラー（x が 0 以下のとき）
x に 0 以下の値を指定すると #NUM! エラーが表示されます。
=BESSELK(0, 1) → #NUM! エラー
=BESSELK(-1, 1) → #NUM! エラー対処法: BESSELK関数では x は正の数のみ有効です。K_n(x) は x=0 で正の無限大に発散するため、0を含む値は指定できません。x の値が 0 以下になっていないか確認してください。
BESSELI との違いとして、BESSELI(0, n) は正常に計算できますが、BESSELK(0, n) はエラーになります。
#NUM! エラー（次数が負のとき）
n に負の整数を指定しても #NUM! エラーになります。
=BESSELK(1.5, -1) → #NUM! エラー対処法: 次数 n は 0 以上の整数を指定してください。
#VALUE! エラー（文字列を渡したとき）
x または n に数値以外の値を指定すると、#VALUE! エラーが表示されます。
=BESSELK("abc", 1) → #VALUE! エラー対処法: 引数に文字列が入っていないか確認してください。セル参照の場合、参照先のセルが空白や文字列になっていることがありますよ。
#NAME? エラー（スペルミス）
関数名のスペルミスで発生します。
対処法: 関数名が BESSELK （スペースなし）になっているか確認してください。スプレッドシートの数式入力時に表示される候補から選ぶと、スペルミスを防げますよ。
まとめ
この記事では、スプレッドシートのBESSELK関数について解説しました。ポイントを振り返っておきましょう。
- BESSELK関数は第2種変形ベッセル関数 K_n(x) の値を返すエンジニアリング関数
- 構文は
=BESSELK(x, n)で、引数は2つとも必須 - x は正の数のみ指定可能（0以下は
#NUM!エラー） - 次数 n は 0 以上の整数（小数は切り捨て）
- ARRAYFORMULA と組み合わせると複数の値を一気に計算できる
- BESSELI（増加）と BESSELK（減衰）は変形ベッセル方程式のペア解
- ベッセル関数は4種類あり、振動系（J/Y）と熱伝導系（I/K）で使い分ける
- x に 0以下で
#NUM!、文字列入力で#VALUE!、スペルミスで#NAME?エラー
ベッセル関数は理工系のニッチな関数ですが、必要になる場面では他の方法では代替できません。手元のスプレッドシートで即座に計算できると、検算や試行錯誤の効率がぐっと上がりますよ。
