sin1.m

利用matlab解海伦方程并画出电流分布图

nmax=1000;
l=1;
xx=[0,l/8,l/4];
a=7.022/1000*l;
dx=l/2/nmax;
z=0;
for n=1:nmax
    r=sqrt((xx-(-l/4+(n-0.5).*dx)).^2+a^2);
    temp=sin(2*2*pi/l.*abs(-l/4+(n-0.5).*dx)).*exp(-i*2*pi/l.*r)./r;
    z=z+temp*dx;
end
z