samples.m

来自「两种数字滤波程序」· M 代码 · 共 36 行

function samples(f,Fs,M) 

N=2^M; % fft点数=取样总点数 

Ts=1/Fs; % 取样时间间隔 

T=N*Ts; % 取样总时间=取样总点数*取样时间间隔 

n=0:N-1; 

t=n*Ts; 

Xn=cos(2*f*pi*t); 

subplot(2,1,1); 

plot(t,Xn,'*'); 

axis([0 T 1.1*min(Xn) 1.1*max(Xn)]); 

xlabel('t -->'); 

ylabel('Xn'); 

Xk=abs(fft(Xn,N)); 

subplot(2,1,2); 

stem(n,Xk); 

axis([0 N 1.1*min(Xk) 1.1*max(Xk)]); 

xlabel('frequency -->'); 

ylabel('!Xk!');