大话频谱分析（三）

信号处理是一门实验性很强的学科

先看图

      x(n)=cos(0.48πn)+cos(0.52πn).我们做DFT分析之后的图x(n)=cos(0.48\pi n)+cos(0.52\pi n).我们做DFT分析之后的图

[code]
clear all;clc;
N=11;
n=[0:1:N-1];k=[0:1:N-1];
xn=cos(0.48*pi*n)+cos(0.52*pi*n);
WN=exp(-j*2*pi/N);
nk=n'*k;
WNnk=WN.^nk;
Xk=xn*WNnk;
plot(k,Xk);grid
xlabel('k');title('谱线');ylabel('X(k)');

前面说过，补零没有效果。只是提高分析分辨率。

[code]%x(n)=cos(0.48pi*n)+cos(0.52*pi*n);
clear all;clc;
N=101;
n1=[0:1:10];
n=[0:1:N-1];k=[0:1:N-1];
xn=cos(0.48*pi*n1)+cos(0.52*pi*n1);
xn=[xn,zeros(1,N-length(n1))];
WN=exp(-j*2*pi/N);
nk=n'*k;
WNnk=WN.^nk;
Xk=xn*WNnk;
plot(k,Xk);
xlabel('k');title('谱线');ylabel('X(k)');

下面是增加截断长度后的图

[code]%x(n)=cos(0.48pi*n)+cos(0.52*pi*n);
clear all;clc;
N=101;
n=[0:1:N-1];k=[0:1:N-1];
xn=cos(0.48*pi*n)+cos(0.52*pi*n);
xn=[xn,zeros(1,N-length(n))];
WN=exp(-j*2*pi/N);
nk=n'*k;
WNnk=WN.^nk;
Xk=xn*WNnk;
plot(k,abs(Xk));grid
xlabel('k');title('谱线');ylabel('X(k)');

有没有恍然大悟的感觉？@@@

（待续）