实验15平面线性映射的迭代.docVIP

电磁场分布 一、实验目的 ① 学 1. 静电场中库仑力的计算 要计算平面上N个电荷之间的库仑力，我们可以先考虑任意两个电荷之间的库仑力. 根据库仑定律，q2对q1的引力公式为 方向与q2到q1的矢径r相同，因此，q1，q2同号时，该力是斥力. 其分量的公式可写成 编写程序sy221.m时，先输入电荷的数目、各电荷的坐标及电荷量，之后选定一个电荷，求其他电荷对它的作用力，叠加求合力和方向；再选下一个电荷，以此类推. clear all; %sy222.m N=input(输入电荷数目N=:); for ic=1:N disp([对电荷#,num2str(ic)]);%输入N个电荷的位置和电量 rc=input(输入电荷位置[x, y](米）：); x(ic)=rc(1); %电荷ic的x坐标 y(ic)=rc(2); %电荷ic的y坐标 q(ic)=input(输入电荷量（库仑）：); end E0=8.85e-12; %真空介电常数 C0=1/(4*pi*E0); %合并常数 for ic=1:N %循环计算每个电荷所受的力 Fx=0.0;Fy=0.0; %先把力初始化为零