当前位置：首页 > news >正文

电动力学问题中的Matlab可视化

news 2025/9/18 9:40:49

电磁场的经典描述

小说一则

电磁场的经典描述就是没有啥玩意量子力学的经典电动力学下对电磁场的描述,以后有空写个科幻小说,写啥呢,就写有天张三遇见了一个外星人,外星人来自这样一个星球,星球上的物质密度特别低,导致外星人的测量会明显的影响物质的运动,外星人不能同时得到动量与位置啊啥的......
外星人很希望有一天能找到一个不那么量子力学的世界(这依然是经典的"量子观")
张三是个优秀的气象工作员,搞出来的全是李雅普诺夫不稳定的数学模型
张三见了外星人表示很好,咱们臭味相投啊
外星人表示,是的,学数学物理,物理数学有啥意思呢,不如享受当下的生活
张三表示,好哇好哇,好个锤子哇,我们要探索宇宙,你不觉得李雅普诺夫的理论很过时吗,不够先进啦
外星人表示:我说一个像是反驳了但是没有成功反驳的话啊:李雅普诺夫是苏俄时期逝世的
张三表示:...........
外星人表示:还是要好好学习哇!

四个方程组

麦克斯韦方程组

$\left\{\begin{matrix} \iint_S\vec{D}\cdot ds=Q_f \\ \iint+S\vec{B}\cdot ds = 0\\ \oint_L\vec{E}\cdot dl=-\frac{d\Phi_B}{dt}\\ \oint_L\vec{H}\cdot dl = I_f + \frac{d\Phi_D}{dt} \end{matrix}\right.$

$\left\{\begin{matrix} \bigtriangledown \cdot \vec{D}=\rho\\ \bigtriangledown \cdot \vec{B} = 0\\ \bigtriangledown \times \vec{E} = -\frac{\partial \vec{B}}{\partial t}\\ \bigtriangledown \times \vec{H} = \vec{J}_f+\frac{\partial \vec{D}}{\partial t}\\ \end{matrix}\right.$

电磁物质方程组

$\begin{matrix} \vec{J}=\sigma \vec{E}\\ \vec{D}=\varepsilon \vec{E}=\varepsilon_0\vec{E}+\vec{P}\\ \vec{B}=\mu \vec{H}=\mu_0(\vec{H}+\vec{M}) \end{matrix}$

电磁边值方程组

$\left\{\begin{matrix} \vec{e_n}\times (\vec{E_2}-\vec{E_1})=0 \\ \vec{e_n}\times (\vec{H_2}-\vec{H_1})=\alpha \\ \vec{e_n}\cdot(\vec{D_2}-\vec{D_1})=\sigma\\ \vec{e_n}\cdot(\vec{B_2}-\vec{B_1})=0\\ \end{matrix}\right.$

势场边值方程组

电势 $\left\{\begin{matrix} \varphi_1=\varphi_2\\ \varepsilon_2\frac{\partial \varphi_2}{\partial n}-\varepsilon_1\frac{\partial \varphi_1}{\partial n}=-\sigma \end{matrix}\right.$

库伦规范下的磁矢势 $\left\{\begin{matrix} \bigtriangledown \cdot \vec{A}=0\\ \vec{A_1}=\vec{A_2} \\ \Delta \vec{A} = -\mu \vec{J} \end{matrix}\right.$

(旋度的旋度=散度的散度-拉普拉斯算子)

磁标势 $\left\{\begin{matrix} \varphi_1=\varphi_2\\ \frac{\partial \varphi_2}{\partial n}-\frac{\partial \varphi_1}{\partial n}=-\alpha \end{matrix}\right.$

数据可视化

重要提示

不论是数据可视化还是数值计算
- 都不可能替代你的理论分析过程
- 永远是你先得到了一个差不多的结果，才去做数值分析，可视化等等。
数据可视化常常受限于你的计算机和人眼
- 不要太较真

常见绘图语言

一些常见符号不再赘述

流线 Streamlines

将矢量平滑得连接起来的不相交的曲线

[x,y] = meshgrid(0:0.1:1,0:0.1:1);
u = sin(x);
v = -cos(y);figure
quiver(x,y,u,v)startx = 0.1:0.1:1;
starty = ones(size(startx));
streamline(x,y,u,v,startx,starty)

matlab 语法

streamline(X,Y,Z,U,V,W,startx,starty,startz)
streamline(U,V,W,startx,starty,startz)
streamline(XYZ)
streamline(X,Y,U,V,startx,starty)
streamline(U,V,startx,starty)

等值线

连接值相等的多个点的不相交的封闭曲线

x = linspace(-2*pi,2*pi);
y = linspace(-2*pi,2*pi);
[X,Y] = meshgrid(x,y);
Z = tan(X)-cos(Y);
contour(X,Y,Z)