lxcaiyang.m

来自「功率谱分析是信号处理的重要内容」· M 代码 · 共 19 行

%连续信号的采样
f=10;
t=0:0.001:0.2;
x=cos(2*pi*f*t);
plot(t,x);
%%%%%%%%%%采样%%%%%%%%%%%%
Fs=50;        
Ts=1/Fs;
n=0:20;
xn=cos(2*pi*f*n*Ts);
figure,stem(n,xn);
% W=(hanning(length(n)))';xn=xn.*W;加窗
%%%%%%%%%%%%
k=-64:64;
X=xn*exp(-j*pi/64).^(n'*k);
magX=abs(X);
w=pi/64*k
figure,stem(w/pi,magX);