examp57.m

来自「ADSP TOOLBOX: Version 2.0 and gui m-file」· M 代码 · 共 18 行

clc,echo on 
%EXAMPLE 57 
df=5; N=256;N1=1024;f0=30;S=128;n=0:N-1;
x=cos(2*pi*f0*n/S)+cos(2*pi*(f0+df)*n/S);
X=abs(fft(x,N1))/N1;f=(0:N1-1)*S/N1;plot(f,X),grid,pause(2)
xb=x.*window('vonh',N);
XB=abs(fft(xb,N1))/N1;plot(f,XB),grid,pause(2)
% Repeat for df=2 and df=1.

df=5; N=512;N1=1024;f0=30;S=128;n=0:N-1;
x2=cos(2*pi*f0*n/S)+cos(2*pi*(f0+df)*n/S);
X2=abs(fft(x2,N1))/N1;f=(0:N1-1)*S/N1;plot(f,X2),grid,pause(2)
xb2=x2.*window('vonh',N);
XB2=abs(fft(xb2,N1))/N1;plot(f,XB2)
% Repeat for df=2 and df=1.
echo off %end of example