您所在位置：主题：[讨论]DFt的问题？
qingzhou
[专家分：0] 发布于 2007-11-22 18:57:00
我求一个DFT的算法，程序如下：
dt=0.04;
t=0:dt:1;
f=30;
M=2.5;
a=2*f^2*log(M);
x=exp(-a*t.*t).*sin(2*pi*f*t);
N=length(t);
W=exp(-j*2*pi/N);
X=zeros(1,N);
for k0=0:N-1;
for n=0:N-1;
X(k0)=X(k0)+x(n).*W^(n.*k0); //错误行
end
end
它提示有这样是错误：??? Attempted to access X(0); index must be a positive integer or logical.
Error in ==> DFT at 13
X(k0)=X(k0)+x(n).*W^(n.*k0);
这是为什么啊？？[em18][em18]
dt=0.04;
t=0:dt:1;
f=30;
M=2.5;
a=2*f^2*log(M);
x=exp(-a*t.*t).*sin(2*pi*f*t);
N=length(t);
W=exp(-j*2*pi/N);
X=zeros(1,N);
for k0=0:N-1;
for n=0:N-1;
X(k0)=X(k0)+x(n).*W^(n.*k0); //错误行
end
end
它提示有这样是错误：??? Attempted to access X(0); index must be a positive integer or logical.
Error in ==> DFT at 13
X(k0)=X(k0)+x(n).*W^(n.*k0);
这是为什么啊？？[em18][em18]
