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

电磁场的经典描述

小说一则

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

四个方程组

麦克斯韦方程组

\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

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

matlab代码----Copyright 2015 The MathWorks, Inc.

[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)

 matlab 语法

contour(X,Y,Z,'线型',等高线表示的值)

数据可视化案例一:匀强电场中的电导体球

数据可视化案例二:匀强电场中的电介质球

  • 0
    点赞
  • 6
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

River Chandler

谢谢,我会更努力学习工作的!!

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值