大家帮忙看下这程序,总是提示未定义
实现通信里抽样定理的matlab程序1)
function =fftseq(m,ts,df)
=fftseq(m,ts,df) ;
=fftseq(m,ts) ;
fs=1/ts;
if nargin == 2
n1=0;
else
n1=fs/df;
end
n2=length(m);
n=2^(max(nextpow2(n1),nextpow2(n2)));
M=fft(m,n);
m=;
df=fs/n;
2)clear
t0=10;
ts=0.001;
ts1=0.1;
ts2=0.05;
fs=1/ts;
fs1=1/ts1;
fs2=1/ts2;
df=0.5;
t=[-t0/2:ts:t0/2];
x=sin(20*t);
m=x./t;
w=t0/(2*ts)+1;
m(w)=20;
m=m.*cos(100*t);
=fftseq(m,ts,df);
M=M/fs;
f=-fs/2;
pause;
subplot(2,1,1);
plot(f,abs(fftshift(M)));
axis([-50,50,0.3]);
t1=[-t0/2:ts1:t0/2];
x1=sin(20*t1);
m1==x1./t1;
w1=t0/(2*ts1)+1;
m1(w1)=20;
m1=m1.*cos(100*t1);
=fftseq(m1,ts1,df);
M1=M1/fs1;
N1=;
f1=[-7*df1*length(mn1):df1:6*df1*length(mn1)-df1]-fs1/2;
subplot(2,2,3);
plot(f1,abs(fftshift(N1)));
title('不满足采样要求');
t2=[-t0/2:ts2:t0/2];
x2=sin(20*t2);
m2=x2./t2;
w2=t0/(2*ts2)+1;
m2(w2)=20;
m2=m2.*cos(100*t2);
=fftseq(m2,ts2,df);
M2=M2/fs2;
N2=;
f2=[-7*df2*length(mn2):df2:6*df2*length(mn2)-df2]-fs2/2;
subplot(2,2,4);
title('满足采样要求');
xlable('频率')
出现错误:??? Maximum recursion limit of 500 reached. Use set(0,'RecursionLimit',N)
to change the limit.Be aware that exceeding your available stack space can
crash MATLAB and/or your computer.
Error in ==> e:\MATLAB6p5\work\fftseq.m
On line 2==> =fftseq(m,ts,df) ;
Error in ==> e:\MATLAB6p5\work\bandpsample.m
On line 16==> =fftseq(m,ts,df); Use set(0,'RecursionLimit',N) to change the limit
上句不是提醒了你了吗 set(0,'RecursionLimit',N)
加在t=[-t0/2:ts:t0/2];后面?
我菜鸟啊,解释下 原帖由 skyhitskaka 于 2007-4-8 16:36 发表
set(0,'RecursionLimit',N)
加在t=;后面?
我菜鸟啊,解释下
set(0,'RecursionLimit',1000) ;
加在程序开头吧，反正在fftseq前面就行 不行
还是一样的提示
用了2个不同的程序,都报一样的错
郁闷 原帖由 skyhitskaka 于 2007-4-8 17:28 发表
不行
还是一样的提示
用了2个不同的程序,都报一样的错
郁闷
这个不是错误，只是提示你最大的迭代次数已经达到，N可以自己调整，不过如果太大容易造成matlab或者系统崩溃 那是程序本身的错误吗
Error in ==> e:\MATLAB6p5\work\fftseq.m
On line 2==> =fftseq(m,ts,df) ;
Error in ==> e:\MATLAB6p5\work\bandpsample.m
On line 16==> =fftseq(m,ts,df);
那个会不会影响程序
在其他网页上看到过,如果0被除,会导致后面的变换出错
是这样吗
请赐教
页:
[1]