MATLABで磁場を計算する　その１

磁場に関するクーロンの法則は、磁荷\(m_1,m_2\)、距離\(r\)とすると

\begin{align}
F=k\frac{m}{r^2} \mathrm{[A/m]}
\end{align}

ここで\(k\)は単位に関する定数で

\begin{align}
k=\frac{1}{4 \pi \mu}
\end{align}

透磁率\(\mu\)は真空の透磁率\(\mu_0=4 \pi \times 10^{-7}\)と物質により定まる比透磁率\(\mu_r\)を用いて

\begin{align}
\mu=\mu_0 \mu_r
\end{align}

で計算できる。以下コード

N=1000;
mu_0=4*pi*10^-7;
mu_r=1;
mu=mu_0*mu_r;
k=1/(4*pi*mu);
m=3*10^-5;
r=linspace(0,0.1,N);

F=k.*m./r.^2;

figure;
semilogy(r,F,'k');
grid on
xlabel("r[m]")
ylabel("H[A/m]")
xlim([0,0.1])