Estude equações de Maxwell 

Derive a equação de onda para um campo magnético a partir das equações de Maxwell.

Insira as equações de Maxwell em unidades de Lorentz-Heaviside naturais.

In[1]:=

X Inactivate[{gaussE = Div[\[ScriptCapitalE][x, y, z, t], {x, y, z}] == \[Rho][x, y, z, t], gaussB = Div[\[ScriptCapitalB][x, y, z, t], {x, y, z}] == 0, faraday = Curl[\[ScriptCapitalE][x, y, z, t], {x, y, z}] == -D[\[ScriptCapitalB][x, y, z, t], t], ampere = Curl[\[ScriptCapitalB][x, y, z, t], {x, y, z}] == j[x, y, z, t] + D[\[ScriptCapitalE][x, y, z, t], t]}, Div | Curl | D];

Crie uma tabela agradavelmente formatada das equações.

mostre o input completo de Wolfram Languageoculte o input

In[2]:=

X Inactivate[{gaussE = Div[\[ScriptCapitalE][x, y, z, t], {x, y, z}] == \[Rho][x, y, z, t], gaussB = Div[\[ScriptCapitalB][x, y, z, t], {x, y, z}] == 0, faraday = Curl[\[ScriptCapitalE][x, y, z, t], {x, y, z}] == -D[\[ScriptCapitalB][x, y, z, t], t], ampere = Curl[\[ScriptCapitalB][x, y, z, t], {x, y, z}] == j[x, y, z, t] + D[\[ScriptCapitalE][x, y, z, t], t]}, Div | Curl | D];; Grid[Transpose[{{"Gauss's Law for Electricity", "Gauss's Law for Magnetism", "Faraday's Law", "Ampere's Law"}, %}], Alignment -> Left, Dividers -> All, Spacings -> {4, 2}] // TraditionalForm

Out[2]//TraditionalForm=

Tome o rotacional da lei de Ampere no vácuo ( e ).

In[3]:=

X Inactive[Curl][#, {x, y, z}] & /@ ampere /. j -> (0 &)

Out[3]=

Inverta a ordem de diferenciação.

In[4]:=

X Inactive[Curl][#, {x, y, z}] & /@ ampere /. j -> (0 &); % /. Inactivate[Curl[D[v_, t_], x_] :> D[Curl[v, x], t]]

Out[4]=

Substitua a lei de Faraday.

In[5]:=

X Inactivate[{gaussE = Div[\[ScriptCapitalE][x, y, z, t], {x, y, z}] == \[Rho][x, y, z, t], gaussB = Div[\[ScriptCapitalB][x, y, z, t], {x, y, z}] == 0, faraday = Curl[\[ScriptCapitalE][x, y, z, t], {x, y, z}] == -D[\[ScriptCapitalB][x, y, z, t], t], ampere = Curl[\[ScriptCapitalB][x, y, z, t], {x, y, z}] == j[x, y, z, t] + D[\[ScriptCapitalE][x, y, z, t], t]}, Div | Curl | D];; % /. Inactivate[Curl[D[v_, t_], x_] :> D[Curl[v, x], t]]; % /. Solve[faraday, Inactive[Curl][\[ScriptCapitalE][x, y, z, t], {x, y, z}]][[1]]

Out[5]=

Ativando os resultados da equação na equação de onda para o campo magnético.

In[6]:=

X Inactivate[{gaussE = Div[\[ScriptCapitalE][x, y, z, t], {x, y, z}] == \[Rho][x, y, z, t], gaussB = Div[\[ScriptCapitalB][x, y, z, t], {x, y, z}] == 0, faraday = Curl[\[ScriptCapitalE][x, y, z, t], {x, y, z}] == -D[\[ScriptCapitalB][x, y, z, t], t], ampere = Curl[\[ScriptCapitalB][x, y, z, t], {x, y, z}] == j[x, y, z, t] + D[\[ScriptCapitalE][x, y, z, t], t]}, Div | Curl | D];; % /. Solve[faraday, Inactive[Curl][\[ScriptCapitalE][x, y, z, t], {x, y, z}]][[1]]; wave = Activate[%]

Out[6]=

Esta é a equação expressa usando notação tradicional.

Uma solução de onda plana pode ser verificada por meio de uma substituição simples.

In[7]:=

X wave /. \[ScriptCapitalB] -> (A Exp[{kx, ky, kz}.{#1, #2, #3} - I Sqrt[kx^2 + ky^2 + kz^2] #4] &) // Simplify

Out[7]=