divergence

シンボリックなベクトル場の発散

ページ内をすべて折りたたむ

構文

d = divergence(V,X)

d = divergence(V)

説明

例

d = divergence(V,X) は、直交座標のベクトル X について、シンボリックなベクトル場 V の発散を返します。ベクトル V と X は、同じ長さでなければなりません。

d = divergence(V) は、ベクトル場 V のシンボリック変数から作成された既定のベクトルについて、V の発散を返します。

例

すべて折りたたむ

ベクトル場の発散を求める

ライブスクリプトを開く

ベクトル $X = (x, y, z)$ に関して、ベクトル場の発散 $V (x, y, z) = (x, 2 y^{2}, 3 z^{3})$ を計算します。

syms x y z
V = [x 2*y^2 3*z^3];
X = [x y z];
div = divergence(V,X)

div =  $9 z^{2} + 4 y + 1$

ベクトル場の回転の発散は 0 になることを示します。

divCurl = divergence(curl(V,X),X)

divCurl =  $0$

スカラー場 $f (x, y, z) = x^{2} + y^{2} + z^{2}$ の勾配の発散を求めます。結果は、スカラー場のラプラシアンになります。

syms x y z
f = x^2 + y^2 + z^2;
divGrad = divergence(gradient(f,X),X)

divGrad =  $6$

電界から電荷密度を求める

ライブスクリプトを開く

微分形式のガウスの法則は、電界の発散が電荷密度に比例することを記述しています。

$\nabla_{}^{\to} \cdot E_{}^{\to} (r_{}^{\to}) = \frac{ρ (r_{}^{\to})}{ϵ_{0}}$

電界 $E_{}^{\to} = x^{2} i_{}^{ˆ} + y^{2} j_{}^{ˆ}$ に対する電荷密度を求めます。

syms x y ep0
E = [x^2 y^2];
rho = divergence(E,[x y])*ep0

rho =  ${ep}_{0} (2 x + 2 y)$

ep0 = 1 として、 $- 2 < x < 2$ および $- 2 < y < 2$ における電界と電荷密度を可視化します。meshgrid を使用して、x と y の値におけるグリッドを作成します。subs を使用してグリッド値を代入することによって、電界と電荷密度の値を求めます。subs への入力として cell 配列を使用して、グリッド値 xPlot および yPlot を同時に、電荷密度 rho に代入します。

rho = subs(rho,ep0,1);
v = -2:0.1:2;
[xPlot,yPlot] = meshgrid(v);
Ex = subs(E(1),x,xPlot);
Ey = subs(E(2),y,yPlot);
rhoPlot = double(subs(rho,{x,y},{xPlot,yPlot}));

quiver を使用して電界をプロットします。contour を使用して電荷密度を重ね合わせます。等高線は電荷の値を示します。

quiver(xPlot,yPlot,Ex,Ey)
hold on
contour(xPlot,yPlot,rhoPlot,"ShowText","on")
title("Contour Plot of Charge Density Over Electric Field")
xlabel("x")
ylabel("y")

Figure contains an axes object. The axes object with title Contour Plot of Charge Density Over Electric Field, xlabel x, ylabel y contains 2 objects of type quiver, contour.

電磁波方程式の導出

R2023a 以降

ライブスクリプトを開く

マクスウェル方程式から、電荷および電流源のない自由空間における電磁波の方程式を導出します。

最初に、真空透磁率と誘電率を表すシンボリックスカラー変数を作成します。直交座標を表すシンボリック行列変数を作成します。電界と磁場を空間と時間の関数として表す 2 つのシンボリック行列関数を作成します。

syms mu_0 epsilon_0
syms X [3 1] matrix
syms E(X,t) B(X,t) [3 1] matrix keepargs

次に、マクスウェル方程式を表す 4 つの方程式を作成します。

Maxwell1 = divergence(E,X) == 0

Maxwell1(X, t) =  $\nabla_{X} \cdot E (X, t) = 0_{1, 1}$

Maxwell2 = curl(E,X) == -diff(B,t)

Maxwell2(X, t) = 
 $\nabla_{X} \times E (X, t) = - \frac{\partial}{\partial t} B (X, t)$

Maxwell3 = divergence(B,X) == 0

Maxwell3(X, t) =  $\nabla_{X} \cdot B (X, t) = 0_{1, 1}$

Maxwell4 = curl(B,X) == mu_0*epsilon_0*diff(E,t)

Maxwell4(X, t) = 
 $\nabla_{X} \times B (X, t) = (ε_{0} μ_{0}) \frac{\partial}{\partial t} E (X, t)$

次に、電界の波動方程式を求めます。2 番目のマクスウェル方程式の回転を計算します。

wave_E = curl(Maxwell2,X)

wave_E(X, t) = 
 $\nabla_{X} (\nabla_{X} \cdot E (X, t)) - Δ_{X} E (X, t) = - \nabla_{X} \times (\frac{\partial}{\partial t} B (X, t))$

