スプレッドシートのBESSELJ関数の使い方｜第1種ベッセル関数 Jn(x) を計算する

スポンサーリンク

スプレッドシートで第1種ベッセル関数 J_n(x) を求めたいけれど、どの関数を使えばいいかわからない。そんな場面はありませんか。

振動や波動の解析、電磁場の計算では、ベッセル関数の値が必要になります。手計算や紙の数表で求めるのは現実的ではないですよね。

GoogleスプレッドシートのBESSELJ関数を使えば、セルに数式を入力するだけで第1種ベッセル関数の値を一発で計算できますよ。この記事では、BESSELJ関数の構文から実務での使いどころ、エラー対処法まで丁寧に解説します。

スプレッドシートのBESSELJ関数とは？第1種ベッセル関数を計算する関数

BESSELJ関数は、第1種ベッセル関数 J_n(x) の値を返すスプレッドシートの関数です。読み方は「ベッセル・ジェイ」です。

ベッセル関数は、円筒座標系の微分方程式を解くときに登場する特殊関数です。工学や物理学の分野で広く使われています。

GoogleスプレッドシートではExcelと同じ仕様で動作します。Excel 2007以降と互換性があるため、Excel経験者にも違和感なく使えますよ。

BESSELJ関数の数学的な意味

BESSELJ関数が扱う「第1種ベッセル関数」は、ベッセルの微分方程式の有界解として定義される関数です。

J_n(x) の特徴は、x の増加とともに振動しながら減衰していくことです。ちょうど sin(x) や cos(x) が減衰していくようなイメージです。これが変形ベッセル関数 I_n(x) と大きく異なる点で、BESSELI は指数関数的に増加しますが、BESSELJ は振動します。

特に重要な性質として、J_0(0) = 1、J_n(0) = 0（n ≥ 1）があります。原点でゼロでない値をもつかどうかが、第2種ベッセル関数 Y_n(x) との使い分けの基準になります。

どんなときに使う？

BESSELJ関数は、主に次のような場面で活躍します。

  • 円形膜（太鼓の膜）の振動モード計算
  • 電磁波の解析（円形導波管の固有モード、アンテナ設計）
  • 音響解析（スピーカーの放射パターン）
  • FM変調のスペクトル解析（搬送波と側波帯の振幅計算）

理工系の研究やエンジニアリング業務で使うことが多い関数ですね。

BESSELJ関数の構文と引数

基本構文

=BESSELJ(x, n)

引数は2つで、どちらも必須です。省略するとエラーになります。

引数の説明

引数必須/省略可説明
x必須関数を評価する値（実数）
n必須ベッセル関数の次数（0以上の整数）

引数 x には、関数に代入したい数値を指定します。正の値・負の値・0のいずれも指定できますよ。セル参照で渡すのが一般的です。

引数 n には、ベッセル関数の次数（じすう）を指定します。0以上の整数を入力してください。小数を渡すと小数点以下が自動的に切り捨てられます。たとえば n に 2.7 を指定すると、2として計算されますよ。

BESSELJ関数の基本的な使い方

実際にスプレッドシートでBESSELJ関数を使ってみましょう。代表的なパターンを3つ紹介します。

数値を直接指定する

セルに次のように入力します。

=BESSELJ(1.5, 1)

この数式は、x=1.5、次数 n=1 のときの第1種ベッセル関数 J_1(1.5) の値を返します。結果は約 0.5579 になります。

数式バーに直接数値を打ち込むだけなので、ちょっと値を確認したいときに便利ですよ。

セル参照を使う

A1セルに x の値、B1セルに次数 n を入力しておけば、次のように書けます。

=BESSELJ(A1, B1)

A1 に 1.5、B1 に 1 を入力すると、先ほどと同じ約 0.5579 という結果になります。パラメータを変えて繰り返し計算するときは、セル参照のほうが圧倒的にラクです。

次数ごとの値を一覧で確認する

x の値を固定して、次数 n を 0 から順に変えていくと、関数の振る舞いがわかりやすくなります。

x の値次数 n数式結果（概算）
1.50=BESSELJ(1.5, 0)0.5118
1.51=BESSELJ(1.5, 1)0.5579
1.52=BESSELJ(1.5, 2)0.2321
1.53=BESSELJ(1.5, 3)0.0610

次数が大きくなるほど x の小さい領域での値が小さくなる特徴がよくわかりますね。この一覧をスプレッドシートで作っておくと、パラメータの傾向をつかみやすいですよ。

x=0 のときの特殊な性質

x=0 を代入した場合、次数によって結果が変わります。

=BESSELJ(0, 0)  → 結果: 1
=BESSELJ(0, 1)  → 結果: 0
=BESSELJ(0, 2)  → 結果: 0

J_0(0) = 1、J_n(0) = 0（n ≥ 1）というのは第1種ベッセル関数の重要な性質です。第2種ベッセル関数 Y_n(x) は x=0 で発散するため BESSELY(0, n) はエラーになります。この挙動の違いが、BESSELJ と BESSELY の使い分けの重要なポイントです。

BESSELJ関数の実務活用パターン

ここからはスプレッドシートならではの便利な使い方を紹介します。表形式での一括計算と、振る舞い確認の2パターンです。

