无限大障板上的活塞辐射

%clear all;
close all;
f=100000;w=f*pi/180; %频率
k=w/343; %波束
a=2; %活塞半径
r=1000; %远场距离

theta=（-90:0.5:90）*pi/180; %方向性

sum=zeros（1length（theta））;
sum1=zeros（1length（theta））;
for i=1:length（theta）
    h=r;
    sum（i）=sum（i）+w/h*exp（-1j*k*h）;
    sum1（i）=sum（i）;
    for R=0.1:0.01:a %循环半径
        l=R*80;
        for n=1:l %循环个数 一圈内的个数
            %不近似
            h=sqrt（r*r+R*R-2*r*R*sin（theta（i））*cos（n*2*pi/l））;
            sum（i）=sum（i）+w/h*exp（-1j*k*h）;
            %近似
            h1=r-R*sin（theta（i））*cos（n*2*pi/l）;
            sum1（i）=sum（i）+w/r*exp（-1j*k*h）;

        end
    end
end

figure
subplot（121）
plot（thetaabs（sum））;title（‘不近似‘）
subplot（122）
plot（thetaabs（sum1））;title（‘近似‘）


figure
subplot（121）
plot（thetaabs（sum））;title（‘叠加模拟‘）


D=zeros（1length（theta