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(0,2*pi);
E=ones(1,100);
subplot(1,2,1);
polar(fai,abs(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(1,2,2);
polar(fai,abs(e));
title('电基本振子向x轴正方向运动速度为10^8m/s时电场赤道面');