使用库仑定律在二维平面上由偶极子产生的电场分布（Matlab代码实现）

  💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

 ⛳️赠与读者

💥1 概述

偶极子电场的计算

二维平面内电场的简化

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现

---

 ⛳️赠与读者

👨‍💻做科研，涉及到一个深在的思想系统，需要科研者逻辑缜密，踏实认真，但是不能只是努力，很多时候借力比努力更重要，然后还要有仰望星空的创新点和启发点。当哲学课上老师问你什么是科学，什么是电的时候，不要觉得这些问题搞笑。哲学是科学之母，哲学就是追究终极问题，寻找那些不言自明只有小孩子会问的但是你却回答不出来的问题。建议读者按目录次序逐一浏览，免得骤然跌入幽暗的迷宫找不到来时的路，它不足为你揭示全部问题的答案，但若能让人胸中升起一朵朵疑云，也未尝不会酿成晚霞斑斓的别一番景致，万一它居然给你带来了一场精神世界的苦雨，那就借机洗刷一下原来存放在那儿的“躺平”上的尘埃吧。

     或许，雨过云收，神驰的天地更清朗.......🔎🔎🔎

💥1 概述

使用库仑定律计算一个偶极子的电场。本文计算在二维平面上由偶极子产生的电场分布。使用库仑定律，在二维平面的每个格点计算电场。在x和y方向各有100个格点，每个格点长度为1米。坐标轴被操纵以使图形的中心看起来是笛卡尔坐标原点，即(x,y)=(0,0)。偶极子电荷为1纳库仑，偶极子长度为10米（相隔10个格点）。

本研究致力于通过应用经典的库仑定律来精确量化二维平面上由一个典型偶极子所产生的电场分布特性。为了深入探究这一物理现象，我们采取了一种细致的空间离散化策略，即在x轴和y轴上均匀布置100个网格点，每个网格的边长精确定义为1米，以此构建了一个详尽的二维坐标系统。通过巧妙地调整坐标轴的布局，我们确保图形的视觉中心精确对应于常规笛卡尔坐标系的原点位置，即坐标点(x,y)=(0,0)，这种安排便于直观理解和分析电场分布的对称性和特性。

在这一框架下，偶极子被理想化为一对等量异号的电荷，总电量设为微小的1纳库仑（1nC），这两点电荷沿垂直于xy平面的方向相距10米，换算成我们的离散化网格，即相当于它们在10个网格点间隔上对置。通过精确计算每一对网格节点上的库仑力贡献，并综合考虑所有网格节点上电场矢量的叠加效应，我们能够重构出整个平面上由偶极子引起的电场强度和方向的精细图像。这种方法不仅深刻揭示了偶极子电场的基本规律，还为后续在更复杂数学模型中模拟和预测电场行为奠定了坚实的基础，尤其是在微观粒子相互作用分析、材料科学以及电磁学领域的相关研究中展现出重要价值。

库仑定律是静电学的基本定律之一，描述了两个静止点电荷之间的相互作用力，其数学表达式为 𝐹=𝑘∣𝑞1𝑞2∣𝑟2F=kr2∣q1​q2​∣​，其中 𝐹F 是作用力的大小，𝑘k 是库仑常数（在真空中为 9×109Nm2/C29×109Nm2/C2），𝑞1q1​ 和 𝑞2q2​ 分别是两个点电荷的电荷量，𝑟r 是它们之间的距离。然而，当涉及到偶极子产生的电场分布时，我们需要考虑的是由一对等量异号电荷组成的系统。

在二维平面上研究偶极子产生的电场分布，我们可以简化问题为考虑两个点电荷：一个正电荷 +𝑞+q 和一个负电荷 −𝑞−q，两者相距 2𝑎2a（偶极矩长度），且假设偶极子沿垂直于平面的方向放置（为了简化，因为我们主要关注平面内的电场分布）。在二维情况下，实际分析仍涉及三维空间的电场，但由于对称性，我们仅关心平面内的电场分量。

偶极子电场的计算

