求助高手指点，有关if的用法，其中 if 判断部分不知道怎么编写

maoihua

程序如下：
c0=340;
f=1000;
w=2*pi*f;
k0=0;
M=41;
N=41;
d=0.05;
num1=0:M-1;
num2=0:N-1;
xnum1=(num1-(M-1)/2)*d;
ynum2=(num2-(N-1)/2)*d;
image=sqrt(-1);
syms kx ky;
km=[kx,ky];
W=0;
if kx^2+ky^2<=(2*w/c0)^2
for n=1:41
    rn=[xnum1(n),ynum2(n)];
    pn=dot(km,rn);
    qn=exp(image.*pn);
    W=W+qn;
end
end
ezcontourf(W)
colorbar

maoihua

已经解决

ChaChing

maoihua 发表于 2012-5-12 09:49

                               
登录/注册后可看大图

已经解决

LZ忘记分享成果!?:@)

maoihua

ChaChing 发表于 2012-5-14 01:09

                               
登录/注册后可看大图

LZ忘记分享成果!?

哦 呵呵 那个程序可以不用if语句。程序如下：
lc;
clear;
c0=340;
f=1000;
w=2*pi*f;
k0=0;
M=11;
N=11;
d=0.05;
num1=0:M-1;
num2=0:N-1;
xnum1=(num1-(M-1)/2)*d;
ynum2=(num2-(N-1)/2)*d;
image=sqrt(-1);
syms kx ky;
km=[kx,ky];
W=0;
for n=1:11
    for m=1:11
    rn=[xnum1(n),ynum2(m)];
    pn=dot(km,rn);
    qn=exp(image.*pn);
    W=W+qn;
    end
end
figure(1)
ezcontourf(W)
colorbar