电偶极子磁场在静态和高速运动状态下的电场方向图

电偶极子磁场在静态和高速运动状态下的电场方向图,其中dcc,为赤道面方向图，dcz为子午面方向图

close all;clear all;  clc
v=10.^8;
lambda=0.001;
c=3*10^8;
beta=v/c;
gamma=1/sqrt（1-beta.^2）;
k=（2*pi）/lambda;
omega=c*k;
epsilon=（1/（36*pi））*（10.^-9）;
I=0.000001;
t=5;
dl=0.00001;
k1=gamma*k-beta*gamma*c*t;
r=1;
sita=pi/2;
fai=linspace（02*pi）;
E=ones（1100）;
subplot（121）;
polar（faiabs（E））;
title（‘电基本振子静态电场赤道面‘）;
E1=（（（I*dl*（k1.^3）.*gamma.*sin（sita）.*cos（sita）.*cos（fai））./（4*pi*omega*epsilon））.*（1i./（k1.*r））.*exp（-1i*k1.*r）-c*beta*gamma*t）.*sin（sita）.*cos（fai）+（（I*dl*（k.^3）*sin（sita）.*cos（sita）.*sin（fai））/（4*pi*omega*epsilon））*（1i/（k*r））*exp（-1i*k*r）.*sin（sita）.*sin（fai）-（（I*dl*（k.^3））/（4*pi*omega*epsilon））.*（1i/（k*r））.*exp（-1i*k*r）.*（（sin（sita））.^2）.*cos（sita）;
E2=（（（I*dl*（k1.^3）.*gamma.*sin（sita）.*cos（sita）.*cos（fai））./（4*pi*omega*epsilon））.*（1i./（k1.*r））.*exp（-1i*k1.*r）-c*beta*gamma*t）.*cos（sita）.*cos（fai）+（（I*dl*（k.^3））/（4*pi*omega*epsilon））.*（1i/（k*r））.*exp（-1i*k*r）.*（（sin（sita））.^2）.*sin（sita）;（（I*dl*（k.^3）*sin（sita）.*cos（sita）.*sin（fai））/（4*pi*omega*epsilon））*（1i/（k*r））*exp（-1i*k*r）.*cos（sita）.*sin（fai）+（（I*dl*（k.^3））/（4*pi*omega*epsilon））.*（1i/（k*r））.*exp（-1i*k*r）.*（（sin（sita））.^2）.*sin（sita）;
E3=（（（I*dl*（k1.^3）.*gamma.*sin（sita）.*cos（sita）.*cos（fai））./（4*pi*omega*epsilon））.*（1i./（k1.*r））.*exp（-1i*k1.*r）-c*beta*gamma*t）.*（-sin（fai））+（（I*dl*（k.^3）*sin（sita）.*cos（sita）.*sin（fai））/（4*pi*omega*epsilon））*（1i/（k*r））*exp（-1i*k*r）.*（cos（fai））;
e1=E1.*conj（E1）;
e2=E2.*conj（E2）;
e3=E3.*conj（E3）;
E4=sqrt（e1+e2+e3）;
e=E4/max（E4）;
subplot（122）;
polar（faiabs（e））;
title（‘电基本振子向x轴正方向运动速度为10^8m/s时电场赤道面‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2018-06-27 10:25  电偶极子磁场在静态和高速运动状态下的电场方向图\
     文件        1599  2018-05-31 20:44  电偶极子磁场在静态和高速运动状态下的电场方向图\ddc.m
     文件        1512  2018-05-31 20:45  电偶极子磁场在静态和高速运动状态下的电场方向图\ddz.m