偶极子在任意点 (𝑥,𝑦)(x,y) 处产生的电场 𝐸E 可以通过对两个点电荷单独产生的电场矢量进行矢量合成得到。对于正负电荷各自产生的电场，我们可以应用点电荷的电场公式 𝐸=𝑘𝑞𝑟2𝑟^E=r2kq​r^，其中 𝑟^r^ 是指向电荷的单位向量，但在二维分析中，我们将使用 𝑥−𝑦x−y 平面的坐标表示。

假设观察点 𝑃(𝑥,𝑦)P(x,y) 到正负电荷的距离分别为 𝑟+r+​ 和 𝑟−r−​，那么：

• 正电荷 +𝑞+q 在 𝑃P 点产生的电场 𝐸+E+​ 为：𝐸+=𝑘𝑞𝑟+2𝑟+^E+​=r+2​kq​r+​^​
• 负电荷 −𝑞−q 在 𝑃P 点产生的电场 𝐸−E−​ 为：𝐸−=−𝑘𝑞𝑟−2𝑟−^E−​=−r−2​kq​r−​^​

其中，𝑟+^r+​^​ 和 𝑟−^r−​^​ 分别指向 +𝑞+q 和 −𝑞−q 到 𝑃P 的方向。

二维平面内电场的简化

在二维平面上，我们通常关注沿 𝑥x 或 𝑦y 方向的电场分量，尤其是当偶极子沿着 𝑧z 轴放置时。对于偶极子在原点且 +𝑞+q 位于 +𝑎+a 而 −𝑞−q 位于 −𝑎−a 时，考虑点 (𝑥,𝑦)(x,y) 的电场，我们仅需要计算沿 𝑥x 轴的电场分量，因为沿 𝑦y 轴的分量因对称性相互抵消。

偶极矩 𝑝=2𝑎𝑞p=2aq，电场强度 𝐸E 可以通过下面的公式计算（特别地，在远离偶极子的区域，即 𝑟≫2𝑎r≫2a 时的近似）：

𝐸𝑥=14𝜋𝜖0𝑝cos⁡𝜃𝑟3Ex​=4πϵ0​1​r3pcosθ​

这里，𝜖0ϵ0​ 是真空介电常数，𝜃θ 是从偶极子指向观察点的方向与 𝑥x 轴之间的夹角。

📚2 运行结果

部分代码：

% Constant 1/(4*pi*epsilon_0) = 9*10^9
k = 9*10^9;

% Enter the Relative permittivity
eps_r = 1;
charge_order = 10^-9; % milli, micro, nano etc..
const = k*charge_order/eps_r;

% Enter the dimensions
Nx = 101; % For 1 meter
Ny = 101; % For 1 meter

% Enter the number of charges.
n = 2;

% Electric fields Initialization
E_f = zeros(Nx,Ny);
Ex = E_f;
Ey = E_f;

% Vectors initialization
ex = E_f;
ey = E_f;
r = E_f;
r_square = E_f;

% Array of charges
Q = [1,-1];

% Array of locations
X = [5,-5];
Y = [0,0];

%-------------------------------------------------------------------------%
%                   COMPUTATION OF ELECTRIC FIELDS
%-------------------------------------------------------------------------%

%  Repeat for all the 'n' charges
for k = 1:n

🎉3 参考文献

文章中一些内容引自网络，会注明出处或引用为参考文献，难免有未尽之处，如有不妥，请随时联系删除。

[1]仲维畅.平面上任意形状范围内的磁荷分布[J].无损检测, 2008(1):4.DOI:10.3969/j.issn.1000-6656.2008.01.008.

[2]解树青.计算偶极子磁场与电场的一种简捷方法[J].滨州学院学报, 2004, 020(002):55-57.

[3]王瑶、刘玉颖、朱世秋.静电场及电荷运动的VPython可视化模拟[J].大学物理, 2020, 39(9):6.DOI:CNKI:SUN:DXWL.0.2020-09-016. 

🌈4 Matlab代码实现

资料获取，更多粉丝福利，MATLAB|Simulink|Python资源获取