ARRAYFORMULAで複数の x 値を一気に計算する

スプレッドシートの強みである ARRAYFORMULA関数（範囲全体に数式を一気に展開する関数）と組み合わせると、複数の入力値に対して BESSELJ を一発で計算できます。

A2:A6 に x の値（0, 1, 2, 3, 4）を、B2:B6 に次数（すべて 0）を入れておきます。C2 に次の数式を入れると、C2:C6 まで一気に展開されますよ。

=ARRAYFORMULA(BESSELJ(A2:A6, B2:B6))
xnJ_n(x)（結果）
001.0000
100.7652
200.2239
30-0.2601
40-0.3971

x が大きくなると J_0(x) が振動していく様子が、表からも確認できます。x=3 付近で初めて負の値になり、振動しながら 0 に収束していく特徴が一目でわかりますね。

FM変調スペクトルを確認する

FM（周波数変調）のスペクトル解析では、変調指数 β のときの n 次側波帯の振幅が J_n(β) になります。たとえば β=2 のとき、各側波帯の振幅を一覧にすると次のようになります。

B列に次数 n（0 から 6）を入れて、C列に =BESSELJ(2, B2) を入力 → 下方向にコピーします。

n（次数）J_n(2)
00.2239
10.5767
20.3528
30.1289
40.0340
50.0070
60.0012

この表をスプレッドシートで作成してグラフ化すると、FMスペクトルの棒グラフが完成します。無線技術の資格学習や信号処理の検算にも使えますよ。

BESSELJ/BESSELY/BESSELI/BESSELK の使い分け

スプレッドシートにはBESSELJのほかに、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: 変形ベッセル関数。熱伝導や拡散など指数的な振る舞いの問題で使う

BESSELJ と BESSELY はペアで使うことが多いので、セットで覚えておくと便利ですよ。

場面別の選び方

実務で「どれを使えばいい？」と迷ったときは、解きたい現象の性質で判断しましょう。

  • 振動・波動現象（弦の振動、円形膜の振動、波動方程式）→ BESSELJ / BESSELY
  • 熱伝導・拡散現象（円柱の温度分布、拡散方程式）→ BESSELI / BESSELK
  • 原点で有限の解が欲しい（例: 円形膜の中心が有限） → BESSELJ / BESSELI
  • 原点で特異（無限大）の解を含む（例: 円筒外側の場）→ BESSELY / BESSELK

どの関数を使うかは、最終的には解いている微分方程式の境界条件で決まります。迷ったら関連現象から逆引きしてみてくださいね。

Excel版の各関数についても詳しく知りたいときは、BESSELJ関数・BESSELY関数・BESSELI関数・BESSELK関数 の各記事も参考になりますよ。

BESSELJ関数のよくあるエラーと対処法

BESSELJ関数で発生しやすいエラーと、その対処法をまとめました。

#VALUE! エラー（文字列を渡したとき）

x または n に数値以外の値を指定すると、#VALUE! エラーが表示されます。

=BESSELJ("abc", 1)   → #VALUE! エラー

対処法: 引数に文字列が入っていないか確認してください。セル参照の場合、参照先のセルが空白や文字列になっていることがありますよ。

#NUM! エラー（次数が負のとき）

n に負の整数を指定すると #NUM! エラーになります。

=BESSELJ(1.5, -1)   → #NUM! エラー

対処法: 次数 n は 0 以上の整数を指定してください。

#NAME? エラー（スペルミス）

関数名のスペルミスで発生します。「BESSEL」と「J」の間にスペースを入れたり、綴りを間違えたりすると起こります。

対処法: 関数名が BESSELJ （スペースなし）になっているか、もう一度確認してみてください。スプレッドシートの数式入力時に表示される候補から選ぶと、スペルミスを防げますよ。

BESSELY との比較で混乱したとき

BESSELJ と BESSELY は似た名前ですが、x=0 付近での挙動が全く異なります。

=BESSELJ(0, 0)   → 1（正常）
=BESSELY(0, 0)   → #NUM! エラー（Y_0(0) は負の無限大に発散）

もし「なぜ BESSELY でエラーが出るのか」と疑問に思ったら、この性質を思い出してください。

まとめ

この記事では、スプレッドシートのBESSELJ関数について解説しました。ポイントを振り返っておきましょう。

  • BESSELJ関数は第1種ベッセル関数 J_n(x) の値を返すエンジニアリング関数
  • 構文は =BESSELJ(x, n) で、引数は2つとも必須
  • x には正・負・0 のどれでも指定でき、n は 0 以上の整数（小数は切り捨て）
  • J_0(0) = 1、J_n(0) = 0（n ≥ 1）という重要な性質がある
  • ARRAYFORMULA と組み合わせると複数の値を一気に計算できる
  • 4種類のベッセル関数は振動系（J/Y）と熱伝導系（I/K）で使い分ける
  • n が負で #NUM!、文字列入力で #VALUE!、スペルミスで #NAME? エラー

ベッセル関数は理工系のニッチな関数ですが、必要になる場面では他の方法では代替できません。手元のスプレッドシートで即座に計算できると、検算や試行錯誤の効率がぐっと上がりますよ。

関連記事

スポンサーリンク
タイトルとURLをコピーしました