exm041332_1.m

matlab mfiles实例分析

M=5;		%做2的幂次用。本例把M设得较小，是为了观察混迭。				<1>
tend=1;	%波形取非零值的时间长度。
T=10;		%窗口化长度应足够大，以减小窗口化引起的泄露“旁瓣”效应。		<3>
N=2^M;		%采样点数，取2的幂是为使FFT运算较快。
dt=T/N;	%以上T、N的取值应使N/T=fs采样频率大于两倍时间波形带宽，以克服
			%采样引起的频谱混迭。
			%在本例中，据理论分析知W(f=7.5)=Sa(7.5*pi)=1/(7.5*pi)<5% 。
         	%因此，可近似认为本例时间信号带宽为7.5Hz 。
n=0:N-1;	%采样序列
t=n*dt;	%采样点时间序列
w=zeros(size(t,2),1);
Tow=find((tend-t)>0);			%产生非零波形时段的相应序列
w(Tow,1)=ones(length(Tow),1);	%在窗口时段内定义的完整波形
plot(t,w,'b','LineWidth',2.5),title('Time Waveform');xlabel('t --- >')
 
[AW,f]=cftbyfft(w,t,0);
ff=f+eps;				%为避免下面指令出现0/0而采取的措施
AWW=abs(sin(pi*ff)./(pi*ff));
plot(f,AW,'b-',ff,AWW,'r:')
title('Aliasing caused by undersampling')
xlabel('f --- >');ylabel('|W(f)|'),legend('by FFT','Theoretical')