请求帮助

日月星辰

:@(郁闷请求帮助Sample Text
  两个程序一样的思想，一样的算法，一样的参数，为什么结果不一样？:'(
第一个：
L1=0.4;
L2=0.4;
m=20;
n=20;
k0=1.25;
k=2*pi*400/340;
a=2/k;
M=m*n;
s=(L1*L2)/M;
S=L1*L2;
for i=1:M
    a1=rem(i,m);
    b1=(i-a1)/m;
    if a1~=0
        x1=(a1-1)*L1/m;
        y1=(n-b1-1)*L2/n;
    else
        x1=(m-1)*L1/m;
        y1=(n-b1)*L2/n;
    end
    u(i)=x1+L1/(2*m);
    w(i)=y1+L2/(2*n);
    for j=1:M
        a2=rem(j,m);
        b2=(j-a2)/m;
        if a2~=0
            x2=(a2-1)*L1/m;
            y2=(n-b2-1)*L2/n;
        else
            x2=(m-1)*L1/m;
            y2=(n-b2)*L2/n;
        end
        if j~=i
            r=sqrt((x2-x1)^2+(y2-y1)^2);
            R(i,j)=k0*k0*s*s*sin(k0*r)/(2*pi*S*k0*r);
        else
            R(i,j)=k0*k0*s*s/(2*pi*S);
        end
    end
end
[v,d]=eig(R);
y=v(:,(397));
for i=1:30
    theta(i)=(i-1)*pi/29;
    x11=0.2;
    y11=0.2-a*cos(pi*(i-1)/29);
    z11=a*sin(pi*(i-1)/29);
    for j=1:M
        r(j)=sqrt((u(j)-x11)^2+(w(j)-y11)^2+(z11)^2);
        g(j)=k*s*(0+1i)*exp((0-1i)*k*r(j))/(2*pi*r(j));
        G(i,j)=g(j);
    end
    F=g*y;
    rho(i)=abs(F);
    Y(i)=angle(F);
end
subplot(1,2,1)
polar(theta,rho,'-o')
subplot(1,2,2)
plot(theta,Y,'-o')

第二个：调用两个M函数
fenbutuxing:
k=2*pi*1500/340;
R=2/k;
for i=1:30
    theta(i)=(i-1)*pi/29;
    x=0.2;
    y=0.2-R*cos(theta(i));
    z=R*sin(theta(i));
    q(i)=fenbumotai(4,x,y,z);
    r(i)=abs(q(i));
    ang(i)=angle(q(i));
end
subplot(1,2,1)
polar(theta,r,'-o')
subplot(1,2,2)
plot(theta,ang,'-o')

fenbumotai:
function F=fenbumotai(I,d,e,f)
L1=0.4;
L2=0.4;
m=20;
n=20;
k=2*pi*400/340;
M=m*n;
s=(L1*L2)/M;
S=L1*L2;
for j=1:M
    a1=rem(j,m);
    b1=(j-a1)/m;
    if a1~=0
        x1=(a1-1+0.5)*L1/m;
        y1=(n-b1-1+0.5)*L2/n;
    else
        x1=(m-1+0.5)*L1/m;
        y1=(n-b1+0.5)*L2/n;
    end
    r1=sqrt((e-x1)^2+(d-y1)^2+(f)^2);
    g(j)=k*s*(1i)*exp((-1i)*k*r1)/(2*pi*r1);
end
g;
F=g*mode(I);

mode:
function y=mode(x)
L1=0.4;
L2=0.4;
m=20;
n=20;
k=2*pi*400/340;
k0=0.5/L1;
M=m*n;
s=(L1*L2)/M;
S=L1*L2;
for i=1:M
        a1=rem(i,m);
        b1=(i-a1)/m;
    if a1~=0
        x1=(a1-1)*L1/m;
        y1=(n-b1-1)*L2/n;
    else x1=(m-1)*L1/m;
        y1=(n-b1)*L2/n;
    end
    for j=1:M
        a2=rem(j,m);
        b2=(j-a2)/m;
    if a2~=0
        x2=(a2-1)*L1/m;
        y2=(n-b2-1)*L2/n;
    else x2=(m-1)*L1/m;
        y2=(n-b2)*L2/n;
    end  
    if j~=i
        r=sqrt((x2-x1)*(x2-x1)+(y2-y1)*(y2-y1));
        R(i,j)=k0*k0*s*s*sin(k0*r)/(2*pi*S*k0*r);
    else
        R(i,j)=k0*k0*s*s/(2*pi*S);
        end
      end
end
[v,d]=eig(R);
y=v(:,(M-x+1));