untitled8.m

来自「这是一段模拟电荷法（scm)的matlab源程序」· M 代码 · 共 54 行

clear
grid on
x=-40:0.5:40;
a=13.72;b=12.19;
y=1.5;
c=8.85e-12;
QR=2*pi*c*[54705 -73523 24287];
Qi=2*pi*c*[46378 19701 -67517];
Q=QR+Qi*j;
A=[(y-b)./sqrt((a-x).^2+(b-1.5).^2)-(y-b)./sqrt((a-x).^2+(-b-1.5).^2);(y-b)./sqrt(x.^2+(b-1.5).^2)-(y-b)./sqrt(x.^2+(-b-1.5).^2);
(y-b)./sqrt((-a-x).^2+(b-1.5).^2)-(y-b)./sqrt((-a-x).^2+(-b-1.5).^2)];
Ey=abs(Q*A./2*pi*c)*1e21;
subplot(3,3,1);
y=1.5;
plot(x,y)
hold on;
plot([0 a -a -a 0 a],[b b b -b -b -b],'o')
title('算例示意图');
ylabel('离地面高度');
xlabel('水平距离');
subplot(3,3,2);
plot(x,Ey);
title('算例一的电场强度曲线');
xlabel('x(m)');
ylabel('E(V/m)');
axis([-40 40 0 12000]);
a1=20;b1=8;c1=25;
QR=2*pi*c*[75864 -40753 -31074];
Qi=2*pi*c*[-7926 70575 -69662];
Q=QR+Qi*j;
A=[(y-b1)./sqrt((a1-x).^2+(b1-1.5).^2)-(y-b1)./sqrt((a1-x).^2+(-b1-1.5).^2);(y-c1)./sqrt(x.^2+(b1-1.5).^2)-(y-c1)./sqrt(x.^2+(-b1-1.5).^2);
(y-b1)./sqrt((-a1-x).^2+(b1-1.5).^2)-(y-b1)./sqrt((-a1-x).^2+(-b1-1.5).^2)];
B=[(x-a1)./sqrt((a1-x).^2+(b1-1.5).^2)-(x-a1)./sqrt((a1-x).^2+(-b1-1.5).^2);(x)./sqrt(x.^2+(b1-1.5).^2)-(x)./sqrt(x.^2+(-b1-1.5).^2);
(x+a1)./sqrt((-a1-x).^2+(b1-1.5).^2)-(x+a1)./sqrt((-a1-x).^2+(-b1-1.5).^2)];
Ey1=abs(Q*A./2*pi*c)*1e20;
Ey2=abs(Q*B./2*pi*c)*1e20;
Ey3=sqrt(Ey1.^2+Ey2.^2);
subplot(3,3,3);
y=1.5;
plot(x,y)
hold on;
plot([0 a1 -a1 -a1 0 a1],[c1 b1 b1 -b1 -c1 -b1],'o')
title('算例示意图');
ylabel('离地面高度');
xlabel('水平距离');
subplot(3,3,4);
plot(x,Ey1);
title('算例二的电场强度曲线');
xlabel('x(m)');
ylabel('E(V/m)');
subplot(3,3,5);
plot(x,Ey2)
subplot(3,3,6);
plot(x,Ey3)