電界の波動方程式に最初のマクスウェル方程式を代入します。lhs と rhs を使用して、最初のマクスウェル方程式の左辺と右辺を取得します。

wave_E = subs(wave_E,lhs(Maxwell1),rhs(Maxwell1))

wave_E(X, t) = 
 $- Δ_{X} E (X, t) = - \nabla_{X} \times (\frac{\partial}{\partial t} B (X, t))$

4 番目のマクスウェル方程式の時間微分を計算します。

dMaxwell4 = diff(Maxwell4,t)

dMaxwell4(X, t) = 
 $\nabla_{X} \times (\frac{\partial}{\partial t} B (X, t)) = (ε_{0} μ_{0}) \frac{\partial}{\partial t} \frac{\partial}{\partial t} E (X, t)$

wave_E の磁場 $\nabla_{X} \times \frac{\partial}{\partial t} B (X, t)$ を含む項を dMaxwell4 の右辺で置き換えます。lhs と rhs を使用して、dMaxwell4 からこれらの項を取得します。

wave_E = subs(wave_E,lhs(dMaxwell4),rhs(dMaxwell4))

wave_E(X, t) = 
 $- Δ_{X} E (X, t) = - (ε_{0} μ_{0}) \frac{\partial}{\partial t} \frac{\partial}{\partial t} E (X, t)$

同様の手順を使用して、磁場の波動方程式も求めることができます。

wave_B = curl(Maxwell4,X)

wave_B(X, t) = 
 $\nabla_{X} (\nabla_{X} \cdot B (X, t)) - Δ_{X} B (X, t) = (ε_{0} μ_{0}) (\nabla_{X} \times (\frac{\partial}{\partial t} E (X, t)))$

wave_B = subs(wave_B,lhs(Maxwell3),rhs(Maxwell3))

wave_B(X, t) = 
 $- Δ_{X} B (X, t) = (ε_{0} μ_{0}) (\nabla_{X} \times (\frac{\partial}{\partial t} E (X, t)))$

dMaxwell2 = diff(Maxwell2,t)

dMaxwell2(X, t) = 
 $\nabla_{X} \times (\frac{\partial}{\partial t} E (X, t)) = - \frac{\partial}{\partial t} \frac{\partial}{\partial t} B (X, t)$

wave_B = subs(wave_B,lhs(dMaxwell2),rhs(dMaxwell2))

wave_B(X, t) = 
 $- Δ_{X} B (X, t) = - (ε_{0} μ_{0}) \frac{\partial}{\partial t} \frac{\partial}{\partial t} B (X, t)$

入力引数

すべて折りたたむ

`V` — シンボリックなベクトル場
シンボリックスカラー変数のベクトル | シンボリック関数 | シンボリック行列変数 | シンボリック行列関数

シンボリックなベクトル場。シンボリックスカラー変数、シンボリック関数、シンボリック行列変数、またはシンボリック行列関数のベクトルとして指定します。V は X と同じ長さでなければなりません。

V がシンボリックスカラー変数の関数であり、V が sym 型または symfun 型である場合、ベクトル X は sym 型または symfun 型でなければならない。
V がシンボリック行列変数の関数であり、V が symmatrix 型または symfunmatrix 型である場合、ベクトル X は symmatrix 型または symfunmatrix 型でなければならない。

データ型: sym | symfun | symmatrix | symfunmatrix

`X` — 発散を求める対象のベクトル
シンボリックスカラー変数のベクトル | シンボリック関数 | シンボリック行列変数 | シンボリック行列関数

発散を求める対象のベクトル。シンボリックスカラー変数、シンボリック関数、シンボリック行列変数、またはシンボリック行列関数のベクトルとして指定します。X は V と同じ長さでなければなりません。

X を指定せず、V がシンボリックスカラー変数の関数である場合、既定では、divergence は、symvar(V) で定義された変数の順序で、V のシンボリックスカラー変数からベクトル X を構成する。
X が symmatrix 型のシンボリック行列変数である場合、X は 1 行 N 列または N 行 1 列のサイズでなければならない。
V と X がスカラーの場合、divergence(V,X) = diff(V,X) となる。

データ型: sym | symfun | symmatrix | symfunmatrix

制限

関数 divergence はテンソル微分をサポートしません。入力 V がベクトルではなくテンソル場または行列である場合、関数 divergence はエラーを返します。
現在、Symbolic Math Toolbox™ は、シンボリック行列変数の関数 dot または関数 cross と、symmatrix 型および symfunmatrix 型の関数をサポートしていません。ベクトル微積分の恒等式がドット積またはクロス積を含む場合、ツールボックスは代わりに他のサポートされている関数でそれらの恒等式を表示します。シンボリック行列の変数と関数をサポートするすべての関数のリストを表示するには、コマンド methods symmatrix と methods symfunmatrix を使用します。