cos1.m

来自「利用matlab解海伦方程并画出电流分布图」· M 代码 · 共 12 行

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=cos(2*pi/l.*abs(-l/4+(n-0.5).*dx)).*exp(-i*2*pi/l.*r)./r;
    z=z+temp*dx;
end
z